Does electron-correlation has any role in the quantitative structure-activity relationships?

Theoretical investigation of substitution effect on the sixth and seventh positions of coumarin derivatives

The linear and non-linear optical properties of 6-donor and 7-donor substituted coumarins were compared using density functional theory (DFT) and time-dependent-DFT (TD-DFT). Charge transfer characteristics were investigated through natural bond order analysis, frontier molecular orbital, and molecular electrostatic potential plots. TD-DFT results suggested that the 6-donor substituted coumarins (PS1, PS3, and PS5) showed red-shifted absorption than the 7-donor substituted coumarins (PS2, PS4, and PS6). The chemical potential (μ) and electrophilicity index (ω) showed direct relation with the band gap and an inverse relation with chemical hardness (η) and hyperhardness (Γ). The global reactivity descriptors μ and ω showed direct and η and Γ showed an inverse correlation with first-order hyperpolarizability (β0) and second-order hyperpolarizability (γ). The β0 and γ for 7-donor substituted coumarin are higher than for 6-donor substituted coumarin.

Machine learning, artificial intelligence, and chemistry: how smart algorithms are reshaping simulation and the laboratory

Machine learning and artificial intelligence are increasingly gaining in prominence through image analysis, language processing, and automation, to name a few applications. Machine learning is also making profound changes in chemistry. From revisiting decades-old analytical techniques for the purpose of creating better calibration curves, to assisting and accelerating traditional in silico simulations, to automating entire scientific workflows, to being used as an approach to deduce underlying physics of unexplained chemical phenomena, machine learning and artificial intelligence are reshaping chemistry, accelerating scientific discovery, and yielding new insights. This review provides an overview of machine learning and artificial intelligence from a chemist’s perspective and focuses on a number of examples of the use of these approaches in computational chemistry and in the laboratory.

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.

Concentration-dependent adsorption of organic contaminants by graphene nanosheets: quantum-mechanical models

Adsorption is the key process in the expression of environmentally relevant physicochemical and toxicological properties of carbon nanomaterials. However, the adsorption of organic contaminants on to nanomaterials is a highly complex phenomenon, owing to the heterogeneity of adsorption sites, for example, on graphene surface as well as due to multiple factors operative during the adsorption, particularly, at the quantum-mechanical level. For predicting the concentration-dependent adsorption coefficients of organic contaminants by carbon nanomaterials, one option has been to rely on the existing linear-solvation energy relationship (LSER) models. The present work on the adsorption of aromatic and aliphatic organic contaminants by graphene nanosheets reveals that the existing LSER models are prone to failure when tested for internal and external validation using an external prediction set of compounds unknown to the model. As an alternative to the LSERs, the present work reports pure quantum-mechanical models developed using computational only quantum-mechanical descriptors. The reliability of the quantum-mechanical models was tested using state-of-the-art validation procedures employing an external prediction set of compounds. The proposed quantum-mechanical models reveal mean polarizability, zero-point vibrational energy, and its electron-correlation contribution to be the key descriptors in the prediction of adsorption coefficients of organic contaminants by graphene nanosheets.

Towards an Understanding of the Mode of Action of Human Aromatase Activity for Azoles through Quantum Chemical Descriptors-Based Regression and Structure Activity Relationship Modeling Analysis

