Fragment-based QSAR: perspectives in drug design

3D-QSAR, Scaffold Hopping, Virtual Screening, and Molecular Dynamics Simulations of Pyridin-2-one as mIDH1 Inhibitors

Isocitrate dehydrogenase 1 (IDH1) is a necessary enzyme for cellular respiration in the tricarboxylic acid cycle. Mutant isocitrate dehydrogenase 1 (mIDH1) has been detected overexpressed in a variety of cancers. mIDH1 inhibitor ivosidenib (AG-120) was only approved by the Food and Drug Administration (FDA) for marketing, nevertheless, a range of resistance has been frequently reported. In this study, several mIDH1 inhibitors with the common backbone pyridin-2-one were explored using the three-dimensional structure-activity relationship (3D-QSAR), scaffold hopping, absorption, distribution, metabolism, excretion (ADME) prediction, and molecular dynamics (MD) simulations. Comparative molecular field analysis (CoMFA, R2 = 0.980, Q2 = 0.765) and comparative molecular similarity index analysis (CoMSIA, R2 = 0.997, Q2 = 0.770) were used to build 3D-QSAR models, which yielded notably decent predictive ability. A series of novel structures was designed through scaffold hopping. The predicted pIC50 values of C3, C6, and C9 were higher in the model of 3D-QSAR. Additionally, MD simulations culminated in the identification of potent mIDH1 inhibitors, exhibiting strong binding interactions, while the analyzed parameters were free energy landscape (FEL), radius of gyration (Rg), solvent accessible surface area (SASA), and polar surface area (PSA). Binding free energy demonstrated that C2 exhibited the highest binding free energy with IDH1, which was -93.25 ± 5.20 kcal/mol. This research offers theoretical guidance for the rational design of novel mIDH1 inhibitors.

A Computational QSAR, Molecular Docking and In Vitro Cytotoxicity Study of Novel Thiouracil-Based Drugs with Anticancer Activity against Human-DNA Topoisomerase II

Computational chemistry, molecular docking, and drug design approaches, combined with the biochemical evaluation of the antitumor activity of selected derivatives of the thiouracil-based dihydroindeno pyrido pyrimidines against topoisomerase I and II. The IC50 of other cell lines including the normal human lung cell line W138, lung cancer cell line, A549, breast cancer cell line, MCF-7, cervical cancer, HeLa, and liver cancer cell line HepG2 was evaluated using biochemical methods. The global reactivity descriptors and physicochemical parameters were computed, showing good agreement with the Lipinski and Veber's rules of the drug criteria. The molecular docking study of the ligands with the topoisomerase protein provides the binding sites, binding energies, and deactivation constant for the inhibition pocket. Various biochemical methods were used to evaluate the IC50 of the cell lines. The QSAR model was developed for colorectal cell line HCT as a case study. Four QSAR statistical models were predicted between the IC50 of the colorectal cell line HCT to correlate the anticancer activity and the computed physicochemical and quantum chemical global reactivity descriptors. The predictive power of the models indicates a good correlation between the observed and the predicted activity.

Using computers to ESKAPE the antibiotic resistance crisis

Since the discovery of penicillin, the development and use of antibiotics have promoted safe and effective control of bacterial infections. However, the number of antibiotic-resistance cases has been ever increasing over time. Thus, the drug discovery process demands fast, efficient and cost-effective alternative approaches for developing lead candidates with outstanding performance. Computational approaches are appealing techniques to develop lead candidates in an in silico fashion. In this review, we provide an overview of the implementation of current in silico state-of-the-art techniques, including machine learning (ML) and deep learning (DL), in drug discovery. We also discuss the development of quantum computing and its potential benefits for antibiotics research and current bottlenecks that limit computational drug discovery advancement.

A partial least squares and artificial neural network study for a series of arylpiperazines as antidepressant agents

Depression affects more than 300 million people around the world and can lead to suicide. About 30% of patients on treatment for depression drop out of therapy due to side effects or to latency time associated to therapeutic effects. 5-HT receptor, known as serotonin, is considered the key in depression treatment. Arylpiperazine compounds are responsible for several pharmacological effects and are considered as ligands in serotonin receptors, such as the subtype 5-HT_2a. Here, in silico studies were developed using partial least squares (PLSs) and artificial neural networks (ANNs) to design new arylpiperazine compounds that could interact with the 5-HT_2a receptor. First, molecular and electronic descriptors were calculated and posteriorly selected from correlation matrixes and genetic algorithm (GA). Then, the selected descriptors were used to construct PLS and ANN models that showed to be robust and predictive. Lastly, new arylpiperazine compounds were designed and their biological activity values were predicted by both PLS and ANN models. It is worth to highlight compounds G5 and G7 (predicted by the PLS model) and G3 and G15 (predicted by the ANN model), whose predicted pIC₅₀ values were as high as the three highest values from the arylpiperazine original set studied here. Therefore, it can be asserted that the two models (PLS and ANN) proposed in this work are promising for the prediction of the biological activity of new arylpiperazine compounds and may significantly contribute to the design of new drugs for the treatment of depression.

