Binding Free Energy (BFE) Calculations and Quantitative Structure-Activity Relationship (QSAR) Analysis of Schistosoma mansoni Histone Deacetylase 8 (smHDAC8) Inhibitors

Computational discovery of dual potential inhibitors of SARS-CoV-2 spike/ACE2 and Mpro: 3D-pharmacophore, docking-based virtual screening, quantum mechanics and molecular dynamics

To find drugs against COVID-19, caused by the SARS-CoV-2, promising targets include the fusion of the viral spike with the human angiotensin-converting enzyme 2 (ACE2) as well as the main protease (Mpro). These proteins are responsible for viral entry and replication, respectively. We combined several state-of-the-art computational methods, including, protein-ligand interaction fingerprint, 3D-pharmacophores, molecular-docking, MM-GBSA, DFT, and MD simulations to explore two databases: ChEMBL and NANPDB to identify molecules that could both block spike/ACE2 fusion and inhibit Mpro. A total of 1,690,649 compounds from the two databases were screened using the pharmacophore model obtained from PLIF analysis. Five recent complexes of Mpro co-crystallized with different ligands were used to generate the pharmacophore model, allowing 4,829 compounds that passed this prefilter. These were then submitted to molecular docking against Mpro. The 5% top-ranked docking hits from docking result having scores < -8.32 kcal mol-1 were selected and then docked against spike/ACE2. Only four compounds: ChEMBL244958, ChEMBL266531, ChEMBL3680003, and 1-methoxy-3-indolymethyl glucosinolate (4) displayed binding energies < - 8.21 kcal mol-1 (for the native ligand) were considered as putative dual-target inhibitors. Furthermore, predictive ADMET, MM-GBSA and DFT/6-311G(d,p) were performed on these compounds and compared with those of well-known antivirals. DFT calculations showed that ChEMBL244958 and compound 4 had significant predicted reactivity values. Molecular dynamics simulations of the docked complexes were run for 100 ns and used to validate the stability docked poses and to confirm that these hits are putative dual binders of the spike/ACE2 and the Mpro.

5-chloro-3-(2-(2,4-dinitrophenyl) hydrazono)indolin-2-one: synthesis, characterization, biochemical and computational screening against SARS-CoV-2

Chemical prototypes with broad-spectrum antiviral activity are important toward developing new therapies that can act on both existing and emerging viruses. Binding of the SARS-CoV-2 spike protein to the host angiotensin-converting enzyme 2 (ACE2) receptor is required for cellular entry of SARS-CoV-2. Toward identifying new chemical leads that can disrupt this interaction, including in the presence of SARS-CoV-2 adaptive mutations found in variants like omicron that can circumvent vaccine, immune, and therapeutic antibody responses, we synthesized 5-chloro-3-(2-(2,4-dinitrophenyl)hydrazono)indolin-2-one (H2L) from the condensation reaction of 5-chloroisatin and 2,4-dinitrophenylhydrazine in good yield. H2L was characterised by elemental and spectral (IR, electronic, Mass) analyses. The NMR spectrum of H2L indicated a keto-enol tautomerism, with the keto form being more abundant in solution. H2L was found to selectively interfere with binding of the SARS-CoV-2 spike receptor-binding domain (RBD) to the host angiotensin-converting enzyme 2 receptor with a 50% inhibitory concentration (IC50) of 0.26 μM, compared to an unrelated PD-1/PD-L1 ligand-receptor-binding pair with an IC50 of 2.06 μM in vitro (Selectivity index = 7.9). Molecular docking studies revealed that the synthesized ligand preferentially binds within the ACE2 receptor-binding site in a region distinct from where spike mutations in SARS-CoV-2 variants occur. Consistent with these models, H2L was able to disrupt ACE2 interactions with the RBDs from beta, delta, lambda, and omicron variants with similar activities. These studies indicate that H2L-derived compounds are potential inhibitors of multiple SARS-CoV-2 variants, including those capable of circumventing vaccine and immune responses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11696-023-03274-5.

Design, synthesis, and biochemical and computational screening of novel oxindole derivatives as inhibitors of Aurora A kinase and SARS-CoV-2 spike/host ACE2 interaction

Isatin (indol-2,3-dione), a secondary metabolite of tryptophan, has been used as the core structure to design several compounds that have been tested and identified as potent inhibitors of apoptosis, potential antitumor agents, anticonvulsants, and antiviral agents. In this work, several analogs of isatin hybrids have been synthesized and characterized, and their activities were established as inhibitors of both Aurora A kinase and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike/host angiotensin-converting enzyme II (ACE2) interactions. Amongst the synthesized isatin hybrids, compounds 6a, 6f, 6g, and 6m exhibited Aurora A kinase inhibitory activities (with IC50 values < 5 μ M), with GScore values of -7.9, -7.6, -8.2 and -7.7 kcal/mol, respectively. Compounds 6g and 6i showed activities in blocking SARS-CoV-2 spike/ACE2 binding (with IC50 values in the range < 30 μ M), with GScore values of -6.4 and -6.6 kcal/mol, respectively. Compounds 6f, 6g, and 6i were both capable of inhibiting spike/ACE2 binding and blocking Aurora A kinase. Pharmacophore profiling indicated that compound 6g tightly fits Aurora A kinase and SARS-CoV-2 pharmacophores, while 6d fits SARS-CoV-2 and 6l fits Aurora A kinase pharmacophore. This work is a proof of concept that some existing cancer drugs may possess antiviral properties. Molecular modeling showed that the active compound for each protein adopted different binding modes, hence interacting with a different set of amino acid residues in the binding site. The weaker activities against spike/ACE2 could be explained by the small sizes of the ligands that fail to address the important interactions for binding to the ACE2 receptor site.

