Pharmacophoric Features Research Articles

Discovery of dual-target natural inhibitors of meprins α and β metalloproteases for inflammation regulation: pharmacophore modelling, molecular docking, ADME prediction, and molecular dynamics studies

ABSTRACT Meprins, zinc-dependent metalloproteinases belonging to the metzincin family, have been associated with various inflammatory diseases due to their abnormal expression and activity. In this study, we utilized pharmacophore modelling to identify crucial features for discovering potential dual inhibitors targeting meprins α and β. We screened four pharmacophoric features against a library of 270,540 natural compounds from the Zinc database, resulting in 84,092 matching compounds. Molecular docking was then performed on these compounds, targeting the active sites of meprins α and β. Docking results revealed six compounds capable of interacting with both isoforms, with binding affinities ranging from −10.0 to −10.5 kcal/mol and −6.9 to −9.9 kcal/mol for meprin α and β, respectively. Among these compounds, ZINC000008790788 and ZINC000095099469 displayed superior docking scores and MM-GBSA binding free energy compared to reference ligands. Furthermore, these two compounds exhibited acceptable predicted pharmacokinetic properties and stable interactions with meprins α and β during molecular dynamics simulations. This study presents a comprehensive approach for identifying potential dual inhibitors of meprin α and β, offering insights into the development of therapeutic interventions for inflammatory diseases associated with meprin dysregulation.

Isatin-indoloquinoxaline click adducts with a potential to overcome platinum-based drug-resistance in ovarian cancer

Herein, a series of isatin tethered indolo[2,3-b]quinoxaline hybrids was synthesized by considering the pharmacophoric features of known DNA intercalators and topoisomerase II inhibitors. The anti-proliferative properties of the synthesized compounds were evaluated against ovarian cancer cell lines (SKOV-3 and Hey A8). Four of the compounds exhibited promising anti-proliferative activities, with one of them being 10-fold more potent than cisplatin against drug-resistant Hey A8 cells. Further investigations were carried out to determine the DNA intercalating affinities of the most active compounds as potential mechanisms for their anti-proliferative activities. ADMET in silico studies were performed to assess the physicochemical, pharmacokinetics, and toxicity parameters of active compounds. This study, to the best of our knowledge, is the first report on the potential of isatin-indoloquinoxaline hybrids as structural blueprints for the development of new DNA intercalators. Additionally, it explores their potential to circumvent platinum-based resistance in ovarian cancer.

3D-QSAR and Pharmacophoric study on 2,6-Disubstituted Thiazolo [4,5-b] Pyridines as H3 Receptor Antagonists

Histamine H3 receptor antagonist (H3RA) is a promising therapeutic for CNS disorders including attention deficit hyperactivity disorder [ADHD], sleep disorders, epilepsy, schizophrenia and obesity. 2,6-Disubstituted thiazolo[4,5-b]pyridines reported for their H3 receptor antagonistic activity were selected for three dimensional quantitative structure activity relationship (3D-QSAR) and pharmacophoric study in order to establish structure activity relationship quantitatively and essential structural features. In the current study, VLife Molecular Design Suite software (VlifeMDS) was used for QSAR and biophore studies. Pharmagist (web based server) was used for pharmacophoric study. Partial least square regression (PLSR) analysis showed r2= 0.7902, q2=0.6449 and pred_r2= 0.6650. In this model steric [S_138, S_826] and electrostatic descriptors [E_243, E_652] are involved to play an important role in eliciting biological activity. It showed good internal and external prediction. The contour plots provided further insight of the relationship between structural features of substituted thiazolo[4,5-b]pyridine derivatives and their activities which should be applicable to design newer potential H3R inhibitors. In addition with these studies, pharmacophoric models were also produced using Molsign (VLifeMDS) and Pharmagist (web based server). The identified pharmacophore features are two aromatic and two hydrogen bond acceptor with Molsign whereas common pharmacophoric features with Pharmagist are two aromatic, two hydrophobic and four hydrogen bond acceptors. The present work may be useful for further lead optimization and designing of potent H3 receptor antagonists.

Theoretical investigation on known renin inhibitors and generation of ligand-based pharmacophore models for hypertension treatment

