Diarylthiazole and diarylimidazole selective COX-1 inhibitor analysis through pharmacophore modeling, virtual screening, and DFT-based approaches

New Thiazole Carboxamide Derivatives as COX Inhibitors: Design, Synthesis, Anticancer Screening, In Silico Molecular Docking, and ADME Profile Studies

The goal of this work was to create and test a new series of thiazole carboxamide derivatives for their cyclooxygenase (COX) suppressor and anticancer effects. The compounds were characterized using 1H, 13C NMR, and HRMS spectrum analysis, and their selectivity toward COX-1 and COX-2 was assessed using an in vitro COX inhibition assay kit. Cytotoxicity was assessed using an MTS assay against a panel of cancer and normal cell lines. The docking studies were aided by the Prime MM-GBSA method for estimating binding affinities. The density functional theory (DFT) analysis was performed to assess compound chemical reactivity, which was calculated by computing the border orbital energy of both HOMO and LUMO orbitals, as well as the HOMO-LUMO energy gap. For ADME-T analysis, the QiKProp module was employed. Furthermore, using human X-ray crystal structures, molecular docking studies were carried out to discover the probable binding patterns of these drugs within both COX-1 and COX-2 isozymes. The results demonstrated that the most effective compound against the COX-1 enzyme was 2b with an IC50 of 0.239 μM. It also showed potent activity against COX-2 with an IC50 value of 0.191 μM and a selectivity ratio of 1.251. The highest selectivity ratio was 2.766 for compound 2a against COX-2 with an IC50 dose of 0.958 μM relating to the celecoxib ratio of 23.8 and its IC50 against COX-2 of 0.002 μM. Compound 2j also showed good selectivity toward COX-2 (1.507) with an IC50 value of 0.957 μM. All compounds showed negligible cytotoxic activity against the evaluated normal cell lines, and the IC50 values were more than 300 μM, except for compound 2b, whose IC50 values were 203.71 ± 1.89 and 116.96 ± 2.05 μM against LX-2 and Hek293t cell lines, respectively. Moreover, compound 2b showed moderate anticancer activity against COLO205 and B16F1 cancer cell lines with IC50 values of 30.79 and 74.15 μM, respectively.

New Insights on the Activity and Selectivity of MAO-B Inhibitors through In Silico Methods

To facilitate the identification of novel MAO-B inhibitors, we elaborated a consolidated computational approach, including a pharmacophoric atom-based 3D quantitative structure-activity relationship (QSAR) model, activity cliffs, fingerprint, and molecular docking analysis on a dataset of 126 molecules. An AAHR.2 hypothesis with two hydrogen bond acceptors (A), one hydrophobic (H), and one aromatic ring (R) supplied a statistically significant 3D QSAR model reflected by the parameters: R² = 0.900 (training set); Q² = 0.774 and Pearson's R = 0.884 (test set), stability s = 0.736. Hydrophobic and electron-withdrawing fields portrayed the relationships between structural characteristics and inhibitory activity. The quinolin-2-one scaffold has a key role in selectivity towards MAO-B with an AUC of 0.962, as retrieved by ECFP4 analysis. Two activity cliffs showing meaningful potency variation in the MAO-B chemical space were observed. The docking study revealed interactions with crucial residues TYR:435, TYR:326, CYS:172, and GLN:206 responsible for MAO-B activity. Molecular docking is in consensus with and complementary to pharmacophoric 3D QSAR, ECFP4, and MM-GBSA analysis. The computational scenario provided here will assist chemists in quickly designing and predicting new potent and selective candidates as MAO-B inhibitors for MAO-B-driven diseases. This approach can also be used to identify MAO-B inhibitors from other libraries or screen top molecules for other targets involved in suitable diseases.

Discovery of novel ALK2 inhibitors of pyrazolo-pyrimidines: A computational study

ALK2 is a serine/threonine kinase, involved in different signaling pathways and associated with cell proliferation and differentiation. The present study includes development of pharmacophore, 3-D QSAR, docking and virtual screening studies on 30 different pyrazolo[1,5-a]pyrimidine derivatives. The pharmacophore study provides ARRR_2 hypothesis with four different features essential for ALK2 kinase inhibitory activity. The 3 D-QSAR study determined the statistically significant model by using partial least-square regression (PLS) method with R² value of 0.9711 and Q² value of 0.6846. Validation of 3 D-QSAR has been performed by LOO cross-validation method where with R²CV value of 0.56. The virtual screening study on ZINC database provides compounds such as ZINC66091638, ZINC43524105, ZINC19458227 and ZINC72441013 involved good binding interactions (docking scores -8.91, -7.40, -8.43, and -9.47, respectively) with ALK2 kinase (PDB ID: 3Q4U). The docking study of pyrazolo-pyrimidines derivatives found potent compounds, 7i, 13r, 13d, and 21 with docking scores -9.83, -9.75, -9.76, and -9.75, respectively. The important interactions with amino acid residues were HID 286, ASN341. ADME properties further assist to provide important structural features of ALK2 kinase. The present study may be help to medicinal scientists in the direction to develop potent inhibitors against ALK2 kinase.Communicated by Ramaswamy H. Sarma.

