3-D QSAutogrid/R: an alternative procedure to build 3-D QSAR models. Methodologies and applications

Molecular Docking Assessment of Cathinones as 5-HT2AR Ligands: Developing of Predictive Structure-Based Bioactive Conformations and Three-Dimensional Structure-Activity Relationships Models for Future Recognition of Abuse Drugs

Commercially available cathinones are drugs of long-term abuse drugs whose pharmacology is fairly well understood. While their psychedelic effects are associated with 5-HT2AR, the enclosed study summarizes efforts to shed light on the pharmacodynamic profiles, not yet known at the receptor level, using molecular docking and three-dimensional quantitative structure-activity relationship (3-D QSAR) studies. The bioactive conformations of cathinones were modeled by AutoDock Vina and were used to build structure-based (SB) 3-D QSAR models using the Open3DQSAR engine. Graphical inspection of the results led to the depiction of a 3-D structure analysis-activity relationship (SAR) scheme that could be used as a guideline for molecular determinants by which any untested cathinone molecule can be predicted as a potential 5-HT2AR binder prior to experimental evaluation. The obtained models, which showed a good agreement with the chemical properties of co-crystallized 5-HT2AR ligands, proved to be valuable for future virtual screening campaigns to recognize unused cathinones and similar compounds, such as 5-HT2AR ligands, minimizing both time and financial resources for the characterization of their psychedelic effects.

Bovine Serum Amine Oxidase and Polyamine Analogues: Chemical Synthesis and Biological Evaluation Integrated with Molecular Docking and 3-D QSAR Studies

Natural polyamines (PAs) are key players in cellular homeostasis by regulating cell growth and proliferation. Several observations highlight that PAs are also implicated in pathways regulating cell death. Indeed, the PA accumulation cytotoxic effect, maximized with the use of bovine serum amine oxidase (BSAO) enzyme, represents a valuable strategy against tumor progression. In the present study, along with the design, synthesis, and biological evaluation of a series of new spermine (Spm) analogues (1-23), a mixed structure-based (SB) and ligand-based (LB) protocol was applied. Binding modes of BSAO-PA modeled complexes led to clarify electrostatic and steric features likely affecting the BSAO-PA biochemical kinetics. LB and SB three-dimensional quantitative structure-activity relationship (Py-CoMFA and Py-ComBinE) models were developed by means of the 3d-qsar.com portal, and their analysis represents a strong basis for future design and synthesis of PA BSAO substrates for potential application in oxidative stress-induced chemotherapy.

Human Estrogen Receptor α Antagonists. Part 1: 3-D QSAR-Driven Rational Design of Innovative Coumarin-Related Antiestrogens as Breast Cancer Suppressants through Structure-Based and Ligand-Based Studies

The estrogen receptor α (ERα) represents a 17β-estradiol-inducible transcriptional regulator that initiates the RNA polymerase II-dependent transcriptional machinery, pointed for breast cancer (BC) development via either genomic direct or genomic indirect (i.e., tethered) pathway. To develop innovative ligands, structure-based (SB) three-dimensional (3-D) quantitative structure-activity relationship (QSAR) studies have been undertaken from structural data taken from partial agonists, mixed agonists/antagonists (selective estrogen receptor modulators (SERMs)), and full antagonists (selective ERα downregulators (SERDs)) correlated with either wild-type or mutated ERα receptors. SB and ligand-based (LB) alignments allow us to rule out guidelines for the SB/LB alignment of untested compounds. 3-D QSAR models for ERα ligands, coupled with SB/LB alignment, were revealed to be useful tools to dissect the chemical determinants for ERα-based anticancer activity as well as to predict their potency. The herein developed protocol procedure was verified through the design and potency prediction of 12 new coumarin-based SERMs, namely, 3DQ-1a to 3DQ-1e, that upon synthesis turned to be potent ERα antagonists by means of either in vitro or in vivo assays (described in the second part of this study).

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.

Rational Drug Design of Antineoplastic Agents Using 3D-QSAR, Cheminformatic, and Virtual Screening Approaches

