Studying the Binding Modes of Novel 2-Aminopyridine Derivatives as Effective and Selective c-Met Kinase Type 1 Inhibitors Using Molecular Modeling Approaches

The mesenchymal epithelial cell transforming factor c-Met, encoded by c-Met proto-oncogene and known as a high-affinity receptor for Hepatocyte Growth Factor (HGF), is one of the receptor tyrosine kinases (RTKs) members. The HGF/c-Met signaling pathway has close correlation with tumor growth, invasion and metastasis. Thus, c-Met kinase has emerged as a prominent therapeutic target for cancer drug discovery. Recently a series of novel 2-aminopyridine derivatives targeting c-Met kinase with high biological activity were reported. In this study, 3D quantitative structure-activity relationship (QSAR), molecular docking and molecular dynamics simulations (MD) were employed to research the binding modes of these inhibitors.The results show that both the atom-based and docking-based CoMFA (Q² = 0.596, R² = 0.950 in atom-based model and Q² = 0.563, R² = 0.985 in docking-based model) and CoMSIA (Q² = 0.646, R² = 0.931 in atom-based model and Q² = 0.568, R² = 0.983 in docking-based model) models own satisfactory performance with good reliabilities and powerful external predictabilities. Molecular docking study suggests that Tyr1230 and Arg1208 might be the key residues, and electrostatic and hydrogen bond interactions were shown to be vital to the activity, concordance with QSAR analysis. Then MD simulation was performed to further explore the binding mode of the most potent inhibitor. The obtained results provide important references for further rational design of c-Met Kinase type I inhibitors.     

Identification of novel nt-MGAM inhibitors for potential treatment of type 2 diabetes: Virtual screening,atom based 3D-QSAR model,docking analysis and ADME study

In this study, a virtual screening procedure was applied to identify new potential nt-MGAM inhibitors as a possible medication for type 2 diabetes. To this aim, a series of salacinol analogues were first investigated by docking analysis for their binding to the X-ray structure of the biological target nt-MGAM. Key interactions for ligand binding into the receptor active site were identified which shared common features to those found for other known inhibitors, which strengthen the results of this study. 3D QSAR model was then built and showed to be statistically significant and with a good predictive power for the training (R² = 0.99, SD = 0.17, F = 555.3 and N = 27) and test set (Q² = 0.81, Pearson(r) = 0.92, RMSE = 0.52, N = 08). The model was then used to virtually screen the ZINC database with the aim of identifying novel chemical scaffolds as potential nt-MGAM inhibitors. Further, in silico predicted ADME properties were investigated for the most promising molecules. The outcome of this investigation sheds light on the molecular characteristics of the binding of salacinol analogues to nt-MGAM enzyme and identifies new possible inhibitors which have the potential to be developed into drugs, thus significantly contributing to the design and optimization of therapeutic strategies against type 2 diabetes.     

In Vitro Evaluation of the Anti-Diabetic Potential of Aqueous Acetone Helichrysum petiolare Extract (AAHPE) with Molecular Docking Relevance in Diabetes Mellitus

Diabetes mellitus (DM) is a chronic metabolic condition that can lead to significant complications and a high fatality rate worldwide. Efforts are ramping up to find and develop novel α-glucosidase and α-amylase inhibitors that are both effective and potentially safe. Traditional methodologies are being replaced with new techniques that are less complicated and less time demanding; yet, both the experimental and computational strategies are viable and complementary in drug discovery and development. As a result, this study was conducted to investigate the in vitro anti-diabetic potential of aqueous acetone Helichrysum petiolare and B.L Burtt extract (AAHPE) using a 2-NBDG, 2-(N-(7-Nitrobenz-2-oxa-1,3-diazol-4-yl) amino)-2-deoxy-d-glucose uptake assay. In addition, we performed molecular docking of the flavonoid constituents identified and quantified by liquid chromatography-mass spectrometry (LC-MS) from AAHPE with the potential to serve as effective and safe α-amylase and α-glucosidase inhibitors, which are important in drug discovery and development. The results showed that AAHPE is a potential inhibitor of both α-amylase and α-glucosidase, with IC₅₀ values of 46.50 ± 6.17 (µg/mL) and 37.81 ± 5.15 (µg/mL), respectively. This is demonstrated by a significant increase in the glucose uptake activity percentage in a concentration-dependent manner compared to the control, with the highest AAHPE concentration of 75 µg/mL of glucose uptake activity being higher than metformin, a standard anti-diabetic drug, in the insulin-resistant HepG2 cell line. The molecular docking results displayed that the constituents strongly bind α-amylase and α-glucosidase while achieving better binding affinities that ranged from ΔG = −7.2 to −9.6 kcal/mol (compared with acarbose ΔG = −6.1 kcal/mol) for α-amylase, and ΔG = −7.3 to −9.0 kcal/mol (compared with acarbose ΔG = −6.3 kcal/mol) for α-glucosidase. This study revealed the potential use of the H. petiolare plant extract and its phytochemicals, which could be explored to develop potent and safe α-amylase and α-glucosidase inhibitors to treat postprandial glycemic levels in diabetic patients.     

