QSAR and 3D QSAR in drug design Part 1: methodology

Pharmacophore mapping, 3D QSAR, molecular docking, and ADME prediction studies of novel Benzothiazinone derivatives

In the quest to combat tuberculosis, DprE1, a challenging target for novel anti-tubercular agents due to its small size and membrane location, has been a focus of research. DprE1 catalyzes the transformation of DPR into Ketoribose DPX, with Benzothiazinone emerging as a potent pharmacophore for inhibiting DprE1. Clinical trial drugs such as BTZ043, BTZ038, PBTZ169, and TMC-207 have shown promising results as DprE1 inhibitors. This study employed pharmacophore mapping of Pyrazolopyridine, Dinitrobenzamide, and Benzothiazinone derivatives to identify crucial features for eliciting a biological response. Benzothiazinone (Ligand code: 73) emerged as a reference ligand with a fitness score of 3.000. ROC analysis validated the pharmacophore with an excellent score of 0.71. To build a 3D QSAR model, a series of Benzothiazinone congeneric derivatives were explored. The model exhibited strong performance, with a standard deviation of 0.1531, a correlation coefficient for the training set (R2) value of 0.9754, and a correlation coefficient for test set Q2 value of 0.7632, indicating robust predictive capabilities. Contour maps guided the design of novel benzothiazinone derivatives, emphasizing steric, electrostatic, hydrophobic, H-bond acceptor, and H-bond donor groups for structure-activity relationships. Docking studies against PDB ID: 4NCR demonstrated favorable scores, with interactions aligning well with the in-built ligand 26 J. Docking validation via RMSD values supported the reliability of the docking results. This comprehensive approach aids in the design of novel benzothiazinone derivatives with potential anti-tubercular properties, contributing to the development of novel anti-tubercular agents which can be pivotal in the eradication of tuberculosis.

Prediction method of pharmacokinetic parameters of small molecule drugs based on GCN network model

CONTEXT

Accurately predicting plasma protein binding rate (PPBR) and oral bioavailability (OBA) helps to better reveal the absorption and distribution of drugs in the human body and subsequent drug design. Although machine learning models have achieved good results in prediction accuracy, they often suffer from insufficient accuracy when dealing with data with irregular topological structures.

METHODS

In view of this, this study proposes a pharmacokinetic parameter prediction framework based on graph convolutional networks (GCN), which predicts the PPBR and OBA of small molecule drugs. In the framework, GCN is first used to extract spatial feature information on the topological structure of drug molecules, in order to better learn node features and association information between nodes. Then, based on the principle of drug similarity, this study calculates the similarity between small molecule drugs, selects different thresholds to construct datasets, and establishes a prediction model centered on the GCN algorithm. The experimental results show that compared with traditional machine learning prediction models, the prediction model constructed based on the GCN method performs best on PPBR and OBA datasets with an inter-molecular similarity threshold of 0.25, with MAE of 0.155 and 0.167, respectively. In addition, in order to further improve the accuracy of the prediction model, GCN is combined with other algorithms. Compared to using a single GCN method, the distribution of the predicted values obtained by the combined model is highly consistent with the true values. In summary, this work provides a new method for improving the rate of early drug screening in the future.

Leishmanicidal activity and 4D quantitative structure-activity relationship and molecular docking studies of vanillin-containing 1,2,3-triazole derivatives

Aim: The assessment of the antileishmanial potential of 22 vanillin-containing 1,2,3-triazole derivatives against Leishmania braziliensis is reported. Materials & methods: Initial screening was performed against the parasite promastigote form. The most active compound, 4b, targeted parasites within amastigotes (IC50 = 4.2 ± 1.0 μmol l-1), presenting low cytotoxicity and a selective index value of 39. 4D quantitative structure-activity relationship and molecular docking studies provided insights into structure-activity and biological effects. Conclusion: A vanillin derivative with significant antileishmanial activity was identified. Enhanced activity was linked to increased electrostatic and Van der Waals interactions near the benzyl ring of the derivatives. Molecular docking indicated the inhibition of the Leishmania amazonensis sterol 14α-demethylase, using Leishmania infantum sterol 14α-demethylase as a model, without affecting the human isoform. Inhibition was active site competition with lanosterol.

Three Decades of Targeting Falcipains to Develop Antiplasmodial Agents: What have we Learned and What can be Done Next?

Malaria is a devastating infectious disease that affects large swathes of human populations across the planet's tropical regions. It is caused by parasites of the genus Plasmodium, with Plasmodium falciparum being responsible for the most lethal form of the disease. During the intraerythrocytic stage in the human hosts, malaria parasites multiply and degrade hemoglobin (Hb) using a battery of proteases, which include two cysteine proteases, falcipains 2 and 3 (FP-2 and FP-3). Due to their role as major hemoglobinases, FP-2 and FP-3 have been targeted in studies aiming to discover new antimalarials and numerous inhibitors with activity against these enzymes, and parasites in culture have been identified. Nonetheless, cross-inhibition of human cysteine cathepsins remains a serious hurdle to overcome for these compounds to be used clinically. In this article, we have reviewed key functional and structural properties of FP-2/3 and described different compound series reported as inhibitors of these proteases during decades of active research in the field. Special attention is also paid to the wide range of computer-aided drug design (CADD) techniques successfully applied to discover new active compounds. Finally, we provide guidelines that, in our understanding, will help advance the rational discovery of new FP-2/3 inhibitors.