The renin enzyme is considered a promising target for hypertension and renal diseases. Over the last three decades, several experimental and theoretical studies have been engaged in the discovery of potent renin inhibitors. The identified inhibitors that undergo clinical trials are still failing to meet the criteria of potency and safety. To date, there is no specific FDA-approved drug for renin inhibition. Our theoretical opinion describes that the most potent compounds identified in experimental studies but lacking safety and overdose issues could be solved by finding similar molecules that are stable, very active, and have no side effects, which will kick start the drug discovery process. Here, we utilized the most potent direct renin inhibitors reported earlier, followed further by our theoretical study reported in 2019. Ligand-based virtual screening, density functional theory, and dynamic simulation studies were employed to explore the identified compounds and co-crystallized molecule in the protein structure. From the diverse databases, we have identified several identical molecules based on their structural features, such as functional groups like hydrophobic (H), aromatic rings (R), hydrogen bond acceptor (A), and donor (D). The HHHPR five-point pharmacophore feature was identified as a template pharmacophore to search the potential compounds from the Enamine and LifeChemical databases and have a good fitness score with known renin inhibitors. Furthermore, theoretical validation was done through several studies that confirmed the activity of the identified molecules. Overall, we propose that these compounds might break the failure in adverse events and improve the potency of hypertension treatment. Communicated by Ramaswamy H. Sarma

Design, synthesis and preclinical evaluation of muscarine receptor antagonists via a scaffold-hopping approach

Our research group recently identified a rearrangement product of pirenzepine as starting point for a comprehensive rational drug design approach towards orthosteric muscarinic acetylcholine receptor ligands. Chemical reduction and bioscaffold hop lead to the development of sixteen promising compounds featuring either a benzimidazole or carbamate moiety, all exhibiting comparable pharmacophoric characteristics. The synthesized compounds were characterized by NMR, HR-MS, and RP-HPLC techniques. Subsequent evaluation encompassed binding affinity assessment on CHO-hM1-5 cells, mode of action determination, and analysis of physico-chemical parameters. The CNS MPO score indicated favorable drug-like attributes and potential CNS activity for the antagonistic ligands. The most promising compounds displayed Ki-values within a desirable low nanomolar range, and their structural features allow for potential carbon-11 radiolabeling. Our optimization efforts resulted in compounds with a remarkable 138-fold increase in binding affinity compared to the previously mentioned rearrangement product towards human M5, suggesting their prospective utility in positron emission tomography applications.

Computational study on natural molecules targeting β-ketoacyl reductase domain of fatty acid synthase for the identification of selective inhibitors

Fatty acid synthase (FASN) enzyme is a lipid metabolism protein that provides the essential nutrients to cancer cells through de novo lipogenesis. Also, it plays a key role in other disease conditions, including obesity and inflammation. Hence, targeting the β-ketoacyl reductase (KR) domain of FASN protein, an in-silico study was performed on some selective bioactive natural molecules following a repurposing strategy to identify FASN inhibitors. A molecular docking study followed by Absorption, Distribution, Metabolism and Excretion (ADME) predictions, binding free energy calculations, and molecular dynamics (MD) simulations were performed against FASN protein (PDB ID:6NNA) using Schrodinger Drug Discovery Software. Compounds rutin, trans-chlorogenic acid, norbergenin, myricetin, quercetin, physalolactone, quercetin-3-O-galactoside, kaempferol, asperulosidic acid, luteolin, curcumin, 12-deoxywithastramonolide, pedunculoside, hernandifoline, and withafastuosin E were identified as hits, presenting better docking scores (–16.2, −14.1, −12.3, −12.1, −12.0, −11.3, −10.3, −9.8, −9.3, −9.2, −9.1, −8.5, −8.4, −8.3, –7.9, respectively) and hydrogen bond interactions with Ser 2021 and Tyr 2034 amino acids of the KR domain of FASN. The MD simulations study of top five hits in complex with protein 6NNA uncovered the significant interactions leading to the stabilization of ligand with Root Mean Square Deviation (RMSD) below 5.00 Å and the stability was further validated by evaluating the root mean square fluctuation, solvent accessible surface area, and radius of gyration graphs. Also, the FASN inhibition effect of top four hits (50 µM) was >50% when corroborated using High Performance Liquid Chromatography HPLC-based estimation of palmitic acid in MCF-7 cells. Thus rutin, trans-chlorogenic acid, norbergenin, and myricetin are proposed as prospective FASN inhibitory leads. Overall, the results provided useful modifications in pharmacophoric features that could improve the inhibitory effect. Communicated by Ramaswamy H. Sarma

