Docking and 3D-QSAR studies of 7-hydroxycoumarin derivatives as CK2 inhibitors

Serine/Threonine Protein Kinases as Attractive Targets for Anti-Cancer Drugs-An Innovative Approach to Ligand Tuning Using Combined Quantum Chemical Calculations, Molecular Docking, Molecular Dynamic Simulations, and Network-like Similarity Graphs

Serine/threonine protein kinases (CK2, PIM-1, RIO1) are constitutively active, highly conserved, pleiotropic, and multifunctional kinases, which control several signaling pathways and regulate many cellular functions, such as cell activity, survival, proliferation, and apoptosis. Over the past decades, they have gained increasing attention as potential therapeutic targets, ranging from various cancers and neurological, inflammation, and autoimmune disorders to viral diseases, including COVID-19. Despite the accumulation of a vast amount of experimental data, there is still no "recipe" that would facilitate the search for new effective kinase inhibitors. The aim of our study was to develop an effective screening method that would be useful for this purpose. A combination of Density Functional Theory calculations and molecular docking, supplemented with newly developed quantitative methods for the comparison of the binding modes, provided deep insight into the set of desirable properties responsible for their inhibition. The mathematical metrics helped assess the distance between the binding modes, while heatmaps revealed the locations in the ligand that should be modified according to binding site requirements. The Structure-Binding Affinity Index and Structural-Binding Affinity Landscape proposed in this paper helped to measure the extent to which binding affinity is gained or lost in response to a relatively small change in the ligand's structure. The combination of the physico-chemical profile with the aforementioned factors enabled the identification of both "dead" and "promising" search directions. Tests carried out on experimental data have validated and demonstrated the high efficiency of the proposed innovative approach. Our method for quantifying differences between the ligands and their binding capabilities holds promise for guiding future research on new anti-cancer agents.

Machine Learning Models for the Classification of CK2 Natural Products Inhibitors with Molecular Fingerprint Descriptors

Casein kinase 2 (CK2) is considered an important target for anti-cancer drugs. Given the structural diversity and broad spectrum of pharmaceutical activities of natural products, numerous studies have been performed to prove them as valuable sources of drugs. However, there has been little study relevant to identifying structural factors responsible for their inhibitory activity against CK2 with machine learning methods. In this study, classification studies were conducted on 115 natural products as CK2 inhibitors. Seven machine learning methods along with six molecular fingerprints were employed to develop qualitative classification models. The performances of all models were evaluated by cross-validation and test set. By taking predictive accuracy(CA), the area under receiver operating characteristic (AUC), and (MCC)as three performance indicators, the optimal models with high reliability and predictive ability were obtained, including the Extended Fingerprint-Logistic Regression model (CA = 0.859, AUC = 0.826, MCC = 0.520) for training test andPubChem fingerprint along with the artificial neural model (CA = 0.826, AUC = 0.933, MCC = 0.628) for test set. Meanwhile, the privileged substructures responsible for their inhibitory activity against CK2 were also identified through a combination of frequency analysis and information gain. The results are expected to provide useful information for the further utilization of natural products and the discovery of novel CK2 inhibitors.

Discovery of Keap1-Nrf2 small-molecule inhibitors from phytochemicals based on molecular docking

Various phytochemicals have been reported to protect against oxidative stress. However, the mechanism underlying has not been systematically evaluated, which limited their application in disease treatment. Nuclear factor erythroid 2−related factor 2 (Nrf2), a central transcription factor in oxidative stress response related to numerous diseases, is activated after dissociating from the cytoskeleton−anchored Kelch−like ECH−associated protein 1 (Keap1). The Keap1–Nrf2 protein–protein interaction has become an important drug target. This study was designed to clarify whether antioxidantive phytochemicals inhibit the Keap1–Nrf2 protein–protein interaction and activate the Nrf2-ARE signaling pathway efficiently. Molecular docking and 3D−QSAR were applied to evaluate the interaction effects between 178 antioxidant phytochemicals and the Nrf2 binding site in Keap1. The Nrf2 activation effect was tested on a H₂O₂−induced oxidative−injured cell model. Results showed that the 178 phytochemicals could be divided into high−, medium−, and low−total−score groups depending on their binding affinity with Keap1, and the high−total−score group consisted of 24 compounds with abundant oxygen or glycosides. Meanwhile, these compounds could bind with key amino acids in the structure of the Keap1−Nrf2 interface. Compounds from high−total−score group show effective activation effects on Nrf2. In conclusion, phytochemicals showed high binding affinity with Keap1 are promising new Nrf2 activators.

•

178 phytochemicals were collected and the Keap1-binding-affinity was estimated.

•

24 compounds with high Keap1-binding-affinity was obtained.

•

Compounds with high Keap1-binding-affinity effectively activated Nrf2.

New molecular insights into the tyrosyl-tRNA synthase inhibitors: CoMFA, CoMSIA analyses and molecular docking studies