Aromatase is an enzyme member of the cytochrome P450 superfamily coded by the CYP19A1 gene. Its main action is the conversion of androgens into estrogens, transforming androstenedione into estrone and testosterone into estradiol. This enzyme is present in several tissues and it has a key role in the maintenance of the balance of androgens and estrogens, and therefore in the regulation of the endocrine system. With regard to chemical safety and human health, azoles, which are used as agrochemicals and pharmaceuticals, are potential endocrine disruptors due to their agonist or antagonist interactions with the human aromatase enzyme. This theoretical study investigated the active agonist and antagonist properties of “chemical classes of azoles” to determine the relationships of azole interaction with CYP19A1, using stereochemical and electronic properties of the molecules through classification and multilinear regression (MLR) modeling. The antagonist activities for the same substituent on diazoles and triazoles vary with its chemical composition and its position and both heterocyclic systems require aromatic substituents. The triazoles require the spherical shape and diazoles have to be in proper proportion of the branching index and the number of ring systems for the inhibition. Considering the electronic aspects, triazole antagonist activity depends on the electrophilicity index that originates from interelectronic exchange interaction (ω_HF) and the LUMO energy (ELUMOPM7), and the diazole antagonist activity originates from the penultimate orbital (EHOMONLPM7) of diazoles. The regression models for agonist activity show that it is opposed by the static charges but favored by the delocalized charges on the diazoles and thiazoles. This study proposes that the electron penetration of azoles toward heme group decides the binding behavior and stereochemistry requirement for antagonist activity against CYP19A1 enzyme.

QSAR/QSPR models based on quantum chemistry for risk assessment of pesticides according to current European legislation

In Europe, agencies and official organizations involved in the pesticide control such as the EFSA, ECHA, JRC and ECETOC or even the OECD are pointing out that the software tools based on quantitative structure relationship models, i.e. QSAR and QSPR, have a huge potential to improve the pesticide risk assessment process. In this sense, these non-animal test methods can promote the competitiveness of agriculture in this region: the consumer safety is increased with them due to the possibility of perform an overall better risk assessment of the degradation products and metabolites from pesticides. However, the use of theses computational-based (in silico) tools must be much more systematised and harmonised, improving their validation and including case studies to test them. To open databases, incorporating critical data in an orderly manner for building the models, becomes also necessary. Moreover, quantum chemistry through the Density Functional Theory should be promoted as tool for calculation of quantum descriptors, especially for the study of similar compounds with the same carbon skeleton but differing substitution patterns, e.g. isomers.

Externally predictive quantum-mechanical models for the adsorption of aromatic organic compounds by graphene-oxide nanomaterials

Graphene oxide is most often chosen as an alternative to graphene in the applications of carbon-based nanomaterials where adsorption is the primary process. However, its adsorption properties are poorly understood. The existing reports on the adsorption mechanism of graphene oxide rely on the linear free-energy/solvation-energy relationship (LFER/LSER) models. This computational work explores the role of quantum mechanical descriptors in the adsorption of aromatic organic compounds by graphene-oxide. For this, externally predictive quantitative models based on quantum-mechanical descriptors are developed and compared with the existing LSERs for the prediction of adsorption coefficients of organic compounds at three different adsorbate concentrations. The predictivity of the models is assessed using an external prediction set of compounds not used for developing the models. Notably, the mean polarizability, but originating from the quantum mechanical exchange interactions (between electrons of parallel spin), is found to be the most significant factor in driving the adsorption on graphene oxide. The present work also proposes quantum-mechanical-LSER models based on a combination of quantum-mechanical and LSER descriptors, which are in fact found to be equally predictive as the existing LSERs. The quantum-mechanical models proposed in this work are further utilized for the prediction of adsorption coefficients of aliphatic compounds.

Quantum-mechanical LSERs for the concentration-dependent adsorption of aromatic organic compounds by activated carbon: Applications and comparison with carbon nanotubes

Carbon nanotubes (CNTs) have taken precedence over activated carbon in various applications where adsorption is the primary process. The adsorption of chemical compounds by CNTs and activated carbon is most often predicted through linear free energy/solvation energy relationships (LFERs/LSERs). This work proposes quantum-mechanical LSER models based on a combination of quantum-mechanical descriptors and solvatochromic descriptors of LSERs for predicting the adsorption of aromatic organic compounds by activated carbon at varying adsorbate concentrations. The models are validated using state-of-the-art procedures employing an external prediction set of compounds. This work reveals the hydrogen bond donating and accepting ability of compounds to be the most influencing - but a negative - factor in the adsorption process of activated carbon. The quantum-mechanical LSERs proposed in this work are analysed and found to be equally reliable as the existing LSERs. These were further used to predict the adsorption of nucleobases, steroid hormones, agrochemicals, endocrine disruptors and pharmaceutical drugs. Notably, agrochemicals and endocrine disruptors are predicted to be adsorbed more strongly by activated carbon when compared with their adsorption by CNTs. However, quantum-mechanical LSERs predict the adsorption strength of biomolecules on activated carbon to be similar to that on the CNTs, which can be used to assess the risk associated with using carbon materials.

