Why QSAR fails: an empirical evaluation using conventional computational approach

QSAR modeling approaches to identify a novel ACE2 inhibitor that selectively bind with the C and N terminals of the ectodomain

Due to the high rates of drug development failure and the massive expenses associated with drug discovery, repurposing existing drugs has become more popular. As a result, we have used QSAR modelling on a large and varied dataset of 657 compounds in an effort to discover both explicit and subtle structural features requisite for ACE2 inhibitory activity, with the goal of identifying novel hit molecules. The QSAR modelling yielded a statistically robust QSAR model with high predictivity (R2tr=0.84, R2ex=0.79), previously undisclosed features, and novel mechanistic interpretations. The developed QSAR model predicted the ACE2 inhibitory activity (PIC50) of 1615 ZINC FDA compounds. This led to the detection of a PIC50 of 8.604 M for the hit molecule (ZINC000027990463). The hit molecule's docking score is -9.67 kcal/mol (RMSD 1.4). The hit molecule revealed 25 interactions with the residue ASP40, which defines the N and C termini of the ectodomain of ACE2. The HIT molecule conducted more than thirty contacts with water molecules and exhibited polar interaction with the ARG522 residue coupled with the second chloride ion, which is 10.4 nm away from the zinc ion. Both molecular docking and QSAR produced comparable findings. Moreover, MD simulation and MMGBSA studies verified docking analysis. The MD simulation showed that the hit molecule-ACE2 receptor complex is stable for 400 ns, suggesting that repurposed hit molecule 3 is a viable ACE2 inhibitor.

Mechanistic QSAR modeling derived virtual screening, drug repurposing, ADMET and in-vitro evaluation to identify anticancer lead as lysine-specific demethylase 5a inhibitor

A lysine-specific demethylase is an enzyme that selectively eliminates methyl groups from lysine residues. KDM5A, also known as JARID1A or RBP2, belongs to the KDM5 Jumonji histone demethylase subfamily. To identify novel molecules that interact with the LSD5A receptor, we created a quantitative structure-activity relationship (QSAR) model. A group of 435 compounds was used in a study of the quantitative relationship between structure and activity to guess the IC50 values for blocking LASD5A. We used a genetic algorithm-multilinear regression-based quantitative structure-activity connection model to forecast the bioactivity (PIC50) of 1615 food and drug administration pharmaceuticals from the zinc database with the goal of repurposing clinically used medications. We used molecular docking, molecular dynamic simulation modelling, and molecular mechanics generalised surface area analysis to investigate the molecule's binding mechanism. A genetic algorithm and multi-linear regression method were used to make six variable-based quantitative structure-activity relationship models that worked well (R2 = 0.8521, Q2LOO = 0.8438, and Q2LMO = 0.8414). ZINC000000538621 was found to be a new hit against LSD5A after a quantitative structure-activity relationship-based virtual screening of 1615 zinc food and drug administration compounds. The docking analysis revealed that the hit molecule 11 in the KDM5A binding pocket adopted a conformation similar to the pdb-6bh1 ligand (docking score: -8.61 kcal/mol). The results from molecular docking and the quantitative structure-activity relationship were complementary and consistent. The most active lead molecule 11, which has shown encouraging results, has good absorption, distribution, metabolism, and excretion (ADME) properties, and its toxicity has been shown to be minimal. In addition, the MTT assay of ZINC000000538621 with MCF-7 cell lines backs up the in silico studies. We used molecular mechanics generalise borne surface area analysis and a 200-ns molecular dynamics simulation to find structural motifs for KDM5A enzyme interactions. Thus, our strategy will likely expand food and drug administration molecule repurposing research to find better anticancer drugs and therapies.Communicated by Ramaswamy H. Sarma.

Pharmacophore Synergism in Diverse Scaffold Clinches in Aurora Kinase B

Aurora kinase B (AKB) is a crucial signaling kinase with an important role in cell division. Therefore, inhibition of AKB is an attractive approach to the treatment of cancer. In the present work, extensive quantitative structure-activity relationships (QSAR) analysis has been performed using a set of 561 structurally diverse aurora kinase B inhibitors. The Organization for Economic Cooperation and Development (OECD) guidelines were used to develop a QSAR model that has high statistical performance (R²_tr = 0.815, Q²_LMO = 0.808, R²_ex = 0.814, CCC_ex = 0.899). The seven-variable-based newly developed QSAR model has an excellent balance of external predictive ability (Predictive QSAR) and mechanistic interpretation (Mechanistic QSAR). The QSAR analysis successfully identifies not only the visible pharmacophoric features but also the hidden features. The analysis indicates that the lipophilic and polar groups-especially the H-bond capable groups-must be present at a specific distance from each other. Moreover, the ring nitrogen and ring carbon atoms play important roles in determining the inhibitory activity for AKB. The analysis effectively captures reported as well as unreported pharmacophoric features. The results of the present analysis are also supported by the reported crystal structures of inhibitors bound to AKB.