Discovery of New Hits as Antitrypanosomal Agents by In Silico and In Vitro Assays Using Neolignan-Inspired Natural Products from Nectandra leucantha

In the present study, the phytochemical study of the n-hexane extract from flowers of Nectandra leucantha (Lauraceae) afforded six known neolignans (1–6) as well as one new metabolite (7), which were characterized by analysis of NMR, IR, UV, and ESI-HRMS data. The new compound 7 exhibited potent activity against the clinically relevant intracellular forms of T. cruzi (amastigotes), with an IC₅₀ value of 4.3 μM and no observed mammalian cytotoxicity in fibroblasts (CC₅₀ > 200 μM). Based on the results obtained and our previous antitrypanosomal data of 50 natural and semi-synthetic related neolignans, 2D and 3D molecular modeling techniques were employed to help the design of new neolignan-based compounds with higher activity. The results obtained from the models were important to understand the main structural features related to the biological response of the neolignans and to aid in the design of new neolignan-based compounds with better biological activity. Therefore, the results acquired from phytochemical, biological, and in silico studies showed that the integration of experimental and computational techniques consists of a powerful tool for the discovery of new prototypes for development of new drugs to treat CD.

Chemical Graph Theory for Property Modeling in QSAR and QSPR—Charming QSAR & QSPR

Quantitative structure-activity relationship (QSAR) and Quantitative structure-property relationship (QSPR) are mathematical models for the prediction of the chemical, physical or biological properties of chemical compounds. Usually, they are based on structural (grounded on fragment contribution) or calculated (centered on QSAR three-dimensional (QSAR-3D) or chemical descriptors) parameters. Hereby, we describe a Graph Theory approach for generating and mining molecular fragments to be used in QSAR or QSPR modeling based exclusively on fragment contributions. Merging of Molecular Graph Theory, Simplified Molecular Input Line Entry Specification (SMILES) notation, and the connection table data allows a precise way to differentiate and count the molecular fragments. Machine learning strategies generated models with outstanding root mean square error (RMSE) and R2 values. We also present the software Charming QSAR & QSPR, written in Python, for the property prediction of chemical compounds while using this approach.

Density Functional Theory, Chemical Reactivity, Pharmacological Potential and Molecular Docking of Dihydrothiouracil-Indenopyridopyrimidines with Human-DNA Topoisomerase II

In this work, three computational methods (Hatree-Fock (HF), Møller-Plesset 2 (MP2), and Density Functional Theory (DFT)) using a variety of basis sets are used to determine the atomic and molecular properties of dihydrothiouracil-based indenopyridopyrimidine (TUDHIPP) derivatives. Reactivity descriptors of this system, including chemical potential (µ), chemical hardness (η), electrophilicity (ω), condensed Fukui function and dual descriptors are calculated at B3LYP/6-311++ G (d,p) to identify reactivity changes of these molecules in both gas and aqueous phases. We determined the molecular electrostatic surface potential (MESP) to determine the most active site in these molecules. Molecular docking study of TUDHIPP with topoisomerase II α and β is performed, predicting binding sites and binding energies with amino acids of both proteins. Docking studies of TUDHIPP versus etoposide suggest their potential as antitumor candidates. We have applied Lipinski, Veber's rules and analysis of the Golden triangle and structure activity/property relationship for a series of TUDHIPP derivatives indicate that the proposed compounds exhibit good oral bioavailability. The comparison of the drug likeness descriptors of TUDHIPP with those of etoposide, which is known to be an antitumor drug, indicates that TUDHIPP can be considered as an antitumor drug. The overall study indicates that TUDHIPP has comparable and even better descriptors than etoposide proposing that it can be as effective antitumor drug, especially 2H, 6H and 7H compounds.

Two- and three-dimensional QSAR studies on hURAT1 inhibitors with flexible linkers: topomer CoMFA and HQSAR

hURAT1 (human urate transporter 1) is a successful target for hyperuricemia. Recently, the modification work on hURAT1 inhibitors showed that the flexible linkers would benefit biological activity. The study aimed to investigate the contribution of the linkers and give modification strategies on this kind of structures based on QSAR models (HQSAR and topomer CoMFA). The most effective HQSAR and topomer CoMFA models were generated by applying the training set containing 63 compounds, with the cross-validated q² values of 0.869/0.818 and the non-cross-validated correlation coefficients r² of 0.951/0.978, respectively. The Y-randomization test was applied to ensure the robustness of the models. The external predictive correlation coefficient (r_pred²) grounded on the external test set (21 compounds) of two models was 0.910 and 0.907, respectively. In addition, the models were validated by Golbraikh-Tropsha and Roy methods, as well as other statistical metrics. The results showed that both models were reliable. Topomer CoMFA steric/electrostatic contours and HQSAR atomic contribution maps illustrated the structural features which governed their inhibitory potency. The dependable results could provide important insights to guide the designing of more potential hURAT1 inhibitors.