Background:

Computer-Aided Drug Design has strongly accelerated the development of novel antineoplastic agents by helping in the hit identification, optimization, and evaluation.

Results:

Computational approaches such as cheminformatic search, virtual screening, pharmacophore modeling, molecular docking and dynamics have been developed and applied to explain the activity of bioactive molecules, design novel agents, increase the success rate of drug research, and decrease the total costs of drug discovery. Similarity, searches and virtual screening are used to identify molecules with an increased probability to interact with drug targets of interest, while the other computational approaches are applied for the design and evaluation of molecules with enhanced activity and improved safety profile.

Conclusion:

In this review are described the main in silico techniques used in rational drug design of antineoplastic agents and presented optimal combinations of computational methods for design of more efficient antineoplastic drugs.

Design of novel amyloid β aggregation inhibitors using QSAR, pharmacophore modeling, molecular docking and ADME prediction

The inhibition of abnormal amyloid β (Aβ) aggregation has been regarded as a good target to control Alzheimer's disease. The present study adopted 2D-QSAR, HQSAR and 3D QSAR (CoMFA & CoMSIA) modeling approaches to identify the structural and physicochemical requirements for the potential Aβ aggregation inhibition. A structure-based molecular docking technique is utilized to approve the features that are obtained from the ligand-based techniques on 30 curcumin derivatives. The combined outputs were then used to screen the modified 10 compounds. The 2D QSAR model on curcumin derivatives gave statistical values R2 = 0.9086 and SEE = 0.1837. The model was further confirmed by Y-randomization test and Applicability domain analysis by the standardization approach. The HQSAR study (Q2 = 0.615, Rncv 2 = 0.931, Rpred 2 = 0.956) illustrated the important molecular fingerprints for inhibition. Contour maps of 3D QSAR models, CoMFA (Q2 = 0.687, Rncv 2 = 0.787, Rpred 2 = 0.731) and CoMSIA (Q2 = 0.743, Rncv 2 = 0.972, Rpred 2 = 0.713), depict that the models are robust and provide explanation of the important features, like steric, electrostatic and hydrogen bond acceptor, which play important role for interaction with the receptor site cavity. The molecular docking study of the curcumin derivatives elucidates the important interactions between the amino acid residues at the catalytic site of the receptor and the ligands, indicating the structural requirements of the inhibitors. The ligand-receptor interactions of top hits were analyzed to explore the pharmacophore features of Aβ aggregation inhibition. The Aβ aggregation inhibitory activities of novel chemical entities were then obtained through inverse QSAR. The newly designed molecules were further screened through machine learning, prediction of toxicity and nature of metabolism to get the proposed six lead compounds.

Disruptor of telomeric silencing 1-like (DOT1L): disclosing a new class of non-nucleoside inhibitors by means of ligand-based and structure-based approaches

Chemical inhibition of chromatin-mediated signaling involved proteins is an established strategy to drive expression networks and alter disease progression. Protein methyltransferases are among the most studied proteins in epigenetics and, in particular, disruptor of telomeric silencing 1-like (DOT1L) lysine methyltransferase plays a key role in MLL-rearranged acute leukemia Selective inhibition of DOT1L is an established attractive strategy to breakdown aberrant H3K79 methylation and thus overexpression of leukemia genes, and leukemogenesis. Although numerous DOT1L inhibitors have been several structural data published no pronounced computational efforts have been yet reported. In these studies a first tentative of multi-stage and LB/SB combined approach is reported in order to maximize the use of available data. Using co-crystallized ligand/DOT1L complexes, predictive 3-D QSAR and COMBINE models were built through a python implementation of previously reported methodologies. The models, validated by either modeled or experimental external test sets, proved to have good predictive abilities. The application of these models to an internal library led to the selection of two unreported compounds that were found able to inhibit DOT1L at micromolar level. To the best of our knowledge this is the first report of quantitative LB and SB DOT1L inhibitors models and their application to disclose new potential epigenetic modulators.

Protein-Ligand Docking in Drug Design: Performance Assessment and Binding-Pose Selection