Discovery of 2-Aminothiazole-4-carboxylic Acids as Broad-Spectrum Metallo-β-lactamase Inhibitors by Mimicking Carbapenem Hydrolysate Binding.

Metallo-β-lactamases (MBLs) are zinc-dependent enzymes capable of hydrolyzing all bicyclic β-lactam antibiotics, posing a great threat to public health. However, there are currently no clinically approved MBL inhibitors. Despite variations in their active sites, MBLs share a common catalytic mechanism with carbapenems, forming similar reaction species and hydrolysates. We here report the development of 2-aminothiazole-4-carboxylic acids (AtCs) as broad-spectrum MBL inhibitors by mimicking the anchor pharmacophore features of carbapenem hydrolysate binding. Several AtCs manifested potent activity against B1, B2, and B3 MBLs. Crystallographic analyses revealed a common binding mode of AtCs with B1, B2, and B3 MBLs, resembling binding observed in the MBL-carbapenem product complexes. AtCs restored Meropenem activity against MBL-producing isolates. In the murine sepsis model, AtCs exhibited favorable synergistic efficacy with Meropenem, along with acceptable pharmacokinetics and safety profiles. This work offers promising lead compounds and a structural basis for the development of potential drug candidates to combat MBL-mediated antimicrobial resistance.

A Theobromine Derivative with Anticancer Properties Targeting VEGFR-2: Semisynthesis, in silico and in vitro Studies.

A computer-assisted drug design (CADD) approach was utilized to design a new acetamido-N-(para-fluorophenyl)benzamide) derivative of the naturally occurring alkaloid, theobromine, (T-1-APFPB), following the pharmacophoric features of VEGFR-2 inhibitors. The stability and reactivity of T-1-AFPB were assessed through density functional theory (DFT) calculations. Molecular docking assessments showed T-1-AFPB's potential to bind with and inhibit VEGFR-2. The precise binding of T-1-AFPB against VEGFR-2 with optimal energy was further confirmed through several molecular dynamics (MD) simulations, PLIP, MM-GBSA, and PCA studies. Then, T-1-AFPB (4-(2-(3,7-Dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-1-yl)acetamido)-N-(4-fluorophenyl)benzamide) was semi-synthesized and the in vitro assays showed its potential to inhibit VEGFR-2 with an IC50 value of 69 nM (sorafenib's IC50 was 56 nM) and to inhibit the growth of HepG2 and MCF-7 cancer cell lines with IC50 values of 2.24±0.02 and 3.26±0.02 μM, respectively. Moreover, T-1-AFPB displayed very high selectivity indices against normal Vero cell lines. Furthermore, T-1-AFPB induced early (from 0.72 to 19.12) and late (from 0.13 to 6.37) apoptosis in HepG2 cell lines. In conclusion, the combined computational and experimental approaches demonstrated the efficacy and safety of T-1-APFPB providing it as a promising lead VEGFR-2 inhibitor for further development aiming at cancer therapy.

Design and synthesis of novel pyrazolopyrimidine candidates as promising EGFR-T790M inhibitors and apoptosis inducers.

Aim: Our objective was to design and synthesize a new range of pyrazolopyrimidines while maintaining the key pharmacophoric features of EGFR tyrosine kinase inhibitors. Materials & methods: Percentage inhibition in 14 human cancer cell lines and IC50 values were recorded. Compounds 6c, 7e and 7f were examined against both wild and mutant (T790M) EGFR subtypes. Apoptosis markers, cell cycle arrest, apoptosis assay and molecular docking were performed. Results: Compounds 6c, 7e and 7f demonstrated superiorinhibitory potentials againstwild and mutant (T790M) EGFR subtypes. A molecular docking study showed that compounds 6c and 7e had the best fit. Conclusion: The designed candidatesdemonstrated superior inhibitory potential as promising EGFR-T790M inhibitors that agrees with the proposed rationale.

In Silico Studies of Piperazinyl-4-Nitroimidazole Derivatives in the Non-Nucleoside Inhibitory Binding Pocket of Human Immunodeficiency Virus-1-Reverse Transcriptase Enzyme