QSAR, Molecular Docking, MD Simulation and MMGBSA Calculations Approaches to Recognize Concealed Pharmacophoric Features Requisite for the Optimization of ALK Tyrosine Kinase Inhibitors as Anticancer Leads

ALK tyrosine kinase ALK TK is an important target in the development of anticancer drugs. In the present work, we have performed a QSAR analysis on a dataset of 224 molecules in order to quickly predict anticancer activity on query compounds. Double cross validation assigns an upward plunge to the genetic algorithm−multi linear regression (GA-MLR) based on robust univariate and multivariate QSAR models with high statistical performance reflected in various parameters like, fitting parameters; R2 = 0.69−0.87, F = 403.46−292.11, etc., internal validation parameters; Q2LOO = 0.69−0.86, Q2LMO = 0.69−0.86, CCCcv = 0.82−0.93, etc., or external validation parameters Q2F1 = 0.64−0.82, Q2F2 = 0.63−0.82, Q2F3 = 0.65−0.81, R2ext = 0.65−0.83 including RMSEtr < RMSEcv. The present QSAR evaluation successfully identified certain distinct structural features responsible for ALK TK inhibitory potency, such as planar Nitrogen within four bonds from the Nitrogen atom, Fluorine atom within five bonds beside the non-ring Oxygen atom, lipophilic atoms within two bonds from the ring Carbon atoms. Molecular docking, MD simulation, and MMGBSA computation results are in consensus with and complementary to the QSAR evaluations. As a result, the current study assists medicinal chemists in prioritizing compounds for experimental detection of anticancer activity, as well as their optimization towards more potent ALK tyrosine kinase inhibitor.

Perceiving the Concealed and Unreported Pharmacophoric Features of the 5-Hydroxytryptamine Receptor Using Balanced QSAR Analysis

The 5-hydroxytryptamine receptor 6 (5-HT6) has gained attention as a target for developing therapeutics for Alzheimer’s disease, schizophrenia, cognitive dysfunctions, anxiety, and depression, to list a few. In the present analysis, a larger and diverse dataset of 1278 molecules covering a broad chemical and activity space was used to identify visual and concealed structural features associated with binding affinity for 5-HT6. For this, quantitative structure–activity relationships (QSAR) and molecular docking analyses were executed. This led to the development of a statistically robust QSAR model with a balance of excellent predictivity (R²_tr = 0.78, R²_ex = 0.77), the identification of unreported aspects of known features, and also novel mechanistic interpretations. Molecular docking and QSAR provided similar as well as complementary results. The present analysis indicates that the partial charges on ring carbons present within four bonds from a sulfur atom, the occurrence of sp3-hybridized carbon atoms bonded with donor atoms, and a conditional occurrence of lipophilic atoms/groups from nitrogen atoms, which are prominent but unreported pharmacophores that should be considered while optimizing a molecule for 5-HT6. Thus, the present analysis led to identification of some novel unreported structural features that govern the binding affinity of a molecule. The results could be beneficial in optimizing the molecules for 5-HT6.

Mechanistic Analysis of Chemically Diverse Bromodomain-4 Inhibitors Using Balanced QSAR Analysis and Supported by X-ray Resolved Crystal Structures

Bromodomain-4 (BRD-4) is a key enzyme in post-translational modifications, transcriptional activation, and many other cellular processes. Its inhibitors find their therapeutic usage in cancer, acute heart failure, and inflammation to name a few. In the present study, a dataset of 980 molecules with a significant diversity of structural scaffolds and composition was selected to develop a balanced QSAR model possessing high predictive capability and mechanistic interpretation. The model was built as per the OECD (Organisation for Economic Co-operation and Development) guidelines and fulfills the endorsed threshold values for different validation parameters (R²_tr = 0.76, Q²_LMO = 0.76, and R²_ex = 0.76). The present QSAR analysis identified that anti-BRD-4 activity is associated with structural characters such as the presence of saturated carbocyclic rings, the occurrence of carbon atoms near the center of mass of a molecule, and a specific combination of planer or aromatic nitrogen with ring carbon, donor, and acceptor atoms. The outcomes of the present analysis are also supported by X-ray-resolved crystal structures of compounds with BRD-4. Thus, the QSAR model effectively captured salient as well as unreported hidden pharmacophoric features. Therefore, the present study successfully identified valuable novel pharmacophoric features, which could be beneficial for the future optimization of lead/hit compounds for anti-BRD-4 activity.

Target Specific Inhibition of Protein Tyrosine Kinase in Conjunction With Cancer and SARS-COV-2 by Olive Nutraceuticals