Complementary Quantitative Structure⁻Activity Relationship Models for the Antitrypanosomal Activity of Sesquiterpene Lactones

Three complementary quantitative structure–activity relationship (QSAR) methodologies, namely, regression modeling based on (i) “classical” molecular descriptors, (ii) 3D pharmacophore features, and (iii) 2D molecular holograms (HQSAR) were employed on the antitrypanosomal activity of sesquiterpene lactones (STLs) toward Trypanosoma brucei rhodesiense (Tbr), the causative agent of the East African form of human African trypanosomiasis. In this study, an extension of a previous QSAR study on 69 STLs, models for a much larger and more diverse set of such natural products, now comprising 130 STLs of various structural subclasses, were established. The extended data set comprises a variety of STLs isolated and tested for antitrypanosomal activity within our group and is furthermore enhanced by 12 compounds obtained from literature, which have been tested in the same laboratory under identical conditions. Detailed QSAR analyses yielded models with comparable and good internal and external predictive ability. For a set of compounds as chemically diverse as the one under study, the models exhibited good coefficients of determination (R²) ranging from 0.71 to 0.85, as well as internal (leave-one-out Q² values ranging from 0.62 to 0.72) and external validation coefficients (P² values ranging from 0.54 to 0.73). The contributions of the various tested descriptors to the generated models are in good agreement with the results of previous QSAR studies and corroborate the fact that the antitrypanosomal activity of STLs is very much dependent on the presence and relative position of reactive enone groups within the molecular structure but is influenced by their hydrophilic/hydrophobic properties and molecular shape.

In Silico Studies in Drug Research Against Neurodegenerative Diseases

Background

Neurodegenerative diseases such as Alzheimer's disease (AD), amyotrophic lateral sclerosis, Parkinson's disease (PD), spinal cerebellar ataxias, and spinal and bulbar muscular atrophy are described by slow and selective degeneration of neurons and axons in the central nervous system (CNS) and constitute one of the major challenges of modern medicine. Computer-aided or in silico drug design methods have matured into powerful tools for reducing the number of ligands that should be screened in experimental assays.

Methods

In the present review, the authors provide a basic background about neurodegenerative diseases and in silico techniques in the drug research. Furthermore, they review the various in silico studies reported against various targets in neurodegenerative diseases, including homology modeling, molecular docking, virtual high-throughput screening, quantitative structure activity relationship (QSAR), hologram quantitative structure activity relationship (HQSAR), 3D pharmacophore mapping, proteochemometrics modeling (PCM), fingerprints, fragment-based drug discovery, Monte Carlo simulation, molecular dynamic (MD) simulation, quantum-mechanical methods for drug design, support vector machines, and machine learning approaches.

Results

Detailed analysis of the recently reported case studies revealed that the majority of them use a sequential combination of ligand and structure-based virtual screening techniques, with particular focus on pharmacophore models and the docking approach.

Conclusion

Neurodegenerative diseases have a multifactorial pathoetiological origin, so scientists have become persuaded that a multi-target therapeutic strategy aimed at the simultaneous targeting of multiple proteins (and therefore etiologies) involved in the development of a disease is recommended in future.

The advancement of multidimensional QSAR for novel drug discovery - where are we headed?

INTRODUCTION

The Multidimensional quantitative structure-activity relationship (multidimensional-QSAR) method is one of the most popular computational methods employed to predict interesting biochemical properties of existing or hypothetical molecules. With continuous progress, the QSAR method has made remarkable success in various fields, such as medicinal chemistry, material science and predictive toxicology. Areas covered: In this review, the authors cover the basic elements of multidimensional -QSAR including model construction, validation and application. It includes and emphasizes the very recent developments of multidimensional -QSAR such as: HQSAR, G-QSAR, MIA-QSAR, multi-target QSAR. The advantages and disadvantages of each method are also discussed and typical examples of their application are detailed. Expert opinion: Although there are defects in multidimensional-QSAR modeling, it is still of enormous help to chemists, biologists and other researchers in various fields. In the authors' opinion, the latest more precise and feasible QSAR models should be further developed by integrating new descriptors, algorithms and other relevant computational techniques. Apart from being applied in traditional fields (e.g. lead optimization and predictive risk assessment), QSAR should be used more widely as a routine method in other emerging research fields including the modeling of nanoparticles(NPs), mixture toxicity and peptides.