QSAR Modeling,Molecular Docking and Cytotoxic Evaluation for Novel Oxidovanadium(IV) Complexes as Colon Anticancer Agents

Four new drug-based oxidovanadium (IV) complexes were synthesized and characterized by various spectral techniques, including molar conductance, magnetic measurements, and thermogravimetric analysis. Moreover, optimal structures geometry for all syntheses was obtained by the Gaussian09 program via the DFT/B3LYP method and showed that all of the metal complexes adopted a square-pyramidal structure. The essential parameters, electrophilicity (ω) value and expression for the maximum charge that an electrophile molecule may accept (ΔN_max) showed the practical biological potency of [VO(CTZ)₂] 2H₂O. The complexes were also evaluated for their propensity to bind to DNA through UV–vis absorption titration. The result revealed a high binding ability of the [VO(CTZ)₂] 2H₂O complex with K_b = 1.40 × 10⁶ M⁻¹. Furthermore, molecular docking was carried out to study the behavior of the VO (II) complexes towards colon cancer cell (3IG7) protein. A quantitative structure–activity relationship (QSAR) study was also implemented for the newly synthesized compounds. The results of validation indicate that the generated QSAR model possessed a high predictive power (R² = 0.97). Within the investigated series, the [VO(CTZ)₂] 2H₂O complex showed the greatest potential the most selective compound comparing to the stander chemotherapy drug.     

Combining ligand-based and structure-based drug design approaches to study the structure-activity relationships of a β-carboline derivative series

In this study, we investigated the structure-activity relationships of a series of β-carboline alkaloid derivatives using the 2D-QSAR and molecular docking, in order to identify the mode of interaction between β-carboline derivatives and the PLK1 kinase, and determine their key substituents responsible for the cytotoxic activity. The obtained QSAR models using multiple linear regression (MLR) and partial least squares (PLS) methods showed a high correlation between the experimental activity and the predicted one by PLS (R²_PLS?=?0.82, q²?=?0.72) and MLR (R²_MLR?=?0.82, q²?=?0.72). An external dataset was used to test the extrapolation power of the models which resulted in an R²_PLS (EV)?=?0.76; RMSE?=?0.39. The 2D-QSAR analysis reveals that lipophilicity plays an important role in the cytotoxic activity of this group of β-carboline derivatives. Indeed, the molecular docking study into the active site of the polo-like kinase (PLK1) revealed that the most active ligand 57 shows higher binding energy and interacts, especially by H-bonds and hydrophobic interactions, with the active site of the PLK1 kinase. Consequently, the results obtained from the 2D-QSAR and docking studies provided a useful tool to design new and potent β-carboline derivatives as cytotoxic agents.     

Insight into the structural requirements of benzimidazole derivatives as interleukin‐2 inducible t‐cell kinase inhibitors by computational explorations

In the present work, a set of ligand‐ and receptor‐based 3D‐QSAR models were developed to explore the structure–activity relationship of 109 benzimidazole‐based interleukin‐2‐inducible T‐cell kinase (ITK) inhibitors. In order to reveal the requisite 3D structural features impacting the biological activities, a variety of in silico modeling approaches including the comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), docking, and molecular dynamics were applied. The results showed that the ligand‐based CoMFA model (Q² = 0.552, R²_ncv = 0.908, R²_pred = 0.787, SEE = 0.252, SEP = 0.558) and CoMSIA model (Q² = 0.579, R²_ncv = 0.914, R²_pred = 0.893, SEE = 0.240, SEP = 0.538) were superior to other models with greater predictive power. In addition, a combined analysis between the 3D contour maps and docking results showed that: (1) Compounds with bulky or hydrophobic substituents near ring D and electropositive or hydrogen acceptor groups around rings C and D could increase the activity. (2) The key amino acids impacting the receptor–ligand interactions in the binding pocket are Met438, Asp500, Lys391, and Glu439. The results obtained from this work may provide helpful guidelines in design of novel benzimidazole analogs as inhibitors of ITK. © 2013 Wiley Periodicals, Inc.     

