Molecular modeling of exquisitely selective c-Met inhibitors through 3D-QSAR and molecular dynamics simulations

Molecular and Biological Investigation of Isolated Marine Fungal Metabolites as Anticancer Agents: A Multi-Target Approach

Cancer is the leading cause of death globally, with an increasing number of cases being annually reported. Nature-derived metabolites have been widely studied for their potential programmed necrosis, cytotoxicity, and anti-proliferation leading to enrichment for the modern medicine, particularly within the last couple of decades. At a more rapid pace, the concept of multi-target agents has evolved from being an innovative approach into a regular drug development procedure for hampering the multi-fashioned pathophysiology and high-resistance nature of cancer cells. With the advent of the Red Sea Penicillium chrysogenum strain S003-isolated indole-based alkaloids, we thoroughly investigated the molecular aspects for three major metabolites: meleagrin (MEL), roquefortine C (ROC), and isoroquefortine C (ISO) against three cancer-associated biological targets Cdc-25A, PTP-1B, and c-Met kinase. The study presented, for the first time, the detailed molecular insights and near-physiological affinity for these marine indole alkaloids against the assign targets through molecular docking-coupled all-atom dynamic simulation analysis. Findings highlighted the superiority of MEL's binding affinity/stability being quite in concordance with the in vitro anticancer activity profile conducted via sulforhodamine B bioassay on different cancerous cell lines reaching down to low micromolar or even nanomolar potencies. The advent of lengthy structural topologies via the metabolites' extended tetracyclic cores and aromatic imidazole arm permitted multi-pocket accommodation addressing the selectivity concerns. Additionally, the presence decorating polar functionalities on the core hydrophobic tetracyclic ring contributed compound's pharmacodynamic preferentiality. Introducing ionizable functionality with more lipophilic characters was highlighted to improve binding affinities which was also in concordance with the conducted drug-likeness/pharmacokinetic profiling for obtaining a balanced pharmacokinetic/dynamic profile. Our study adds to the knowledge regarding drug development and optimization of marine-isolated indole-based alkaloids for future iterative synthesis and pre-clinical investigations as multi-target anticancer agents.

Docking-Based 3D-QSAR Studies for 1,3,4-oxadiazol-2-one Derivatives as FAAH Inhibitors

This work aimed to construct 3D-QSAR CoMFA and CoMSIA models for a series of 31 FAAH inhibitors, containing the 1,3,4-oxadiazol-2-one moiety. The obtained models were characterized by good statistical parameters: CoMFA Q2 = 0.61, R2 = 0.98; CoMSIA Q2 = 0.64, R2 = 0.93. The CoMFA model field contributions were 54.1% and 45.9% for steric and electrostatic fields, respectively. In the CoMSIA model, electrostatic, steric, hydrogen bond donor, and hydrogen acceptor properties were equal to 34.6%, 23.9%, 23.4%, and 18.0%, respectively. These models were validated by applying the leave-one-out technique, the seven-element test set (CoMFA r2test-set = 0.91; CoMSIA r2test-set = 0.91), a progressive scrambling test, and external validation criteria developed by Golbraikh and Tropsha (CoMFA r20 = 0.98, k = 0.95; CoMSIA r20 = 0.98, k = 0.89). As the statistical significance of the obtained model was confirmed, the results of the CoMFA and CoMSIA field calculation were mapped onto the enzyme binding site. It gave us the opportunity to discuss the structure-activity relationship based on the ligand-enzyme interactions. In particular, examination of the electrostatic properties of the established CoMFA model revealed fields that correspond to the regions where electropositive substituents are not desired, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one moiety. This highlights the importance of heterocycle, a highly electronegative moiety in this area of each ligand. Examination of hydrogen bond donor and acceptor properties contour maps revealed several spots where the implementation of another hydrogen-bond-donating moiety will positively impact molecules' binding affinity, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one ring. On the other hand, there is a large isopleth that refers to the favorable H-bond properties close to the terminal phenoxy group of a ligand, which means that, generally speaking, H-bond acceptors are desired in this area.

Design, synthesis and biological evaluation of 7H-pyrrolo[2,3-d]pyrimidine derivatives containing 1,8-naphthyridine-4-one fragment

