Conceptual density functional theory: status, prospects, issues

Excited-State Polarizabilities: A Combined Density Functional Theory and Information-Theoretic Approach Study

Accurate and efficient determination of excited-state polarizabilities (α) is an open problem both experimentally and computationally. Following our previous work, (Phys. Chem. Chem. Phys. 2023, 25, 2131−2141), in which we employed simple ground-state (S0) density-related functions from the information-theoretic approach (ITA) to accurately and efficiently evaluate the macromolecular polarizabilities, in this work we aimed to predict the lowest excited-state (S1) polarizabilities. The philosophy is to use density-based functions to depict excited-state polarizabilities. As a proof-of-principle application, employing 2-(2′-hydroxyphenyl)benzimidazole (HBI), its substituents, and some other commonly used ESIPT (excited-state intramolecular proton transfer) fluorophores as model systems, we verified that either with S0 or S1 densities as an input, ITA quantities can be strongly correlated with the excited-state polarizabilities. When transition densities are considered, both S0 and S1 polarizabilities are in good relationships with some ITA quantities. The transferability of the linear regression model is further verified for a series of molecules with little or no similarity to those molecules in the training set. Furthermore, the excitation energies can be predicted based on multivariant linear regression equations of ITA quantities. This study also found that the nature of both the ground-state and excited-state polarizabilities of these species are due to the spatial delocalization of the electron density.

Synthesis of Disubstituted Carboxonium Derivatives of Closo-Decaborate Anion [2,6-B10H8O2CC6H5]-: Theoretical and Experimental Study

A comprehensive study focused on the preparation of disubstituted carboxonium derivatives of closo-decaborate anion [2,6-B₁₀H₈O₂CC₆H₅]^- was carried out. The proposed synthesis of the target product was based on the interaction between the anion [B₁₀H₁₁]^- and benzoic acid C₆H₅COOH. It was shown that the formation of this product proceeds stepwise through the formation of a mono-substituted product [B₁₀H₉OC(OH)C₆H₅]^-. In addition, an alternative one-step approach for obtaining the target derivative is postulated. The structure of tetrabutylammonium salts of carboxonium derivative ((C₄H₉)₄N)[2,6-B₁₀H₈O₂CC₆H₅] was established with the help of X-ray structure analysis. The reaction pathway for the formation of [2,6-B₁₀H₈O₂CC₆H₅]^- was investigated with the help of density functional theory (DFT) calculations. This process has an electrophile induced nucleophilic substitution (EINS) mechanism, and intermediate anionic species play a key role. Such intermediates have a structure in which one boron atom coordinates two hydrogen atoms. The regioselectivity for the process of formation for the 2,6-isomer was also proved by theoretical calculations. Generally, in the experimental part, the simple and available approach for producing disubstituted carboxonium derivative was introduced, and the mechanism of this process was investigated with the help of theoretical calculations. The proposed approach can be applicable for the preparation of a wide range of disubstituted derivatives of closo-borate anions.

External fields in conceptual density functional theory

The necessity of the recent incorporation of new external variables in the context of conceptual DFT (CDFT) is discussed based on the ever-increasing portfolio of experimental reaction conditions in the endeavor of experimentalists to synthesize new molecules with unprecedented properties. Electric and magnetic fields (ε and B), mechanical forces (F), and confinement are proposed as valuable new variables, extending conventional CDFT and its associated response functions. A finite field approach is used to calculate the evolution of both global and local descriptors in a selected series of atomic and molecular applications, and from it derive new response function involving, with one exception, the first derivative to the field considered. The electric field results, displaying, for example, a case of a field-induced enantioselectivity in the Fukui function, may be instrumental in the recent upsurge of chemistry in oriented external electric fields. The study of atomic electronegativity and hardness in magnetic fields displays a piecewise behavior, associated to configurational jumps upon increasing field strength and reveals an overall compression of their ranges for stronger fields, which may be guiding upon investigating chemistry in extremely high fields like in white dwarfs. The evolution of the electronegativity and hardness of diatomics under mechanical force can elegantly be traced back to differences in their equilibrium distance in the neutral, cationic, and anionic state. The well-known reduction of the polarizability under confinement can be seen as a fore-runner of the increasing hardness of atoms under pressure, presently under investigation. Periodicity showing up in a spontaneous way in the variety of properties is a leitmotiv in this study, as well as the interconnections/analogies between the different response functions.

Efficient and accurate density-based prediction of macromolecular polarizabilities

Accurately and efficiently predicting macromolecules' polarizabilities is an open problem. In this work, we employ a few simple density-based quantities from the information-theoretic approach (ITA) to predict polarizability of proteins. We first build quantitative structure/property relationships between molecular polarizabilities and ITA quantities. We then verify the broad applicability of ITA quantities for polarizability prediction for inorganic, organic, and biological systems with both localized and delocalized electronic structure. As a proof-of-concept application, we predict the molecular polarizabilities of complex proteins. Based on the linear regression equations for 20 natural amino acid residues, 400 dipeptides, and 8000 tripeptides, one then predicts the molecular polarizability of a larger peptide or even a protein once the molecular wavefunction is obtained. Because it is extremely costly to determine the wavefunction for a macromolecule like a protein, we propose to combine the ITA with the linear-scaling generalized energy-based fragmentation (GEBF) method to predict the macromolecular polarizability. In GEBF, the total molecular polarizability is obtained as a linear combination of the corresponding quantities from a series of small subsystems. We can predict them based on the subsystem wavefunction and linear regression equations rather than compute them from the nearly-intractable coupled-perturbed Hartree-Fock or Kohn-Sham equations for the whole macromolecule. Computational results showcase that the GEBF-ITA protocol should be an inexpensive yet accurate theoretical tool for predicting macromolecular polarizabilities.