The fact that viruses cause human cancer dates back to the early 1980s. By reprogramming cellular signaling pathways, viruses encoded protein that can regulate altered control of cell cycle events. Viruses can interact with a superfamily of membrane bound protein, receptor tyrosine kinase to modulate their activity in order to increase virus entrance into cells and promotion of viral replication within the host. Therefore, our study aimed at screening of inhibitors of tyrosine kinase using natural compounds from olive. Protein tyrosine kinase (PTK) is an important factor for cancer progression and can be linked to coronavirus. It is evident that over expression of Protein tyrosine kinase (PTK) enhance viral endocytosis and proliferation and the use of tyrosine kinase inhibitors reduced the period of infection period. Functional network studies were carried out using two major PTKs viz. Anaplastic lymphoma kinase (ALK) and B-lymphocytic kinase (BTK). They are associated with coronavirus in regulation of cell signaling proteins for cellular processes. We virtually screened for 161 library of natural compounds from olive found overexpressed in ALK and BTK in metastatic as well as virus host cells. We have employed both ligand and target-based approach for drug designing by high throughput screening using Multilinear regression model based QSAR and docking. The QSAR based virtual screening of 161 olive nutraceutical compounds has successfully identified certain new hit; Wedelosin, in which, the descriptor rsa (ratio of molecular surface area to the solvent accessible surface area) plays crucial role in deciding Wedelosin’s inhibitory potency. The best-docked olive nutraceuticals further investigated for the stability and effectivity of the BTK and ALK during in 150 ns molecular dynamics and simulation. Post simulation analysis and binding energy estimation in MMGBSA further revealed the intensive potential of the olive nutraceuticals in PTK inhibition. This study is therefore expected to widen the use of nutraceuticals from olive in cancer as well as SARS-CoV2 alternative therapy.

An AI-based prediction model for drug-drug interactions in osteoporosis and Paget's diseases from SMILES

Referring to common skeletal-related diseases, osteoporosis and Paget's are two of the most frequently found diseases in the elderly. Nowadays, the combination of multiple drugs is the optimal therapy to decelerate osteoporosis and Paget's pathologic process, which contains various underlying adverse effects due to drug-drug interactions (DDIs). Artificial intelligence (AI) has the potential to evaluate the interaction, pharmacodynamics, and possible side effects between drugs. In this research, we created an AI-based machine-learning model to predict the outcomes of interactions between drugs used for osteoporosis and Paget's treatment, furthermore, to mitigate cost and time in implementing the best combination of medications in clinical practice. Our dataset was collected from the DrugBank database, and we then extracted a variety of chemical features from the simplified molecular-input line-entry system (SMILES) of defined drug pairs that interact with each other. Finally, machine-learning algorithms have been implemented to learn the extracted features. Our stack ensemble model from Random Forest and XGBoost reached an average accuracy of 74% in predicting DDIs. It was superior to individual models and previous methods in most measurement metrics. This study showed the potential of AI models in predicting DDIs of Osteoporosis-Paget's disease in particular, and other diseases in general.

QSAR based virtual screening derived identification of a novel hit as a SARS CoV-229E 3CLpro Inhibitor: GA-MLR QSAR modeling supported by molecular Docking, molecular dynamics simulation and MMGBSA calculation approaches

Congruous coronavirus drug targets and analogous lead molecules must be identified as quickly as possible to produce antiviral therapeutics against human coronavirus (HCoV SARS 3CLpro) infections. In the present communication, we bear recognized a HIT candidate for HCoV SARS 3CLpro inhibition. Four Parametric GA-MLR primarily based QSAR model (R2:0.84, R2adj:0.82, Q2loo: 0.78) was once promoted using a dataset over 37 structurally diverse molecules along QSAR based virtual screening (QSAR-VS), molecular docking (MD) then molecular dynamic simulation (MDS) analysis and MMGBSA calculations. The QSAR-based virtual screening was utilized to find novel lead molecules from an in-house database of 100 molecules. The QSAR-vS successfully offered a hit molecule with an improved _PEC₅₀ value from 5.88 to 6.08. The benzene ring, phenyl ring, amide oxygen and nitrogen, and other important pharmacophoric sites are revealed via MD and MDS studies. Ile164, Pro188, Leu190, Thr25, His41, Asn46, Thr47, Ser49, Asn189, Gln191, Thr47, and Asn141 are among the key amino acid residues in the S1 and S2 pocket. A stable complex of a lead molecule with the HCoV SARS 3CLpro was discovered using MDS. MM-GBSA calculations resulted from MD simulation results well supported with the binding energies calculated from the docking results. The results of this study can be exploited to develop a novel antiviral target, such as an HCoV SARS 3CLpro Inhibitor.

Topological Distance-Based Electron Interaction Tensor to Apply a Convolutional Neural Network on Drug-like Compounds