Pharmacophore-based virtual screening, 3D QSAR, Docking, ADMET, and MD simulation studies: An in silico perspective for the identification of new potential HDAC3 inhibitors

Histone deacetylase 3 (HDAC3) is an epigenetic regulator that involves gene expression, apoptosis, and cell cycle progression, and the overexpression of HDAC3 is accountable for several cancers, neurodegeneracy, and many other diseases. Therefore, HDAC3 emerged as a promising drug target for the novel drug design. Here, we carried out the pharmacophore modeling using 50 benzamide-based HDAC3 selective inhibitors and utilized it for PHASE ligand screening to retrieve the hits with similar pharmacophore features. The dataset inhibitors of best hypotheses used to build the 3D QSAR model and the generated 3D QSAR model resulted in good PLS statistics with a regression coefficient (R2) of 0.89, predictive coefficient (Q2) of 0.88, and Pearson-R factor of 0.94 indicating its excellent predictive ability. The hits retrieved from pharmacophore-based virtual screening were subjected to docking against HDAC3 for the identification of potential inhibitors. A total of 10 hitsM1 to M10 were ranked using their scoring functions and further subject to lead optimization. The Prime MM/GBSA, AutoDock binding free energies, and ADMET studies were implemented for the selection of lead candidates. The four ligand molecules M1, M2, M3, and M4 were identified as potential leads against HDAC3 after lead optimization. The top two leads M1 and M2 were subjected to MD simulations for their stability evaluation with HDAC3. The newly designed leads M11 and M12 were identified as HDAC3 potential inhibitors from MD simulations studies. Therefore, the outcomes of the present study could provide insights into the discovery of new potential HDAC3 inhibitors with improved selectivity and activity against a variety of cancers and neurodegenerative diseases.

Steroids Bearing Heteroatom as Potential Drugs for Medicine

Heteroatom steroids, a diverse class of organic compounds, have attracted significant attention in the field of medicinal chemistry and drug discovery. The biological profiles of heteroatom steroids are of considerable interest to chemists, biologists, pharmacologists, and the pharmaceutical industry. These compounds have shown promise as potential therapeutic agents in the treatment of various diseases, such as cancer, infectious diseases, cardiovascular disorders, and neurodegenerative conditions. Moreover, the incorporation of heteroatoms has led to the development of targeted drug delivery systems, prodrugs, and other innovative pharmaceutical approaches. Heteroatom steroids represent a fascinating area of research, bridging the fields of organic chemistry, medicinal chemistry, and pharmacology. The exploration of their chemical diversity and biological activities holds promise for the discovery of novel drug candidates and the development of more effective and targeted treatments.

Using the local symmetry in amino acids sequences of polypeptides to improve the predictive potential of models of their inhibitor activity

The minimal inhibitory concentrations (pMIC) are a valuable measure of the biological activity of polypeptides. Numerical data on the pMIC are necessary to systematize knowledge on polypeptides' biochemical behaviour. The model of negative decimal logarithm of pMIC of polypeptides in the form of a mathematical function of a sequence of amino acids is suggested. The suggested model is based on the so-called correlation weights of amino acids together with the correlation weights of fragments of local symmetry (FLS). Three kinds of the FLS are considered: (i) three-symbol fragments '…xyx…', (ii) four-symbol fragments '…xyyx…', and (iii) five-symbol fragments '…xyzyx…'. The models built using the Monte Carlo technique improved by applying the index of ideality of correlation (IIC) and correlation intensity index (CII).

Extended study on atomic featurization in graph neural networks for molecular property prediction

Graph neural networks have recently become a standard method for analyzing chemical compounds. In the field of molecular property prediction, the emphasis is now on designing new model architectures, and the importance of atom featurization is oftentimes belittled. When contrasting two graph neural networks, the use of different representations possibly leads to incorrect attribution of the results solely to the network architecture. To better understand this issue, we compare multiple atom representations by evaluating them on the prediction of free energy, solubility, and metabolic stability using graph convolutional networks. We discover that the choice of atom representation has a significant impact on model performance and that the optimal subset of features is task-specific. Additional experiments involving more sophisticated architectures, including graph transformers, support these findings. Moreover, we demonstrate that some commonly used atom features, such as the number of neighbors or the number of hydrogens, can be easily predicted using only information about bonds and atom type, yet their explicit inclusion in the representation has a positive impact on model performance. Finally, we explain the predictions of the best-performing models to better understand how they utilize the available atomic features.

Discovery of ANO1 Inhibitors based on Machine learning and molecule docking simulation approaches

Calcium-activated chloride channels (CaCCs) are chloride channels that are regulated according to intracellular calcium ion concentrations. The channel protein ANO1 is widely present in cells and is involved in physiological activities including cellular secretion, signaling, cell proliferation and vasoconstriction and diastole. In this study, the ANO1 inhibitors were investigated with machine learning and molecular simulation. Two-dimensional structure-activity relationship (2D-SAR) and three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed for the qualitative and quantitative prediction of ANO1 inhibitors. The results showed that the prediction accuracies of the model were 85.9% and 87.8% for the training and test sets, respectively, and 85.9% and 87.8% for the rotating forest (RF) in the 2D-SAR model. The CoMFA and CoMSIA methods were then used for 3D QSAR modeling of ANO1 inhibitors, respectively. The q² coefficients for model cross-validation were all greater than 0.5, implying that we were able to obtain a stable model for drug activity prediction. Molecular docking was further used to simulate the interactions between the five most promising compounds predicted by the model and the ANO1 protein. The total score for the docking results between all five compounds and the target protein was greater than 6, indicating that they interacted strongly in the form of hydrogen bonds. Finally, simulations of amino acid mutations around the docking cavity of the target proteins showed that each molecule had two or more sites of reduced affinity following a single mutation, indicating outstanding specificity of the screened drug molecules and their protein ligands.