In Vitro Study of the Fibrinolytic Activity via Single Chain Urokinase-Type Plasminogen Activator and Molecular Docking of FGFC1

Fungi fibrinolytic compound 1 (FGFC1) is a rare marine-derived compound that can enhance fibrinolysis both in vitro and in vivo. The fibrinolytic activity characterization of FGFC1 mediated by plasminogen (Glu-/Lys-) and a single-chain urokinase-type plasminogen activator (pro-uPA) was further evaluated. The binding sites and mode of binding between FGFC1 and plasminogen were investigated by means of a combination of in vitro experiments and molecular docking. A 2.2-fold enhancement of fibrinolytic activity was achieved at 0.096 mM FGFC1, whereas the inhibition of fibrinolytic activity occurred when the FGFC1 concentration was above 0.24 mM. The inhibition of fibrinolytic activity of FGFC1 by 6-aminohexanoic acid (EACA) and tranexamic acid (TXA) together with the docking results revealed that the lysine-binding sites (LBSs) play a crucial role in the process of FGFC1 binding to plasminogen. The action mechanism of FGFC1 binding to plasminogen was inferred, and FGFC1 was able to induce plasminogen to exhibit an open conformation by binding through the LBSs. The molecular docking results showed that docking of ligands (EACA, FGFC1) with receptors (KR1–KR5) mainly occurred through hydrophilic and hydrophobic interactions. In addition, the binding affinity values of EACA to KR1–KR5 were −5.2, −4.3, −3.7, −4.5, and −4.3 kcal/moL, respectively, and those of FGFC1 to KR1–KR5 were −7.4, −9.0, −6.3, −8.3, and −6.7 kcal/moL, respectively. The findings demonstrate that both EACA and FGFC1 bound to KR1–KR5 with moderately high affinity. This study could provide a theoretical basis for the clinical pharmacology of FGFC1 and establish a foundation for practical applications of FGFC1.     

Targeting Cytotoxin-Associated Antigen A,a Virulent Factor of Helicobacter pylori-Associated Gastric Cancer: Structure-Based In Silico Screening of Natural Compounds

Gastric cancer is the fifth most frequent cancer and the third major cause of mortality worldwide. Helicobacter pylori, a bacterial infection linked with GC, injects the cytotoxin-associated antigen A (CagA; an oncoprotein) into host cells. When the phosphorylated CagA protein enters the cell, it attaches to other cellular components, interfering with normal cellular signaling pathways. CagA plays an important role in the progression of GC by interacting with phosphatidylserine of the host cell membrane. Therefore, disrupting the CagA–phosphatidylserine connection using small molecules appears to be a promising therapeutic approach. In this report, we screened the natural compounds from ZINC database against the CagA protein using the bioinformatics tools. Hits were initially chosen based on their physicochemical, absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics, as well as other drug-like characteristics. To locate safe and effective hits, the PAINS filter, binding affinities estimation, and interaction analysis were used. Three compounds with high binding affinity and specificity for the CagA binding pocket were discovered. The final hits, ZINC153731, ZINC69482055, and ZINC164387, were found to bind strongly with CagA protein, with binding energies of −11.53, −10.67, and −9.21 kcal/mol, respectively, which were higher than that of the control compound (−7.25 kcal/mol). Further, based on binding affinity and interaction pattern, two leads (ZINC153731, ZINC69482055) were chosen for molecular dynamics (MD) simulation analysis. MD results showed that they displayed stability in their vicinity at 100 ns. This study suggested that these compounds could be used as possible inhibitors of CagA protein in the fight against GC. However, additional benchwork tests are required to validate them as CagA protein inhibitors.     

Identification of novel Plasmodium falciparum PI4KB inhibitors as potential anti-malarial drugs: Homology modeling,molecular docking and molecular dynamics simulations

The current study was set to discover selective Plasmodium falciparum phosphatidylinositol-4-OH kinase type III beta (pfPI4KB) inhibitors as potential antimalarial agents using combined structure-based and ligand-based drug discovery approach. A comparative model of pfPI4KB was first constructed and validated using molecular docking techniques. Performance of Autodock4.2 and Vina4 software in predicting the inhibitor-PI4KB binding mode and energy was assessed based on two Test Sets: Test Set I contained five ligands with resolved crystal structures with PI4KB, while Test Set II considered eleven compounds with known IC50 value towards PI4KB. The outperformance of Autodock as compared to Vina was reported, giving a correlation coefficient (R²) value of 0.87 and 0.90 for Test Set I and Test Set II, respectively. Pharmacophore-based screening was then conducted to identify drug-like molecules from ZINC database with physicochemical similarity to two potent pfPI4KB inhibitors –namely cpa and cpb. For each query inhibitor, the best 1000 hits in terms of TanimotoCombo scores were selected and subjected to molecular docking and molecular dynamics (MD) calculations. Binding energy was then estimated using molecular mechanics–generalized Born surface area (MM-GBSA) approach over 50 ns MD simulations of the inhibitor-pfPI4KB complexes. According to the calculated MM-GBSA binding energies, ZINC78988474 and ZINC20564116 were identified as potent pfPI4KB inhibitors with binding energies better than those of cpa and cpb, with ΔG_binding ≥ −34.56 kcal/mol. The inhibitor-pfPI4KB interaction and stability were examined over 50 ns MD simulation; as well the selectivity of the identified inhibitors towards pfPI4KB over PI4KB was reported.     

