Multitemplate Alignment Method for the Development of a Reliable 3D-QSAR Model for the Analysis of MMP3 Inhibitors

From Anti-infective Agents to Cancer Therapy: a Drug Repositioning Study Revealed a New Use for Nitrofuran Derivatives

BACKGROUND

The progression of ovarian cancer seems to be related to HDAC1, HDAC3 and HDAC6 activity. A possible strategy for improving therapies for treating ovarian carcinoma, minimizing the preclinical screenings, is the repurposing of already approved pharmaceutical products as inhibitors of these enzymes.

OBJECTIVE

This work was aimed to implement a computational strategy for identifying new HDAC inhibitors for ovarian carcinoma treatment among approved drugs.

METHOD

The CHEMBL database was used to construct training, test and decoys sets for performing and validating HDAC1, HDAC3 and HDAC6 3D-QSAR models obtained by using FLAP program. Docking and MD simulations were used in combination with the generated models to identify novel potential HDAC inhibitors. Cell viability assays and Western blot analyses were performed on normal and cancer cells for a direct evaluation of the anti-proliferative activity and an in vitro estimation of HDAC inhibition of the compounds selected through in silico screening.

RESULT

The best quantitative prediction was obtained for the HDAC6 3D-QSAR model. The screening of approved drugs highlighted a new potential use as HDAC inhibitors for some compounds, in particular nitrofuran derivatives, usually known for their antibacterial activity, and frequently used as antimicrobial adjuvant therapy in cancer treatment. Experimental evaluation of these derivatives highlighted a significant antiproliferative activity against cancer cell lines overexpressing HDAC6, and an increase in acetylated alpha-tubulin levels.

CONCLUSION

Experimental results support the hypothesis of a potential direct interaction of nitrofuran derivatives with HDACs. In addition to the possible repurposing of already approved drugs, this work suggests the nitro group as a new zinc binding group, able to interact with the catalytic zinc ion of HDACs.

Query-guided protein-protein interaction inhibitor discovery

Protein-protein interactions (PPIs) are central to biological mechanisms, and can serve as compelling targets for drug discovery. Yet, the discovery of small molecule inhibitors of PPIs remains challenging given the large and typically shallow topography of the interacting protein surfaces. Here, we describe a general approach to the discovery of orthosteric PPI inhibitors that mimic specific secondary protein structures. Initially, hot residues at protein-protein interfaces are identified in silico or from experimental data, and incorporated into secondary structure-based queries. Virtual libraries of small molecules are then shape-matched against the queries, and promising ligands docked to target proteins. The approach is exemplified experimentally using two unrelated PPIs that are mediated by an α-helix (p53/hDM2) and a β-strand (GKAP/SHANK1-PDZ). In each case, selective PPI inhibitors are discovered with low μM activity as determined by a combination of fluorescence anisotropy and 1H-15N HSQC experiments. In addition, hit expansion yields a series of PPI inhibitors with defined structure-activity relationships. It is envisaged that the generality of the approach will enable discovery of inhibitors of a wide range of unrelated secondary structure-mediated PPIs.

Molecular docking and 4D-QSAR studies of metastatic cancer inhibitor thiazoles

By using the molecular docking and 4D-QSAR analysis, it is aimed to find the interaction points in the receptor binding site of transforming growth factor-beta (TGF-beta) used to inhibit invasion and metastasis. To elucidate the interaction points of receptor, different types of local reactive descriptor (LRD) of ligands have been used. Activity values related to interaction energy between the ligand-receptor (L-R) were determined by nonlinear least squares (NLLS) using the Levenberg-Marquardt (LM) algorithm. Using the Molecule Comparative Electron Topology (MCET) method, the 3D pharmacophore model (3D-PhaM) was obtained after alignment and superimposition of the molecules, and also confirmed by molecular docking method. With the leave one out-cross validation (LOO-CV) method, the best predictions are q² or r_CV² = 0.789 for the 51 compounds in the internal training set and r² = 0.785 for the 13 compounds in the external test set. Furthermore, the predictive capability of the advanced QSAR model is more precisely calculated with the r_m² metric (r_m² = 0.769).

In silico ligand-based modeling of hBACE-1 inhibitors