Human immunodeficiency virus-1 (HIV-1) reverse transcriptase enzyme catalyzes the conversion of viral RNA to DNA. This viral DNA infects a healthy host body, it leads to the production of viral DNA in higher concentrations and leads to a disease called acquired immunodeficiency syndrome; which is a very fatal disease and leads to the development of several other diseases. Therefore, the development of effective inhibitors is required to execute the therapeutic effects against this disease. In this context, in the present study, the docking studies of compounds showing HIV-1 inhibition were performed in the non-nucleoside reverse transcriptase inhibitory binding pocket of the enzyme HIV-1 reverse transcriptase (protein data bank [PDB] ID-1RT2), using AutoDock 4.2.6 and to identify the important pharmacophoric features required for the development of effective non-nucleoside reverse transcriptase inhibitors, pharmacophore analysis was performed. The findings of this study can be used as a valuable tool for identifying and developing potent and effective inhibitors.

Unraveling the mechanistic details of rho-associated coiled-coil kinase1 (ROCK1) inhibition: Implications from pharmacophore mapping and free energy analysis

Rho-associated coiled-coil kinases (ROCKs) are essential for various cellular processes and play a crucial role in neurological and inflammatory disorders. The two ROCK isoforms (ROCK1 and ROCK2) are involved in the pathogenesis of multiple diseases. Increased levels of ROCK1 are observed in patients with Alzheimer's disease (AD), making it a promising drug target. This work presents the mechanism of inhibitor binding at the ROCK1 active site using a combination of ligand and structure-based drug design techniques. Ligand-based pharmacophore studies were performed on previously reported ROCK1 inhibitors to unravel the 3D steric and electronic features essential for ROCK1 inhibition. Molecular docking and molecular dynamics-based free energy analysis studies were performed to recognize the amino acids assisting in inhibitor binding by stabilizing the identified pharmacophore features. The crucial 3D pharmacophore features of ROCK1 inhibitors included hydrophobic (HYD), hydrogen-bond acceptor (HBA), and hydrogen-bond donor (HBD) characteristics. These features were oriented next to the complementary amino acids in the ROCK1 active site, as evident from the detailed interaction analysis. The free energy calculations provided a quantifiable estimate of the favorable energetic contribution of Gly83, Val90, Lys105, Asp160, Leu205, and Asp216 to the total binding energy. These amino acids significantly contributed to the ROCK1-inhibitor complexes' overall binding free energy, indicating their importance in inhibitor binding.

Identification of potential inhibitors against E.coli via novel approaches based on deep learning and quantum mechanics-based atomistic investigations

Currently, drug resistance to commercially available antibiotics is imparting negative consequences to global health, and the development of novel antibiotics in a timely manner is a prime need of the hour. In the current study, an e-pharmacophore model was built using the 3D structure of DNA gyrase in complex with a standard inhibitor. The generated model was subjected to a pharmacophore based virtual screening against 45,257,086 molecules having 223,460,579 conformers available in MCULE database. Pharmacophore based screening retrieved eight molecules as top hit based on pharmacophoric features in comparison to standard inhibitors. Afterward, all eight compounds were subjected molecular docking based on deep learning algorithm. The molecular docking revealed that compound MCULE-6042843173 and MCULE-2362244223 had significant binding orientation inside active pocket of targeted protein with binding affinity of −9.52 and −9.24 kcal/mol respectively. In addition, density functional theory studies (DFT) were performed to evaluate quantum mechanics of top ranked compounds which were investigated through quantum mechanics (QM) computations which strongly assisted the findings of other in-silico investigations. Consequently, the MCULE-6042843173 and MCULE-2362244223 were subjected to MD simulation studies for evaluation of stability, hydrogen bond analysis, van der Waals interactions, and the contact profile of compounds with targeted amino acid residues. Findings of current study suggested MCULE-6042843173 and MCULE-2362244223 as potential and novel inhibitor of DNA Gyrase enzyme.

Exploring natural products as multi-target-directed drugs for Parkinson’s disease: an in-silico approach integrating QSAR, pharmacophore modeling, and molecular dynamics simulations