Anticancer Activity of Thiophene Carboxamide Derivatives as CA-4 Biomimetics: Synthesis, Biological Potency, 3D Spheroid Model, and Molecular Dynamics Simulation

The present study aimed to synthesize thiophene carboxamide derivatives, which are considered biomimetics of the anticancer medication Combretastatin A-4 (CA-4), and compare the similarity in the polar surface area (PSA) between the novel series and CA-4. Our results showed that the PSA of the most synthesized structures was biomimetic to CA-4, and similar chemical and biological properties were observed against Hep3B cancer cell line. Among the synthesized series 2b and 2e compounds were the most active molecules on Hep3B (IC50 = 5.46 and 12.58 µM, respectively). The 3D results revealed that both 2b and 2e structures confuse the surface of Hep3B cancer cell lines' spheroid formation and force these cells to aggregate into a globular-shaped spheroid. The 2b and 2e showed a comparable interaction pattern to that observed for CA-4 and colchicine within the tubulin-colchicine-binding pocket. The thiophene ring, due to holding a high aromaticity character, participated critically in that observed interaction profile and showed additional advanced interactions over CA-4. The 2b and 2e tubulin complexes showed optimal dynamics trajectories within a time scale of 100 ns at 300 K temperature, which asserts their high stability and compactness. Together, these findings revealed the biomimetic role of 2b and 2e compounds in CA-4 in preventing cancer progression.

Natural Compounds as DPP-4 Inhibitors: 3D-Similarity Search, ADME Toxicity, and Molecular Docking Approaches

Type 2 diabetes mellitus is one of the most common diseases of the 21st century, caused by a sedentary lifestyle, poor diet, high blood pressure, family history, and obesity. To date, there are no known complete cures for type 2 diabetes. To identify bioactive natural products (NPs) to manage type 2 diabetes, the NPs from the ZINC15 database (ZINC-NPs DB) were screened using a 3D shape similarity search, molecular docking approaches, and ADMETox approaches. Frequently, in silico studies result in asymmetric structures as “hit” molecules. Therefore, the asymmetrical FDA-approved diabetes drugs linagliptin (8-[(3R)-3-aminopiperidin-1-yl]-7-but-2-ynyl-3-methyl-1-[(4-methylquinazolin-2-yl)methyl]purine-2,6-dione), sitagliptin ((3R)-3-amino-1-[3-(trifluoromethyl)-6,8-dihydro-5H-[1,2,4]triazolo[4,3-a]pyrazin-7-yl]-4-(2,4,5-trifluorophenyl)butan-1-one), and alogliptin (2-[[6-[(3R)-3-aminopiperidin-1-yl]-3-methyl-2,4-dioxopyrimidin-1-yl]methyl]benzonitrile) were used as queries to virtually screen the ZINC-NPs DB and detect novel potential dipeptidyl peptidase-4 (DPP-4) inhibitors. The most promising NPs, characterized by the best sets of similarity and ADMETox features, were used during the molecular docking stage. The results highlight that 11 asymmetrical NPs out of 224,205 NPs are potential DPP-4 candidates from natural sources and deserve consideration for further in vitro/in vivo tests.

Homology Modeling and Molecular Docking Approaches for the Proposal of Novel Insecticides against the African Malaria Mosquito (Anopheles gambiae)

Vector-borne infectious diseases are responsible for the deaths of over 700,000 people annually, than 400,000 of them resulting from malaria. The mosquito Anopheles gambiae is one of the dominant vector species of human malaria transmission. A significant issue of the conventional insecticides which target the arthropod borne infectious diseases is their induced resistance. To overcome this inconvenience, insecticides with new modes of action are required. One of the most promising targets for the development of new potential insecticides as evidenced by current studies is the D1-like dopamine receptor (DAR). To get a deeper understanding of the structural information of this receptor, the 3D homology model was built. The possible sites within the protein were identified and the most probable binding site was highlighted. The homology model along with a series of DAR antagonists with known activity against Anopheles gambiae larvae were used in docking experiments to gain insight into their intermolecular interactions. Furthermore, virtual screening of the natural compounds from the SPECS database led to the prediction of toxicity and environmental hazards for one potential new insecticide against the Anopheles gambiae mosquito.