Synthesis of fluorescent chalcones, photophysical properties, quantitative structure-activity relationship and their biological application

The development of fluorescent pigments is an area of interest in several research fields due to their high sensitivity. In the current study-eight known and three new N,N-dimethylamino-chalcones (12a-k) were synthesized with good yields using the Claisen-Schmidt reaction. For each molecular system, the photophysical properties, including the maximum absorption wavelength (λ_Absorption), molar absorption coefficient (ε), maximum excitation wavelength (λ_Excitation), maximum emission wavelength (λ_Emission), Stokes Shift (Δλ), fluorescence quantum yield (Φ_fl), fluorescence lifetime (τ_fl), radiative and non-radiative rate constants (k_R and k_NR, respectively) were evaluated. Variations in each of these properties were analyzed depending on the substituents present on each compound. To relate the chemical structures of the synthesized compounds to their photophysical properties, Hansch analysis (2D-QSPR) was applied. As a result of Hansch analysis, we found different photophysical properties related to molecular orbitals and the energy of their derivatives (Highest Occupied Molecular Orbital-HOMO, Lowest Unoccupied Molecular Orbital-LUMO, Difference between LUMO-HOMO-ΔLH, Chemical potential-µ, Hardness-η, Softness-S, and electrophilic global index-ω) as well as to the atomic charges on atoms C₅, C_α, C_β, and CO. The application of this type of analysis has made it possible to understand and subsequently design new molecules with defined photophysical properties. Finally, the compounds were use as fluorescent pigment to get living cell imaging on breast cancer cells, obtaining the compound 12a as promissory alternative.

Synthesis of vanillin derivatives with 1,2,3-triazole fragments and evaluation of their fungicide and fungistatic activities

Vanillin is the main component of natural vanilla extract and is responsible for its flavoring properties. Besides its well-known applications as an additive in food and cosmetics, it has also been reported that vanillin can inhibit fungi of clinical interest, such as Candida spp., Cryptococcus spp., Aspergillus spp., as well as dermatophytes. Thus, the present work approaches the synthesis of a series of vanillin derivatives with 1,2,3-triazole fragments and the evaluation of their antifungal activities against Candida albicans, Candida glabrata, Candida parapsilosis, Candida tropicalis, Cryptococcus neoformans, Cryptococcus gattii, Trichophyton rubrum, and Trichophyton interdigitale strains. Twenty-two vanillin derivatives were obtained, with yields in the range of 60%-91%, from copper(I)-catalyzed alkyne-azide cycloaddition (CuAAC) click reaction between two terminal alkynes prepared from vanillin and different benzyl azides. In general, the evaluated compounds showed moderate activity against the microorganisms tested, with minimum inhibitory concentration (MIC) values ranging from 32 to >512 µg mL^-1 . Except for compound 3b against the C. gattii R265 strain, all vanillin derivatives showed fungicidal activity for the yeasts tested. The predicted physicochemical and ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties for the compounds indicated favorable profiles for drug development. In addition, a four-dimensional structure-activity relationship (4D-SAR) analysis was carried out and provided useful insights concerning the structures of the compounds and their biological profile. Finally, molecular docking calculations showed that all compounds bind favorably at the lanosterol 14α-demethylase enzyme active site with binding energies ranging from -9.1 to -12.2 kcal/mol.

Density Functional Theory and Machine Learning-Based Quantitative Structure-Activity Relationship Models Enabling Prediction of Contaminant Degradation Performance with Heterogeneous Peroxymonosulfate Treatments

Heterogeneous peroxymonosulfate (PMS) treatment is recognized as an effective advanced oxidation process (AOP) for the treatment of organic contaminants. Quantitative structure-activity relationship (QSAR) models have been applied to predict the oxidation reaction rates of contaminants in homogeneous PMS treatment systems but are seldom applied in heterogeneous treatment systems. Herein, we established QSAR models updated with density functional theory (DFT) and machine learning approaches to predict the degradation performance for a series of contaminants in heterogeneous PMS systems. We imported the characteristics of organic molecules calculated using constrained DFT as input descriptors and predicted the apparent degradation rate constants of contaminants as the output. The genetic algorithm and deep neural networks were used to improve the predictive accuracy. The qualitative and quantitative results from the QSAR model for the degradation of contaminants can be used to select the most appropriate treatment system. A strategy for selection of the optimum catalyst for PMS treatment of specific contaminants was also established according to the QSAR models. This work not only increases our understanding of contaminant degradation in PMS treatment systems but also highlights a novel QSAR model to predict the degradation performance in complicated heterogeneous AOPs.

Hansch-Type QSAR Models for the Rational Design of MAO Inhibitors: Basic Principles and Methodology