Parkinson’s disease is a neurodegenerative disorder characterized by the progressive loss of dopaminergic neurons in the midbrain. Current treatments provide limited symptomatic relief without halting disease progression. A multi-targeting approach has shown potential benefits in treating neurodegenerative diseases. In this study, we employed in silico approaches to explore the COCONUT natural products database and identify novel drug candidates with multi-target potential against relevant Parkinson’s disease targets. QSAR models were developed to screen for potential bioactive molecules, followed by a hybrid virtual screening approach involving pharmacophore modeling and molecular docking against MAO-B, AA2AR, and NMDAR. ADME evaluation was performed to assess drug-like properties. Our findings revealed 22 candidates that exhibited the desired pharmacophoric features. Particularly, two compounds: CNP0121426 and CNP0242698 exhibited remarkable binding affinities, with energies lower than −10 kcal/mol and promising interaction profiles with the chosen targets. Furthermore, all the ligands displayed desirable pharmacokinetic properties for brain-targeted drugs. Lastly, molecular dynamics simulations were conducted on the lead candidates, belonging to the dihydrochalcone and curcuminoid class, to evaluate their stability over a 100 ns timeframe and compare their dynamics with reference complexes. Our findings revealed the curcuminoid CNP0242698 to have an overall better stability with the three targets compared to the dihydrochalcone, despite the high ligand RMSD, the curcuminoid CNP0242698 showed better protein stability, implying ligand exploration of different orientations. Similarly, AA2AR exhibited higher stability with CNP0242698 compared to the reference complex, despite the high initial ligand RMSD due to the bulkier active site. In NMDAR, CNP0242698 displayed good stability and less fluctuations implying a more restricted conformation within the smaller active site of NMDAR. These results may serve as lead compounds for the development and optimization of natural products as multi-target disease-modifying natural remedies for Parkinson’s disease patients. However, experimental assays remain necessary to validate these findings. Communicated by Ramaswamy H. Sarma.

Pharmacophore-based virtual screening, 3D QSAR, Docking, ADMET, and MD simulation studies: An in silico perspective for the identification of new potential HDAC3 inhibitors

Histone deacetylase 3 (HDAC3) is an epigenetic regulator that involves gene expression, apoptosis, and cell cycle progression, and the overexpression of HDAC3 is accountable for several cancers, neurodegeneracy, and many other diseases. Therefore, HDAC3 emerged as a promising drug target for the novel drug design. Here, we carried out the pharmacophore modeling using 50 benzamide-based HDAC3 selective inhibitors and utilized it for PHASE ligand screening to retrieve the hits with similar pharmacophore features. The dataset inhibitors of best hypotheses used to build the 3D QSAR model and the generated 3D QSAR model resulted in good PLS statistics with a regression coefficient (R2) of 0.89, predictive coefficient (Q2) of 0.88, and Pearson-R factor of 0.94 indicating its excellent predictive ability. The hits retrieved from pharmacophore-based virtual screening were subjected to docking against HDAC3 for the identification of potential inhibitors. A total of 10 hitsM1 to M10 were ranked using their scoring functions and further subject to lead optimization. The Prime MM/GBSA, AutoDock binding free energies, and ADMET studies were implemented for the selection of lead candidates. The four ligand molecules M1, M2, M3, and M4 were identified as potential leads against HDAC3 after lead optimization. The top two leads M1 and M2 were subjected to MD simulations for their stability evaluation with HDAC3. The newly designed leads M11 and M12 were identified as HDAC3 potential inhibitors from MD simulations studies. Therefore, the outcomes of the present study could provide insights into the discovery of new potential HDAC3 inhibitors with improved selectivity and activity against a variety of cancers and neurodegenerative diseases.

Pharmacophore derived 3D-QSAR, molecular docking, and simulation studies of quinoxaline derivatives as ALR2 inhibitors

Aldose Reductase 2 (ALR2), a key enzyme of the polyol pathway, plays a crucial role in the pathogenesis of diabetic complications. Quinoxaline scaffold-based compounds have been identified as potential ALR2 inhibitors for the management of diabetic complications. In the present work, molecular dynamic simulation studies in conjugation with pharmacophore mapping and atom-based 3D-QSAR were performed on a dataset of 99 molecules in comparison with Epalrestat (reference) to mark the desirable structural features of quinoxaline analogs to generate a probable template for designing novel and effective ALR2 inhibitors. The most potent compound 81 was subjected to MD simulation studies and found to be stable, with better interactions with the binding pocket as compared to Epalrestat. The MM-GBSA and MM-PBSA calculations showed that compound 81 possessed binding free energies of −35.96 and −4.92 kcal/mol, respectively. Atom-based 3D-QSAR yielded various pharmacophoric features with excellent statistical measures, such as correlation coefficient (R 2 value), F-value (Fischer ratio), Q 2 value (cross-validated correlation coefficient), and Pearson’s R-value for training and test sets. Furthermore, the pharmacophore mapping provided a five-point hypothesis (AADRR) and docking analysis revealed the active ligand-binding orientations on the active site’s amino acid residues TYR 48, HIE 110, TRP 111, and TRP 219. The results of this study will help in designing potent inhibitors of ALR2 for the management of diabetic complications. Communicated by Ramaswamy H. Sarma