Impacts of external fields on aromaticity and acidity of benzoic acid: a density functional theory, conceptual density functional theory and information-theoretic approach study

The impact of external fields on the molecular structure and reactivity properties has been of considerable interest in the recent literature. Benzoic acid as one of the most widely used compounds in medicinal and materials sciences is known for its dual propensity in aromaticity and acidity. In this work, we systematically investigate the impact of a uniform external electric field on these properties. We apply density functional theory, conceptual density functional theory, and an information-theoretic approach to appreciate the change pattern of aromaticity and acidity properties in external fields with different strengths. Our results show that they possess different change patterns under external fields, which can be satisfactorily rationalized by variations in reactivity descriptors and partial charges. The surprising yet novel results from this study should enrich the body of our knowledge about the impact of external fields for different kinds of electronic properties and provide guidance and foundation for future studies of this phenomenon in other molecular systems.

Integrated approach for studying bioactive compounds from Cladosporium spp. against estrogen receptor alpha as breast cancer drug target

Cladosporium spp. have been reported for their great diversity of secondary metabolites which represent as a prominent base material for verifying the biological activities. Several bioactive compounds which have antimicrobial, cytotoxic, quorum sensing inhibitory and phytotoxic activities have been isolated from Cladosporium species. Most of them are still needed to be explored for their anticancer properties. Therefore, the present study is focused on screening and identifying the bioactive compounds of Cladosporium spp. for their anticancer activity via the integrated approaches of Molecular Docking (MD), Molecular Dynamics Simulation (MDS) and Density Functional Theory (DFT) studies. A total of 123 bioactive compounds of Cladosporium spp. were explored for their binding affinity with the selected breast cancer drug target receptor such as estrogen receptor alpha (PDB:6CBZ). The Molecular Docking studies revealed that amongst the bioactive compounds screened, Altertoxin X and Cladosporol H showed a good binding affinity of - 10.5 kcal/mol and - 10.3 kcal/mol, respectively, with the estrogen receptor alpha when compared to the reference compound (17[Formula: see text]-Estradiol: - 10.2 kcal/mol). The MDS study indicated the stable binding patterns and conformation of the estrogen receptor alpha-Altertoxin X complex in a stimulating environment. In addition, in silico absorption, distribution, metabolism, excretion and toxicity (ADMET) study suggested that Altertoxin X has a good oral bioavailability with a high LD[Formula: see text] value of 2.375 mol/kg and did not cause any hepatotoxicity and skin sensitization. In summary, the integrated approaches revealed that Altertoxin X possesses a promising anticancer activity and could serve as a new therapeutic drug for breast cancer treatment.

Development and Applications of the Density-Based Theory of Chemical Reactivity

Density functional theory, which is well-recognized for its accuracy and efficiency, has become the workhorse for modeling the electronic structure of molecules and extended materials in recent decades. Nevertheless, establishing a density-based conceptual framework to appreciate bonding, stability, function, reactivity, and other physicochemical properties is still an unaccomplished task. In this Perspective, we at first provide an overview of the four pathways currently available in the literature to tackle the matter, including orbital-free density functional theory, conceptual density functional theory, direct use of density-associated quantities, and the information-theoretic approach. Then, we highlight several recent advances of employing these approaches to realize new understandings for chemical concepts such as covalent bonding, noncovalent interactions, cooperation, frustration, homochirality, chirality hierarchy, electrophilicity, nucleophilicity, regioselectivity, and stereoselectivity. Finally, we provide a few possibilities for the future development of this relatively uncharted territory. Opportunities are abundant, and they are all ours for the taking.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Electrophilicity of Hoveyda-Grubbs Olefin Metathesis Catalysts as the Driving Force that Controls Initiation Rates

The dissociative mechanism of initiation for a series of Hoveyda-Grubbs type metathesis catalysts modified at the para and meta positions in the isopropoxybenzylidene ligand is investigated by means of DFT calculations. The electron donating/withdrawing capacity of the ligand was screened through the incorporation of various substituents such as halogens, nitro, alkoxides, ketones, esters, amines, and amides. Variations in structural parameters, energy barriers for the Ru-O bond dissociation, and Ru-O bond strength were examined as a function of the Hammett constant. It was found that electronic properties of the catalysts such as chemical potential, hardness, and electrophilicity correlate linearly with the dissociative energy barriers. These findings enable a systematic rationalization and prediction of rate of precatalyst initiation through the calculation of only the HOMO-LUMO gap of catalysts, as the faster the initiation, the more electrophilic the catalyst.

Exploration of Anti-HIV Phytocompounds against SARS-CoV-2 Main Protease: Structure-Based Screening, Molecular Simulation, ADME Analysis and Conceptual DFT Studies