Hansch-type regression analysis enables the derivation of quantitative structure-activity relationship (QSAR) equations correlating bioactivity data with physicochemical parameters accounting for hydrophobicity, electronic properties, and steric effects of molecules or functional groups (substituents). Two datasets of MAO A and B inhibitors were enrolled in prototypical workflows employing multiparametric stepwise regression analysis, which includes linear and nonlinear (generally quadratic) terms. The optimal choice of variables (and/or combinations thereof) along with statistical validation yielded two robust equations describing MAO B potency and B/A selectivity, which included three and one parameter(s), respectively, and explained more than 80% of y-variance (r²) with low standard deviation (s) and good statistical significance (F, Fisher value).

An Optimized Cloud Computing Method for Extracting Molecular Descriptors

Extracting molecular descriptors from chemical compounds is an essential preprocessing phase for developing accurate classification models. Supervised machine learning algorithms offer the capability to detect "hidden" patterns that may exist in a large dataset of compounds, which are represented by their molecular descriptors. Assuming that molecules with similar structure tend to share similar physicochemical properties, large chemical libraries can be screened by applying similarity sourcing techniques in order to detect potential bioactive compounds against a molecular target. However, the process of generating these compound features is time-consuming. Our proposed methodology not only employs cloud computing to accelerate the process of extracting molecular descriptors but also introduces an optimized approach to utilize the computational resources in the most efficient way.

QSPR Modeling and Experimental Determination of the Antioxidant Activity of Some Polycyclic Compounds in the Radical-Chain Oxidation Reaction of Organic Substrates

The present work addresses the quantitative structure−antioxidant activity relationship in a series of 148 sulfur-containing alkylphenols, natural phenols, chromane, betulonic and betulinic acids, and 20-hydroxyecdysone using GUSAR2019 software. Statistically significant valid models were constructed to predict the parameter logk7, where k7 is the rate constant for the oxidation chain termination by the antioxidant molecule. These results can be used to search for new potentially effective antioxidants in virtual libraries and databases and adequately predict logk7 for test samples. A combination of MNA- and QNA-descriptors with three whole molecule descriptors (topological length, topological volume, and lipophilicity) was used to develop six statistically significant valid consensus QSPR models, which have a satisfactory accuracy in predicting logk7 for training and test set structures: R2TR > 0.6; Q2TR > 0.5; R2TS > 0.5. Our theoretical prediction of logk7 for antioxidants AO1 and AO2, based on consensus models agrees well with the experimental value of the measure in this paper. Thus, the descriptor calculation algorithms implemented in the GUSAR2019 software allowed us to model the kinetic parameters of the reactions underlying the liquid-phase oxidation of organic hydrocarbons.

Fundamental considerations in drug design

The drug discovery paradigm has been very time-consuming, challenging, and expensive; however, the disease conditions originating from bacteria, virus, protozoa, fungus and other microorganisms are steadily shooting up. For instance, COVID-19 is the latest viral infection that affects millions of people and the world’s economy very severely. Therefore, the quest for discovery of novel and potent drug compounds against deadly pathogens is crucial at the moment. Despite a lot of drawbacks in drug discovery and development and its pertaining technology, the advancement must be taken into account so the time duration and cost would be minimized. In this chapter, basic principles in drug design and discovery have been discussed together with advances in drug development.

Application of quasi-SMILES to the model of gold-nanoparticles uptake in A549 cells

Cell death is critical to human health and is associated with a variety of medical conditions. Therefore, new controllers of cell death are needed for the treatment of diverse diseases. In particular, nanoparticles (NP) are now regularly used in various applications, including a variety of products and medicines. Gold nanoparticles (GNPs) are widely used in the medical field against A549 lung carcinoma cells. The present study is devoted to developing computational models of the cellular uptake potentials by A549 cells of gold nanoparticles (GNPs) under various conditions. Simplified molecular input-line entry system (SMILES) is an efficient tool to represent the molecular structure by a sequence of symbols. Quasi-SMILES represents an extended version of SMILES where symbols to denote physicochemical and/or biochemical conditions are added. In other words, the quasi-SMILES represents a biochemical (medical) phenomenon related to the whole matter (not only molecular structure). We developed models for the cellular utpake potential of gold nanoparticles (GNPs) in A549 [10-11 g Au/Cell] under various conditions based on quasi-SMILES using the Monte Carlo method. The statistical quality of these models is quite good.

Leishmanicidal and cytotoxic activities and 4D-QSAR of 2-arylidene indan-1,3-diones

The indan-1,3-dione and its derivatives are important building blocks in organic synthesis and present important biological activities. Herein, the leishmanicidal and cytotoxicity evaluation of 16 2-arylidene indan-1,3-diones is described. The compounds were evaluated against the leukemia cell lines HL60 and Nalm6, and the most effective ones were 2-(4-nitrobenzylidene)-1H-indene-1,3(2H)-dione (4) and 4-[(1,3-dioxo-1H-inden-2(3H)-ylidene)methyl]benzonitrile (10), presenting IC₅₀ values of around 30 µmol/L against Nalm6. The leishmanicidal activity was assessed on Leishmania amazonensis, with derivative 4 (IC₅₀ = 16.6 µmol/L) being the most active. A four-dimensional quantitative structure-activity analysis (4D-QSAR) was applied to the indandione derivatives, through partial least-squares regression. The statistics presented by the regression models built with the selected field descriptors of Coulomb (C) and Lennard-Jones (L) nature, considering the activities against L. amazonensis, HL60, and Nalm6 leukemia cells, were, respectively, R² = 0.88, 0.92, and 0.98; Q² = 0.83, 0.88, and 0.97. The presence of positive Coulomb descriptors near the carbonyl groups indicates that these polar groups are related to the activities. Besides, the presence of positive Lennard-Jones descriptors close to substituents R³ or R¹ indicates that bulky nonpolar substituents in these positions tend to increase the activities. This study provides useful insights into the mode of action of indandione derivatives for each biological activity involved.