Main goal in drug discovery is the identification of drug-like compounds capable to modulate specific biological targets. Thus, the prediction of reliable binding poses of candidate ligands, through molecular docking simulations, represents a key step to be pursued in structure-based drug design (SBDD). Since the increasing number of resolved three-dimensional ligand-protein structures, together with the expansion of computational power and software development, the comprehensive and systematic use of experimental data can be proficiently employed to validate the docking performance. This allows to select and refine the protocol to adopt when predicting the binding pose of trial compounds in a target. Given the availability of multiple docking software, a comparative docking assessment in an early research stage represents a must-use step to minimize fails in molecular modeling. This chapter describes how to perform a docking assessment, using freely available tools, in a semiautomated fashion.

Synthesis, biological evaluation and quantitative structure-active relationships of 1,3-thiazolidin-4-one derivatives. A promising chemical scaffold endowed with high antifungal potency and low cytotoxicity

With reference to recent studies reporting on the various biological properties of the thiazolidinone scaffold, we synthesized more than a hundred compounds characterized by a 1,3-thiazolidin-4-one nucleus derivatised at the C2 with a hydrazine bridge linked to (cyclo)aliphatic or hetero(aryl) moieties, and their N-benzylated derivatives. These molecules were assayed as potential anti-Candida agents and they were shown to possess comparable, and in some cases higher biological activity than well-established topical and systemic antimycotic drugs (i.e. clotrimazole, fluconazole, ketoconazole, miconazole, tioconazole, amphotericin B). Compounds endowed with the lowest MICs underwent further testing in order to assess their cytotoxic effect (CC₅₀) on Hep2 cells, which demonstrated their relative safety. Finally, QSAR and 3-D QSAR models were used to predict putative chemical modifications of the 1,3-thiazolidin-4-one scaffold in order to design new and potential more active compounds against Candida spp.

New Inhibitors of Indoleamine 2,3-Dioxygenase 1: Molecular Modeling Studies, Synthesis, and Biological Evaluation

Indoleamine 2,3-dioxygenase 1 (IDO1) is an attractive target for anticancer therapy. Herein, we report a virtual screening study which led to the identification of compound 5 as a new IDO1 inhibitor. In order to improve the biological activity of the identified hit, arylthioindoles 6-30 were synthesized and tested. Among these, derivative 21 exhibited an IC₅₀ value of 7 μM, being the most active compound of the series. Furthermore, compounds 5 and 21 induced a dose-dependent growth inhibition in IDO1 expressing cancer cell lines HTC116 and HT29. Three-dimensional quantitative structure-activity relationship studies were carried out in order to rationalize obtained results and suggest new chemical modifications.

Chemical Structure-Biological Activity Models for Pharmacophores' 3D-Interactions

Within medicinal chemistry nowadays, the so-called pharmaco-dynamics seeks for qualitative (for understanding) and quantitative (for predicting) mechanisms/models by which given chemical structure or series of congeners actively act on biological sites either by focused interaction/therapy or by diffuse/hazardous influence. To this aim, the present review exposes three of the fertile directions in approaching the biological activity by chemical structural causes: the special computing trace of the algebraic structure-activity relationship (SPECTRAL-SAR) offering the full analytical counterpart for multi-variate computational regression, the minimal topological difference (MTD) as the revived precursor for comparative molecular field analyses (CoMFA) and comparative molecular similarity indices analysis (CoMSIA); all of these methods and algorithms were presented, discussed and exemplified on relevant chemical medicinal systems as proton pump inhibitors belonging to the 4-indolyl,2-guanidinothiazole class of derivatives blocking the acid secretion from parietal cells in the stomach, the 1-[(2-hydroxyethoxy)-methyl]-6-(phenylthio)thymine congeners’ (HEPT ligands) antiviral activity against Human Immunodeficiency Virus of first type (HIV-1) and new pharmacophores in treating severe genetic disorders (like depression and psychosis), respectively, all involving 3D pharmacophore interactions.