Identification of novel scaffolds targeting SIRT3 through molecular modeling techniques for the treatment of Hepatocellular carcinoma

Hepatocellular carcinoma is one of the top causes of cancer-related death globally. SIRT3 belongs to the Sirtuin family of proteins, a collection of NAD+-dependent enzymes that play a role in controlling several cellular functions, including metabolism, aging, and stress response. SIRT3 expression has been discovered to be often downregulated in HCC tissues relative to normal liver tissues. Hence, SIRT3 may function as a tumor suppressor in HCC. In the present study, pharmacophore-based virtual screening of a small molecule database was performed initially, and then the screened hits were docked to the active site of SIRT3 to choose the best binding modes. One co-crystal ligand (PDB name: 1NQ) was utilized as a template to generate pharmacophore model query. A total of 0.2 million compounds from the VITAS-M Lab database were downloaded and prepared for virtual screening. Following database preparation, ligand-based virtual screening was performed using the pharmacophore query model generated in the previous phase. The compounds with the same pharmacophoric characteristics as the query at the same distance were screened. There were a total of 74 hits that matched the query model. These compounds were then docked to the SIRT3 using the standard precision protocol of the glide tool. To select hits with high binding affinities, a threshold of −8 kcal/mol was used. Based on the glide gscore, two hits were chosen. These two hits were selected to investigate the stability of the protein-ligand complex by molecular dynamics simulation. All of these findings indicate that the selected hit compounds C1 and C2 can serve as lead compounds in inhibiting the biological activity of SIRT3 requiring further detailed investigations. Communicated by Ramaswamy H. Sarma

In Silico Modeling and Structural Analysis of Soluble Epoxide Hydrolase Inhibitors for Enhanced Therapeutic Design.

Human soluble epoxide hydrolase (sEH), a dual-functioning homodimeric enzyme with hydrolase and phosphatase activities, is known for its pivotal role in the hydrolysis of epoxyeicosatrienoic acids. Inhibitors targeting sEH have shown promising potential in the treatment of various life-threatening diseases. In this study, we employed a range of in silico modeling approaches to investigate a diverse dataset of structurally distinct sEH inhibitors. Our primary aim was to develop predictive and validated models while gaining insights into the structural requirements necessary for achieving higher inhibitory potential. To accomplish this, we initially calculated molecular descriptors using nine different descriptor-calculating tools, coupled with stochastic and non-stochastic feature selection strategies, to identify the most statistically significant linear 2D-QSAR model. The resulting model highlighted the critical roles played by topological characteristics, 2D pharmacophore features, and specific physicochemical properties in enhancing inhibitory potential. In addition to conventional 2D-QSAR modeling, we implemented the Transformer-CNN methodology to develop QSAR models, enabling us to obtain structural interpretations based on the Layer-wise Relevance Propagation (LRP) algorithm. Moreover, a comprehensive 3D-QSAR analysis provided additional insights into the structural requirements of these compounds as potent sEH inhibitors. To validate the findings from the QSAR modeling studies, we performed molecular dynamics (MD) simulations using selected compounds from the dataset. The simulation results offered crucial insights into receptor-ligand interactions, supporting the predictions obtained from the QSAR models. Collectively, our work serves as an essential guideline for the rational design of novel sEH inhibitors with enhanced therapeutic potential. Importantly, all the in silico studies were performed using open-access tools to ensure reproducibility and accessibility.

Design, semi-synthesis, anti-cancer assessment, docking, MD simulation, and DFT studies of novel theobromine-based derivatives as VEGFR-2 inhibitors and apoptosis inducers