Concentration dependent adsorption of aromatic organic compounds by SWCNTs: Quantum-mechanical descriptors for nano-toxicological studies of biomolecules and agrochemicals

For evaluating the environmental risk associated using carbon nanotubes (CNTs), a successful prediction is desired for the adsorption of organic compounds (OCs) by CNTs at different adsorbate concentrations. This is most often achieved through poly-parameter linear solvation energy relationships (LSERs) based on solvatochromic descriptors. This study examines the real predictivity of the existing LSERs for predicting the adsorption of OCs by single-walled CNTs (SWCNTs) while comparing it with that of the models developed in the present work using quantum-mechanical descriptors. The real predictivity of the quantum-mechanical models and existing LSERs is compared using state-of-the-art statistical procedures employing an external prediction set of compounds not used in the model development. The quantum-mechanically computed mean polarizability, but originating from the interactions between electrons of parallel spin, is found to play an essential role in the adsorption of OCs by SWCNTs. Besides the solvatochromic descriptors (McGowan volume and molar excess refractivity), the instantaneous inter-electronic interactions, captured through electron-correlation based quantum-mechanical descriptors, are found to significantly affect the adsorption at varying adsorbate concentration. The models developed using a combination of quantum-mechanical and solvatochromic descriptors are found to be quite reliable. The models proposed were further employed to predict the adsorption of agrochemicals such as insecticides, pesticides, herbicides, as well as adsorption of endocrine disruptors and biomolecules such as nucleobases and steroid hormones. These are predicted to be strongly adsorbed by SWCNTs with Progesterone and Guanine exhibiting maximal interaction with the SWCNTs among biomolecules. The quantum-mechanical descriptors proposed in this work can be used for the risk assessment of SWCNTs in systems where adsorption is the primary process.

Quantum-mechanical parameters for the risk assessment of multi-walled carbon-nanotubes: A study using adsorption of probe compounds and its application to biomolecules

This work forwards new insights into the risk-assessment of multi-walled carbon-nanotubes (MWCNTs) while analysing the role of quantum-mechanical interactions between the electrons in the adsorption of probe compounds and biomolecules by MWCNTs. For this, the quantitative models are developed using quantum-chemical descriptors and their electron-correlation contribution. The major quantum-chemical factors contributing to the adsorption are found to be mean polarizability, electron-correlation energy, and electron-correlation contribution to the absolute electronegativity and LUMO energy. The proposed models, based on only three quantum-chemical factors, are found to be even more robust and predictive than the previously known five or four factors based linear free-energy and solvation-energy relationships. The proposed models are employed to predict the adsorption of biomolecules including steroid hormones and DNA bases. The steroid hormones are predicted to be strongly adsorbed by the MWCNTs, with the order: hydrocortisone > aldosterone > progesterone > ethinyl-oestradiol > testosterone > oestradiol, whereas the DNA bases are found to be relatively less adsorbed but follow the order as: guanine > adenine > thymine > cytosine > uracil. Besides these, the developed electron-correlation based models predict several insecticides, pesticides, herbicides, fungicides, plasticizers and antimicrobial agents in cosmetics, to be strongly adsorbed by the carbon-nanotubes. The present study proposes that the instantaneous inter-electronic interactions may be quite significant in various physico-chemical processes involving MWCNTs, and can be used as a reliable predictor for their risk assessment.