On the use of electronegativity and electron affinity based pseudo-molecular field descriptors in developing correlations for quantitative structure-activity relationship modeling of drug activities

For quantitative structure-activity relationship (QSAR) modeling in ligand-based drug discovery programs, pseudo-molecular field (PMF) descriptors using intrinsic atomic properties, namely, electronegativity and electron affinity are studied. In combination with partial least squares analysis and Procrustes transformation, these PMF descriptors were employed successfully to develop correlations that predict the activities of target protein inhibitors involved in various diseases (cancer, neurodegenerative disorders, HIV, and malaria). The results show that the present QSAR approach is competitive to existing QSAR models. In order to demonstrate the use of this algorithm, we present results of screening naturally occurring molecules with unknown bioactivities. The pIC₅₀ predictions can screen molecules that have desirable activity before assessment by docking studies.

QSAR Modeling, Molecular Docking and Molecular Dynamics Simulations Studies of Lysine-Specific Demethylase 1 (LSD1) Inhibitors as Anticancer Agents

BACKGROUND

Histone Lysine Demetylases1 (LSD1) is a promising medication to treat cancer, which plays a crucial role in epigenetic modulation of gene expression. Inhibition of LSD1with small molecules has emerged as a vital mechanism to treat cancer.

OBJECTIVE

In the present research, molecular modeling investigations, such as CoMFA, CoMFA-RF, CoMSIA and HQSAR, molecular docking and Molecular Dynamics (MD) simulations were carried out on some tranylcypromine derivatives as LSD1 inhibitors.

METHODS

The QSAR models were carried out on a series of Tranylcypromine derivatives as data set via the SYBYL-X2.1.1 program. Molecular docking and MD simulations were carried out by the MOE software and the SYBYL program, respectively. The internal and external predictability performances related to the generated models for these LSD1 inhibitors were justified by evaluating cross-validated correlation coefficient (q²), noncross- validated correlation coefficient (r²_ncv) and predicted correlation coefficient (r²_pred) of the training and test set molecules, respectively.

RESULTS

The CoMFA (q², 0.670; r²_ncv, 0.930; r²_pred, 0.968), CoMFA-RF (q², 0.694; r²_ncr, 0.926; r²_pred, 0.927), CoMSIA (q², 0.834; r²_ncv, 0.956; r²_pred, 0.958) and HQSAR models (q², 0.854; r²_ncv, 0.900; r²_pred, 0.728) for training as well as the test set of LSD1 inhibition resulted in significant findings.

CONCLUSION

These QSAR models were found to be perfect and strong with better predictability. Contour maps of all models were generated and it was proven by molecular docking studies and molecular dynamics simulation that the hydrophobic, electrostatic and hydrogen bonding fields are crucial in these models for improving the binding affinity and determining the structure-activity relationship. These theoretical results are possibly beneficial to design new strong LSD1 inhibitors with enhanced activity to treat cancer.

A machine-learning-assisted study of the permeability of small drug-like molecules across lipid membranes

Study of the permeability of small organic molecules across lipid membranes plays a significant role in designing potential drugs in the field of drug discovery. Approaches to design promising drug molecules have gone through many stages, from experiment-based trail-and-error approaches, to the well-established avenue of the quantitative structure-activity relationship, and currently to the stage guided by machine learning (ML) and artificial intelligence techniques. In this work, we present a study of the permeability of small drug-like molecules across lipid membranes by two types of ML models, namely the least absolute shrinkage and selection operator (LASSO) and deep neural network (DNN) models. Molecular descriptors and fingerprints are used for featurization of organic molecules. Using molecular descriptors, the LASSO model uncovers that the electro-topological, electrostatic, polarizability, and hydrophobicity/hydrophilicity properties are the most important physical properties to determine the membrane permeability of small drug-like molecules. Additionally, with molecular fingerprints, the LASSO model suggests that certain chemical substructures can significantly affect the permeability of organic molecules, which closely connects to the identified main physical properties. Moreover, the DNN model using molecular fingerprints can help develop a more accurate mapping between molecular structures and their membrane permeability than LASSO models. Our results provide deep understanding of drug-membrane interactions and useful guidance for the inverse molecular design of drug-like molecules. Last but not least, while the current focus is on the permeability of drug-like molecules, the methodology of this work is general and can be applied for other complex physical chemistry problems to gain molecular insights.

Machine Learning in Computational Surface Science and Catalysis: Case Studies on Water and Metal-Oxide Interfaces