Deep learning (DL) models in quantitative structure-activity relationship fed the molecular structure directly to the network without using human-designed descriptors by representing molecule as a graph or string (e.g., SMILES code). However, these two representations were oversimplification of real molecules to reflect chemical properties of molecular structures. Given that the choice of molecular representation determines the architecture of the DL model to apply, a novel way of molecular representation can open a way to apply diverse DL networks developed and used in other fields. A topological distance-based electron interaction (TDEi) tensor has been developed in this study inspired by the quantum mechanical model of the molecule, which defines a molecule with electrons and protons. In the TDEi tensor, the atomic orbital (AO) of each atom is represented by an electron configuration (EC) vector, which is a bit string based on the presence and absence of electrons in each AO according to spin indicated by positive and negative signs. Interactions between EC vectors were calculated based on the topological distance between atoms in a molecule. As a molecular structure was translated into 3D array, CNN models (modified VGGNet) were applied using a TDEi tensor to predict four physicochemical properties of drug-like compound datasets: MP (275,131), Lipop (4193), Esol (1127), and Freesolv (639). Models achieved good prediction accuracy. PCA showed that a stronger correlation was observed between the extracted features and the target endpoint as features were extracted from the deeper layer.

QSAR, pharmacophore modeling and molecular docking studies to identify structural alerts for some nitrogen heterocycles as dual inhibitor of telomerase reverse transcriptase and human telomeric G-quadruplex DNA

Background

Telomerase reverse transcriptase (TERT) and human telomeric G-quadruplex DNA are amongst the favorable target for researchers to discover novel and more effective anticancer agents. To understand and elucidate structure activity relationship and mechanism of inhibition of telomerase reverse transcriptase (TERT) and human telomeric G-quadruplex DNA, a QSAR modeling and molecular docking were conducted.

Results

Two robust QSAR model were obtained which consist of full set QSAR model (R2: 0.8174, CCCtr: 0.8995, Q2loo: 0.7881, Q2LMO: 0.7814) and divided set QSAR model (R2: 0.8217, CCCtr: 0.9021, Q2loo: 0.7886, Q2LMO: 0.7783, Q2-F1: 0.7078, Q2-F2: 0.6865, Q2-F3: 0.7346) for envisaging the inhibitory activity of telomerase reverse transcriptase (TERT) and human telomeric G-quadruplex DNA. The analysis reveals that carbon atom exactly at 3 bonds from aromatic carbon atom, nitrogen atom exactly at six bonds from planer nitrogen atom, aromatic carbon atom within 2 A0 from the center of mass of molecule and occurrence of element hydrogen within 2 A0 from donar atom are the key pharmacophoric features important for dual inhibition of TERT and human telomeric G-quadruplex DNA. To validate this analysis, pharmacophore modeling and the molecular docking is performed. Molecular docking analysis support QSAR analysis and revealed that, dual inhibition of TERT and human telomeric DNA is mainly contributed from hydrophobic and hydrogen bonding interactions.

Conclusion

The findings of molecular docking, pharmacophore modelling, and QSAR are all consistent and in strong agreement. The validated QSAR analyses can detect structural alerts, pharmacophore modelling can classify a molecule's consensus pharmacophore involving hydrophobic and acceptor regions, whereas docking analysis can reveal the mechanism of dual inhibition of telomerase reverse transcriptase (TERT) and human telomeric G-quadruplex DNA. The combination of QSAR, pharmacophore modeling and molecular docking may be useful for the future drug design of dual inhibitors to combat the devastating issue of resistance.

Graphical abstract

QSAR analysis of sodium glucose co-transporter 2 (SGLT2) inhibitors for anti-hyperglycaemic lead development

QSAR (Quantitative Structure Activity Relationship) modelling was performed on a dataset of 90 sodium-dependent glucose cotransporter 2 (SGLT2) inhibitors. The quantitative and explicative evaluations revealed some of the subtle and distinguished structural features that are responsible for the inhibitory potency of these compounds against SGLT2, such as less possible number of ring carbons at 8 Å from the lipophilic atoms in the molecule (fringClipo8A) and more possible value for the sum of the partial charges of the lipophilic atoms present within seven bonds from the donor atoms (lipo_don_7Bc). Multivariate GA-MLR (genetic algorithm-multi linear regression) and thorough validation methodology out-turned a statistically robust QSAR model with a very high predictability shown from various statistical parameters. A QSAR model with r² = 0.83, F = 51.54, Q²_LOO = 0.79, Q²_LMO = 0.79, CCC_cv = 0.88, Q²Fⁿ = 0.76-0.81, r²_ext = 0.77, CCC_ext = 0.85, and with RMSE_tr < RMSE_cv was proposed. This QSAR model will assist synthetic chemists in the development of the SGLT2 inhibitors as the antidiabetic leads.

Quantitative Structure-Activity Relationship Evaluation of MDA-MB-231 Cell Anti-Proliferative Leads

In the present endeavor, for the dataset of 219 in vitro MDA-MB-231 TNBC cell antagonists, a (QSAR) quantitative structure–activity relationships model has been carried out. The quantitative and explicative assessments were performed to identify inconspicuous yet pre-eminent structural features that govern the anti-tumor activity of these compounds. GA-MLR (genetic algorithm multi-linear regression) methodology was employed to build statistically robust and highly predictive multiple QSAR models, abiding by the OECD guidelines. Thoroughly validated QSAR models attained values for various statistical parameters well above the threshold values (i.e., R² = 0.79, Q²_LOO = 0.77, Q²_LMO = 0.76–0.77, Q²-Fⁿ = 0.72–0.76). Both de novo QSAR models have a sound balance of descriptive and statistical approaches. Decidedly, these QSAR models are serviceable in the development of MDA-MB-231 TNBC cell antagonists.

