Conceptual DFT properties-based 3D QSAR: Analysis of inhibitors of the nicotine metabolizing CYP2A6 enzyme

Atom and Bond Fukui Functions and Matrices: A Hirshfeld-I Atoms-in-Molecule Approach

The Fukui function is often used in its atom-condensed form by isolating it from the molecular Fukui function using a chosen weight function for the atom in the molecule. Recently, Fukui functions and matrices for both atoms and bonds separately were introduced for semiempirical and ab initio levels of theory using Hückel and Mulliken atoms-in-molecule models. In this work, a double partitioning method of the Fukui matrix is proposed within the Hirshfeld-I atoms-in-molecule framework. Diagonalizing the resulting atomic and bond matrices gives eigenvalues and eigenvectors (Fukui orbitals) describing the reactivity of atoms and bonds. The Fukui function is the diagonal element of the Fukui matrix and may be resolved in atom and bond contributions. The extra information contained in the atom and bond resolution of the Fukui matrices and functions is highlighted. The effect of the choice of weight function arising from the Hirshfeld-I approach to obtain atom- and bond-condensed Fukui functions is studied. A comparison of the results with those generated by using the Mulliken atoms-in-molecule approach shows low correlation between the two partitioning schemes.

Analyzing the substitution effect on the CoMFA results within the framework of density functional theory (DFT)

Though QSAR was originally developed in the context of physical organic chemistry, it has been applied very extensively to chemicals (drugs) which act on biological systems, in this idea one of the most important QSAR methods is the 3D QSAR model. However, due to the complexity of understanding the results it is necessary to postulate new methodologies to highlight their physical-chemical meaning. In this sense, this work postulates new insights to understand the CoMFA results using molecular quantum similarity and chemical reactivity descriptors within the framework of density functional theory. To obtain these insights a simple theoretical scheme involving quantum similarity (overlap, coulomb operators, their euclidean distances) and chemical reactivity descriptors such as chemical potential (μ), hardness (ɳ), softness (S), electrophilicity (ω), and the Fukui functions, was used to understand the substitution effect. In this sense, this methodology can be applied to analyze the biological activity and the stabilization process in the non-covalent interactions on a particular molecular set taking a reference compound.

Insights on cytochrome p450 enzymes and inhibitors obtained through QSAR studies

The cytochrome P450 (CYP) superfamily of heme enzymes play an important role in the metabolism of a large number of endogenous and exogenous compounds, including most of the drugs currently on the market. Inhibitors of CYP enzymes have important roles in the treatment of several disease conditions such as numerous cancers and fungal infections in addition to their critical role in drug-drug interactions. Structure activity relationships (SAR), and three-dimensional quantitative structure activity relationships (3D-QSAR) represent important tools in understanding the interactions of the inhibitors with the active sites of the CYP enzymes. A comprehensive account of the QSAR studies on the major human CYPs 1A1, 1A2, 1B1, 2A6, 2B6, 2C9, 2C19, 2D6, 2E1, 3A4 and a few other CYPs are detailed in this review which will provide us with an insight into the individual/common characteristics of the active sites of these enzymes and the enzyme-inhibitor interactions.

Identification of novel CYP2A6 inhibitors by virtual screening

The human CYP2A6 enzyme metabolises several xenobiotics including nicotine, the addictive component in tobacco. Reduced activity of CYP2A6, either for genetic reasons or by administering inhibitors of CYP2A6, reduces tobacco smoking. The aim was to design novel inhibitors of CYP2A6 using 3D-QSAR analysis combined with virtual screening. A 3D-QSAR model was utilised to identify the most important features of the inhibitors, and this knowledge was used to design inhibitors for CYP2A6. Chemical database screening yielded several potent inhibitor candidates such as alkylamine derivatives (compound no. 5, IC(50)=0.1 μM) and 1-benzothiophene-3-carbaldehyde that can be used as lead compounds in the development of drugs for smoking reduction therapy.

Integrating molecular docking, DFT and CoMFA/CoMSIA approaches for a series of naphthoquinone fused cyclic α-aminophosphonates that act as novel topoisomerase II inhibitors

Since they are potential topoisomerase II (Topo II) inhibitors, naphthoquinone fused cyclic α-aminophosphonates display anticancer activity. In order to explore the inhibitory mechanisms of these compounds, they were docked into the active site of Topo II structure, which allowed their probable binding modes to be predicted. Some meaningful results concerning their structure-activity relationships were obtained from density functional theory calculations. Models based on quantitative comparative molecular field analysis and comparative molecular similarity index analysis were derived for the steric, electrostatic, hydrophobic and H-bonding features of the compounds. The present study provides valuable results that enhance our understanding of the anticancer activities of these inhibitors and will aid the rational drug design of novel Topo II inhibitors in the future.

Quantitative Prediction of Regioselectivity Toward Cytochrome P450/3A4 Using Machine Learning Approaches

In the drug discovery process, it is important to know the properties of both drug candidates and their metabolites. Fast and precise prediction of metabolites is essential. However, it has been difficult to predict metabolites because of the complexity of the mechanism of cytochrome P450/3A4 (CYP 3A4), which is the main metabolite enzyme of drugs. In this study, we focus on the regioselectivity of CYP 3A4, i.e., the selectivity of metabolic sites. We have developed a model to predict the regioselectivity of drug candidates by using machine learning (ML) approaches.

QSAR of cytochrome inhibitors

Cytochrome P450 (CYP450) enzymes are predominantly involved in the Phase I metabolism of xenobiotics. Metabolic inhibition and induction can give rise to clinically important drug-drug interactions. Metabolic stability is a prerequisite for sustaining the therapeutically relevant concentrations, and very often drug candidates are sacrificed due to poor metabolic profiles. Computational tools such as quantitative structure-activity relationships are widely used to study different metabolic end points successfully to accelerate the drug discovery process. There are a lot of computational studies on clinically important CYPs already reported in recent years. But other clinically significant families are to yet be explored computationally. Powerfulness of quantitative structure-activity relationship will drive computational chemists to develop new potent and selective inhibitors of different classes of CYPs for the treatment of different diseases with least drug-drug interactions. Furthermore, there is a need to enhance the accuracy, interpretability and confidence in the computational models in accelerating the drug discovery pathways.