The ever-expanding pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection has gained attention as COVID-19 and caused an emergency in public health to an unmatched level to date. However, the treatments used are the only options; currently, no effective and licensed medications are available to combat disease transmission, necessitating further research. In the present study, an in silico-based virtual screening of anti-HIV bioactive compounds from medicinal plants was carried out through molecular docking against the main protease (M^pro) (PDB: 6LU7) of SARS-CoV-2, which is a key enzyme responsible for virus replication. A total of 16 anti-HIV compounds were found to have a binding affinity greater than -8.9 kcal/mol out of 150 compounds screened. Pseudohypericin had a high affinity with the energy of -10.2 kcal/mol, demonstrating amino acid residual interactions with LEU141, GLU166, ARG188, and GLN192, followed by Hypericin (-10.1 kcal/mol). Moreover, the ADME (Absorption, Distribution, Metabolism and Excretion) analysis of Pseudohypericin and Hypericin recorded a low bioavailability (BA) score of 0.17 and violated Lipinski's rule of drug-likeness. The docking and molecular simulations indicated that the quinone compound, Pseudohypericin, could be tested in vitro and in vivo as potent molecules against COVID-19 disease prior to clinical trials.This was also supported by the theoretical and computational studies conducted. The global and local descriptors, which are the underpinnings of Conceptual Density FunctionalTheory (CDFT) have beenpredicted through successful model chemistry, hoping that they could be of help in the comprehension of the chemical reactivity properties of the molecular systems considered in this study.

Quantitative Structure-Toxicity Relationship in Bioactive Molecules from a Conceptual DFT Perspective

The preclinical drug discovery stage often requires a large amount of costly and time-consuming experiments using huge sets of chemical compounds. In the last few decades, this process has undergone significant improvements by the introduction of quantitative structure-activity relationship (QSAR) modelling that uses a certain percentage of experimental data to predict the biological activity/property of compounds with similar structural skeleton and/or containing a particular functional group(s). The use of machine learning tools along with it has made life even easier for pharmaceutical researchers. Here, we discuss the toxicity of certain sets of bioactive compounds towards Pimephales promelas and Tetrahymena pyriformis in terms of the global conceptual density functional theory (CDFT)-based descriptor, electrophilicity index (ω). We have compared the results with those obtained by using the commonly used hydrophobicity parameter, logP (where P is the n-octanol/water partition coefficient), considering the greater ease of computing the ω descriptor. The Human African trypanosomiasis (HAT) curing activity of 32 pyridyl benzamide derivatives is also studied against Tryphanosoma brucei. In this review article, we summarize these multiple linear regression (MLR)-based QSAR studies in terms of electrophilicity (ω, ω2) and hydrophobicity (logP, (logP)2) parameters.

Elbehairi S, Achar RR, Silina E, Stupin V, Manturova N, Glossman-Mitnik D, Shivamallu C, Kollur SP. Anticholinesterase activity of Areca Catechu: In Vitro and in silico green synthesis approach in search for therapeutic agents against Alzheimer's disease

For many years, the primary focus has been on finding effective treatments for Alzheimer's disease (AD), which has led to the identification of promising therapeutic targets. The necessity for AD stage-dependent optimal settings necessitated a herbal therapy strategy. The plant species Areca Catechu L. (AC) was selected based on the traditional uses against CNS-related diseases. AC leaf extract were prepared using a Soxhlet extraction method and hydroxyapatite nanoparticles (HAp-NPs) were synthesized from the same (AC-HAp-NPs). Powder X-ray diffractometer (XRD), scanning electron microscopy (SEM), transmission electron microscopy (TEM), selected area electron diffraction (SAED) and fourier transform infrared spectroscopy (FTIR) were used to confirm the structure and morphology of the as-prepared AC-HAp-NPs. The crystalline character of the AC-HAp-NPs was visible in the XRD pattern. The synthesized material was found to be nanoflake, with an average diameter of 15-20 nm, according to SEM analysis. The TEM and SAED pictures also revealed the form and size of AC-HAp-NPs. In vitro anti-acetylcholinesterase and butyrylcholinesterase (AChE and BChE) activities of hydroxyapatite nanoparticles produced from an AC leaf extract was tested in this study. When compared to control, AC-HAp-NPs had higher anti-AChE and BChE activity. The anti-acetylcholinesterase action of phytoconstituents generated from AC leaf extract was mediated by 4AQD and 4EY7, according to a mechanistic study conducted utilizing in silico research. The global and local descriptors, which are the underpinnings of Conceptual Density Functional Theory (CDFT), have been predicted through the MN12SX/Def2TZVP/H2O model chemistry to help in the comprehension of the chemical reactivity properties of the five ligands considered in this study. The CDFT experiments are supplemented by the calculation of several useful calculated pharmacokinetics indices, their expected biological targets connected to the bioavailability of the five ligands in order to further the goal of studying their bioactivity.

Application of Fundamental Chemical Principles for Solvation Effects: A Unified Perspective for Interaction Patterns in Solution