The goal of many computational physicists and chemists is the ability to bridge the gap between atomistic length scales of about a few multiples of an Ångström (Å), i. e., 10⁻¹⁰ m, and meso- or macroscopic length scales by virtue of simulations. The same applies to timescales. Machine learning techniques appear to bring this goal into reach. This work applies the recently published on-the-fly machine-learned force field techniques using a variant of the Gaussian approximation potentials combined with Bayesian regression and molecular dynamics as efficiently implemented in the Vienna ab initio simulation package, VASP. The generation of these force fields follows active-learning schemes. We apply these force fields to simple oxides such as MgO and more complex reducible oxides such as iron oxide, examine their generalizability, and further increase complexity by studying water adsorption on these metal oxide surfaces. We successfully examined surface properties of pristine and reconstructed MgO and Fe₃O₄ surfaces. However, the accurate description of water–oxide interfaces by machine-learned force fields, especially for iron oxides, remains a field offering plenty of research opportunities.

Antiangiogenic Activity and in Silico Cereblon Binding Analysis of Novel Thalidomide Analogs

Due to its antiangiogenic and anti-immunomodulatory activity, thalidomide continues to be of clinical interest despite its teratogenic actions, and efforts to synthesize safer, clinically active thalidomide analogs are continually underway. In this study, a cohort of 27 chemically diverse thalidomide analogs was evaluated for antiangiogenic activity in an ex vivo rat aorta ring assay. The protein cereblon has been identified as the target for thalidomide, and in silico pharmacophore analysis and molecular docking with a crystal structure of human cereblon were used to investigate the cereblon binding abilities of the thalidomide analogs. The results suggest that not all antiangiogenic thalidomide analogs can bind cereblon, and multiple targets and mechanisms of action may be involved.

Flexi-pharma: a molecule-ranking strategy for virtual screening using pharmacophores from ligand-free conformational ensembles

Computer-aided strategies are useful for reducing the costs and increasing the success-rate in drug discovery. Among these strategies, methods based on pharmacophores (an ensemble of electronic and steric features representing the target active site) are efficient to implement over large compound libraries. However, traditional pharmacophore-based methods require knowledge of active compounds or ligand-receptor structures, and only few ones account for target flexibility. Here, we developed a pharmacophore-based virtual screening protocol, Flexi-pharma, that overcomes these limitations. The protocol uses molecular dynamics (MD) simulations to explore receptor flexibility, and performs a pharmacophore-based virtual screening over a set of MD conformations without requiring prior knowledge about known ligands or ligand-receptor structures for building the pharmacophores. The results from the different receptor conformations are combined using a "voting" approach, where a vote is given to each molecule that matches at least one pharmacophore from each MD conformation. Contrarily to other approaches that reduce the pharmacophore ensemble to some representative models and score according to the matching models or molecule conformers, the Flexi-pharma approach takes directly into account the receptor flexibility by scoring in regards to the receptor conformations. We tested the method over twenty systems, finding an enrichment of the dataset for 19 of them. Flexi-pharma is computationally efficient allowing for the screening of thousands of compounds in minutes on a single CPU core. Moreover, the ranking of molecules by vote is a general strategy that can be applied with any pharmacophore-filtering program.

Teaching and Learning Computational Drug Design: Student Investigations of 3D Quantitative Structure-Activity Relationships through Web Applications

The increasing use of information technology in the discovery of new molecular entities encourages the use of modern molecular-modeling tools to help teach important concepts of drug design to chemistry and pharmacy undergraduate students. In particular, statistical models such as quantitative structure-activity relationships (QSAR)-often as its 3D QSAR variant-are commonly used in the development and optimization of a leading compound. We describe how these drug discovery methods can be taught and learned by means of free and open-source web applications, specifically the online platform www.3d-qsar.com. This new suite of web applications has been integrated into a drug design teaching course, one that provides both theoretical and practical perspectives. We include the teaching protocol by which pharmaceutical biotechnology master students at Pharmacy Faculty of Sapienza Rome University are introduced to drug design. Starting with a choice among recent articles describing the potencies of a series of molecules tested against a biological target, each student is expected to build a 3D QSAR ligand-based model from their chosen publication, proceeding as follows: creating the initial data set (Py-MolEdit); generating the global minimum conformations (Py-ConfSearch); proposing a promising mutual alignment (Py-Align); and finally, building, and optimizing a robust 3D QSAR models (Py-CoMFA). These student activities also help validate these new molecular modeling tools, especially for their usability by inexperienced hands. To more fully demonstrate the effectiveness of this protocol and its tools, we include the work performed by four of these students (four of the coauthors), detailing the satisfactory 3D QSAR models they obtained. Such scientifically complete experiences by undergraduates, made possible by the efficiency of the 3D QSAR methodology, provide exposure to computational tools in the same spirit as traditional laboratory exercises. With the obsolescence of the classic Comparative Molecular Field Analysis Sybyl host, the 3dqsar web portal offers one of the few available means of performing this well-established 3D QSAR method.

Antioxidant Activity of Selected Phenolic Acids-Ferric Reducing Antioxidant Power Assay and QSAR Analysis of the Structural Features