A group of theobromine derivatives was designed based on the key pharmacophoric characteristics of VEGFR-2 inhibitors. HepG2 and MCF-7 cancer cell lines were used to test the obtained compounds for their in vitro anti-proliferative activities. Compound 15 (2-(3,7-Dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-1-yl)-N-(4-(1-(2-(4-hydroxybenzoyl)hydrazono)ethyl) phenyl)acetamide) was the most potent cytotoxic member against MCF-7 (IC50 = 0.42 µM) and HepG2 (IC50 = 0.22 µM). The effectiveness of VEGFR-2 inhibition was assessed for compound 15, and its IC50 value was calculated to be 0.067 µM. Additional cellular mechanistic investigations showed that compound 15 dramatically increased the population of apoptotic HepG2 cells in both early and late apoptosis. The investigation of apoptotic markers confirmed that compound 15 upregulated the levels of BAX (2.26-fold) and downregulated the levels of Bcl-2 (4.4-fold). The molecular docking investigations, MM-GPSA, PLIP studies, and MD simulations validated the potential of compound 15 to be a VEGFR-2 inhibitor. DFT calculations have been completed to comprehend how the electrical charge is distributed within compound 15 and to predict how it would bond to VEGFR-2. Lastly, ADMET prediction showed that the designed members have drug-like characteristics and minimal levels of toxicity. In conclusion, our in vitro and in silico investigations showed that compound 15 exhibited promising apoptotic anticancer potential through the suppression of VEGFR-2.

Novel VEGFR2 inhibitors with thiazoloquinoxaline scaffold targeting hepatocellular carcinoma with lower cardiotoxic impact

Hepatocellular carcinoma (HCC) is a fatal tumor which is usually diagnosed at advanced stage. Molecular targeted drugs were used recently to treat HCC, however, due to serious side effects, mainly cardiotoxicity and emergence of resistance, there is demanding to explore new chemotherapeutics. 10 novel thiazoloquinoxaline derivatives coupled with different sulfonamide moieties 4(a–j) were designed and synthesized fulfilling pharmacophoric features of VEGFR-2 inhibition. Structures of all new compounds were verified via spectral and microanalytical data. After carrying in-vitro VEGFR-2 assay for compounds 4(a–j); sulfapyridine and sulfamethoxazole derivatives 4d and 4f showed potential inhibitory effect [61.04 and 83.35 nM], respectively, comparable to standard sorafenib [51.41 nM]. Both were then further evaluated for their cytocidal activity against HepG2 cell-line and against myocardium cells using H9C2 cell-line. As a result, only sulfapyridine derivative 4d exhibited a significant inhibition of HepG2 cells viability [IC50 = 4.31 μM]. Furthermore, it showed relatively lower cytotoxic impact against normal H9C2 myocardium cells [IC50, 33.47 μM] compared to that of sorafenib [IC50, 98.07 μM]. In-vivo study was carried out to determine myocardium safety of compound 4d on irradiated mice (8 Gy). In-vivo results of sulfapyridine derivative 4d showed normal cardiac enzyme function (CK) and serum catalase activity with significant reductions in LDH, cardiac TNF-α and caspase-9 levels, alongside with its efficacy in suppressing the expression of hepatic VEGF. In conclusion, sulfapyridine derivative 4d could be considered a promising candidate as VEGFR-2 inhibitor with less myocardium side effect.

Insilico molecular modelling to identify PDK-1 targeting agents based on its protein–protein docking interaction

PDK1, an attractive cancer target that downstreams 23 other kinases towards cell growth, survival and metabolism has gaining attention due to allosteric effect of ligands bound to it. Generally, the drug design strategy using pharmacophores is either a single protein structure or ensemble or ligand-based. Apart from these methods, yet another new approach of protein–protein docking with state of art computational tool like Schrodinger Suite to generate pharmacophores based on the interacting partners of the protein is proposed in this work. The structure-based pharmacophoric features were picked up from docking the ten interacting partners of PDK1 and screened against the Enamine libraries containing protein–protein interacting compound collection, advanced, protein mimetic and allosteric compounds. High throughput virtual screening against the PIF pocket of PDK1 yields an indole scaffold. The identified indole derivative is proposed to be a strong activator that binds in the protein–protein interaction site of PDK1 which was further confirmed by molecular metadynamics simulations, free energy surface analysis and MM-GBSA calculations. Thus, the pharmacophores generated by the interacting proteins for PPI can facilitate the virtual screening in structure-based drug discovery of similar therapeutic targets. Communicated by Ramaswamy H. Sarma