Drug resistance caused by excessive and indiscriminate antibiotic usage has become a serious public health problem. The need of finding new antibacterial drugs is more urgent than ever before. Tyrosyl-tRNA synthase was proved to be a potent target in combating drug-resistant bacteria. In silico methodologies including molecular docking and 3D-QSAR were employed to investigate a series of newly reported tyrosyl-tRNA synthase inhibitors of furanone derivatives. Both internal and external cross-validation were conducted to obtain high predictive and satisfactory CoMFA model (q 2 = 0.611, r 2pred  = 0.933, r 2m  = 0.954) and CoMSIA model (q 2 = 0.546, r 2pred  = 0.959, r 2m  = 0.923). Docking results, which correspond with CoMFA/CoMSIA contour maps, gave the information for interactive mode exploration. Ten new molecules designed on the basis of QSAR and docking models have been predicted more potent than the most active compound 3-(4-hydroxyphenyl)-4-(2-morpholinoethoxy)furan-2(5H)-one (15) in the literatures. The results expand our understanding of furanones as inhibitors of tyrosyl-tRNA synthase and could be helpful in rationally designing of new analogs with more potent inhibitory activities.

The Development of CK2 Inhibitors: From Traditional Pharmacology to in Silico Rational Drug Design

Casein kinase II (CK2) is an ubiquitous and pleiotropic serine/threonine protein kinase able to phosphorylate hundreds of substrates. Being implicated in several human diseases, from neurodegeneration to cancer, the biological roles of CK2 have been intensively studied. Upregulation of CK2 has been shown to be critical to tumor progression, making this kinase an attractive target for cancer therapy. Several CK2 inhibitors have been developed so far, the first being discovered by "trial and error testing". In the last decade, the development of in silico rational drug design has prompted the discovery, de novo design and optimization of several CK2 inhibitors, active in the low nanomolar range. The screening of big chemical libraries and the optimization of hit compounds by Structure Based Drug Design (SBDD) provide telling examples of a fruitful application of rational drug design to the development of CK2 inhibitors. Ligand Based Drug Design (LBDD) models have been also applied to CK2 drug discovery, however they were mainly focused on methodology improvements rather than being critical for de novo design and optimization. This manuscript provides detailed description of in silico methodologies whose applications to the design and development of CK2 inhibitors proved successful and promising.

Prediction and evaluation of the lipase inhibitory activities of tea polyphenols with 3D-QSAR models

The extraordinary hypolipidemic effects of polyphenolic compounds from tea have been confirmed in our previous study. To gain compounds with more potent activities, using the conformations of the most active compound revealed by molecular docking, a 3D-QSAR pancreatic lipase inhibitor model with good predictive ability was established and validated by CoMFA and CoMISA methods. With good statistical significance in CoMFA (r²_cv = 0.622, r² = 0.956, F = 261.463, SEE = 0.096) and CoMISA (r²_cv = 0.631, r² = 0.932, F = 75.408, SEE = 0.212) model, we summarized the structure-activity relationship between polyphenolic compounds and pancreatic lipase inhibitory activities and find the bulky substituents in R₂, R₄ and R₅, hydrophilic substituents in R₁ and electron withdrawing groups in R₂ are the key factors to enhance the lipase inhibitory activities. Under the guidance of the 3D-QSAR results, (2R,3R,2′R,3′R)-desgalloyloolongtheanin-3,3′-O-digallate (DOTD), a potent lipase inhibitor with an IC50 of 0.08 μg/ml, was obtained from EGCG oxidative polymerization catalyzed by crude polyphenol oxidase. Furthermore, DOTD was found to inhibit lipid absorption in olive oil-loaded rats, which was related with inhibiting the activities of lipase in the intestinal mucosa and contents.

Structural insights into flavones as protein kinase CK2 inhibitors derived from a combined computational study

Binding conformation of flavone inhibitors to protein kinase CK2.

Insight into the interactions between novel isoquinolin-1,3-dione derivatives and cyclin-dependent kinase 4 combining QSAR and molecular docking

Several small-molecule CDK inhibitors have been identified, but none have been approved for clinical use in the past few years. A new series of 4-[(3-hydroxybenzylamino)-methylene]-4H-isoquinoline-1,3-diones were reported as highly potent and selective CDK4 inhibitors. In order to find more potent CDK4 inhibitors, the interactions between these novel isoquinoline-1,3-diones and cyclin-dependent kinase 4 was explored via in silico methodologies such as 3D-QSAR and docking on eighty-one compounds displaying potent selective activities against cyclin-dependent kinase 4. Internal and external cross-validation techniques were investigated as well as region focusing, bootstraping and leave-group-out. A training set of 66 compounds gave the satisfactory CoMFA model (q² = 0.695, r² = 0.947) and CoMSIA model (q² = 0.641, r² = 0.933). The remaining 15 compounds as a test set also gave good external predictive abilities with r²_pred values of 0.875 and 0.769 for CoMFA and CoMSIA, respectively. The 3D-QSAR models generated here predicted that all five parameters are important for activity toward CDK4. Surflex-dock results, coincident with CoMFA/CoMSIA contour maps, gave the path for binding mode exploration between the inhibitors and CDK4 protein. Based on the QSAR and docking models, twenty new potent molecules have been designed and predicted better than the most active compound 12 in the literatures. The QSAR, docking and interactions analysis expand the structure-activity relationships of constrained isoquinoline-1,3-diones and contribute towards the development of more active CDK4 subtype-selective inhibitors.

