Conceptual density functional theory based electronic structure principles

Adiabatic models for the quantum dynamics of surface scattering with lattice effects

In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.

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.

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.

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.

Phenylamine/Amide Grafted in Silica as Sensing Nanocomposites for the Removal of Carbamazepine: A DFT Approach

This study aimed to remove carbamazepine from aqueous solutions, using functional silica phenylamine (SiBN), which is characterized and showed excellent chemical and thermal stability. Adsorbents based on silica were developed due to their unusually large surface area, homogenous pore structure, and well-modified surface properties, as silica sparked tremendous interest. It was determined to develop a novel silica adsorbent including phenylamine and amide (SiBCON). The adsorbents obtained were analyzed by various spectroscopy devices, including SEM, FT-IR and TGA analysis. The maximum removal rates for carbamazepine were 98.37% and 98.22% for SiBN and SiBCON, respectively, when optimized at room temperature, pH 9.0, initial concentration of 10 mg·L−1 and contact time of 15 min. Theoretical tools are widely used in the prediction of the power of interactions between chemical systems. The computed data showed that new amine modified silica is quite effective in terms of the removal of carbamazepine from aqueous solution. Calculation binding energies and DFT data showed that there is a powerful interaction between amine-modified silica and carbamazepine.

Harvesting the fragment-based nature of bifunctional organocatalysts to enhance their activity

Enhancing the activity of bifunctional organocatalysts: a fragment-based approach coupled with activity maps helps identifying better-performing catalytic motifs.

A straightforward method to quantify the electron-delocalizing ability of pi-conjugated molecules

Electronic delocalization is essential to the properties of pi-conjugated molecules. We introduce the inter-fragment delocalization index (IFDI) as an easy-to-use computational method for quantifying the electronic delocalization in pi-conjugated oligomers...

XNgNSi (X = HCC, F; Ng = Kr, Xe, Rn): A New Class of Metastable Insertion Compounds Containing Ng-C/F and Ng-N Bonds and Possible Isomerization therein

Recently, astronomically important silaisocyanoacetylene (HCCNSi) possessing a large dipole moment has been detected for the first time with the help of crossed molecular beam experiments. Quantum chemical computations at higher levels of theory have also been performed to characterize the transient species. In this study, we have analyzed the equilibrium geometry, stability, reactivity, and energetics as well as the nature of bonding in the noble gas (Ng) inserted HCCNSi compound. We have also considered its F analogue to understand the influence of the most electronegative atom in the compound. Metastable behavior of the XNgNSi compounds (X = HCC, F; Ng = Kr-Rn) is examined by calculating thermochemical parameters like free energy change (ΔG) and zero-point-energy-corrected dissociation energy (D₀) at 298 K for all possible two-body (2B) and three-body (3B) (both neutral as well as ionic) dissociation channels using coupled-cluster theory [CCSD(T)] in addition to density functional theory (DFT) as well as second order Møller-Plesset perturbation theory (MP2). The set of predicted compounds is found to be endergonic in nature, having high positive free energy change suggesting the thermochemical stability of the compounds except for the 2B Ng-release paths. Though thermodynamically feasible, they are kinetically protected with very high activation free energy barriers. Interestingly, the release of Ng from the parent moiety XNgNSi produces the XSiN isomer, by 180° flipping of the NSi moiety. This can also be seen in the dynamical simulation carried out with the help of atom-centered density matrix propagation (ADMP) technique at 2000K for 1 ps. The bonding in Ng-C, Ng-F, and Ng-N bonds of the studied compounds is analyzed and described with the aid of natural bond orbital (NBO), topological parameters computed using atoms-in-molecules theory (AIM), energy decomposition analysis (EDA), and adaptive natural density partitioning (AdNDP) methods. The natural charge distribution on the constituent atoms suggests that the compounds can be partitioned into both ways of representations, viz., neutral radical as well as ionic fragments. Lastly, the reactivity of the compounds is scrutinized using certain reactivity descriptors calculated within the domain of conceptual density functional theory (CDFT).

Virtual Screening for Potential Phytobioactives as Therapeutic Leads to Inhibit NQO1 for Selective Anticancer Therapy