In this study, a series of pyrrolo [2,3-d]pyrimidine derivatives containing 1,8-naphthyridine-4-one fragment were synthesized and their biological activity were tested. Most of the target compounds displayed moderate to excellent activity against one or more cancer cell lines and low activity against human normal cell LO2 in vitro. The most promising compound 51, of which the IC₅₀ values were 0.66 μM, 0.38 μM and 0.44 μM against cell lines A549, Hela and MCF-7, shown more remarkable activity and better apoptosis effect than the positive control Cabozantinib. The structure-activity relationships (SARs) indicated that double-EWGs (such as R₃ = 2-Cl-4-CF₃) on the terminal phenyl rings was a key factor in improving the biological activity. In addition, the further research on compound 51 mainly included c-Met kinase activity and selectivity, concentration dependence, and molecular docking.

Discovery of Ubiquitin-Specific Protease 7 (USP7) Inhibitors with Novel Scaffold Structures by Virtual Screening, Molecular Dynamics Simulation, and Biological Evaluation

USP7 has been regarded as a potential therapeutic target for cancer. In this study, virtual screening, molecular dynamics (MD) simulation, and biological evaluation have been applied for the discovery of novel USP7 inhibitors targeting the catalytic active site. Among the obtained compounds, compound 12 with a novel scaffold structure exhibited certain USP7 inhibitory activity (Ub-AMC assay IC₅₀ = 18.40 ± 1.75 μM, Ub-Rho assay IC₅₀ = 7.75 μM). The binding affinity between USP7_CD (USP7 catalytic domain) and this hit compound was confirmed with a K_D value of 4.46 ± 0.86 μM. Preliminary in vitro studies disclosed its antiproliferative activity on human prostate cancer cell line LNCaP with an IC₅₀ value of 15.43 ± 3.49 μM. MD simulation revealed the detailed differences of protein-ligand interactions between USP7_CD and the ligands, including the reference compound ALM4 and compound 12, providing some important information for improving the bioactivity of 12. This hit compound will serve as a promising starting point for facilitating the further discovery of novel USP7 inhibitors.

Structural insights into ligand binding features of dual FABP4/5 inhibitors by molecular dynamics simulations

The fatty acid binding protein (FABP) 4 and 5 have been considered as potential targets for the treatment of metabolic diseases. A compensatory upregulation of FABP5 due to the gene ablation of FABP4 in adipocytes indicated the importance of dual FABP4/5 inhibitors. A few compounds have been discovered as dual FABP4/5 inhibitors. However, none exhibited equivalent inhibitory activity against both FABP4 and FABP5, and almost all compounds showed weaker inhibition against FABP5. To provide a better structural understanding for the design of potent dual FABP4/5 inhibitors, molecular dynamics simulations have been performed for 100 ns to disclose the ligand binding features in FABP4 and FABP5 using Amber14, respectively. Key residues were identified by analysis of close contact, hydrogen bond occupancy, binding free energy and alanine scanning mutagenesis. In addition, induced-fit effects have been observed upon ligand binding in the process of simulations. The shifted alkyl chain of ligand in FABP4 was significantly different from that in FABP5 due to the corresponding residues (Phe58^FABP4 and Leu60^FABP5). Thus, to avoid different steric effects made by these two residues, hydrophobic groups of suitable size should be taken into account. Besides, electrostatic and steric effects with Arg107^FABP4 and Arg109^FABP5 should be paid more attention to. The results will facilitate the rational design of dual FABP4/5 inhibitors.

Comparative molecular field analysis and molecular dynamics studies of the dopamine D2 receptor antagonists without a protonatable nitrogen atom

The dopaminergic hypothesis of schizophrenia is the main concept explaining the direct reasons of schizophrenia and the effectiveness of current antipsychotics. All antipsychotics present on the market are potent dopamine D₂ receptor antagonists or partial agonists. In this work we investigate a series of dopamine D₂ receptor antagonists which do not fulfill the criteria of the classical pharmacophore model as they do not possess a protonatable nitrogen atom necessary to interact with the conserved Asp(3.32). Such compounds are interesting, inter alia, due to possible better pharmacokinetic profile when compared to basic, ionizable molecules. By means of homology modeling, molecular docking and molecular dynamics we determined that the compounds investigated interact with Asp(3.32) via their amide nitrogen atom. It was found that the studied compounds stabilize the receptor inactive conformation through the effect on the ionic lock, which is typical for GPCR antagonists. We constructed a CoMFA model for the studied compounds with the following statistics: R² = 0.95, Q² = 0.63. The quality of the CoMFA model was confirmed by high value of R² of the test set, equal 0.96. The CoMFA model indicated two regions where bulky substituents are favored and two regions where bulky substituents are not beneficial. Two red contour regions near carbonyl groups were identified meaning that negative charge would be favored here. Furthermore, the S-oxide group is connected with blue contour region meaning that positive charge is favored in this position. These findings may be applied for further optimization of the studied compound series.