Alzheimer's disease is a chronic neurodegenerative disease affecting more than 30 million people worldwide. Development of small molecule inhibitors of human β-secretase 1 (hBACE-1) is being the focus of pharmaceutical industry for the past 15-20 years. Here, we successfully applied multiple ligand-based in silico modeling techniques to understand the inhibitory activities of a diverse set of small molecule hBACE-1 inhibitors reported in the scientific literature. Strikingly, the use of only a small subset of 230 (13%) molecules allowed us to develop quality models that performed reasonably well on the validation set of 1,476 (87%) inhibitors. Varying the descriptor sets and the complexity of the modeling techniques resulted in only minor improvements to the model's performance. The current results demonstrate that predictive models can be built by choosing appropriate modeling techniques in spite of using small datasets consisting of diverse chemical classes, a scenario typical in triaging of high-throughput screening results to identify false negatives. We hope that these encouraging results will help the community to develop more predictive models that would support research efforts for the debilitating Alzheimer's disease. Additionally, the integrated diversity of the techniques employed will stimulate scientists in the field to use in silico statistical modeling techniques like these to derive better models to help advance the drug discovery projects faster.

Alignment-independent technique for 3D QSAR analysis

Molecular biochemistry is controlled by 3D phenomena but structure-activity models based on 3D descriptors are infrequently used for large data sets because of the computational overhead for determining molecular conformations. A diverse dataset of 146 androgen receptor binders was used to investigate how different methods for defining molecular conformations affect the performance of 3D-quantitative spectral data activity relationship models. Molecular conformations tested: (1) global minimum of molecules' potential energy surface; (2) alignment-to-templates using equal electronic and steric force field contributions; (3) alignment using contributions "Best-for-Each" template; (4) non-energy optimized, non-aligned (2D > 3D). Aggregate predictions from models were compared. Highest average coefficients of determination ranged from R Test (2) = 0.56 to 0.61. The best model using 2D > 3D (imported directly from ChemSpider) produced R Test (2) = 0.61. It was superior to energy-minimized and conformation-aligned models and was achieved in only 3-7 % of the time required using the other conformation strategies. Predictions averaged from models built on different conformations achieved a consensus R Test (2) = 0.65. The best 2D > 3D model was analyzed for underlying structure-activity relationships. For the compound strongest binding to the androgen receptor, 10 substructural features contributing to binding were flagged. Utility of 2D > 3D was compared for two other activity endpoints, each modeling a medium sized data set. Results suggested that large scale, accurate predictions using 2D > 3D SDAR descriptors may be produced for interactions involving endocrine system nuclear receptors and other data sets in which strongest activities are produced by fairly inflexible substrates.

POVME 2.0: An Enhanced Tool for Determining Pocket Shape and Volume Characteristics

Analysis of macromolecular/small-molecule binding pockets can provide important insights into molecular recognition and receptor dynamics. Since its release in 2011, the POVME (POcket Volume MEasurer) algorithm has been widely adopted as a simple-to-use tool for measuring and characterizing pocket volumes and shapes. We here present POVME 2.0, which is an order of magnitude faster, has improved accuracy, includes a graphical user interface, and can produce volumetric density maps for improved pocket analysis. To demonstrate the utility of the algorithm, we use it to analyze the binding pocket of RNA editing ligase 1 from the unicellular parasite Trypanosoma brucei, the etiological agent of African sleeping sickness. The POVME analysis characterizes the full dynamics of a potentially druggable transient binding pocket and so may guide future antitrypanosomal drug-discovery efforts. We are hopeful that this new version will be a useful tool for the computational- and medicinal-chemist community.

Challenging the gold standard for 3D-QSAR: template CoMFA versus X-ray alignment

X-ray-based alignments of bioactive compounds are commonly used to correlate structural changes with changes in potencies, ultimately leading to three-dimensional quantitative structure-activity relationships such as CoMFA or CoMSIA models that can provide further guidance for the design of new compounds. We have analyzed data sets where the alignment of the compounds is entirely based on experimentally derived ligand poses from X-ray-crystallography. We developed CoMFA and CoMSIA models from these X-ray-determined receptor-bound conformations and compared the results with models generated from ligand-centric Template CoMFA, finding that the fluctuations in the positions and conformations of compounds dominate X-ray-based alignments can yield poorer predictions than those from the self-consistent template CoMFA alignments. Also, when there exist multiple different binding modes, structural interpretation in terms of binding site constraints can often be simpler with template-based alignments than with X-ray-based alignments.