Quantum-chemically-calculated mechanistically interpretable molecular descriptors for drug-action mechanism study – a case study of anthracycline anticancer antibiotics

Identification of drug-action mechanisms and understanding of chemical substituents affecting the anticancer activity of drugs are important for drug development.

Externally predictive single-descriptor based QSPRs for physico-chemical properties of polychlorinated-naphthalenes: Exploring relationships of logS(W), logK(OA), and logK(OW) with electron-correlation

Quantitative structure-property relationships (QSPRs), based only on a single-parameter, are proposed for the prediction of physico-chemical properties, namely, aqueous solubility (logSW), octanol-water partition coefficient (logKOW) and octanol-air partition coefficient (logKOA) of polychloronaphthalenes (PCNs) including all the 75 chloronaphthalene congeners. The QSPR models are developed using molecular descriptors computed through quantum mechanical methods including ab-initio as well as advanced semi-empirical methods. The predictivity of the developed models is tested through state-of-the-art external validation procedures employing an external prediction set of compounds. To analyse the role of instantaneous interactions between electrons (the electron-correlation), the models are also compared with those developed using only the electron-correlation contribution of the quantum chemical descriptor. The electron-correlation contribution towards the chemical hardness and the LUMO energy are observed to be the best predictors for octanol-water partition coefficient, whereas for the octanol-air partition coefficient, the total electronic energy and electron-correlation energy are found to be reliable descriptors, in fact, even better than the polarisability. For aqueous solubility of PCNs, the absolute electronegativity is observed to be the best predictor. This work suggests that the electron-correlation contribution of a quantum-chemical descriptor can be used as a reliable indicator for physico-chemical properties, particularly the partition coefficients.

Exploring the role of quantum chemical descriptors in modeling acute toxicity of diverse chemicals to Daphnia magna

Various quantum-mechanically computed molecular and thermodynamic descriptors along with physico-chemical, electrostatic and topological descriptors are compared while developing quantitative structure-activity relationships (QSARs) for the acute toxicity of 252 diverse organic chemicals towards Daphnia magna. QSAR models based on the quantum-chemical descriptors, computed with routinely employed advanced semi-empirical and ab-initio methods, along with the electron-correlation contribution (CORR) of the descriptors, are analyzed for the external predictivity of the acute toxicity. The models with reliable internal stability and external predictivity are found to be based on the HOMO energy along with the physico-chemical, electrostatic and topological descriptors. Besides this, the total energy and electron-correlation energy are also observed as highly reliable descriptors, suggesting that the intra-molecular interactions between the electrons play an important role in the origin of the acute toxicity, which is in fact an unexplored phenomenon. The models based on quantum-chemical descriptors such as chemical hardness, absolute electronegativity, standard Gibbs free energy and enthalpy are also observed to be reliable. A comparison of the robust models based on the quantum-chemical descriptors computed with various quantum-mechanical methods suggests that the advanced semi-empirical methods such as PM7 can be more reliable than the ab-initio methods which are computationally more expensive.

Single-descriptor based quantum-chemical QSPRs for physico-chemical properties of polychlorinated-dibenzo-p-dioxins and -dibenzo-furans (PCDD/Fs): exploring the role of electron-correlation

In this work, we propose and analyze single-descriptor based quantitative structure-property relationships (QSPRs) developed using quantum-chemical descriptors and their electron-correlation contribution, for various physico-chemical properties, namely, aqueous solubility, subcooled liquid vapour pressure, n-octanol/water and n-octanol/air partition coefficients of polychlorinated-dibenzo-p-dioxins (PCDDs) and -dibenzo-furans (PCDFs). The predictivity of the developed QSPRs is examined through external validation procedures employing an external prediction set of chemicals not used in the QSPR model development. From the comparison of different models, it is observed that the total energy and the mean polarizability of a molecule are highly significant in determining the predictivity of a QSPR, whereas dipole moment is found to be a poor descriptor in the models developed for all the properties analyzed. Besides these, the electron density based descriptors, namely, absolute electronegativity and electrophilicity index, were also observed to be highly influential in determining the external prediction for partition coefficients.