Combined 3D-QSAR and Docking Modelling Study on Indolocarbazole Series Compounds as Tie-2 Inhibitors

Tie-2, a kind of endothelial cell tyrosine kinase receptor, is required for embryonic blood vessel development and tumor angiogenesis. Several compounds that showed potent activity toward this attractive anticancer drug target in the assay have been reported. In order to investigate the structure-activity correlation of indolocarbazole series compounds and modify them to improve their selectivity and activity, 3D-QSAR models were built using CoMFA and CoMSIA methods and molecular docking was used to check the results. Based on the common sketch align, two good QSAR models with high predictabilities (CoMFA model: q² = 0.823, r² = 0.979; CoMSIA model: q² = 0.804, r² = 0.967) were obtained and the contour maps obtained from both models were applied to identify the influence on the biological activity. Molecular docking was then used to confirm the results. Combined with the molecular docking results, the detail binding mode between the ligands and Tie-2 was elucidated, which enabled us to interpret the structure-activity relationship. These satisf actory results not only offered help to comprehend the action mechanism of indolocarbazole series compounds, but also provide new information for the design of new potent inhibitors.

A 3D-QSAR model based screen for dihydropyridine-like compound library to identify inhibitors of amyloid beta (Aβ) production

Abnormal accumulation of amyloid beta peptide (Aβ) is one of the hallmarks of Alzheimer's disease progression. Practical limitations such as cost , poor hit rates and a lack of well characterized targets are a major bottle neck in the in vitro screening of a large number of chemical libraries and profiling them to identify Aβ inhibitors. We used a limited set of 1,4 dihydropyridine (DHP)-like compounds from our model set (MS) of 24 compounds which inhibit Aβ as a training set and built 3D-QSAR (Three-dimensional Quantitative Structure-Activity Relationship) models using the Phase program (SchrÖdinger, USA). We developed a 3D-QSAR model that showed the best prediction for Aβ inhibition in the test set of compounds and used this model to screen a 1,043 DHP-like small library set of (LS) compounds. We found that our model can effectively predict potent hits at a very high rate and result in significant cost savings when screening larger libraries. We describe here our in silico model building strategy, model selection parameters and the chemical features that are useful for successful screening of DHP and DHP-like chemical libraries for Aβ inhibitors.

3D-QSAR and molecular docking studies on fused pyrazoles as p38α mitogen-activated protein kinase inhibitors

The p38α mitogen-activated protein kinase (MAPK) has become an attractive target for the treatment of many diseases such as rheumatoid arthritis, inflammatory bowel disease and Crohn’s disease. In this paper, 3D-QSAR and molecular docking studies were performed on 59 p38α MAPK inhibitors. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were applied to determine the structural requirements for potency in inhibiting p38α MAPK. The resulting model of CoMFA and CoMSIA exhibited good r²_cv values of 0.725 and 0.609, and r² values of 0.961 and 0.905, respectively. Molecular docking was used to explore the binding mode between the inhibitors and p38α MAPK. We have accordingly designed a series of novel p38α MAPK inhibitors by utilizing the structure-activity relationship (SAR) results revealed in the present study, which were predicted with excellent potencies in the developed models. The results provided a useful guide to design new compounds for p38α MAPK inhibitors.

3D-QSAR and molecular docking studies of azaindole derivatives as Aurora B kinase inhibitors

The Aurora kinases have been regarded as attractive targets for the development of new anticancer agents. Recently a series of azaindole derivatives with Aurora B inhibitory activities were reported. To explore the relationship between the structures of substituted azaindole derivatives and their inhibition of Aurora B, 3D-QSAR and molecular docking studies were performed on a dataset of 41 compounds. 3D-QSAR, including CoMFA and CoMSIA, were applied to identify the key structures impacting their inhibitory potencies. The CoMSIA model showed better results than CoMFA, with r(2)(cv) value of 0.575 and r(2) value of 0.987. 3D contour maps generated from CoMFA and CoMSIA along with the docking binding structures provided enough information about the structural requirements for better activity. Based on the structure-activity relationship revealed by the present study, we have designed a set of novel Aurora B inhibitors that showed excellent potencies in the developed models. Thus, our results allowed us to design new derivatives with desired activities.