Mechanistic and Predictive QSAR Analysis of Diverse Molecules to Capture Salient and Hidden Pharmacophores for Anti-Thrombotic Activity

Thrombosis is a life-threatening disease with a high mortality rate in many countries. Even though anti-thrombotic drugs are available, their serious side effects compel the search for safer drugs. In search of a safer anti-thrombotic drug, Quantitative Structure-Activity Relationship (QSAR) could be useful to identify crucial pharmacophoric features. The present work is based on a larger data set comprising 1121 diverse compounds to develop a QSAR model having a balance of acceptable predictive ability (Predictive QSAR) and mechanistic interpretation (Mechanistic QSAR). The developed six parametric model fulfils the recommended values for internal and external validation along with Y-randomization parameters such as R²_tr = 0.831, Q²_LMO = 0.828, R²_ex = 0.783. The present analysis reveals that anti-thrombotic activity is found to be correlated with concealed structural traits such as positively charged ring carbon atoms, specific combination of aromatic Nitrogen and sp2-hybridized carbon atoms, etc. Thus, the model captured reported as well as novel pharmacophoric features. The results of QSAR analysis are further vindicated by reported crystal structures of compounds with factor Xa. The analysis led to the identification of useful novel pharmacophoric features, which could be used for future optimization of lead compounds.

Dapsone is not a Pharmacodynamic Lead Compound for its Aryl Derivatives

BACKGROUND

The relatedness between the linear equations of thermodynamics and QSAR was studied thanks to the recently elucidated crystal structure complexes between sulfonamide pterin conjugates and dihydropteroate synthase (DHPS) together with a published set of thirty- six synthetic dapsone derivatives with their reported entropy-driven activity data. Only a few congeners were slightly better than dapsone.

OBJECTIVE

Our study aimed at demonstrating the applicability of thermodynamic QSAR and to shed light on the mechanistic aspects of sulfone binding to DHPS.

METHODS

To this end ligand docking to DHPS, quantum mechanical properties, 2D- and 3D-QSAR as well as Principle Component Analysis (PCA) were carried out.

RESULTS

The short aryl substituents of the docked pterin-sulfa conjugates were outward oriented into the solvent space without interacting with target residues which explains why binding enthalpy (ΔH) did not correlate with potency. PCA revealed how chemically informative descriptors are evenly loaded on the first three PCs (interpreted as ΔG, ΔH and ΔS), while chemically cryptic ones reflected higher dimensional (complex) loadings.

CONCLUSION

It is safe to utter that synthesis efforts to introduce short side chains for aryl derivatization of the dapsone scaffold have failed in the past. On theoretical grounds we provide computed evidence why dapsone is not a pharmacodynamic lead for drug profiling because enthalpic terms do not change significantly at the moment of ligand binding to target.

QSAR and Pharmacophore Modeling of Nitrogen Heterocycles as Potent Human N-Myristoyltransferase (Hs-NMT) Inhibitors

N-myristoyltransferase (NMT) is an important eukaryotic monomeric enzyme which has emerged as an attractive target for developing a drug for cancer, leishmaniasis, ischemia-reperfusion injury, malaria, inflammation, etc. In the present work, statistically robust machine leaning models (QSAR (Quantitative Structure–Activity Relationship) approach) for Human NMT (Hs-NMT) inhibitory has been performed for a dataset of 309 Nitrogen heterocycles screened for NMT inhibitory activity. Hundreds of QSAR models were derived. Of these, the model 1 and 2 were chosen as they not only fulfil the recommended values for a good number of validation parameters (e.g., R² = 0.77–0.79, Q²_LMO = 0.75–0.76, CCC_ex = 0.86–0.87, Q²-F³ = 0.74–0.76, etc.) but also provide useful insights into the structural features that sway the Hs-NMT inhibitory activity of Nitrogen heterocycles. That is, they have an acceptable equipoise of descriptive and predictive qualities as per Organisation for Economic Co-operation and Development (OECD) guidelines. The developed QSAR models identified a good number of molecular descriptors like solvent accessible surface area of all atoms having specific partial charge, absolute surface area of Carbon atoms, etc. as important features to be considered in future optimizations. In addition, pharmacophore modeling has been performed to get additional insight into the pharmacophoric features, which provided additional results.

Method Optimisation in Hydrophilic-Interaction Liquid Chromatography by Design of Experiments Combined with Quantitative Structure–Retention Relationships