Combined HQSAR, topomer CoMFA, homology modeling and docking studies on triazole derivatives as SGLT2 inhibitors

Aim: Sodium–glucose cotransporter 2 (SGLT2) is a promising target for diabetes therapy. We aimed to develop computational approaches to identify structural features for more potential SGLT2 inhibitors. Materials & methods: In this work, 46 triazole derivatives as SGLT2 inhibitors were studied using a combination of several approaches, including hologram quantitative structure–activity relationships (HQSAR), topomer comparative molecular field analysis (CoMFA), homology modeling, and molecular docking. HQSAR and topomer CoMFA were used to construct models. Molecular docking was conducted to investigate the interaction of triazole derivatives and homology modeling of SGLT2, as well as to validate the results of the HQSAR and topomer CoMFA models. Results: The most effective HQSAR and topomer CoMFA models exhibited noncross-validated correlation coefficients of 0.928 and 0.891 for the training set, respectively. External predictions were made successfully on a test set and then compared with previously reported models. The graphical results of HQSAR and topomer CoMFA were proven to be consistent with the binding mode of the inhibitors and SGLT2 from molecular docking. Conclusion: The models and docking provided important insights into the design of potent inhibitors for SGLT2.

Studies of Staphylococcus aureus FabI inhibitors: fragment-based approach based on holographic structure-activity relationship analyses

AIM

FabI is a key enzyme in the fatty acid metabolism of Gram-positive bacteria such as Staphylococcus aureus and is an established drug target for known antibiotics such as triclosan. However, due to increasing antibacterial resistance, there is an urgent demand for new drug discovery. Recently, aminopyridine derivatives have been proposed as promising competitive inhibitors of FabI.

METHODS

In the present study, holographic structure-activity relationship (HQSAR) analyses were employed for determining structural contributions of a series containing 105 FabI inhibitors.

RESULTS & CONCLUSION

The final HQSAR model was robust and predictive according to statistical validation (q² and r²_pred equal to 0.696 and 0.854, respectively) and could be further employed to generate fragment contribution maps. Then, final HQSAR model together with FabI active site information can be useful for designing novel bioactive ligands.

Computational strategies to explore antimalarial thiazine alkaloid lead compounds based on an Australian marine sponge Plakortis Lita

In this work, an attempt was made to propose new leads based on the natural scaffold Thiaplakortone-A active against malaria. The 2D QSAR studies suggested that three descriptors correlate with the anti-malarial activity with an R² value of 0.814. Robustness, reliability, and predictive power of the model were tested by internal validation, external validation, Y-scrambling, and applicability domain analysis. HQSAR studies were carried out as an additional tool to find the sub-structural fingerprints. The CoMFA and CoMSIA models gave Q² values of 0.813 and 0.647, and [Formula: see text] values of 0.994 and 0.984, respectively. Using the 2D-QSAR equation, the activity values of the seven modified compounds were calculated and it was found that three molecules showed good anti-malarial activity. Molecular docking of the 42 Thiaplakortone-A derivatives with Plasmodium falciparum calcium-dependent protein kinase 1 (PfCDPK1) was carried out to find out protein-ligand interactions. Data mining of the bioassay data-set AID: 504850 using the classifier based on Random Forest of Weka suggested that all of the eight molecules selected and three out of the seven virtual molecules were anti-malarial active. Both the virtual molecules and drug molecules were docked with CYP3A4, indicating that the virtual molecules could metabolize easily. Toxicity studies using Osiris shows that three molecules showed no toxic characters.

QSPR models for predicting generator-column-derived octanol/water and octanol/air partition coefficients of polychlorinated biphenyls

Octanol/water (K(OW)) and octanol/air (K(OA)) partition coefficients are two important physicochemical properties of organic substances. In current practice, K(OW) and K(OA) values of some polychlorinated biphenyls (PCBs) are measured using generator column method. Quantitative structure-property relationship (QSPR) models can serve as a valuable alternative method of replacing or reducing experimental steps in the determination of K(OW) and K(OA). In this paper, two different methods, i.e., multiple linear regression based on dragon descriptors and hologram quantitative structure-activity relationship, were used to predict generator-column-derived log K(OW) and log K(OA) values of PCBs. The predictive ability of the developed models was validated using a test set, and the performances of all generated models were compared with those of three previously reported models. All results indicated that the proposed models were robust and satisfactory and can thus be used as alternative models for the rapid assessment of the K(OW) and K(OA) of PCBs.

Convergent QSAR studies on a series of NK₃ receptor antagonists for schizophrenia treatment