Phenolic acids are naturally occurring compounds that are known for their antioxidant and antiradical activity. We present experimental and theoretical studies on the antioxidant potential of the set of 22 phenolic acids with different models of hydroxylation and methoxylation of aromatic rings. Ferric reducing antioxidant power assay was used to evaluate this property. 2,3-dihydroxybenzoic acid was found to be the strongest antioxidant, while mono hydroxylated and methoxylated structures had the lowest activities. A comprehensive structure-activity investigation with density functional theory methods elucidated the influence of compounds topology, resonance stabilization, and intramolecular hydrogen bonding on the exhibited activity. The key factor was found to be a presence of two or more hydroxyl groups being located in ortho or para position to each other. Finally, the quantitative structure-activity relationship approach was used to build a multiple linear regression model describing the dependence of antioxidant activity on structure of compounds, using features exclusively related to their topology. Coefficients of determination for training set and for the test set equaled 0.9918 and 0.9993 respectively, and Q2 value for leave-one-out was 0.9716. In addition, the presented model was used to predict activities of phenolic acids that haven't been tested here experimentally.

QSAR without borders

Prediction of chemical bioactivity and physical properties has been one of the most important applications of statistical and more recently, machine learning and artificial intelligence methods in chemical sciences. This field of research, broadly known as quantitative structure-activity relationships (QSAR) modeling, has developed many important algorithms and has found a broad range of applications in physical organic and medicinal chemistry in the past 55+ years. This Perspective summarizes recent technological advances in QSAR modeling but it also highlights the applicability of algorithms, modeling methods, and validation practices developed in QSAR to a wide range of research areas outside of traditional QSAR boundaries including synthesis planning, nanotechnology, materials science, biomaterials, and clinical informatics. As modern research methods generate rapidly increasing amounts of data, the knowledge of robust data-driven modelling methods professed within the QSAR field can become essential for scientists working both within and outside of chemical research. We hope that this contribution highlighting the generalizable components of QSAR modeling will serve to address this challenge.

Docking and 3D-QSAR Studies of Hydrazone and Triazole Derivatives for Selective Inhibition of GRK2 over ROCK2

Introduction:

G protein-coupled receptor kinase 2 (GRK2) is known to be implicated in heart failure, and therefore serves as an important drug target. GRK2 belongs to the protein kinase A, G, and C family and shares high sequence similarity with its closely related protein, the Rhoassociated coiled-coil protein kinase 2 (ROCK2). Therefore, selective inhibition of GRK2 over ROCK2 is considered crucial for heart failure therapy.

Objective:

To understand the structural factors for enhancing the inhibitory activity for GRK2 and selectivity over ROCK2, we analyzed and compared molecular interactions using the same set of ligands against both receptors.

Methods:

We have performed molecular docking and three-dimensional quantitative structure activity relationship (3D-QSAR) studies on a series of hydrazone and triazole derivatives.

Results:

The presence of hydrophobic substituents at the triazole ring, electronegative substituents between the pyridine and triazole ring and hydrophobic substituents near the benzene ring increases the activity of both kinases. Whereas, having non-bulky substituents near the triazole ring, bulky and hydrophobic substations at the benzene ring and electronegative and H-bond acceptor substituents at the triazole ring showed a higher inhibitory preference for GRK2 over ROCK2.

Conclusion:

The outcome of this study may be used in the future development of potent GRK2 inhibitors having ROCK2 selectivity.

A Simple Approach to the Toxicity Prediction of Anilines and Phenols Towards Aquatic Organisms

Chemicals pollution in the environment has attracted attention all over the world, and the toxicity prediction of chemical pollutants has become quite important. In this paper, we introduce a simple approach to predict the toxicity of some chemical components, in which the Tchebichef image moment (TM) method was employed to extract useful chemical information from the images of molecular structures to establish quantitative structure-activity relationship (QSAR) prediction models. The proposed approach was applied to predict the toxicity of anilines and phenols for the aquatic organisms of P. subcapitata and V. fischeri, in which the obtained TMs were defined as the independent variables, while the biological toxicity (pEC₅₀) was regarded to be the dependent variable. Then, the predictive models were established by stepwise regression, respectively. The obtained squared correlation coefficients of leave-one-out cross-validation (Q²) for training sets and the predictive squared correlation coefficients (R_p²) for test sets of the two groups of data were higher than 0.79 and 0.75, respectively, which indicated that the obtained models possessed satisfactory accuracy and reliability. Compared with several reported methods, the proposed approach was more convenient and has a higher predictive capability. Our study provides another perspective in QSAR research.

Towards reproducible computational drug discovery

The reproducibility of experiments has been a long standing impediment for further scientific progress. Computational methods have been instrumental in drug discovery efforts owing to its multifaceted utilization for data collection, pre-processing, analysis and inference. This article provides an in-depth coverage on the reproducibility of computational drug discovery. This review explores the following topics: (1) the current state-of-the-art on reproducible research, (2) research documentation (e.g. electronic laboratory notebook, Jupyter notebook, etc.), (3) science of reproducible research (i.e. comparison and contrast with related concepts as replicability, reusability and reliability), (4) model development in computational drug discovery, (5) computational issues on model development and deployment, (6) use case scenarios for streamlining the computational drug discovery protocol. In computational disciplines, it has become common practice to share data and programming codes used for numerical calculations as to not only facilitate reproducibility, but also to foster collaborations (i.e. to drive the project further by introducing new ideas, growing the data, augmenting the code, etc.). It is therefore inevitable that the field of computational drug design would adopt an open approach towards the collection, curation and sharing of data/code.

Prediction of Lipophilicity of some Quinolone Derivatives by using Quantitative Structure-Activity Relationship