We demonstrate the utility of basic chemical principles like the "|Δμ| big is good" (DMB) rule for the study of solvation interactions between distinct solutes such as ions and solvents. The corresponding approach allows us to define relevant criteria for maximum solvation energies of ion pairs in different solvents in terms of electronegativities and chemical hardnesses. Our findings reveal that the DMB principle culminates into the strong and weak acids and bases concept as recently derived for specific ion effects in various solvents. The further application of the DMB approach highlights a similar condition for the chemical hardnesses with a reminiscence to the hard/soft acids and bases principle. Comparable conclusions can also be drawn with regard to the change of the solvent. We show that favorable solvent interactions are mainly driven by low chemical hardnesses as well as high electronegativity differences between the ions and the solvent. Our findings highlight that solvation interactions are governed by basic chemical principles, which demonstrates the close similarity between solvation mechanisms and chemical reactions.

Phytoconstituents of Withania somnifera unveiled Ashwagandhanolide as a potential drug targeting breast cancer: Investigations through computational, molecular docking and conceptual DFT studies

Breast cancer is the second most common malignancy in females worldwide and poses a great challenge that necessitates the identification of novel therapeutic agents from several sources. This research aimed to study the molecular docking and molecular dynamics simulations of four proteins (such as PDB: 6CBZ, 1FDW, 5GWK and 2WTT) with the selected phytochemicals from Withania somnifera to identify the potential inhibitors for breast cancer. The molecular docking result showed that among 44 compounds, two of them, Ashwagandhanolide and Withanolide sulfoxide have the potential to inhibit estrogen receptor alpha (ERα), 17-beta-hydroxysteroid -dehydrogenase type 1 (17β-HSD1), topoisomerase II alpha (TOP2A) and p73 tetramerization domain that are expressed during breast cancer. The molecular dynamics (MD) simulations results suggested that Ashwagandhanolide remained inside the binding cavity of four targeted proteins and contributed favorably towards forming a stable protein-ligand complex throughout the simulation. Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties confirmed that Ashwagandhanolide is hydrophobic and has moderate intestinal permeability, good intestinal absorption, and poor skin permeability. The compound has a relatively low VDss value (-1.652) and can be transported across ABC transporter and good central nervous system (CNS) permeability but did not easily cross the blood-brain barrier (BBB). This compound does not possess any mutagenicity, hepatotoxicity and skin sensitization. Based on the results obtained, the present study highlights the anticancer potential of Ashwagandhanolide, a compound from W. somnifera. Furthermore, in vitro and in vivo studies are necessary to perform before clinical trials to prove the potentiality of Ashwagandhanolide.

On the Periodicity of the Information Theory and Conceptual DFT-Based Reactivity Descriptors

The periodic trends in conceptual density functional and information theory-based reactivity descriptors are reported for the atoms H to Ba (Z = 1 to 56). Ionization potential, electron affinity, electronegativity, and hardness show periodic behavior following the Aufbau principle and popular electronic structure principles. They are in agreement with those reported in standard chemistry textbooks. The trend in the electrophilicity index, however, shows an interesting behavior, where it contradicts earlier reports. Our calculation reveals that the noble gas elements correspond to minimum ω values in each period which obey the minimum electrophilicity principle as well as reflect their low reactivity. Periodic trends in electroaccepting and electrodonating powers, along with that of net electrophilicity, are as expected. The behavior of information theory-based Shannon and GBP entropies, along with the Shannon entropy of shape function are also explored across the periodic table.

Reactive orbital energy theory serving a theoretical foundation for the electronic theory of organic chemistry

It is established that the reactive orbital energy theory (ROET) theoretically reproduces the rule-based electronic theory diagrams of organic chemistry by a comparative study on the charge transfer natures of typical organic carbon-carbon and carbon-heteroatom bond formation reactions: aldol, Mannich, α-aminooxylation, and isogyric reactions. The ROET, which is an expansion of the reaction electronic theories (e.g., the frontier orbital theory) in terms of orbital energies, elucidates the reactive orbitals driving reactions and the charge transferability indices of the reactions. Performing the ROET analyses of these reactions shows that the charge transfer directions given in the rule-based diagrams of the electronic theory are reproduced even for the functional groups of charge transfer destinations in all but only two processes for 38 reaction processes. The ROET analyses also make clear the detailed orbital-based pictures of these bond formation reactions: that is, the use of the out-of-plane antibonding π orbitals in acidic conditions (enol-mode) and in-plane antibonding π orbitals in basic conditions (enolate-mode), which explain the experimentally assumed mechanisms such as the π-bond formations in acidic conditions and σ-bond formations at α-carbons in basic conditions. Furthermore, the ROET analyses explicate that the methyl group initially accepts electrons and then donates them to the bond formations in the target reactions. It is, consequently, suggested that the ROET serves a theoretical foundation for the electronic theory of organic chemistry.

In Vitro Anticancer Screening, Molecular Docking and Antimicrobial Studies of Triazole-Based Nickel(II) Metal Complexes