Design and Synthesis of Simplified Largazole Analogues as Isoform-Selective Human Lysine Deacetylase Inhibitors

Selective inhibition of KDAC isoforms while maintaining potency remains a challenge. Using the largazole macrocyclic depsipeptide structure as a starting point for developing new KDACIs with increased selectivity, a combination of four different simplified largazole analogue (SLA) scaffolds with diverse zinc-binding groups (for a total of 60 compounds) were designed, synthesized, and evaluated against class I KDACs 1, 3, and 8, and class II KDAC6. Experimental evidence as well as molecular docking poses converged to establish the cyclic tetrapeptides (CTPs) as the primary determinant of both potency and selectivity by influencing the correct alignment of the zinc-binding group in the KDAC active site, providing a further basis for developing new KDACIs of higher isoform selectivity and potency.

An Automated Strategy for Binding-Pose Selection and Docking Assessment in Structure-Based Drug Design

Molecular docking is a widely used technique in drug design to predict the binding pose of a candidate compound in a defined therapeutic target. Numerous docking protocols are available, each characterized by different search methods and scoring functions, thus providing variable predictive capability on a same ligand-protein system. To validate a docking protocol, it is necessary to determine a priori the ability to reproduce the experimental binding pose (i.e., by determining the docking accuracy (DA)) in order to select the most appropriate docking procedure and thus estimate the rate of success in docking novel compounds. As common docking programs use generally different root-mean-square deviation (RMSD) formulas, scoring functions, and format results, it is both difficult and time-consuming to consistently determine and compare their predictive capabilities in order to identify the best protocol to use for the target of interest and to extrapolate the binding poses (i.e., best-docked (BD), best-cluster (BC), and best-fit (BF) poses) when applying a given docking program over thousands/millions of molecules during virtual screening. To reduce this difficulty, two new procedures called Clusterizer and DockAccessor have been developed and implemented for use with some common and "free-for-academics" programs such as AutoDock4, AutoDock4(Zn), AutoDock Vina, DOCK, MpSDockZn, PLANTS, and Surflex-Dock to automatically extrapolate BD, BC, and BF poses as well as to perform consistent cluster and DA analyses. Clusterizer and DockAccessor (code available over the Internet) represent two novel tools to collect computationally determined poses and detect the most predictive docking approach. Herein an application to human lysine deacetylase (hKDAC) inhibitors is illustrated.

Vascular endothelial growth factor receptor-2 (VEGFR-2) inhibitors: development and validation of predictive 3-D QSAR models through extensive ligand- and structure-based approaches

Vascular endothelial growth factor receptor-2, (VEGFR-2), is a key element in angiogenesis, the process by which new blood vessels are formed, and is thus an important pharmaceutical target. Here, 3-D quantitative structure-activity relationship (3-D QSAR) were used to build a quantitative screening and pharmacophore model of the VEGFR-2 receptors for design of inhibitors with improved activities. Most of available experimental data information has been used as training set to derive optimized and fully cross-validated eight mono-probe and a multi-probe quantitative models. Notable is the use of 262 molecules, aligned following both structure-based and ligand-based protocols, as external test set confirming the 3-D QSAR models' predictive capability and their usefulness in design new VEGFR-2 inhibitors. From a survey on literature, this is the first generation of a wide-ranging computational medicinal chemistry application on VEGFR2 inhibitors.

Hsp90 inhibitors, part 1: definition of 3-D QSAutogrid/R models as a tool for virtual screening

The multichaperone heat shock protein (Hsp) 90 complex mediates the maturation and stability of a variety of oncogenic signaling proteins. For this reason, Hsp90 has emerged as a promising target for anticancer drug development. Herein, we describe a complete computational procedure for building several 3-D QSAR models used as a ligand-based (LB) component of a comprehensive ligand-based (LB) and structure-based (SB) virtual screening (VS) protocol to identify novel molecular scaffolds of Hsp90 inhibitors. By the application of the 3-D QSAutogrid/R method, eight SB PLS 3-D QSAR models were generated, leading to a final multiprobe (MP) 3-D QSAR pharmacophoric model capable of recognizing the most significant chemical features for Hsp90 inhibition. Both the monoprobe and multiprobe models were optimized, cross-validated, and tested against an external test set. The obtained statistical results confirmed the models as robust and predictive to be used in a subsequent VS.