NAD(P)H:quinone acceptor oxidoreductase-1 (NQO1) is a ubiquitous flavin adenine dinucleotide-dependent flavoprotein that promotes obligatory two-electron reductions of quinones, quinonimines, nitroaromatics, and azo dyes. NQO1 is a multifunctional antioxidant enzyme whose expression and deletion are linked to reduced and increased oxidative stress susceptibilities. NQO1 acts as both a tumor suppressor and tumor promoter; thus, the inhibition of NQO1 results in less tumor burden. In addition, the high expression of NQO1 is associated with a shorter survival time of cancer patients. Inhibiting NQO1 also enables certain anticancer agents to evade the detoxification process. In this study, a series of phytobioactives were screened based on their chemical classes such as coumarins, flavonoids, and triterpenoids for their action on NQO1. The in silico evaluations were conducted using PyRx virtual screening tools, where the flavone compound, Orientin showed a better binding affinity score of -8.18 when compared with standard inhibitor Dicumarol with favorable ADME properties. An MD simulation study found that the Orientin binding to NQO1 away from the substrate-binding site induces a potential conformational change in the substrate-binding site, thereby inhibiting substrate accessibility towards the FAD-binding domain. Furthermore, with this computational approach we are offering a scope for validation of the new therapeutic components for their in vitro and in vivo efficacy against NQO1.

Synthesis, Computational Pharmacokinetics Report, Conceptual DFT-Based Calculations and Anti-Acetylcholinesterase Activity of Hydroxyapatite Nanoparticles Derived From Acorus Calamus Plant Extract

Over the years, Alzheimer's disease (AD) treatments have been a major focus, culminating in the identification of promising therapeutic targets. A herbal therapy approach has been required by the demand of AD stage-dependent optimal settings. Present study describes the evaluation of anti-acetylcholinesterase (AChE) activity of hydroxyapatite nanoparticles derived from an Acorus calamus rhizome extract (AC-HAp NPs). The structure and morphology of as-prepared (AC-HAp NPs) was confirmed using powder X-ray diffractometer (XRD), scanning electron microscopy (SEM), transmission electron microscopy (TEM) and high-resolution transmission electron microscopy (HR-TEM). The crystalline nature of as-prepared AC-HAp NPs was evident from XRD pattern. The SEM analysis suggested the spherical nature of the synthesized material with an average diameter between 30 and 50 nm. Further, the TEM and HR-TEM images revealed the shape and size of as-prepared (AC-HAp NPs). The interplanar distance between two lattice fringes was found to be 0.342 nm, which further supported the crystalline nature of the material synthesized. The anti-acetylcholinesterase activity of AC-HAp NPs was greater as compared to that of pure HAp NPs. The mechanistic evaluation of such an activity carried out using in silico studies suggested that the anti-acetylcholinesterase activity of phytoconstituents derived from Acorus calamus rhizome extract was mediated by BNDF, APOE4, PKC-γ, BACE1 and γ-secretase proteins. 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. With the further objective of analyzing their bioactivity, the CDFT studies are complemented with the estimation of some useful computed pharmacokinetics indices, their predicted biological targets, and the ADMET parameters related to the bioavailability of the five ligands are also reported.

An integrated molecular modeling protocol for drug screening based on conceptual density functional theory and chemoinformatics for the study of marine cyclopeptides

An integrated molecular modeling protocol resulting from the combination of conceptual density functional theory (CDFT) chemical reactivity descriptors with several chemoinformatics tools has been used for the study of the chemical reactivity and bioactivity properties of a group of marine cyclic peptides. CP-CDFT is a branch of computational chemistry and molecular modeling dedicated to the study of peptides. The protocol allowed the estimation of the CDFT-based reactivity indices together with the associated physicochemical parameters that can help to identify the ability of the studied peptides to behave as potential useful drugs. This was complemented with an analysis of the bioactivity and pharmacokinetics parameters related to the ADMET (absorption, distribution, metabolism, excretion, and toxicity) features. Some examples related to the ability of the CDFT-based chemical reactivity descriptors for the prediction of the pKas of the peptides as well as their potential as AGE inhibitors are also presented.

In Silico Studies on Selected Neutral Molecules, CGa2Ge2, CAlGaGe2, and CSiGa2Ge Containing Planar Tetracoordinate Carbon

Density functional theory (DFT) was used to study the structure, stability, and bonding in some selected neutral pentaatomic systems, viz., CGa2Ge2, CAlGaGe2, and CSiGa2Ge containing planar tetracoordinate carbon. The systems are kinetically stable, as predicted from the ab initio molecular dynamics simulations. The natural bond orbital (NBO) analysis showed that strong electron donation occurs to the central planar carbon atom by the peripheral atoms in all the studied systems. From the nucleus independent chemical shift (NICS) analysis, it is shown that the systems possess both σ- and π- aromaticity. The presence of 18 valence electrons in these systems, in their neutral form, appears to be important for their stability with planar geometries rather than tetrahedral structures. The nature of bonding is understood through the adaptive natural density partitioning analysis (AdNDP), quantum theory of atoms in molecules (QTAIM) analysis, and also via Wiberg bond index (WBI) and electron localization function (ELF).