Development of antidiabetic drugs from benzamide derivatives as glucokinase activator: A computational approach

Hyperglycemia is a condition known for the impairment of insulin secretion and is responsible for diabetes mellitus. Various small molecule inhibitors have been discovered as glucokinase activators. Recent studies on benzamide derivatives showed their importance in the treatment of diabetes as glucokinase activator. The present manuscript showed a computation study on benzamide derivatives to help in the production of potent glucokinase activators. In the present study, pharmacophore development, 3D-QSAR, and docking studies were performed on benzamide derivatives to find out the important features required for the development of a potential glucokinase activator. The generated pharmacophore hypothesis ADRR_1 consisted of essential features required for the activity. The resultant statistical data showed high significant values with R² > 0.99; 0.98 for the training set and Q² > 0.52; 0.71 for test set based on atom-based and field-based models, respectively. The potent compound 15b of the series showed a good docking score via binding with different amino acid residues such as (NH…ARG63), (SO₂…ARG250, THR65), and π-π staking with (phenyl……TYR214). The virtual screening study used 3563 compounds from ZINC database and screened hit compound ZINC08974524, binds with similar amino acids as shown by compound 15b and crystal ligand with docking scores SP (-11.17 kcal/mol) and XP (-8.43 kcal/mol). Compounds were further evaluated by ADME and MMGBSA parameters. Ligands and ZINC hits showed no violation of Lipinski rules. All the screened compounds showed good synthetic accessibility. The present study may be used by researchers for the development of novel benzamide derivatives as glucokinase activator.

Computational Approaches for the Design of Novel Anticancer Compounds Based on Pyrazolo[3,4-d]pyrimidine Derivatives as TRAP1 Inhibitor

In the present in-silico study, various computational techniques were applied to determine potent compounds against TRAP1 kinase. The pharmacophore hypothesis DHHRR_1 consists of important features required for activity. The 3D QSAR study showed a statistically significant model with R² = 0.96 and Q² = 0.57. Leave one out (LOO) cross-validation (R² CV = 0.58) was used to validate the QSAR model. The molecular docking study showed maximum XP docking scores (−11.265, −10.532, −10.422, −10.827, −10.753 kcal/mol) for potent pyrazole analogs (42, 46, 49, 56, 43), respectively, with significant interactions with amino acid residues (ASP 594, CYS 532, PHE 583, SER 536) against TRAP1 kinase receptors (PDB ID: 5Y3N). Furthermore, the docking results were validated using the 100 ns MD simulations performed for the selected five docked complexes. The selected inhibitors showed relatively higher binding affinities than the TRAP1 inhibitor molecules present in the literature. The ZINC database was used for a virtual screening study that screened ZINC05297837, ZINC05434822, and ZINC72286418, which showed similar binding interactions to those shown by potent ligands. Absorption, distribution, metabolism, and excretion (ADME) analysis showed noticeable results. The results of the study may be helpful for the further development of potent TRAP1 inhibitors

In silico guided designing of 4-(1H-benzo[d]imidazol-2-yl)phenol-based mutual-prodrugs of NSAIDs: synthesis and biological evaluation

The free COOH group of conventional NSAIDs is a structural feature for non-selective cyclooxygenase (COX) inhibition and the molecular cause of their gastrointestinal (GI) toxicity. In this context, an in house database of synthesizable ester prodrugs of some well-known NSAIDs was developed by combining their -COOH group with -OH of a newly identified antioxidant 4-(1H-benzo[d]imidazol-2-yl)phenol (BZ). The antioxidant potential of BZ was unveiled through in silico PASS prediction and in vitro/in vivo evaluation. The in house database of NSAIDs-BZ prodrugs was first subjected to screening with our previously reported pharmacophore models of hCES1 (AAHRR.430) and hCES2 (AHHR.21) for determining hydrolytic susceptibility. Biotransformation behaviour of screened prodrugs was then assessed by using QM/MM and sterimol parameterization, followed by ADMET calculations to predict the drug likeness. On the basis of in silico results, five prodrugs were duly synthesized and the best three were subject to the in vivo evaluation for their anti-inflammatory, analgesic, antioxidant activities, and ulcerogenic index. Among these prodrugs, BN₂ and BN₅ displayed better anti-inflammatory and analgesics potential in comparison to their parent drugs. All the prodrugs were found to be gastro sparing in the rat model and significantly improved the levels of oxidative stress biomarkers in both blood plasma as well as gastric homogenate.