Herein we describe the synthesis of a series of nickel(II) complexes (C1–C3) with Schiff bases (HL1–HL3) derived from 4-amino-5-mercapto-3-methyl-1,2,4-triazole and ortho/meta/para-nitrobenzaldehyde having composition [Ni(L)₂(H₂O)₂]. The obtained ligands and their complexes were characterized using physico-chemical techniques viz., elemental analysis, magnetic moment study, spectral (electronic, FT-IR, 1H-NMR) and thermal analysis. The elemental analysis and spectral analysis revealed that Schiff bases behave as monoanionic bidentate ligands towards the Ni(II) ion. Whereas, the magnetic moment study suggested the octahedral geometry of all the Ni(II) complexes. The thermal behavior of the complexes has been studied by thermogravimetric analysis and agrees well with the composition of complexes. Further, the biological activities such as antimicrobial and antifungal studies of the Schiff bases and Ni(II) complexes have been screened against bacterial species (Staphylococcus aureus and Pseudomonas aeruginosa) and fungal species (Aspergillus niger and Candida albicans) activity by MIC method, the results of which revealed that metal complexes exhibited significant antimicrobial activities than their respective ligands against the tested microbial species. Furthermore, the molecular docking technique was employed to investigate the active sites of the selected protein, which indeed helped us to screen the potential anticancer agents among the synthesized ligand and complexes. Further, these compounds have been screened for their in vitro anticancer activity using OVCAR-3 cell line. The results revealed that the complexes are more active than the ligands.

Electronic properties of amino acids and nucleobases: similarity classes and pairing principles from chemical reactivity indices

We present a new classification scheme for amino acids and nucleobases based on the electronic properties of the individual molecules. Using chemical reactivity indices such as electronegativity, electrophilicity, and chemical hardness, we can identify similarities and differences between each class of amino acids and nucleobases. Notable differences emerge in particular with regard to high, neutral or low electronegativity as well as different combinations of chemical hardness. Our approach allows us to relate these insights to the properties of the side groups in terms of a unique reference scheme. We further show that hydrophobic differences between amino acids are rather negligible in the context of electronic properties. Our classification scheme also rationalizes the occurrence of distinct stable nucleobase pairs and clearly emphasizes certain differences between individual molecules. The stability and abundant occurrence of Watson-Crick nucleobase pairs is further discussed in the context of the minimum electrophilicity principle.

Properties of the density functional response kernels and its implications on chemistry

An overview of mathematical properties of the non-local second order derivatives of the canonical, grand canonical, isomorphic, and grand isomorphic ensembles is given. The significance of their positive or negative semidefiniteness and the implications of these properties for atoms and molecules are discussed. Based on this property, many other interesting properties can be derived, such as the expansion in eigenfunctions, bounds on the diagonal and off-diagonal elements, and the eigenvalues of these kernels. We also prove Kato's theorem for the softness kernel and linear response and the dissociation limit of the linear responses as the sum of the linear responses of the individual fragments when dissociating a system into two non-interacting molecular fragments. Finally, strategies for the practical calculation of these kernels, their eigenfunctions, and their eigenvalues are discussed.

Rationalizing the Substituent Effects in Diels-Alder Reactions of Triazolinediones with Anthracene

In this work we tackle the problem of the substituent effects in the Diels-Alder cycloadditions between triazolinediones (TADs) and anthracene. Experiments showed that aryl TADs substituted with electron-withdrawing groups (EWG) are more reactive than those substituted with electron-donating (EDG) or alkyl groups. However, the molecular origin of this preference is not yet understood. By a combination of methods including the activation strain model (ASM), energy decomposition analysis (EDA), molecular orbital (MO) theory, and conceptual density functional theory (CDFT), we disclosed the substituent effects of TADs. First, ASM/EDA analysis revealed that the reactivity of alkyl and aryl-substituted TADs is controlled by interaction energies, ΔE_int, which are ultimately defined by orbital interactions between frontier molecular orbitals. Moreover, alkyl-TADs are also controlled by the extent of strain at the transition state. The MO analysis suggested that the rate acceleration for EWG-substituted TADs is due to a more favorable orbital interaction between the HOMO of anthracene and the LUMO of the TADs, which is corroborated by calculations of charge transfer at the transition states. From CDFT, the chemical potential of anthracene is higher than those of TADs, indicating a flow of electron density from anthracene to TADs, in agreement with the results from the electrophilicity index.

From Density Functional Theory to Conceptual Density Functional Theory and Biosystems

The position of conceptual density functional theory (CDFT) in the history of density functional theory (DFT) is sketched followed by a chronological report on the introduction of the various DFT descriptors such as the electronegativity, hardness, softness, Fukui function, local version of softness and hardness, dual descriptor, linear response function, and softness kernel. Through a perturbational approach they can all be characterized as response functions, reflecting the intrinsic reactivity of an atom or molecule upon perturbation by a different system, including recent extensions by external fields. Derived descriptors such as the electrophilicity or generalized philicity, derived from the nature of the energy vs. N behavior, complete this picture. These descriptors can be used as such or in the context of principles such as Sanderson’s electronegativity equalization principle, Pearson’s hard and soft acids and bases principle, the maximum hardness, and more recently, the minimum electrophilicity principle. CDFT has known an ever-growing use in various subdisciplines of chemistry: from organic to inorganic chemistry, from polymer to materials chemistry, and from catalysis to nanotechnology. The increasing size of the systems under study has been coped with thanks to methodological evolutions but also through the impressive evolution in software and hardware. In this flow, biosystems entered the application portfolio in the past twenty years with studies varying (among others) from enzymatic catalysis to biological activity and/or the toxicity of organic molecules and to computational peptidology. On the basis of this evolution, one can expect that “the best is yet to come”.

The roles of charge transfer and polarization in non-covalent interactions: a perspective from ab initio valence bond methods