Electron-correlation based externally predictive QSARs for mutagenicity of nitrated-PAHs in Salmonella typhimurium TA100

In quantitative modeling, there are two major aspects that decide reliability and real external predictivity of a structure-activity relationship (SAR) based on quantum chemical descriptors. First, the information encoded in employed molecular descriptors, computed through a quantum-mechanical method, should be precisely estimated. The accuracy of the quantum-mechanical method, however, is dependent upon the amount of electron-correlation it incorporates. Second, the real external predictivity of a developed quantitative SAR (QSAR) should be validated employing an external prediction set. In this work, to analyze the role of electron-correlation, QSAR models are developed for a set of 51 ubiquitous pollutants, namely, nitrated monocyclic and polycyclic aromatic hydrocarbons (nitrated-AHs and PAHs) having mutagenic activity in TA100 strain of Salmonella typhimurium. The quality of the models, through state-of-the-art external validation procedures employing an external prediction set, is compared to the best models known in the literature for mutagenicity. The molecular descriptors whose electron-correlation contribution is analyzed include total energy, energy of HOMO and LUMO, and commonly employed electron-density based descriptors such as chemical hardness, chemical softness, absolute electronegativity and electrophilicity index. The electron-correlation based QSARs are also compared with those developed using quantum-mechanical descriptors computed with advanced semi-empirical (SE) methods such as PM6, PM7, RM1, and ab initio methods, namely, the Hartree-Fock (HF) and the density functional theory (DFT). The models, developed using electron-correlation contribution of the quantum-mechanical descriptors, are found to be not only reliable but also satisfactorily predictive when compared to the existing robust models. The robustness of the models based on descriptors computed through advanced SE methods, is also observed to be comparable to those developed with the electron-correlation based descriptors. The work emphasizes that the correlation-energy can serve as a reliable descriptor to explore the origin of biological activities at the level of electron-dynamics.

Externally predictive quantitative modeling of supercooled liquid vapor pressure of polychlorinated-naphthalenes through electron-correlation based quantum-mechanical descriptors

For predicting physico-chemical properties related to environmental fate of molecules, quantitative structure-property relationships (QSPRs) are valuable tools in environmental chemistry. For developing a QSPR, molecular descriptors computed through quantum-mechanical methods are generally employed. The accuracy of a quantum-mechanical method, however, rests on the amount of electron-correlation estimated by the method. In this work, single-descriptor QSPRs for supercooled liquid vapor pressure of chloronaphthalenes and polychlorinated-naphthalenes are developed using molecular descriptors based on the electron-correlation contribution of the quantum-mechanical descriptor. The quantum-mechanical descriptors for which the electron-correlation contribution is analyzed include total-energy, mean polarizability, dipole moment, frontier orbital (HOMO/LUMO) energy, and density-functional theory (DFT) based descriptors, namely, absolute electronegativity, chemical hardness, and electrophilicity index. A total of 40 single-descriptor QSPRs were developed using molecular descriptors computed with advanced semi-empirical (SE) methods, namely, RM1, PM7, and ab intio methods, namely, Hartree-Fock and DFT. The developed QSPRs are validated using state-of-the-art external validation procedures employing an external prediction set. From the comparison of external predictivity of the models, it is observed that the single-descriptor QSPRs developed using total energy and correlation energy are found to be far more robust and predictive than those developed using commonly employed descriptors such as HOMO/LUMO energy and dipole moment. The work proposes that if real external predictivity of a QSPR model is desired to be explored, particularly, in terms of intra-molecular interactions, correlation-energy serves as a more appropriate descriptor than the polarizability. However, for developing QSPRs, computationally inexpensive advanced SE methods such as PM7 can be more reliable than the expensive ab inito methods.