Accurate prediction of the separation conditions for a set of target analytes with no retention data available is fundamental for routine analytical assays but remains a very challenging task. In this paper, a quality by design (QbD) optimisation workflow capable of discovering the optimal chromatographic conditions for separation of new compounds in hydrophilic-interaction liquid chromatography (HILIC) is introduced. This workflow features the application of quantitative structure−retention relationship (QSRR) methodology in conjunction with design of experiments (DoE) principles and was used to carry out a two-level full factorial DoE optimisation for a mixture of pharmaceutical analytes on zwitterionic, amide, amine, and bare silica HILIC stationary phases, with mobile phases containing varying acetonitrile content, mobile phase pH, and salt concentration. A dual-filtering approach that considers both retention time (tR) and structural similarity was used to identify the optimal set of analytes to train the QSRR in order to maximise prediction accuracy. Highly predictive retention models (average R2 of 0.98) were obtained and statistical analysis of the prediction performance of the QSRR models demonstrated their ability to predict the retention times of new compounds based solely on their molecular structures, with root-mean-square errors of prediction in the range 7.6–11.0 %. Further, the obtained retention data for pharmaceutical test compounds were used to compute their separation selectivity, which was used as input into a DoE optimiser in order to select the optimal separation conditions. Experimental separations performed under the chosen optimal working conditions showed good agreement with the theoretical predictions. To the best of our knowledge, this is the first study of a QbD optimisation workflow assisted with dual-filtering-based retention modelling to facilitate the method development process in HILIC.

Structure features of peptide-type SARS-CoV main protease inhibitors: Quantitative structure activity relationship study

In the present work, an extensive QSAR (Quantitative Structure Activity Relationships) analysis of a series of peptide-type SARS-CoV main protease (MPro) inhibitors following the OECD guidelines has been accomplished. The analysis was aimed to identify salient and concealed structural features that govern the MPro inhibitory activity of peptide-type compounds. The QSAR analysis is based on a dataset of sixty-two peptide-type compounds which resulted in the generation of statistically robust and highly predictive multiple models. All the developed models were validated extensively and satisfy the threshold values for many statistical parameters (for e.g. R² ​= ​0.80-0.82, Q² _loo ​= ​0.74-0.77, Q ² _LMO  ​= ​0.66-0.67). The developed QSAR models identified number of sp² hybridized Oxygen atoms within seven bonds from aromatic Carbon atoms, the presence of Carbon and Nitrogen atoms at a topological distance of 3 and other interrelations of atom pairs as important pharmacophoric features. Hence, the present QSAR models have a good balance of Qualitative (Descriptive QSARs) and Quantitative (Predictive QSARs) approaches, therefore useful for future modifications of peptide-type compounds for anti- SARS-CoV activity.

Computational Approaches in Preclinical Studies on Drug Discovery and Development

Because undesirable pharmacokinetics and toxicity are significant reasons for the failure of drug development in the costly late stage, it has been widely recognized that drug ADMET properties should be considered as early as possible to reduce failure rates in the clinical phase of drug discovery. Concurrently, drug recalls have become increasingly common in recent years, prompting pharmaceutical companies to increase attention toward the safety evaluation of preclinical drugs. In vitro and in vivo drug evaluation techniques are currently more mature in preclinical applications, but these technologies are costly. In recent years, with the rapid development of computer science, in silico technology has been widely used to evaluate the relevant properties of drugs in the preclinical stage and has produced many software programs and in silico models, further promoting the study of ADMET in vitro. In this review, we first introduce the two ADMET prediction categories (molecular modeling and data modeling). Then, we perform a systematic classification and description of the databases and software commonly used for ADMET prediction. We focus on some widely studied ADMT properties as well as PBPK simulation, and we list some applications that are related to the prediction categories and web tools. Finally, we discuss challenges and limitations in the preclinical area and propose some suggestions and prospects for the future.

Quantum Mechanical Investigation of Three-Dimensional Activity Cliffs Using the Molecules-in-Molecules Fragmentation-Based Method

The concept of activity cliff (AC) (i.e., a small structural modification resulting in a substantial bioactivity change) is widely encountered in medicinal chemistry during compound design. Whereas the study of ACs is of high interest as it provides a wealth of opportunities for effective drug design, its practical application in the actual drug development process has been difficult because of significant computational challenges. To provide some understanding of the ACs, we have carried out a rigorous quantum-mechanical investigation of the electronic interactions of a wide range of ACs (205 cliffs formed by 261 protein-ligand complexes covering 37 different receptor types) using multilayer molecules-in-molecules (MIM) fragmentation-based methodology. The MIM methodology enables performing accurate high-level quantum mechanical (QM) calculations at a substantially lower computational cost, while allowing for a quantitative decomposition of the protein-ligand binding energy into the contributions from individual residues, solvation, and entropy. Our investigation in this study is mainly focused on whether the QM binding energy calculation can correctly identify the higher potency cliff partner for a given ligand pair having a sufficiently high activity difference. We have also analyzed the effect of including crystal water molecules as a part of the receptor as well as the impact of ligand desolvation energy on the correct identification of the more potent ligand in a cliff pair. Our analysis reveals that, in the majority of the cases, the AC prediction could be significantly improved by carefully identifying the critical crystal water molecules, whereas the contribution from the ligand desolvation also remains essential. Additionally, we have exploited the residue-specific interaction energies provided by MIM to identify the key residues and interaction hot-spots that are responsible for the experimentally observed drastic activity changes. The results show that our MIM fragmentation-based protocol provides comprehensive interaction energy profiles that can be employed to understand the distinctiveness of ligand modifications, for potential applications in structure-based drug design.