Noncovalent interactions are ubiquitous and have been well recognized in chemistry, biology and material science. Yet, there are still recurring controversies over their natures, due to the wide range of noncovalent interaction terms. In this Essay, we employed the Valence Bond (VB) methods to address two types of interactions which recently have drawn intensive attention, i.e., the halogen bonding and the CH‧‧‧HC dihydrogen bonding. The VB methods have the advantage of interpreting molecular structures and properties in the term of electron-localized Lewis (resonance) states (structures), which thereby shed specific light on the alteration of the bonding patterns. Due to the electron localization nature of Lewis states, it is possible to define individually and measure both polarization and charge transfer effects which have different physical origins. We demonstrated that both the ab initio VB method and the block-localized wavefunction (BLW) method can provide consistent pictures for halogen bonding systems, where strong Lewis bases NH₃, H₂O and NMe₃ partake as the halogen bond acceptors, and the halogen bond donors include dihalogen molecules and XNO₂ (X = Cl, Br, I). Based on the structural, spectral, and energetic changes, we confirm the remarkable roles of charge transfer in these halogen bonding complexes. Although the weak C-H∙∙∙H-C interactions in alkane dimers and graphene sheets are thought to involve dispersion only, we show that this term embeds delicate yet important charge transfer, bond reorganization and polarization interactions.

Conceptual density functional theory under pressure: Part I. XP-PCM method applied to atoms

High pressure chemistry offers the chemical community a range of possibilities to control chemical reactivity, develop new materials and fine-tune chemical properties. Despite the large changes that extreme pressure brings to the table, the field has mainly been restricted to the effects of volume changes and thermodynamics with less attention devoted to electronic effects at the molecular scale. This paper combines the conceptual DFT framework for analyzing chemical reactivity with the XP-PCM method for simulating pressures in the GPa range. Starting from the new derivatives of the energy with respect to external pressure, an electronic atomic volume and an atomic compressibility are found, comparable to their enthalpy analogues, respectively. The corresponding radii correlate well with major known sets of this quantity. The ionization potential and electron affinity are both found to decrease with pressure using two different methods. For the electronegativity and chemical hardness, a decreasing and increasing trend is obtained, respectively, and an electronic volume-based argument is proposed to rationalize the observed periodic trends. The cube of the softness is found to correlate well with the polarizability, both decreasing under pressure, while the interpretation of the electrophilicity becomes ambiguous at extreme pressures. Regarding the electron density, the radial distribution function shows a clear concentration of the electron density towards the inner region of the atom and periodic trends can be found in the density using the Carbó quantum similarity index and the Kullback-Leibler information deficiency. Overall, the extension of the CDFT framework with pressure yields clear periodic patterns.

[Fe4S4] cubane in sulfite reductases: new insights into bonding properties and reactivity

The dissimilatory sulfite reductase enzyme has very characteristic active site where the substrate binds to an iron site, ligated by a siroheme macrocycle and a thiol directly connected to a [Fe₄S₄] cluster. This arrangement gives the enzyme remarkable efficiency in reducing sulfite and nitrite all the way to hydrogen sulfide and ammonia. For the first time we present a theoretical study where substrate binding modalities and activation are elucidated using active site models containing proton supply side chains and the [Fe₄S₄] cluster. Density functional theory (DFT) was deployed in conjunction with the energy decomposition scheme (as implemented in AMS), the quantum theory of atoms in molecules (QTAIM), and conceptual DFT (cDFT) descriptors. We quantified the role of the electrostatic interactions inside the active site created by the side chains as well as the influence of the [Fe₄S₄] cluster on the substrate binding. Furthermore, using conceptual DFT results we shed light of the activation process, thus, laying foundation for further mechanistic studies. We found that the bonding of the ligands to the iron complex is dominated by electrostatic interactions, but the presence of the [Fe₄S₄] cubane leads to substantial changes in electronic interaction. The spin state of the cubane, however, affects the binding energy only marginally. The conceptual DFT results show that the presence of the [Fe₄S₄] cubane affects the reactivity of the active site as it is involved in electron transfer. This is corroborated by an increase in the electrophilicity index, thus making the active site more prone to react with the ligands. The interaction energies between the ligand and the siroheme group are also increased upon the presence of the cubane group, thus, suggesting that the siroheme group is not an innocent spectator but plays an active role in the reactivity of the dSIR active site.

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

An insight into reactivity and bioactivity properties of quorum sensing peptides against PDE10A: a computational peptidology approach

Peptides are currently the most promising lead molecules. Quorum sensing peptides have a variety of structural features and are regularly exposed to post-translational modifications. Antiparkinsonian drugs lose their efficacy after a long period of use, and patients develop motor problems such as drug-induced dyskinesia (DIDs). The interaction between PDE10A and cAMP is necessary for dopamine neurotransmission and may play a role in Parkinson's disease pathogenesis. cAMP and cGMP are cyclic nucleotides that act as secondary messengers in the signal transduction pathway, influencing a range of CNS activities. PDE enzymes hydrolyze phosphodiester bonds to break down cAMP and cGMP, allowing them to control intracellular levels of these second messengers effectively. PDE expression, and hence cyclic nucleotide levels and their downstream targets, may change with age and in numerous age-related illnesses, including Parkinson's disease, according to mounting evidence. At the peak of dyskinesias, cyclic nucleotide levels were lower, and using phosphodiesterase inhibitors before antiparkinsonian medicines reduced the severity of dyskinesias. In a recent study, PapRIV was found to have the ability to activate BV-2 microglia cells, indicating that this quorum sensing peptide may play a role in gut-brain contact. As a result of the current in silico work, mainly focused on QSPs as a lead molecule for inhibiting PDE10A, the SRNAT QSP sequence has been a potent molecule in molecular docking and molecular dynamics simulations. Furthermore, we can test the efficiency of therapeutic components in vitro and in vivo utilizing this computational approach against PDE10A.