Risedronate and Methotrexate Are High-Affinity Inhibitors of New Delhi Metallo-β-Lactamase-1 (NDM-1): A Drug Repurposing Approach

Bacteria expressing New Delhi metallo-β-lactamase-1 (NDM-1) can hydrolyze β-lactam antibiotics (penicillins, cephalosporins, and carbapenems) and, thus, mediate multidrug resistance. The worldwide dissemination of NDM-1 poses a serious threat to public health, imposing a huge economic burden in the development of new antibiotics. Thus, there is an urgent need for the identification of novel NDM-1 inhibitors from a pool of already-known drug molecules. Here, we screened a library of FDA-approved drugs to identify novel non-β-lactam ring-containing inhibitors of NDM-1 by applying computational as well as in vitro experimental approaches. Different steps of high-throughput virtual screening, molecular docking, molecular dynamics simulation, and enzyme kinetics were performed to identify risedronate and methotrexate as the inhibitors with the most potential. The molecular mechanics/generalized Born surface area (MM/GBSA) and molecular dynamics (MD) simulations showed that both of the compounds (risedronate and methotrexate) formed a stable complex with NDM-1. Furthermore, analyses of the binding pose revealed that risedronate formed two hydrogen bonds and three electrostatic interactions with the catalytic residues of NDM-1. Similarly, methotrexate formed four hydrogen bonds and one electrostatic interaction with NDM-1’s active site residues. The docking scores of risedronate and methotrexate for NDM-1 were –10.543 kcal mol⁻¹ and −10.189 kcal mol⁻¹, respectively. Steady-state enzyme kinetics in the presence of risedronate and methotrexate showed a decreased catalytic efficiency (i.e., kcat/Km) of NDM-1 on various antibiotics, owing to poor catalytic proficiency and affinity. The results were further validated by determining the MICs of imipenem and meropenem in the presence of risedronate and methotrexate. The IC₅₀ values of the identified inhibitors were in the micromolar range. The findings of this study should be helpful in further characterizing the potential of risedronate and methotrexate to treat bacterial infections.     

Theoretical studies on QSAR and mechanism of 2‐indolinone derivatives as tubulin inhibitors

The theoretical studies on three‐dimensional quantitative structure activity relationship (3D‐QSAR) and action mechanism of a series of 2‐indolinone derivatives as tubulin inhibitors against human breast cancer cell line MDA‐MB‐231 have been carried out. The established 3D‐QSAR model from the comparative molecular field analysis (CoMFA) shows not only significant statistical quality but also predictive ability, with high correlation coefficient (R² = 0.986) and cross‐validation coefficient (q² = 0.683). In particular, the appropriate binding orientations and conformations of these 2‐indolinone derivatives interacting with tubulin are located by docking study, and it is very interesting to find that the plot of the energy scores of these compounds in DOCK versus the corresponding experimental pIC₅₀ values exhibits a considerable linear correlation. Therefore, the inhibition mechanism that 2‐indolinone derivatives were regarded as tubulin inhibitors can be theoretically confirmed. Based on such an inhibition mechanism along with 3D‐QSAR results, some important factors improving the activities of these compounds were discussed in detail. These factors can be summarized as follows: the H atom adopted as substituent R₁, the substituent R₂ with higher electropositivity and smaller bulk, the substituents R₄–R₆ (on the phenyl ring) with higher electropositivity and larger bulk, and so on. These results can offer useful theoretical references for understanding the action mechanism, designing more potent inhibitors, and predicting their activities prior to synthesis. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

3D-QSAR study,docking molecular and simulation dynamic on series of benzimidazole derivatives as anti-cancer agents