The dopamine hypothesis states that decreased dopaminergic neurotransmission reduces schizophrenia symptoms. Neurokinin-3 receptor (NK3) antagonists reduce dopamine release and have shown positive effects in pre-clinical and clinical trials. We employed 2D and 3D-QSAR analysis on a series of 40 non-peptide NK3 antagonists. Multivariate statistical analysis, PCA and HCA, were performed to rational training/test set splitting and PLS regression was employed to construct all QSAR models. We constructed one highly predictive CoMFA model (q(2)= 0.810 and r(2)= 0.929) and acceptable HQSAR and CoMSIA models (HQSAR q(2)= 0.644 and r(2)= 0.910; CoMSIA q(2)= 0.691, r(2)= 0.911). The three different techniques provided convergent physicochemical results. All models indicate cyclopropane, piperidine and di-chloro-phenyl ring attached to cyclopropane ring and also the amide group attached to the piperidine ring could play an important role in ligand-receptor interactions. These findings may contribute to develop potential NK3 receptor antagonists for schizophrenia.

Identification of anthranilamide derivatives as potential factor Xa inhibitors: drug design, synthesis and biological evaluation

The coagulation enzyme factor Xa (fXa) plays a crucial role in the blood coagulation cascade. In this study, three-dimensional fragment based drug design (FBDD) combined with structure-based pharmacophore (SBP) model and structural consensus docking were employed to identify novel fXa inhibitors. After a multi-stage virtual screening (VS) workflow, two hit compounds 3780 and 319 having persistent high performance were identified. Then, these two hit compounds and several analogs were synthesized and screened for in-vitro inhibition of fXa. The experimental data showed that most of the designed compounds displayed significant in vitro potency against fXa. Among them, compound 9b displayed the greatest in vitro potency against fXa with the IC50 value of 23 nM and excellent selectivity versus thrombin (IC50 = 40 μM). Moreover, the prolongation of the prothrombin time (PT) was measured for compound 9b to evaluate its in vitro anticoagulant activity. As a result, compound 9b exhibited pronounced anticoagulant activity with the 2 × PT value of 8.7 μM.

Ligand biological activity predictions using fingerprint-based artificial neural networks (FANN-QSAR)

This chapter focuses on the fingerprint-based artificial neural networks QSAR (FANN-QSAR) approach to predict biological activities of structurally diverse compounds. Three types of fingerprints, namely ECFP6, FP2, and MACCS, were used as inputs to train the FANN-QSAR models. The results were benchmarked against known 2D and 3D QSAR methods, and the derived models were used to predict cannabinoid (CB) ligand binding activities as a case study. In addition, the FANN-QSAR model was used as a virtual screening tool to search a large NCI compound database for lead cannabinoid compounds. We discovered several compounds with good CB2 binding affinities ranging from 6.70 nM to 3.75 μM. The studies proved that the FANN-QSAR method is a useful approach to predict bioactivities or properties of ligands and to find novel lead compounds for drug discovery research.

Hologram QSAR studies of antiprotozoal activities of sesquiterpene lactones

Infectious diseases such as trypanosomiasis and leishmaniasis are considered neglected tropical diseases due the lack for many years of research and development into new drug treatments besides the high incidence of mortality and the lack of current safe and effective drug therapies. Natural products such as sesquiterpene lactones have shown activity against T. brucei and L. donovani, the parasites responsible for these neglected diseases. To evaluate structure activity relationships, HQSAR models were constructed to relate a series of 40 sesquiterpene lactones (STLs) with activity against T. brucei, T. cruzi, L. donovani and P. falciparum and also with their cytotoxicity. All constructed models showed good internal (leave-one-out q² values ranging from 0.637 to 0.775) and external validation coefficients (r²_test values ranging from 0.653 to 0.944). From HQSAR contribution maps, several differences between the most and least potent compounds were found. The fragment contribution of PLS-generated models confirmed the results of previous QSAR studies that the presence of α,β-unsatured carbonyl groups is fundamental to biological activity. QSAR models for the activity of these compounds against T. cruzi, L. donovani and P. falciparum are reported here for the first time. The constructed HQSAR models are suitable to predict the activity of untested STLs.

Flavonoids promoting HaCaT migration: I. Hologram quantitative structure-activity relationships

Cell migration plays an important role in multicellular development and preservation. Because wound healing requires cell migration, compounds promoting cell migration can be used for wound repair therapy. Several plant-derived polyphenols are known to promote cell migration, which improves wound healing. Previous studies of flavonoids on cell lines have focused on their inhibitory effects and not on wound healing. In addition, studies of flavonoids on wound healing have been performed using mixtures. In this study, individual flavonoids were used for cellular migration measurements. Relationships between the cell migration effects of flavonoids and their structural properties have never been reported. Here, we investigated the migration of keratinocytes caused by 100 flavonoids and examined their relationships using hologram quantitative structure-activity relationships. The structural conditions responsible for efficient cell migration on keratinocyte cell lines determined from the current study will facilitate the design of flavonoids with improved activity.