Quinoxalinones Based Aldose Reductase Inhibitors: 2D and 3D-QSAR Analysis

In the present work, 2D- and 3D-quantitative structure-activity relationship (QSAR) analysis has been employed for a diverse set of eighty-nine quinoxalinones to identify the pharmacophoric features with significant correlation with the aldose reductase inhibitory activity. Using genetic algorithm (GA) as a variable selection method, multivariate linear regression (MLR) models were derived using a pool of molecular descriptors. All the six-descriptor based GA-MLR QSAR models are statistically robust with coefficient of determination (R² )>0.80 and cross-validated R² >0.77. The derived GA-MLR models were thoroughly validated using internal and external and Y-scrambling techniques. The CoMFA like model, which is based on a combination of steric and electrostatic effects and graphically inferred using contour plots, is highly robust with R² >0.93 and cross-validated R² >0.73. The established QSAR and CoMFA like models are proficient in identify key pharmacophoric features that govern the aldose reductase inhibitory activity of quinoxalinones.

Best Practices for QSAR Model Reporting: Physical and Chemical Properties, Ecotoxicity, Environmental Fate, Human Health, and Toxicokinetics Endpoints

BACKGROUND

Quantitative and qualitative structure–activity relationships (QSARs) have been used to understand chemical behavior for almost a century. The main source of QSAR models is the scientific literature, but the open question is how well these models are documented.

OBJECTIVES

The main aim of this study was to critically analyze the publication practices of QSARs with regard to transparency, potential reproducibility, and independent verification. The focus was on the level of technical completeness of the published QSARs.

METHODS

A total of 1,533 QSAR articles reporting 79 individual endpoints, mostly in environmental and health science, were reviewed. The QSAR parameters required for technical completeness were grouped into five categories: chemical structures, experimental endpoint values, descriptor values, mathematical representation of the model, and predicted endpoint values. The data were summarized and discussed using Circos plots.

RESULTS

Altogether, 42.5% of the reviewed articles were found to be potentially reproducible. The potential reproducibility for different endpoint groups varied; the respective rates were 39% for physical and chemical properties, 52% for ecotoxicity, 56% for environmental fate, 30% for human health, and 32% for toxicokinetics. The reproducibility of QSARs is discussed and placed in the context of the reproducibility of the experimental methods. Included are 65 references to open QSAR datasets as examples of models restored from scientific articles.

DISCUSSION

Strikingly poor documentation of QSARs was observed, which reduces the transparency, availability, and consequently, the application of research results in scientific, industrial, and regulatory areas. A list of the components needed to ensure the best practices for QSAR reporting is provided, allowing long-term use and preservation of the models. This list also allows an assessment of the reproducibility of models by interested parties such as journal editors, reviewers, regulators, evaluators, and potential users. https://doi.org/10.1289/EHP3264.

Antiprotozoal Nitazoxanide Derivatives: Synthesis, Bioassays and QSAR Study Combined with Docking for Mechanistic Insight

In view of the serious health problems concerning infectious diseases in heavily populated areas, we followed the strategy of lead compound diversification to evaluate the near-by chemical space for new organic compounds. To this end, twenty derivatives of nitazoxanide (NTZ) were synthesized and tested for activity against Entamoeba histolytica parasites. To ensure drug-likeliness and activity relatedness of the new compounds, the synthetic work was assisted by a quantitative structure-activity relationships study (QSAR). Many of the inherent downsides – well-known to QSAR practitioners – we circumvented thanks to workarounds which we proposed in prior QSAR publication. To gain further mechanistic insight on a molecular level, ligand-enzyme docking simulations were carried out since NTZ is known to inhibit the protozoal pyruvate ferredoxin oxidoreductase (PFOR) enzyme as its biomolecular target.

AutoWeka: toward an automated data mining software for QSAR and QSPR studies

In biology and chemistry, a key goal is to discover novel compounds affording potent biological activity or chemical properties. This could be achieved through a chemical intuition-driven trial-and-error process or via data-driven predictive modeling. The latter is based on the concept of quantitative structure-activity/property relationship (QSAR/QSPR) when applied in modeling the biological activity and chemical properties, respectively, of compounds. Data mining is a powerful technology underlying QSAR/QSPR as it harnesses knowledge from large volumes of high-dimensional data via multivariate analysis. Although extremely useful, the technicalities of data mining may overwhelm potential users, especially those in the life sciences. Herein, we aim to lower the barriers to access and utilization of data mining software for QSAR/QSPR studies. AutoWeka is an automated data mining software tool that is powered by the widely used machine learning package Weka. The software provides a user-friendly graphical interface along with an automated parameter search capability. It employs two robust and popular machine learning methods: artificial neural networks and support vector machines. This chapter describes the practical usage of AutoWeka and relevant tools in the development of predictive QSAR/QSPR models.