In Silico Pharmacokinetics, ADMET Study and Conceptual DFT Analysis of Two Plant Cyclopeptides Isolated From Rosaceae as a Computational Peptidology Approach

This research presents the outcomes of a computational determination of the chemical reactivity and bioactivity properties of two plant cyclopeptides isolated from Rosaceae through the consideration of Computational Peptidology (CP), a protocol employed previously in the research of similar molecular systems. CP allows the prediction of the global and local descriptors that are the integral foundations of Conceptual Density Functional Theory (CDFT) and which could help in getting in the understanding of the chemical reactivity properties of the two plant cyclopeptides under study, hoping that they could be related to their bioactivity. The methodology based on the Koopmans in DFT (KID) approach and the MN12SX/Def2TZVP/H2O model chemistry has been successfully validated. Various Chemoinformatics tools have been used to improve the process of virtual screening, thus identifying some additional properties of these two plant cyclopeptides connected to their ability to behave as potentially useful drugs. With the further objective of analyzing their bioactivity, the CP protocol is complemented with the estimation of some useful parameters related to pharmacokinetics, their predicted biological targets, and the Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) parameters related to the bioavailability of the two plant cyclopeptides under study are also reported.

Atomic Clusters: Structure, Reactivity, Bonding, and Dynamics

Atomic clusters lie somewhere in between isolated atoms and extended solids with distinctly different reactivity patterns. They are known to be useful as catalysts facilitating several reactions of industrial importance. Various machine learning based techniques have been adopted in generating their global minimum energy structures. Bond-stretch isomerism, aromatic stabilization, Rener-Teller effect, improved superhalogen/superalkali properties, and electride characteristics are some of the hallmarks of these clusters. Different all-metal and nonmetal clusters exhibit a variety of aromatic characteristics. Some of these clusters are dynamically stable as exemplified through their fluxional behavior. Several of these cluster cavitands are found to be agents for effective confinement. The confined media cause drastic changes in bonding, reactivity, and other properties, for example, bonding between two noble gas atoms, and remarkable acceleration in the rate of a chemical reaction under confinement. They have potential to be good hydrogen storage materials and also to activate small molecules for various purposes. Many atomic clusters show exceptional opto-electronic, magnetic, and nonlinear optical properties. In this Review article, we intend to highlight all these aspects.

Substituent Effects on Electride Characteristics of Mg2(η5-C5H5)2: A Theoretical Study

An ab initio study has been carried out on the substituted binuclear sandwich complexes of Mg₂(η⁵-C₅H₅)₂. We have checked whether the substitution destroys the electride properties of a complex, as it needs to satisfy several stringent criteria to obtain the status of an electride. The thermochemical results show that the complexes are stable at room temperature and 1 atm pressure. From the analysis of the various electron density descriptors and the natural bond orbital (NBO) for all the complexes, it is confirmed that the Mg-Mg bonds are covalent and the metal-ligand bonds are ionic in nature. The charges on each Mg atom in the studied complexes are +1. Analysis of the electron density descriptors shows the presence of a non-nuclear attractor (NNA) at the middle of the bond formed by the two Mg atoms when attached to the ligands. The electride characteristics are exhibited by all of the designed complexes. We also report the aromaticity behavior and reactivity descriptors of these complexes. The electride characteristics of Mg₂(η⁵-C₅H₅)₂ complex get affected on substitution, as both the NNA population and the nonlinear optical properties (NLO) of the complexes are changed.

Local Temperature as a Chemical Reactivity Descriptor

Using the electron density and its associated quantities in a molecular system to quantify chemical reactivity in density functional theory is of considerable recent interest. Local temperature based on the kinetic energy density is an intrinsic property of a molecular system, which can be employed for this purpose. In this work, we explore such a possibility. To this end, we examine the local behavior of local temperature with a few choices of the kinetic energy density, apply it to determine regioselectivity of nucleophilic and electrophilic compounds, and then investigate its performance in appreciating reactions along the intrinsic reaction pathway for exothermic, endothermic, and thermoneutral transformations. Our results confirm that local temperature can be used as an effective descriptor of molecular reactivity.