Combined 2D and 3D-QSAR, molecular modelling and docking studies of pyrazolodiazepinones as novel phosphodiesterase 2 inhibitors

Selective inhibition of phosphodiesterase 2 (PDE2) in cells where it is located elevates cyclic guanosine monophosphate (cGMP) and acts as novel analgesic with antinociceptive activity. Three-dimensional quantitative structure-activity relationship (QSAR) studies for pyrazolodiazepinone inhibitors exhibiting PDE2 inhibition were performed using comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) and Topomer CoMFA, and two-dimensional QSAR study was performed using a Hologram QSAR (HQSAR) method. QSAR models were generated using training set of 23 compounds and were validated using test set of nine compounds. The optimum partial least squares (PLS) for CoMFA-Focusing, CoMSIA-SDH, Topomer CoMFA and HQSAR models exhibited good 'leave-one-out' cross validated correlation coefficient (q(2)) of 0.790, 0.769, 0.840 and 0.787, coefficient of determination (r(2)) of 0.999, 0.964, 0.979 and 0.980, and high predictive power (r(2)(pred)) of 0.796, 0.833, 0.820 and 0.803 respectively. Docking studies revealed that those inhibitors able to bind to amino acid Gln859 by cGMP binding orientation called 'glutamine-switch', and also bind to the hydrophobic clamp of PDE2 isoform, could possess high selectivity for PDE2. From the results of all the studies, structure-activity relationships and structural requirements for binding to active site of PDE2 were established which provide useful guidance for the design and future synthesis of potent PDE2 inhibitors.

Holographic quantitative structure-activity relationships of tryptamine derivatives at NMDA, 5HT(1A) and 5HT(2A) receptors

Tryptamine derivatives (Ts) were found to inhibit the binding of [³H]MK-801, [³H]ketanserin and [³H]8-OH-DPAT to rat brain membranes. [³H]MK-801 labels the NMDA (N-methyl-D-aspartate) receptor, a ionotropic glutamate receptor which controls synaptic plasticity and memory function in the brain, whereas [³H]ketanserin and [³H]8-OH-DPAT label 5HT_2A and 5HT_1A receptors, respectively. The inhibitory potencies of 64 Ts (as given by IC₅₀ values) were correlated with their structural properties by using the Holographic QSAR procedure (HQSAR). This method uses structural fragments and connectivities as descriptors which were encoded in a hologram thus avoiding the usual problems with conformation and alignment of the structures. Four correlation equations with high predictive ability and appropriate statistical test values could be established. The results are visualized by generation of maps reflecting the contribution of individual structural parts to the biological activities.

Computational studies of TGF-βRI (ALK-5) inhibitors: analysis of the binding interactions between ligand-receptor using 2D and 3D techniques

ALK-5 (Activin-Like Kinase 5) is a biological receptor involved in a variety of pathological processes such as cancer and fibrosis. ALK-5 receptor propagates an intracellular signaling that forms a protein complex capable of reaching the nucleus and modulating the gene transcription. In the present study, comparative molecular field analysis (CoMFA) and hologram quantitative structure-activity relationship (HQSAR) studies were conducted on a series of potent ALK-5 inhibitors. Significant correlation coefficients (CoMFA, r(2)=0.99 and q(2)=0.85; HQSAR, r(2)=0.92 and q(2)=0.72) were obtained, indicating the predictive potential of the 2D and 3D models for untested compounds. The models were then used to predict the potency of a test set, and the predicted values from the HQSAR and CoMFA models were in good agreement with the experimental results. The final QSAR models, along with the information obtained from 3D (steric and electrostatic) contour maps and 2D contribution maps, can be useful for the design of novel bioactive ligands.

Molecular fingerprint-based artificial neural networks QSAR for ligand biological activity predictions

In this manuscript, we have reported a novel 2D fingerprint-based artificial neural network QSAR (FANN-QSAR) method in order to effectively predict biological activities of structurally diverse chemical ligands. Three different types of fingerprints, namely, ECFP6, FP2 and MACCS, were used in FANN-QSAR algorithm development, and FANN-QSAR models were compared to known 3D and 2D QSAR methods using five data sets previously reported. In addition, the derived models were used to predict GPCR cannabinoid ligand binding affinities using our manually curated cannabinoid ligand database containing 1699 structurally diverse compounds with reported cannabinoid receptor subtype CB(2) activities. To demonstrate its useful applications, the established FANN-QSAR algorithm was used as a virtual screening tool to search a large NCI compound database for lead cannabinoid compounds, and we have discovered several compounds with good CB(2) binding affinities ranging from 6.70 nM to 3.75 μM. To the best of our knowledge, this is the first report for a fingerprint-based neural network approach validated with a successful virtual screening application in identifying lead compounds. The studies proved that the FANN-QSAR method is a useful approach to predict bioactivities or properties of ligands and to find novel lead compounds for drug discovery research.