Cheminfomatic-based Drug Discovery of Human Tyrosine Kinase Inhibitors

Receptor Tyrosine Kinases (RTKs) are essential components for regulating cell-cell signaling and communication events in cell growth, proliferation, differentiation, survival and metabolism. Deregulation of RTKs and their associated signaling pathways can lead to a wide variety of human diseases such as immunodeficiency, diabetes, arterosclerosis, psoriasis and cancer. Thus RTKs have become one of the most important drug targets families in recent decade. Pharmaceutical companies have dedicated their research efforts towards the discovery of small-molecule inhibitors of RTKs, many of which had been approved by the U.S. Food and Drug Administration (US FDA) or are currently in clinical trials. The great successes in the development of small-molecule inhibitors of RTKs are largely attributed to the use of modern cheminformatic approaches to identifying lead scaffolds. Those include the quantitative structure-activity relationship (QSAR) modeling, as well as the structure-, and ligand-based pharmacophore modeling techniques in this case. Herein we inspected the literature thoroughly in an effort to conduct a comparative analysis of major findings regarding the essential structure-activity relationships (SARs)/pharmacophore features of known active RTK inhibitors, most of which were collected from cheminformatic modeling approaches.

Comparative molecular field analysis and molecular dynamics studies of α/β hydrolase domain containing 6 (ABHD6) inhibitors

The endocannabinoid system remains an attractive molecular target for pharmacological intervention due to its roles in the central nervous system in learning, thinking, emotional function, regulation of food intake or pain sensation, as well as in the peripheral nervous system, where it modulates the action of cardiovascular, immune, metabolic or reproductive function. α/β hydrolase domain containing 6 (ABHD6)—an enzyme forming part of the endocannabinoid system—is a newly discovered post-genomic protein acting as a 2-AG (2-arachidonoylglycerol) serine hydrolase. We have recently reported a series of 1,2,5-thiadiazole carbamates as potent and selective ABHD6 inhibitors. Here, we present comparative molecular field analysis (CoMFA) and molecular dynamics studies of these compounds. First, we performed a homology modeling study of ABHD6 based on the assumption that the catalytic triad of ABHD6 comprises Ser148–His306–Asp 278 and the oxyanion hole is formed by Met149 and Phe80. A total of 42 compounds was docked to the homology model using the Glide module from the Schrödinger suite of software and the selected docking poses were used for CoMFA alignment. A model with the following statistics was obtained: R² = 0.98, Q² = 0.55. In order to study the molecular interactions of the inhibitors with ABHD6 in detail, molecular dynamics was performed with the Desmond program. It was found that, during the simulations, the hydrogen bond between the inhibitor carbonyl group and the main chain of Phe80 is weakened, whereas a new hydrogen bond with the side chain of Ser148 is formed, facilitating the possible formation of a covalent bond.

Graphical Abstract

Quantitative structure-activity relationship study on BTK inhibitors by modified multivariate adaptive regression spline and CoMSIA methods

Bruton's tyrosine kinase (BTK) plays a crucial role in B-cell activation and development, and has emerged as a new molecular target for the treatment of autoimmune diseases and B-cell malignancies. In this study, two- and three-dimensional quantitative structure-activity relationship (2D and 3D-QSAR) analyses were performed on a series of pyridine and pyrimidine-based BTK inhibitors by means of genetic algorithm optimized multivariate adaptive regression spline (GA-MARS) and comparative molecular similarity index analysis (CoMSIA) methods. Here, we propose a modified MARS algorithm to develop 2D-QSAR models. The top ranked models showed satisfactory statistical results (2D-QSAR: Q(2) = 0.884, r(2) = 0.929, r(2)pred = 0.878; 3D-QSAR: q(2) = 0.616, r(2) = 0.987, r(2)pred = 0.905). Key descriptors selected by 2D-QSAR were in good agreement with the conclusions of 3D-QSAR, and the 3D-CoMSIA contour maps facilitated interpretation of the structure-activity relationship. A new molecular database was generated by molecular fragment replacement (MFR) and further evaluated with GA-MARS and CoMSIA prediction. Twenty-five pyridine and pyrimidine derivatives as novel potential BTK inhibitors were finally selected for further study. These results also demonstrated that our method can be a very efficient tool for the discovery of novel potent BTK inhibitors.