Discovery of AcrAB-TolC pump inhibitors: Virtual screening and molecular dynamics simulation approach

AcrAB-TolC tripartite efflux pump, which belongs to the RND superfamily, is a main multi-drug efflux system of Escherichia coli (E. coli) because of the broad resistance on various antibiotics. With the discovering of efflux pump inhibitors (EPIs), a combination between these and antibiotics is one of the most promising therapies. Therefore, building a virtual screening model with prediction capacities for the efflux pump inhibitory activities of candidates from DrugBank and ZINC15 dataset, is one of the key goals of this project. Based on the database of 170 diverse chemical structures collected from 28 research journals, two 2D-QSAR models and a 3D-pharmacophore model have been performed. On the AcrB protein (PDB 4DX7), two binding sites have been discovered that match to the hydrophobic trap in the distal pocket and the switch loop in the proximal pocket. After virtual screening processes, twenty candidate AcrAB-TolC inhibitors have been subjected to molecular dynamics simulations, binding free energy calculations and ADMET predictions. The results indicate that three compounds namely DB09233, DB02581, and DB15224 are potential inhibitors with ΔGbind of -42.30 ± 4.58, -40.76 ± 7.30 and -31.06 ± 7.63 kcal.mol-1, respectively.Communicated by Ramaswamy H. Sarma.

QSAR, simulation techniques, and ADMET/pharmacokinetics assessment of a set of compounds that target MAO-B as anti-Alzheimer agent

Background

Alzheimer’s disease (AD), the most common cause of dementia in the elderly, is a progressive neurodegenerative disorder that gradually affects cognitive function and eventually causes death. Most approved drugs can only treat the disease alleviating the disease symptoms; therefore, there is a need to develop drugs that can treat this illness holistically. The medical community is searching for new drugs and new drug targets to cure this disease. In this study, QSAR, molecular docking evaluation, and ADMET/pharmacokinetics assessment were used as modeling methods to identify the compounds with outstanding physicochemical properties.

Results

The 37 MAO-B compounds were screened using the aforementioned methods and yielded a model with the following molecular properties: AATS1v, AATS3v, GATS4m, and GATS6e. Good statistical values were R2train = 0.69, R2adj = 0.63, R2pred = 0.57, LOF = 0.23, and RMSE = 0.38. The model was validated using an evaluation set that confirmed its robustness. The molecular docking was also utilized using crystal structure of human monoamine oxidase B in complex with chlorophenylchromone-carboxamide with ID code of 6FW0, and three compounds were identified with outstanding high binding affinity (13 = − 30.51 kcal mol−1, 31 = − 31.85 kcal mol−1, and 33 = − 33.70 kcal mol−1), and better than the Eldepryl (referenced) drug (− 11.40 kcal mol−1).

Conclusions

These three compounds (13, 31, and 33) were analyzed for ADMET/pharmacokinetics evaluation and found worthy of further analysis as promising drug candidates to cure AD and could also serve as a template to design several monoamine oxidase B inhibitors in the future to cure AD.

Docking, Binding Free Energy Calculations and In Vitro Characterization of Pyrazine Linked 2-Aminobenzamides as Novel Class I Histone Deacetylase (HDAC) Inhibitors

Class I histone deacetylases, HDAC1, HDAC2, and HDAC3, represent potential targets for cancer treatment. However, the development of isoform-selective drugs for these enzymes remains challenging due to their high sequence and structural similarity. In the current study, we applied a computational approach to predict the selectivity profile of developed inhibitors. Molecular docking followed by MD simulation and calculation of binding free energy was performed for a dataset of 2-aminobenzamides comprising 30 previously developed inhibitors. For each HDAC isoform, a significant correlation was found between the binding free energy values and in vitro inhibitory activities. The predictive accuracy and reliability of the best preforming models were assessed on an external test set of newly designed and synthesized inhibitors. The developed binding free-energy models are cost-effective methods and help to reduce the time required to prioritize compounds for further studies.

Histone Deacetylase (HDAC) Inhibitors for the Treatment of Schistosomiasis

Schistosomiasis is a major neglected parasitic disease that affects more than 240 million people worldwide and for which the control strategy consists of mass treatment with the only available drug, praziquantel. Schistosomes display morphologically distinct stages during their life cycle and the transformations between stages are controlled by epigenetic mechanisms. The targeting of epigenetic actors might therefore represent the parasites’ Achilles’ heel. Specifically, histone deacetylases have been recently characterized as drug targets for the treatment of schistosomiasis. This review focuses on the recent development of inhibitors for schistosome histone deacetylases. In particular, advances in the development of inhibitors of Schistosoma mansoni histone deacetylase 8 have indicated that targeting this enzyme is a promising approach for the treatment of this infection.

Synthesis, structure-activity relationships, cocrystallization and cellular characterization of novel smHDAC8 inhibitors for the treatment of schistosomiasis

Schistosomiasis is a major neglected parasitic disease that affects more than 265 million people worldwide and for which the control strategy consists of mass treatment with the only available drug, praziquantel. In this study, we chemically optimized our previously reported benzhydroxamate-based inhibitors of Schistosoma mansoni histone deacetylase 8 (smHDAC8). Crystallographic analysis provided insights into the inhibition mode of smHDAC8 activity by the highly potent inhibitor 5o. Structure-based optimization of the novel inhibitors was carried out using the available crystal structures as well as docking studies on smHDAC8. The compounds were evaluated in screens for inhibitory activity against schistosome and human HDACs (hHDAC). The in vitro and docking results were used for detailed structure activity relationships. The synthesized compounds were further investigated for their lethality against the schistosome larval stage using a fluorescence-based assay. The most promising inhibitor 5o showed significant dose-dependent killing of the schistosome larvae and markedly impaired egg laying of adult worm pairs maintained in culture.