Toward rational fragment-based lead design without 3D structures

Fragment-based lead discovery (FBLD) has become a prime component of the armamentarium of modern drug design programs. FBLD identifies low molecular weight ligands that weakly bind to important biological targets. Three-dimensional structural information about the binding mode is provided by X-ray crystallography or NMR spectroscopy and is subsequently used to improve the lead compounds. Despite tremendous success rates, FBLD relies on the availability of high-resolution structural information, still a bottleneck in drug discovery programs. To overcome these limitations, we recently demonstrated that the meta-structure approach provides an alternative route to rational lead identification in cases where no 3D structure information about the biological target is available. Combined with information-rich NMR data, this strategy provides valuable information for lead development programs. We demonstrate with several examples the feasibility of the combined NMR and meta-structure approach to devise a rational strategy for fragment evolution without resorting to highly resolved protein complex structures.

Structure- and ligand-based drug design approaches for neglected tropical diseases

Drug discovery has moved toward more rational strategies based on our increasing understanding of the fundamental principles of protein–ligand interactions. Structure- (SBDD) and ligand-based drug design (LBDD) approaches bring together the most powerful concepts in modern chemistry and biology, linking medicinal chemistry with structural biology. The definition and assessment of both chemical and biological space have revitalized the importance of exploring the intrinsic complementary nature of experimental and computational methods in drug design. Major challenges in this field include the identification of promising hits and the development of high-quality leads for further development into clinical candidates. It becomes particularly important in the case of neglected tropical diseases (NTDs) that affect disproportionately poor people living in rural and remote regions worldwide, and for which there is an insufficient number of new chemical entities being evaluated owing to the lack of innovation and R&D investment by the pharmaceutical industry. This perspective paper outlines the utility and applications of SBDD and LBDD approaches for the identification and design of new small-molecule agents for NTDs.

Fragment-similarity-based QSAR (FS-QSAR) algorithm for ligand biological activity predictions

Quantitative structure-activity relationship (QSAR) studies are useful computational tools often used in drug discovery research and in many scientific disciplines. In this study, a robust fragment-similarity-based QSAR (FS-QSAR) algorithm was developed to correlate structures with biological activities by integrating fragment-based drug design concept and a multiple linear regression method. Similarity between any pair of training and testing fragments was determined by calculating the difference of lowest or highest eigenvalues of the chemistry space BCUT matrices of corresponding fragments. In addition to the BCUT-similarity function, molecular fingerprint Tanimoto coefficient (Tc) similarity function was also used as an alternative for comparison. For validation studies, the FS-QSAR algorithm was applied to several case studies, including a dataset of COX2 inhibitors and a dataset of cannabinoid CB2 triaryl bis-sulfone antagonist analogues, to build predictive models achieving average coefficient of determination (r(2)) of 0.62 and 0.68, respectively. The developed FS-QSAR method is proved to give more accurate predictions than the traditional and one-nearest-neighbour QSAR methods and can be a useful tool in the fragment-based drug discovery for ligand activity prediction.

4-aminoquinolines active against chloroquine-resistant Plasmodium falciparum: basis of antiparasite activity and quantitative structure-activity relationship analyses

Chloroquine (CQ) is a safe and economical 4-aminoquinoline (AQ) antimalarial. However, its value has been severely compromised by the increasing prevalence of CQ resistance. This study examined 108 AQs, including 68 newly synthesized compounds. Of these 108 AQs, 32 (30%) were active only against CQ-susceptible Plasmodium falciparum strains and 59 (55%) were active against both CQ-susceptible and CQ-resistant P. falciparum strains (50% inhibitory concentrations [IC50s], ≤25 nM). All AQs active against both CQ-susceptible and CQ-resistant P. falciparum strains shared four structural features: (i) an AQ ring without alkyl substitution, (ii) a halogen at position 7 (Cl, Br, or I but not F), (iii) a protonatable nitrogen at position 1, and (iv) a second protonatable nitrogen at the end of the side chain distal from the point of attachment to the AQ ring via the nitrogen at position 4. For activity against CQ-resistant parasites, side chain lengths of ≤3 or ≥10 carbons were necessary but not sufficient; they were identified as essential factors by visual comparison of 2-dimensional (2-D) structures in relation to the antiparasite activities of the AQs and were confirmed by computer-based 3-D comparisons and differential contour plots of activity against P. falciparum. The advantage of the method reported here (refinement of quantitative structure-activity relationship [QSAR] descriptors by random assignment of compounds to multiple training and test sets) is that it retains QSAR descriptors according to their abilities to predict the activities of unknown test compounds rather than according to how well they fit the activities of the compounds in the training sets.

Computational medicinal chemistry in fragment-based drug discovery: what, how and when