Pharmacophore assessment through 3-D QSAR: evaluation of the predictive ability on new derivatives by the application on a series of antitubercular agents

Pharmacophoric mapping is a useful procedure to frame, especially when crystallographic receptor structures are unavailable as in ligand-based studies, the hypothetical site of interaction. In this study, 71 pyrrole derivatives active against M. tuberculosis were used to derive through a recent new 3-D QSAR protocol, 3-D QSAutogrid/R, several predictive 3-D QSAR models on compounds aligned by a previously reported pharmacophoric application. A final multiprobe (MP) 3-D QSAR model was then obtained configuring itself as a tool to derive pharmacophoric quantitative models. To stress the applicability of the described models, an external test set of unrelated and newly synthesized series of R-4-amino-3-isoxazolidinone derivatives found to be active at micromolar level against M. tuberculosis was used, and the predicted bioactivities were in good agreement with the experimental values. The 3-D QSAutogrid/R procedure proved to be able to correlate by a single multi-informative scenario the different activity molecular profiles thus confirming its usefulness in the rational drug design approach.

Pharmacophore and 3D-QSAR characterization of 6-arylquinazolin-4-amines as Cdc2-like kinase 4 (Clk4) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (Dyrk1A) inhibitors

Cdc2-like kinase 4 (Clk4) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (Dyrk1A) are protein kinases that are promising targets for treatment of diseases caused by abnormal gene splicing. 6-Arylquinazolin-4-amines have been recently identified as potent Clk4 and Dyrk1A inhibitors. In order to understand the structure-activity correlation of these analogs, we have applied ligand-based pharmacophore and 3D-QSAR modeling combined with structure-based homology modeling and docking. The high R(2) and Q(2) (0.88 and 0.79 for Clk4, 0.85 and 0.82 for Dyrk1A, respectively) based on validation with training and test set compounds suggested that the generated 3D-QSAR models are reliable in predicting novel ligand activities against Clk4 and Dyrk1A. The binding mode identified through docking ligands to the ATP binding domain of Clk4 was consistent with the structural properties and energy field contour maps characterized by pharmacophore and 3D-QSAR models and gave valuable insights into the structure-activity profile of 6-arylquinazolin-4-amine analogs. The obtained 3D-QSAR and pharmacophore models in combination with the binding mode between inhibitor and residues of Clk4 will be helpful for future lead compound identification and optimization to design potent and selective Clk4 and Dyrk1A inhibitors.

Pharmacophore and 3D-QSAR characterization of 6-arylquinazolin-4-amines as Cdc2-like kinase 4 (Clk4) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (Dyrk1A) inhibitors

Cdc2-like kinase 4 (Clk4) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (Dyrk1A) are protein kinases that are promising targets for treatment of diseases caused by abnormal gene splicing. 6-Arylquinazolin-4-amines have been recently identified as potent Clk4 and Dyrk1A inhibitors. In order to understand the structure-activity correlation of these analogs, we have applied ligand-based pharmacophore and 3D-QSAR modeling combined with structure-based homology modeling and docking. The high R(2) and Q(2) (0.88 and 0.79 for Clk4, 0.85 and 0.82 for Dyrk1A, respectively) based on validation with training and test set compounds suggested that the generated 3D-QSAR models are reliable in predicting novel ligand activities against Clk4 and Dyrk1A. The binding mode identified through docking ligands to the ATP binding domain of Clk4 was consistent with the structural properties and energy field contour maps characterized by pharmacophore and 3D-QSAR models and gave valuable insights into the structure-activity profile of 6-arylquinazolin-4-amine analogs. The obtained 3D-QSAR and pharmacophore models in combination with the binding mode between inhibitor and residues of Clk4 will be helpful for future lead compound identification and optimization to design potent and selective Clk4 and Dyrk1A inhibitors.