AVAILABILITY

The software is freely available at http://www.mt.mahidol.ac.th/autoweka.

Chemical predictive modelling to improve compound quality

The 'quality' of small-molecule drug candidates, encompassing aspects including their potency, selectivity and ADMET (absorption, distribution, metabolism, excretion and toxicity) characteristics, is a key factor influencing the chances of success in clinical trials. Importantly, such characteristics are under the control of chemists during the identification and optimization of lead compounds. Here, we discuss the application of computational methods, particularly quantitative structure-activity relationships (QSARs), in guiding the selection of higher-quality drug candidates, as well as cultural factors that may have affected their use and impact.

Enhanced QSAR model performance by integrating structural and gene expression information

Despite decades of intensive research and a number of demonstrable successes, quantitative structure-activity relationship (QSAR) models still fail to yield predictions with reasonable accuracy in some circumstances, especially when the QSAR paradox occurs. In this study, to avoid the QSAR paradox, we proposed a novel integrated approach to improve the model performance through using both structural and biological information from compounds. As a proof-of-concept, the integrated models were built on a toxicological dataset to predict non-genotoxic carcinogenicity of compounds, using not only the conventional molecular descriptors but also expression profiles of significant genes selected from microarray data. For test set data, our results demonstrated that the prediction accuracy of QSAR model was dramatically increased from 0.57 to 0.67 with incorporation of expression data of just one selected signature gene. Our successful integration of biological information into classic QSAR model provided a new insight and methodology for building predictive models especially when QSAR paradox occurred.

Prediction of acute aquatic toxicity in Tetrahymena pyriformis--'Eco-Derek', a knowledge-based system approach

A 'proof-of-concept' version of a software tool for making transparent predictions of acute aquatic toxicity has been developed. It is primarily limited to semi-quantitative predictions in one species, the ciliated protozoan, Tetrahymena pyriformis. A freely available system, 'Eco-Derek', was derived by adapting a well-established, knowledge-based structure-activity and reasoning platform (Derek for Windows, Lhasa Limited). The Derek reasoning code was modified to express potency rather than confidence. Structure-activity relationship (SAR) development utilised a curated version of a published dataset, supplemented with the CADASTER Challenge datasets. Forty-five structural alerts were produced. The dependence on log P was examined for each alert and entered into the system as qualitative reasoning rules specifying the predicted potency as Very Low, Low, Moderate, High or Very High. Evaluation studies showed: (a) moderate accuracy for the training set but low accuracy for an external test set; (b) non-linearity in the toxicity-log P relationship for chemicals without identified structural alerts; (c) insufficient differentiation of substituent effects in some of the reactivity-based structural alerts resulting in too few chemicals predicted with Very High toxicity; and (d) the need for additional structural alerts covering polar narcosis and less common reactive or metabolically activated chemical functionality.

Health and environment regulation: ethical standards of evidence

The regulation of health and environmental hazards is coercive to the point of imposing substantial fines and even detention to transgressors. In free societies, those regulations have ethical standing if grounded on the scientific evidence of physical measurements relevant to what is being protected including humans, or on transparent use and benefit tradeoffs. These requirements are commonly observed for risks and hazards that are present and measurable, and they sustain much of current regulations that remain uncontroversial. By contrast, most supposable distant risks that cannot be measured are currently regulated on the basis of arbitrary and very costly assumptions of evidence, simulated as if relevant, objective and scientific. Such illusory assumptions are enforced to impart an ersatz moral authority to regulations, which in reality end up as the opaque condensations of political and bureaucratic pressures. The common outcome of such pressures is regulations that whimsically sanction the most minimal exposures that still allow practical applications. Such official misrepresentations of evidence are intolerable and should be removed, because only relevant and factually measured evidence provides moral standing to regulations. When that evidence is not possible, precautionary minimal exposures compatible with necessary applications should be ethically defined by reasoned and transparent use and benefits tradeoffs.

Computational tools and resources for metabolism-related property predictions. 2. Application to prediction of half-life time in human liver microsomes

BACKGROUND

The most important factor affecting metabolic excretion of compounds from the body is their half-life time. This provides an indication of compound stability of, for example, drug molecules. We report on our efforts to develop QSAR models for metabolic stability of compounds, based on in vitro half-life assay data measured in human liver microsomes.

METHOD

A variety of QSAR models generated using different statistical methods and descriptor sets implemented in both open-source and commercial programs (KNIME, GUSAR and StarDrop) were analyzed. The models obtained were compared using four different external validation sets from public and commercial data sources, including two smaller sets of in vivo half-life data in humans.

CONCLUSION

In many cases, the accuracy of prediction achieved on one external test set did not correspond to the results achieved with another test set. The most predictive models were used for predicting the metabolic stability of compounds from the open NCI database, the results of which are publicly available on the NCI/CADD Group web server ( http://cactus.nci.nih.gov ).