The use of fragment-based drug discovery (FBDD) has increased in the last decade due to the encouraging results obtained to date. In this scenario, computational approaches, together with experimental information, play an important role to guide and speed up the process. By default, FBDD is generally considered as a constructive approach. However, such additive behavior is not always present, therefore, simple fragment maturation will not always deliver the expected results. In this review, computational approaches utilized in FBDD are reported together with real case studies, where applicability domains are exemplified, in order to analyze them, and then, maximize their performance and reliability. Thus, a proper use of these computational tools can minimize misleading conclusions, keeping the credit on FBDD strategy, as well as achieve higher impact in the drug-discovery process. FBDD goes one step beyond a simple constructive approach. A broad set of computational tools: docking, R group quantitative structure–activity relationship, fragmentation tools, fragments management tools, patents analysis and fragment-hopping, for example, can be utilized in FBDD, providing a clear positive impact if they are utilized in the proper scenario – what, how and when. An initial assessment of additive/non-additive behavior is a critical point to define the most convenient approach for fragments elaboration.

Ligand-based virtual screening procedure for the prediction and the identification of novel β-amyloid aggregation inhibitors using Kohonen maps and Counterpropagation Artificial Neural Networks

In this work we have developed an in silico model to predict the inhibition of β-amyloid aggregation by small organic molecules. In particular we have explored the inhibitory activity of a series of 62 N-phenylanthranilic acids using Kohonen maps and Counterpropagation Artificial Neural Networks. The effects of various structural modifications on biological activity are investigated and novel structures are designed using the developed in silico model. More specifically a search for optimized pharmacophore patterns by insertions, substitutions, and ring fusions of pharmacophoric substituents of the main building block scaffolds is described. The detection of the domain of applicability defines compounds whose estimations can be accepted with confidence.

Recent advances in fragment-based QSAR and multi-dimensional QSAR methods

This paper provides an overview of recently developed two dimensional (2D) fragment-based QSAR methods as well as other multi-dimensional approaches. In particular, we present recent fragment-based QSAR methods such as fragment-similarity-based QSAR (FS-QSAR), fragment-based QSAR (FB-QSAR), Hologram QSAR (HQSAR), and top priority fragment QSAR in addition to 3D- and nD-QSAR methods such as comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA), Topomer CoMFA, self-organizing molecular field analysis (SOMFA), comparative molecular moment analysis (COMMA), autocorrelation of molecular surfaces properties (AMSP), weighted holistic invariant molecular (WHIM) descriptor-based QSAR (WHIM), grid-independent descriptors (GRIND)-based QSAR, 4D-QSAR, 5D-QSAR and 6D-QSAR methods.

4D-QSAR: perspectives in drug design

Drug design is a process driven by innovation and technological breakthroughs involving a combination of advanced experimental and computational methods. A broad variety of medicinal chemistry approaches can be used for the identification of hits, generation of leads, as well as to accelerate the optimization of leads into drug candidates. The quantitative structure–activity relationship (QSAR) formalisms are among the most important strategies that can be applied for the successful design new molecules. This review provides a comprehensive review on the evolution and current status of 4D-QSAR, highlighting present challenges and new opportunities in drug design.

Identification of Metabotropic Glutamate Receptor Subtype 5 Potentiators Using Virtual High-Throughput Screening

Selective potentiators of glutamate response at metabotropic glutamate receptor subtype 5 (mGluR5) have exciting potential for the development of novel treatment strategies for schizophrenia. A total of 1,382 compounds with positive allosteric modulation (PAM) of the mGluR5 glutamate response were identified through high-throughput screening (HTS) of a diverse library of 144,475 substances utilizing a functional assay measuring receptor-induced intracellular release of calcium. Primary hits were tested for concentration-dependent activity, and potency data (EC₅₀ values) were used for training artificial neural network (ANN) quantitative structure−activity relationship (QSAR) models that predict biological potency from the chemical structure. While all models were trained to predict EC₅₀, the quality of the models was assessed by using both continuous measures and binary classification. Numerical descriptors of chemical structure were used as input for the machine learning procedure and optimized in an iterative protocol. The ANN models achieved theoretical enrichment ratios of up to 38 for an independent data set not used in training the model. A database of ∼450,000 commercially available drug-like compounds was targeted in a virtual screen. A set of 824 compounds was obtained for testing based on the highest predicted potency values. Biological testing found 28.2% (232/824) of these compounds with various activities at mGluR5 including 177 pure potentiators and 55 partial agonists. These results represent an enrichment factor of 23 for pure potentiation of the mGluR5 glutamate response and 30 for overall mGluR5 modulation activity when compared with those of the original mGluR5 experimental screening data (0.94% hit rate). The active compounds identified contained 72% close derivatives of previously identified PAMs as well as 28% nontrivial derivatives of known active compounds.