Benzimidazole is an important heterocyclic organic compound which has a structural analogy to nucleotides found in human body and hence is an important pharmacophore in medicinal chemistry. The anti-cancer activities for a diverse set of benzimidazole as anti-cancer agents against breast cancer cell line (MCF7) assay have been subjected to 3D-QSAR (3-Dimensional Quantitative Structural-Activity Relationship) studies. Both CoMFA and CoMSIA models exhibit significant results in terms of statistical parameters as determination coefficients R² > 0.9 and Leave One Out cross-validation determination coefficients Q²> 6. The predictive quality of both 3D QSAR models have been assessed by external validation and Y-randomization test. Five new compounds have been designed and predicted by in silico ADMET method. In the second part, we have used the docking molecular and simulation dynamics (MD) to investigate the bonding interactions and stability of the designed compounds into the Pin1. Then, we have compared them to Trastuzumab and Tamoxifen as a standard inhibitors drug of breast cancer. The designed compounds form stable hydrogen and hydrophobic bonding interactions with the residues Lys63, Gln131, Ser154, Arg 68 and Arg69 of Pin1 receptor during 100 ns as a time of the simulation. The obtained results showed that the new benzimidazole are useful as a template for future design of more potent inhibitors against breast cancer cell lines (MCF7).     

Chalcone Scaffolds Exhibiting Acetylcholinesterase Enzyme Inhibition: Mechanistic and Computational Investigations

This study was aimed to perform the mechanistic investigations of chalcone scaffold as inhibitors of acetylcholinesterase (AChE) enzyme using molecular docking and molecular dynamics simulation tools. Basic chalcones (C1–C5) were synthesized and their in vitro AChE inhibition was tested. Binding interactions were studied using AutoDock and Surflex-Dock programs, whereas the molecular dynamics simulation studies were performed to check the stability of the ligand–protein complex. Good AChE inhibition (IC₅₀ = 22 ± 2.8 to 37.6 ± 0.75 μM) in correlation with the in silico results (binding energies = −8.55 to −8.14 Kcal/mol) were obtained. The mechanistic studies showed that all of the functionalities present in the chalcone scaffold were involved in binding with the amino acid residues at the binding site through hydrogen bonding, π–π, π–cation, π–sigma, and hydrophobic interactions. Molecular dynamics simulation studies showed the formation of stable complex between the AChE enzyme and C4 ligand.     

Computational Screening of Phenylamino-Phenoxy-Quinoline Derivatives against the Main Protease of SARS-CoV-2 Using Molecular Docking and the ONIOM Method

In the search for new anti-HIV-1 agents, two forms of phenylamino-phenoxy-quinoline derivatives have been synthesized, namely, 2-phenylamino-4-phenoxy-quinoline and 6-phenylamino-4-phenoxy-quinoline. In this study, the binding interactions of phenylamino-phenoxy-quinoline derivatives and six commercially available drugs (hydroxychloroquine, ritonavir, remdesivir, S-217622, N3, and PF-07321332) with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (M^pro) were investigated using molecular docking and the ONIOM method. The molecular docking showed the hydrogen bonding and hydrophobic interactions of all the compounds in the pocket of SARS-CoV-2 main protease (M^pro), which plays an important role for the division and proliferation of the virus into the cell. The binding free energy values between the ligands and M^pro ranged from −7.06 to −10.61 kcal/mol. The molecular docking and ONIOM results suggested that 4-(2′,6′-dimethyl-4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline and 4-(4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline have low binding energy values and appropriate molecular properties; moreover, both compounds could bind to M^pro via hydrogen bonding and Pi-Pi stacking interactions with amino acid residues, namely, HIS41, GLU166, and GLN192. These amino acids are related to the proteolytic cleavage process of the catalytic triad mechanisms. Therefore, this study provides important information for further studies on synthetic quinoline derivatives as antiviral candidates in the treatment of SARS-CoV-2.     

3D QSAR Pharmacophore Based Virtual Screening for Identification of Potential Inhibitors for CDC25B

Owing to its fundamental roles in cell cycle phases, the cell division cycle 25B (CDC25B) was broadly considered as potent clinical drug target for cancers. In this study, 3D QSAR pharmacophore models for CDC25B inhibitors were developed by the module of Hypogen. Three methods (cost analysis, test set prediction, and Fisher’s test) were applied to validate that the models could be used to predict the biological activities of compounds. Subsequently, 26 compounds satisfied Lipinski’s rule of five were obtained by the virtual screening of the Hypo-1-CDC25B against ZINC databases. It was then discovered that 9 identified molecules had better binding affinity than a known CDC25B inhibitors-compound 1 using docking studies. The molecular dynamics simulations showed that the compound had favorable conformations for binding to the CDC25B. Thus, our findings here would be helpful to discover potent lead compounds for the treatment of cancers.