The structure-activity relationship of aromatic compounds in advanced oxidation processes：a review

Advanced oxidation processes (AOPs) are widely used as efficient technologies to treat highly toxic and harmful substances in wastewater. Taking the most representative aromatic compounds (monosubstituted benzenes, substituted phenols and heterocyclic compounds) as examples, this paper firstly introduces their structures and the structural descriptors studied in AOPs before, and the influence of structural differences in AOPs with different reactive oxygen species (ROS) on the degradation rate was discussed in detail. The structure-activity relationship of pollutants has been previously analyzed through quantitative structure-activity relationship (QSAR) model, in which ROS is a very important influencing factor. When electrophilic oxidative species attacks pollutants, aromatic compounds with electron donating groups are more favorable for degradation than aromatic compounds with electron donating groups. While nucleophilic oxidative species comes to the opposite conclusion. The choice of advanced oxidation processes, the synergistic effect of various active oxygen species and the used catalysts will also change the degradation mechanism. This makes the structure-dependent activity relationship uncertain, and different conclusions are obtained under the influence of various experimental factors.

Investigation of water substitution at RuII complexes by conceptual density function theory approach

In this paper, we investigated water exchange reactions and substitution of aqua Ru^II complexes of general formula [Ru(terpy)(N^N)(H₂ O)]²⁺ (where N^N = ethylenediamine (en), 1,2-(aminomethyl)pyridine (ampy) and 2,2'-bipyridine (bipy)) by ammonia and thioformaldehyde. These reactions were studied in detail by applying conceptual density functional theory. This approach enabled us to gain further insight into the underlying reaction mechanism at the microscopic level (involving only direct participants of the reaction, without the influence of the solvent) and to put the concept of reaction mechanism on a quantitative basis. The course of the chemical reaction along the reaction coordinate ξ, is rationalized in terms of reaction energy, force, dipole moment, and reaction electronic flux (REF). The results yield and characterize the significant influence of an intermolecular hydrogen bond formed between the entering and the spectator ligand to the overall energy barrier of the reactions.

Meta-Hybrid Density Functional Theory Prediction of the Reactivity, Stability, and IGM of Azepane, Oxepane, Thiepane, and Halogenated Cycloheptane

The application of plain cycloalkanes and heterocyclic derivatives in the synthesis of valuable natural products and pharmacologically active intermediates has increased tremendously in recent times with much attention being paid to the lower cycloalkane members. The structural and molecular properties of higher seven-membered and nonaromatic heterocyclic derivatives are less known despite their stable nature and vast application; thus, an insight into their structural and electronic properties is still needed. Appropriate quantum chemical calculations utilizing the ab initio (MP2) method, meta-hybrid (M06-2X) functional, and long-range-separated functionals (ωB97XD) have been utilized in this work to investigate the structural reactivity, stability, and behavior of substituents on cycloheptane (CHP) and its derivatives: azepane, oxepane, thiepane, fluorocycloheptane (FCHP), bromocycloheptane (BrCHP), and chlorocycloheptane (ClCHP). Molecular global reactivity descriptors such as Fukui function, frontier molecular orbitals (FMOs), and molecular electrostatic potential were computed and compared with lower members. The results of two population methods CHELPG and Atomic Dipole Corrected Hirshfeld Charges (ADCH) were equally compared to scrutinize the charge distribution in the molecules. The susceptibility of intramolecular interactions between the substituents and cycloalkane ring is revealed by natural bond orbital analysis and intramolecular weak interactions by the independent gradient model (IGM). Other properties such as atomic density of states, intrinsic bond strength index (IBSI), and dipole moments are considered. It is acclaimed that the strain effect is a major determinant effect in the energy balance of cyclic molecules; thus, the ring strain energies and validation of spectroscopic specificities with reference to the X-ray crystallographic data are also considered.

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.

Assessing Alkene Reactivity toward Cytochrome P450-Mediated Epoxidation through Localized Descriptors and Regression Modeling

The prediction of sites of epoxidation by cytochrome P450s during metabolism is particularly important in drug design, as epoxides are capable of alkylating biological macromolecules. Reliable methods are needed to quantitatively predict P450-mediated epoxidation barriers for inclusion in high-throughput screening campaigns alongside protein-ligand docking. Utilizing the fractional occupation number weighted density (FOD) and orbital-weighted Fukui index (f_w⁺) as descriptors of local reactivity and a data set of 36 alkene epoxidation barriers computed with density functional theory (DFT), we developed and validated a multiple linear regression model for the reliable estimation of epoxidation barriers using only substrate structures as input. Using our recommended level of theory (GFN2-xTB//GFN-FF), mean absolute errors in the training and test sets were found to be 0.66 and 0.70 kcal/mol, respectively, with coefficients of determination of ca. 0.80. We demonstrate the utility of this approach on three known substrates of CYP101A1 and further show that this approach is inappropriate for particularly electron-rich alkenes. By employing a modern semiempirical method on force-field-generated geometries, the required descriptors can be calculated on the millisecond timescale per structure, making the approach well suited for incorporation into high-throughput methodologies alongside docking.

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.

From the Electron Density Gradient to the Quantitative Reactivity Indicators: Local Softness and the Fukui Function

Important reactivity measures such as the local softness, the Fukui function, and the global hardness have been calculated directly from first principles with the use of the electron density function, beyond the finite difference approximation. Our recently derived density gradient theorem and the principle of nearsightedness of the electronic matter have been instrumental in obtaining the original, albeit approximate, result on the local softness of an atom. By integration of the local softness s(r), we obtain the global softness S and the Fukui function f(r) = s(r)/S. Local and global softness values have also been calculated analytically for the basic hydrogenic orbitals; the general relation to the atomic number S = σZ ^-2 has been demonstrated, with constants σ characteristic for each orbital type. Global hardness η = 1/S calculated for atoms and ions has been favorably tested against its conventional measure given by the finite difference approximation: (I - A). Calculated test results for atoms and ions in rows 1-4 of the periodic table have been presented.

Reactivity and a Charge-Transfer Model Analysis in Aminopolycarboxylic-Metal Complexes

In the present work, we have calculated several density functional theory (DFT) reactivity descriptors for the aminopolycarboxylate (APC) acids at the B3LYP/6311++G (d,p) levels of theory, aiming to analyze their reactivity. Reactivity descriptors such as ionization energy, molecular hardness, electrophilicity, and condensed Fukui function local indices have been determined to predict the reactivity of APCs. The influence of the solvent was taken into account by employing the CPCM model. The results indicate that the solvation slightly modifies the tendency of the reactivity of the APCs studied. On the other hand, we applied a global and local charge-transfer partitioning model, which introduces two charge-transfer channels [one for accepting electrons (electrophilic) and another for donating one (nucleophilic)] to the complexation reaction of a set of APC acids with transition metals (Mn, Co, and Ni targets enlarged by Fe, Cu, and Zn). The correlation between the charges obtained for the interaction between APC acids and transition metal stability constants provides support for their interpretation as measures of the electrophilicity and nucleophilicity of a chemical species and, at the same time, allows one to describe the donation and back-donation processes in terms of the DFT of chemical reactivity. Also, the application of dual descriptors for these acids provides valuable information concerning the atoms in the reactants playing the most important roles in the reaction, thus helping to improve our understanding of the reaction under study.

N-((1H-Pyrrol-2-yl)methylene)-6-methoxypyridin-3-amine and Its Co(II) and Cu(II) Complexes as Antimicrobial Agents: Chemical Preparation, In Vitro Antimicrobial Evaluation, In Silico Analysis and Computational and Theoretical Chemistry Investigations

Researchers are interested in Schiff bases and their metal complexes because they offer a wide range of applications. The chemistry of Schiff bases of heterocompounds has got a lot of attention because of the metal's ability to coordinate with Schiff base ligands. In the current study, a new bidentate Schiff base ligand, N-((1H-pyrrol-2-yl)methylene)-6-methoxypyridin-3-amine (MPM) has been synthesized by condensing 6-methoxypyridine-3-amine with pyrrole-2-carbaldehyde. Further, MPM is used to prepare Cu(II) and Co(II) metal complexes. Analytical and spectroscopic techniques are used for the structural elucidation of the synthesized compounds. Both MPM and its metal complexes were screened against Escherichia coli, Bacillus subtilis, Staphylococcus aureus and Klebsiella pneumoniae species for antimicrobial studies. Furthermore, these compounds were subjected to in silico studies against bacterial proteins to comprehend their best non-bonded interactions. The results confirmed that the Schiff base ligand show considerably higher binding affinity with good hydrogen bonding and hydrophobic interactions against various tested microbial species. These results were complemented with a report of the Conceptual DFT global reactivity descriptors of the studied compounds together with their biological scores and their ADMET computed parameters.

Exploration of plant-derived natural polyphenols toward COVID-19 main protease inhibitors: DFT, molecular docking approach, and molecular dynamics simulations

Recent outbreaks of coronavirus have brought serious challenges to public health around the world, and it is essential to find effective treatments. In this study, the 3C-like proteinase (3CLpro) of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has been considered as an important drug target because of its role in viral replication. We initially optimized 251 compounds at the PM7 level of theory for docking with 3CLpro, and then we selected the top 12 compounds for further optimization with the B3LYP-D3/6-311G** method and obtained the top four compounds by further molecular docking. Quantum chemistry calculations were performed to predict molecular properties, such as the electrostatic potential and some CDFT descriptors. We also performed molecular dynamics simulations and free energy calculations to determine the relative stability of the selected four potential compounds. We have identified key residues controlling the 3CLpro/ligand binding from per-residue based decomposition of the binding free energy. Convincingly, the comprehensive results support the conclusion that the compounds have the potential to become a candidate for anti-coronavirus treatment.