OBJECTIVES

Quantitative structure activity relationship (QSAR) was used to study the partition coefficient of some quinolones and their derivatives.

METHODS

These molecules are broad-spectrum antibiotic pharmaceutics. First, data were divided into two categories of train and test (validation) sets using a random selection method. Second, three approaches, including stepwise selection (STS) (forward), genetic algorithm (GA), and simulated annealing (SA) were used to select the descriptors, to examine the effect feature selection methods. To find the relation between descriptors and partition coefficient, multiple linear regression (MLR), principal component regression (PCR) and partial least squares (PLS) were used.

RESULTS

QSAR study showed that both regression and descriptor selection methods have a vital role in the results. Different statistical metrics showed that the MLR-SA approach with (r2=0.96, q2=0.91, pred_r2=0.95) gives the best outcome.

CONCLUSION

The proposed expression by the MLR-SA approach can be used in the better design of novel quinolones and their derivatives.

Insights into the EGFR SAR of N-phenylquinazolin-4-amine-derivatives using quantum mechanical pairwise-interaction energies

Protein kinases are an important class of enzymes that play an essential role in virtually all major disease areas. In addition, they account for approximately 50% of the current targets pursued in drug discovery research. In this work, we explore the generation of structure-based quantum mechanical (QM) quantitative structure-activity relationship models (QSAR) as a means to facilitate structure-guided optimization of protein kinase inhibitors. We explore whether more accurate, interpretable QSAR models can be generated for a series of 76 N-phenylquinazolin-4-amine inhibitors of epidermal growth factor receptor (EGFR) kinase by comparing and contrasting them to other standard QSAR methodologies. The QM-based method involved molecular docking of inhibitors followed by their QM optimization within a ~ 300 atom cluster model of the EGFR active site at the M062X/6-31G(d,p) level. Pairwise computations of the interaction energies with each active site residue were performed. QSAR models were generated by splitting the datasets 75:25 into a training and test set followed by modelling using partial least squares (PLS). Additional QSAR models were generated using alignment dependent CoMFA and CoMSIA methods as well as alignment independent physicochemical, e-state indices and fingerprint descriptors. The structure-based QM-QSAR model displayed good performance on the training and test sets (r² ~ 0.7) and was demonstrably more predictive than the QSAR models built using other methods. The descriptor coefficients from the QM-QSAR models allowed for a detailed rationalization of the active site SAR, which has implications for subsequent design iterations.

Value of p-Value

The goal of this manuscript is to discuss important aspects of external validation of classification and category Quantitative Structure - Activity/Property/Toxicity Relationship QS/A/P/T/R models that to the best of author's knowledge are not addressed in publications. Statistical significance (in terms of p-value) and accuracy of prediction (in terms of Correct Classification Rate (CCR)) of external validation set compounds are among most important characteristics of the models. We assert that in most cases the models built for classification or category response variable should be statistically significant and predictive for each class or category. We show that three thresholds of the number of compounds in each class or category of the external validation sets should be satisfied. 1) The p-value criterion can never be satisfied, if the number of compounds is below the first threshold. 2) If the number of compounds is between the first and the second thresholds, p-value criterion should be used. 3) If it is higher than the third threshold, classification or category accuracy criterion should be used. 4) If the number of compounds is between second and third thresholds, either one or the other criterion should be used depending on the value of p-value. 5) When the number of compounds in the class approaches infinity, the maximum relative error of prediction approaches the relative expected error. The results are of interest in other areas of multidimensional data analysis.

Quantitative structure-activity relationship models for compounds with anticonvulsant activity

Introduction: Third-generation antiepileptic drugs have seemingly failed to improve the global figures of seizure control and can still be regarded as symptomatic treatments. Quantitative structure-activity relationships (QSAR) can be used to guide hit-to-lead and lead optimization projects and applied to the large-scale virtual screening of chemical libraries. Areas covered: In this review, the authors cover reports on QSAR models related to antiepileptic drugs and drug targets in epilepsy, analyzing whether they refer to classic or non-classic QSAR and if they apply QSAR as a descriptive or predictive approach, among other considerations. The article finally focuses on a more detailed discussion of those predictive studies which include some sort of experimental validation, i.e. papers in which the reported models have been used to identify novel active compounds which have been tested in vitro and/or in vivo. Expert opinion: There are significant opportunities to apply the QSAR methodology to assist the discovery of more efficacious antiepileptic drugs. Considering the intrinsic complexity of the disorder, such applications should focus on state-of-the-art approximations (e.g. systemic, multi-target and multi-scale QSAR as well as ensemble and deep learning) and modeling the effects on novel drug targets and modern screening tools.

Design and synthesis of anticancer 1-hydroxynaphthalene-2-carboxanilides with a p53 independent mechanism of action

A series of 116 small-molecule 1-hydroxynaphthalene-2-carboxanilides was designed based on the fragment-based approach and was synthesized according to the microwave-assisted protocol. The biological activity of all of the compounds was tested on human colon carcinoma cell lines including a deleted TP53 tumor suppressor gene. The mechanism of activity was studied according to the p53 status in the cell. Several compounds revealed a good to excellent activity that was similar to or better than the standard anticancer drugs. Some of these appeared to be more active against the p53 null cells than their wild-type counterparts. Intercalating the properties of these compounds could be responsible for their mechanism of action.