Structure-based ensemble-QSAR model: a novel approach to the study of the EGFR tyrosine kinase and its inhibitors

AIM

To develop a novel 3D-QSAR approach for study of the epidermal growth factor receptor tyrosine kinase (EGFR TK) and its inhibitors.

METHODS

One hundred thirty nine EGFR TK inhibitors were classified into 3 clusters. Ensemble docking of these inhibitors with 19 EGFR TK crystal structures was performed. Three protein structures that showed the best recognition of each cluster were selected based on the docking results. Then, a novel QSAR (ensemble-QSAR) building method was developed based on the ligand conformations determined by the corresponding protein structures.

RESULTS

Compared with the 3D-QSAR model, in which the ligand conformations were determined by a single protein structure, ensemble-QSAR exhibited higher R2 (0.87) and Q2 (0.78) values and thus appeared to be a more reliable and better predictive model. Ensemble-QSAR was also able to more accurately describe the interactions between the target and the ligands.

CONCLUSION

The novel ensemble-QSAR model built in this study outperforms the traditional 3D-QSAR model in rationality, and provides a good example of selecting suitable protein structures for docking prediction and for building structure-based QSAR using available protein structures.

PLS and shape-based similarity analysis of maleimides--GSK-3 inhibitors

CONTEXT

Glycogen synthase kinase-3 (GSK-3) overactivity was correlated with several pathologies including type 2 diabetes mellitus, Alzheimer's disease, cancer, inflammation, obesity, etc.

OBJECTIVE

The aim of the current investigation was to model the inhibitory activity of maleimide derivatives--inhibitors of GSK-3, to evaluate the impact of alignment on statistical performances of the Quantitative Structure-Activity Relationship (QSAR) and the effect of the template on shape-similarity--binding affinity relationship.

MATERIALS AND METHODS

Dragon descriptors were used to generate Projection to Latent Structures (PLS) models in order to identify the structural prerequisites of maleimides to inhibit GSK-3. Additionally, shape/volume structural analysis of binding site interactions was evaluated.

RESULTS

Reliable statistics R(2)(Y(CUM)) = 0.938/0.920, Q((2)(Y)(CUM)) = 0.866/0.838 for aligned and alignment free QSAR models and significant (Pearson, Kendall and Spearman) correlations between shape/volume similarity and affinities were obtained.

DISCUSSION AND CONCLUSIONS

The crucial structural features modulating the activity of maleimides include topology, charge, geometry, 2D autocorrelations, 3D-MoRSE as well as shape/volume and molecular flexibility.

Structure- and ligand-based drug design of novel p38-alpha MAPK inhibitors in the fight against the Alzheimer's disease

Alzheimer's disease (AD) is characterized microscopically by the presence of amyloid plaques, which are accumulations of beta-amyloid protein inter-neurons, and neurofibrillary tangles formed predominantly by highly phosphorylated forms of the microtubule-associated protein, tau, which form tangled masses that consume neuronal cell body, possibly leading to neuronal dysfunction and ultimately death. p38α mitogen-activated protein kinase (MAPK) has been implicated in both events associated with AD, tau phosphorylation and inflammation. p38α MAPK pathway is activated by a dual phosphorylation at Thr180 and Tyr182 residues. Drug design of p38α MAPK inhibitors is mainly focused on small molecules that compete for Adenosine triphosphate in the catalytic site. Here, we used different approaches of structure- and ligand-based drug design and medicinal chemistry strategies based on a selected p38α MAPK structure deposited in the Protein Data Bank in complex with inhibitor, as well as others reported in literature. As a result of the virtual screening experiments performed here, as well as molecular dynamics, molecular interaction fields studies, shape and electrostatic similarities, activity and toxicity predictions, and pharmacokinetic and physicochemical properties, we have selected 13 compounds that meet the criteria of low or no toxicity potential, good pharmacotherapeutic profile, predicted activities, and calculated values ​​comparable with those obtained for the reference compounds, while maintaining the main interactions observed for the most potent inhibitors.