Exploiting ChEMBL database to identify indole analogs as HCV replication inhibitors

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR, molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure-activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q² = 0.743,  r² = 0.980) and the CoMSIA model (q² = 0.813,  r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.

In silico design of novel benzohydroxamate-based compounds as inhibitors of histone deacetylase 6 based on 3D-QSAR, molecular docking, and molecular dynamics simulations

In silico design of benzohydroxamate-based selective HDAC6 inhibitors.

The Acid Strength of the Lewis-Brønsted Superacids - A QSPR Study

The acidity of Lewis-Brønsted superacids can be derived from the theoretical calculations as the Gibbs free energy of the deprotonation reaction (ΔG_acid ), which describes the tendency of a studied compound to donate a proton. This paper presents the first Quantitative Structure - Property Relationship (QSPR) model that correlates the ΔG_acid of superacid (HF/MeX₃ formula (X=F, Cl, Br)) with their structure. Developed model is well fitted, roubustness, has good predictive abilities, fulfills all OECD recommendation for good model. Obtained results provide the insight into the relation of structural features of superacids, which are responsible for their acid strength - the structures characterized by strong F-Me dative bond (with relatively large vibrational frequency), small positive partial atomic charge on Me central atom, possibly large polarity exhibit large acid strength. Such assumption can be used in the future as valuable information in the process of the designing new, stronger, more effective superacids.

Combined comparative molecular field analysis, comparative molecular similarity indices analysis, molecular docking and molecular dynamics studies of histone deacetylase 6 inhibitors

Human histone deacetylase isoform 6 (HDAC6) has been shown to have an immense role in cell motility and aggresome formation and is being an attractive selective target for the treatment of multiple tumour types and neurodegenerative conditions. The discovery of selective HDAC6 inhibitors with new chemical functionalities is therefore of utmost interest to researchers. In order to examine the structural requirements for HDAC6-specific inhibitors and to derive predictive model which can be used for designing new selective HDAC6 inhibitors, a three-dimensional quantitative structure-activity relationship study was carried out on a diverse set of ligands using common feature-based pharmacophore alignment followed by employing comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. The models displayed high correlation of 0.978 and 0.991 for best CoMFA and CoMSIA models, respectively, and a good statistical significance. The model could be used for predicting activities of the test set compounds as well as for deriving useful information regarding steric, electrostatic, hydrophobic properties of the molecules used in this study. Further, the training and test set molecules were docked into the HDAC6 binding site and molecular dynamics was carried out to suggest structural requirements for design of new inhibitors.

A nanoinformatics decision support tool for the virtual screening of gold nanoparticle cellular association using protein corona fingerprints

The increasing use of nanoparticles (NPs) in a wide range of consumer and industrial applications has necessitated significant effort to address the challenge of characterizing and quantifying the underlying nanostructure - biological response relationships to ensure that these novel materials can be exploited responsibly and safely. Such efforts demand reliable experimental data not only in terms of the biological dose-response, but also regarding the physicochemical properties of the NPs and their interaction with the biological environment. The latter has not been extensively studied, as a large surface to bind biological macromolecules is a unique feature of NPs that is not relevant for chemicals or pharmaceuticals, and thus only limited data have been reported in the literature quantifying the protein corona formed when NPs interact with a biological medium and linking this with NP cellular association/uptake. In this work we report the development of a predictive model for the assessment of the biological response (cellular association, which can include both internalized NPs and those attached to the cell surface) of surface-modified gold NPs, based on their physicochemical properties and protein corona fingerprints, utilizing a dataset of 105 unique NPs. Cellular association was chosen as the end-point for the original experimental study due to its relevance to inflammatory responses, biodistribution, and toxicity in vivo. The validated predictive model is freely available online through the Enalos Cloud Platform ( http://enalos.insilicotox.com/NanoProteinCorona/ ) to be used as part of a regulatory or NP safe-by-design decision support system. This online tool will allow the virtual screening of NPs, based on a list of the significant NP descriptors, identifying those NPs that would warrant further toxicity testing on the basis of predicted NP cellular association.

HIVprotI: an integrated web based platform for prediction and design of HIV proteins inhibitors

A number of anti-retroviral drugs are being used for treating Human Immunodeficiency Virus (HIV) infection. Due to emergence of drug resistant strains, there is a constant quest to discover more effective anti-HIV compounds. In this endeavor, computational tools have proven useful in accelerating drug discovery. Although methods were published to design a class of compounds against a specific HIV protein, but an integrated web server for the same is lacking. Therefore, we have developed support vector machine based regression models using experimentally validated data from ChEMBL repository. Quantitative structure activity relationship based features were selected for predicting inhibition activity of a compound against HIV proteins namely protease (PR), reverse transcriptase (RT) and integrase (IN). The models presented a maximum Pearson correlation coefficient of 0.78, 0.76, 0.74 and 0.76, 0.68, 0.72 during tenfold cross-validation on IC₅₀ and percent inhibition datasets of PR, RT, IN respectively. These models performed equally well on the independent datasets. Chemical space mapping, applicability domain analyses and other statistical tests further support robustness of the predictive models. Currently, we have identified a number of chemical descriptors that are imperative in predicting the compound inhibition potential. HIVprotI platform (http://bioinfo.imtech.res.in/manojk/hivproti) would be useful in virtual screening of inhibitors as well as designing of new molecules against the important HIV proteins for therapeutics development.

Quantum molecular modeling of hepatitis C virus inhibition through non-structural protein 5B polymerase receptor binding of C5-arylidene rhodanines

We have carried out high-level quantum chemical computations followed by molecular docking studies on a set of 17C₅-arylidene rhodanine isomers to provide insights into the binding modes with different reported binding pockets of the nonstructural protein 5B (NS5B) polymerase that contribute to the hepatitis C virus (HCV) inhibition. We optimized the multi-target profile of the selected rhodanine analogs to investigate potential non-nucleotide inhibitors (NNIs) by quantum chemical optimization of the 18 isomers followed by docking with quantum chemically optimized structures of each isomer with NS5B polymerase at multiple binding pockets. The binding affinities of the PP-I, PP-II and TP-II pockets of NS5B polymerase were analyzed for all the 17 isomers of 2-[(5Z)-5-(2,4-dichlorobenzylidene)-4-oxo-2-thioxo-1,3-thiazolidin-3-yl]-3-phenylpropanoic acid. On the basis of binding propensity at the different pockets and inhibitor constants, we ranked these isomers as potential candidates for the HCV inhibition. We have identified four isomers as promising NNIs of NS5B polymerase with comparable binding and inhibition to the standard (1,3) dichloro substituted isomer that exhibits in vitro activity and several other isomers as candidates in a "multi-targeted drug" approach.

Enalos Suite: New Cheminformatics Platform for Drug Discovery and Computational Toxicology

In this chapter we present and discuss, with the aid of several representative case studies from drug discovery and computational toxicology, a new cheminformatics platform, Enalos Suite, that was developed with open source and freely available software. Enalos Suite ( http://enalossuite.novamechanics.com/ ) was designed and developed as a useful tool to address a variety of cheminformatics problems, given that it expedites tasks performed in predictive modeling and allows access, data mining and manipulation for multiple chemical databases (PubChem, UniChem, etc.). Enalos Suite was carefully designed to permit its extension and adjustment to the special field of interest of each user, including, for instance, nanoinformatics, biomedical, and other applications. To demonstrate the functionalities of Enalos Suite that are useful in different cheminformatics applications, we present indicative case studies that include the exploitation of chemical databases within a drug discovery project, the calculation of molecular descriptors, and finally the development of a predictive QSAR model validated according to OECD principles. We aspire that at the end of this chapter, the reader will capture the effectiveness of different functionalities included in the Enalos Suite that could be of significant value in a multitude of cheminformatics applications.

Integrating docking scores and key interaction profiles to improve the accuracy of molecular docking: towards novel B-RafV600E inhibitors

A set of ninety-eight B-RafV600E inhibitors was used for the development of a molecular docking based QSAR model using linear and non-linear regression models. The integration of docking scores and key interaction profiles significantly improved the accuracy of the QSAR models, providing reasonable statistical parameters (Rtrain2 = 0.935, Rtest2 = 0.728 and QCV2 = 0.905). The established MD-SVR (molecular docking based SMV regression) model as well as model screening of a natural product database was carried out and two natural products (quercetin and myricetin) with good prediction activities were biologically evaluated. Both compounds exhibited promising B-RafV600E inhibitory activities (ICQuercetin50 = 7.59 μM and ICMyricetin50 = 1.56 μM), suggesting a high reliability and good applicability of the established MD-SVR model in the future development of B-RafV600E inhibitors with high efficacy.

Multi-Layer Identification of Highly-Potent ABCA1 Up-Regulators Targeting LXRβ Using Multiple QSAR Modeling, Structural Similarity Analysis, and Molecular Docking

In this study, in silico approaches, including multiple QSAR modeling, structural similarity analysis, and molecular docking, were applied to develop QSAR classification models as a fast screening tool for identifying highly-potent ABCA1 up-regulators targeting LXRβ based on a series of new flavonoids. Initially, four modeling approaches, including linear discriminant analysis, support vector machine, radial basis function neural network, and classification and regression trees, were applied to construct different QSAR classification models. The statistics results indicated that these four kinds of QSAR models were powerful tools for screening highly potent ABCA1 up-regulators. Then, a consensus QSAR model was developed by combining the predictions from these four models. To discover new ABCA1 up-regulators at maximum accuracy, the compounds in the ZINC database that fulfilled the requirement of structural similarity of 0.7 compared to known potent ABCA1 up-regulator were subjected to the consensus QSAR model, which led to the discovery of 50 compounds. Finally, they were docked into the LXRβ binding site to understand their role in up-regulating ABCA1 expression. The excellent binding modes and docking scores of 10 hit compounds suggested they were highly-potent ABCA1 up-regulators targeting LXRβ. Overall, this study provided an effective strategy to discover highly potent ABCA1 up-regulators.

Structural Investigation for Optimization of Anthranilic Acid Derivatives as Partial FXR Agonists by in Silico Approaches

In this paper, a three level in silico approach was applied to investigate some important structural and physicochemical aspects of a series of anthranilic acid derivatives (AAD) newly identified as potent partial farnesoid X receptor (FXR) agonists. Initially, both two and three-dimensional quantitative structure activity relationship (2D- and 3D-QSAR) studies were performed based on such AAD by a stepwise technology combined with multiple linear regression and comparative molecular field analysis. The obtained 2D-QSAR model gave a high predictive ability (R²(train) = 0.935, R²(test) = 0.902, Q²(LOO) = 0.899). It also uncovered that number of rotatable single bonds (b_rotN), relative negative partial charges (RPC(-)), oprea's lead-like (opr_leadlike), subdivided van der Waal's surface area (SlogP_VSA2) and accessible surface area (ASA) were important features in defining activity. Additionally, the derived3D-QSAR model presented a higher predictive ability (R²(train) = 0.944, R²(test) = 0.892, Q²(LOO) = 0.802). Meanwhile, the derived contour maps from the 3D-QSAR model revealed the significant structural features (steric and electronic effects) required for improving FXR agonist activity. Finally, nine newly designed AAD with higher predicted EC50 values than the known template compound were docked into the FXR active site. The excellent molecular binding patterns of these molecules also suggested that they can be robust and potent partial FXR agonists in agreement with the QSAR results. Overall, these derived models may help to identify and design novel AAD with better FXR agonist activity.

Study of novel pyrazolo[3,4-d]pyrimidine derivatives as selective TgCDPK1 inhibitors: molecular docking, structure-based 3D-QSAR and molecular dynamics simulation

We explored the structural features that have an impact on TgCDPK1 activity and TgCDPK1/Src selectivity by multi-computational methods with different statistical models.

The KNIME based classification models for yellow fever virus inhibition

The Naïve Bayes method as implemented in KNIME platform for classification of YFV inhibition. The best classification model is able to correctly discriminate >90% of inhibitors and non-inhibitors.

Enalos InSilicoNano platform: an online decision support tool for the design and virtual screening of nanoparticles

A QNAR model, available online through Enalos InSilicoNano platform, has been developed and validated for the risk assessment of nanoparticles (NPs).