Optimizing Bedaquiline for cardiotoxicity by structure based virtual screening, DFT analysis and molecular dynamic simulation studies to identify selective MDR-TB inhibitors

Screening of Phytocompounds for Identification of Prospective Histone Deacetylase 1 (HDAC1) Inhibitor: An In Silico Molecular Docking, Molecular Dynamics Simulation, and MM-GBSA Approach

The upregulation of HDAC1 facilitate the induction of epigenetic repression of genes responsible for suppressing tumourigenesis, thereby triggering the development of cancer. HDAC1 inhibitors have thus emerged as possible therapeutic approaches against a variety of human malignancies, as they can inhibit the activity of certain HDACs, repair the overexpression of tumour suppressor genes, and induce cell differentiation, cell cycle arrest, and apoptosis. In this study, among 810 virtually screened compounds, Pinocembrin (PHUB000396) had a significant binding affinity (-7.99 kcal/mol). In molecular dynamics simulation (MD) studies for 200 ns time scale, the compound Pinocembrin effectively undergoes conformational optimization, thereby enabling its accommodation within the active site of the receptor. This outcome serves as a rational for the observed binding affinity. The optimal binding free energy calculations using the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) (-35.86 ± 7.52 kcal/mol) showed the significant role of van der Waals forces and Coulomb interactions in the stability of the respective complex. The pharmacokinetic study showed its potential as a lead compound. The in-silico cytotoxicity prediction also confirmed its potential as an active anticancer phytocompound in lung and brain cancer. Therefore, it can be predicted that Pinocembrin could be a useful bioactive compound as an HDAC1 inhibitor and could be used in developing epigenetic therapy in cancer such as brain cancer and lung cancer to regulate gene expression.

Newly discovered clouting interplay between matrix metalloproteinases structures and novel quaternary Ammonium K21: computational and in-vivo testing

AIMS AND OBJECTIVES

To analyze anti-MMP mode of action of Quaternary Ammonium Silane (QAS, codenamed as k21) by binding onto specific MMP site using computational molecular simulation and Anti-Sortase A (SrtA) mode of action by binding onto specific site using computational molecular simulation.

MATERIALS AND METHODS

In silico Molecular Dynamics (MD) was used to determine the interactions of K21 inside the pocket of the targeted protein (crystal structure of fibroblast collagenase-1 complexed to a diphenyl-ether sulphone based hydroxamic acid; PDB ID: 966C; Crystal structure of MMP-2 active site mutant in complex with APP-derived decapeptide inhibitor. MD simulations were accomplished with the Desmond package in Schrödinger Drug Discovery Suite. Blood samples (~ 0.5 mL) collected into K2EDTA were immediately transferred for further processing using the Litron MicroFlow® PLUS micronucleus analysis kit for mouse blood according to the manufacturer's instructions. Bacterial Reverse Mutation Test of K21 Molecule was performed to evaluate K21 and any possible metabolites for their potential to induce point mutations in amino acid-requiring strains of Escherichia coli (E. coli) (WP2 uvrA (tryptophan-deficient)).

RESULTS

Molecular Simulation depicted that K21 has a specific pocket binding on various MMPs and SrtA surfaces producing a classical clouting effect. K21 did not induce micronuclei, which are the result of chromosomal damage or damage to the mitotic apparatus, in the peripheral blood reticulocytes of male and female CD-1 mice when administered by oral gavage up to the maximum recommended dose of 2000 mg/kg. The test item, K21, was not mutagenic to Salmonella typhimurium (S. typhimurium) strains TA98, TA100, TA1535 and TA1537 and E. coli strain WP2 uvrA in the absence and presence of metabolic activation when tested up to the limit of cytotoxicity or solubility under the conditions of the test.

CONCLUSION

K21 could serve as a potent protease inhibitor maintaining the physical and biochemical properties of dental structures.

Application of temperature-dependent and steered molecular dynamics simulation to screen anti-dengue compounds against Marburg virus

Marburg virus infections are extremely fatal with a fatality range of 23% to 90%, therefore there is an urgent requirement to design and develop efficient therapeutic molecules. Here, a comprehensive temperature-dependent molecular dynamics (MD) simulation method was implemented to identify the potential molecule from the anti-dengue compound library that can inhibit the function of the VP24 protein of Marburg. Virtual high throughput screening identified five effective binders of VP24 after screening 484 anti-dengue compounds. These compounds were treated in MD simulation at four different temperatures: 300, 340, 380, and 420 K. Higher temperatures showed dissociation of hit compounds from the protein. Further, triplicates of 100 ns MD simulation were conducted which showed that compounds ID = 118717693, and ID = 5361 showed strong stability with the protein molecule. These compounds were further validated using Δ G binding free energies and they showed: -30.38 kcal/mol, and -67.83 kcal/mol binding free energies, respectively. Later, these two compounds were used in steered MD simulation to detect its dissociation. Compound ID = 5361 showed the maximum pulling force of 199.02 kcal/mol/nm to dissociate the protein-ligand complex while ID = 118717693 had a pulling force of 101.11 kcal/mol/nm, respectively. This ligand highest number of hydrogen bonds with varying occupancies at 89.93%, 69.80%, 57.93%, 52.33%, and 50.63%. This study showed that ID = 5361 can bind with the VP24 strongly and has the potential to inhibit its function which can be validated in the in-vitro experiment.Communicated by Ramaswamy H. Sarma.

Synthesis and Evaluation of Biological Activities for a Novel 1,2,3,4-Tetrahydroisoquinoline Conjugate with Dipeptide Derivatives: Insights from Molecular Docking and Molecular Dynamics Simulations

Peptide synthesis has opened new frontiers in the quest for bioactive molecules with limitless biological applications. This study presents the synthesis of a series of novel isoquinoline dipeptides using advanced spectroscopic techniques for characterization. These compounds were designed with the goal of discovering unexplored biological activities that could contribute to the development of novel pharmaceuticals. We evaluated the biological activities of novel compounds including their antimicrobial, antibacterial, and antifungal properties. The results show promising activity against Escherichia coli and potent antibacterial activity against MTCC 443 and MTCC 1688. Furthermore, these compounds demonstrate strong antifungal activity, outperforming existing standard drugs. Computational binding affinity studies of tetrahydroisoquinoline-conjugated dipeptides against E. coli DNA gyrase displayed significant binding interactions and binding affinity, which are reflected in antimicrobial activities of compounds. Our integrative significant molecular findings from both wet and dry laboratories would help pave a path for the development of antimicrobial therapeutics. The findings suggest that these isoquinoline-conjugated dipeptides could be excellent candidates for drug development, with potential applications in the fight against bacterial and fungal infections. This research represents an exciting step forward in the field of peptide synthesis and its potential to discover novel bioactive molecules with significant implications for human health.

Unveiling the antiviral potential of Plant compounds from the Meliaceae family against the Zika virus through QSAR modeling and MD simulation analysis

Zika virus (ZIKV) is a flavivirus transmitted by mosquitoes, causing neurological disorders and congenital malformations. RNA-dependent RNA polymerase (RdRp) is one of its essential enzymes and a promising drug target for antiviral therapy due to its involvement in the growth and multiplication of the virus. In this study, we conducted a QSAR-based chemical library screening from the Meliaceae family to identify potential RdRp inhibitors. The QSAR model was built using the known inhibitors of RdRp NS5 of ZIKV and their biological activity (EC50), along with the structural and chemical characteristics of the compounds. The top two hit compounds were selected from QSAR screening for further analysis using molecular docking to evaluate their binding energies and intermolecular interactions with RdRp, including the critical residue Trp485. Furthermore, molecular dynamics (MD) simulations were performed to evaluate their binding stability and flexibility upon binding to RdRp. The MD results showed that the selected compounds formed stable complexes with RdRp, and their binding interactions were similar to those observed for the native ligand. The binding energies of the top two hits (-8.6 and -7.7 kcal/mole) were comparable to those of previously reported ZIKV RdRp inhibitors (-8.9 kcal/mole). The compound IMPHY009135 showed the strongest binding affinity with RdRp, forming multiple hydrogen bonds and hydrophobic interactions with key residues. However, compound IMPHY009276 showed the most stable and consistent RMSD, which was similar to the native ligand. Our findings suggest that IMPHY009135 and IMPHY009276 are potential lead compounds for developing novel antiviral agents against ZIKV.Communicated by Ramaswamy H. Sarma.

Investigating the Mechanism of Action of Anti-Dengue Compounds as Potential Binders of Zika Virus RNA-Dependent RNA Polymerase

The World Health Organization (WHO) has designated the Zika virus (ZIKV) as a significant risk to the general public's health. Currently, there are no vaccinations or medications available to treat or prevent infection with the Zika virus. Thus, it is urgently required to develop a highly efficient therapeutic molecule. In the presented study, a computationally intensive search was carried out to identify potent compounds that have the potential to bind and block the activity of ZIKV NS5 RNA-dependent RNA polymerase (RdRp). The anti-dengue chemical library was subjected to high-throughput virtual screening and MM/GBSA analysis in order to rate the potential candidates. The top three compounds were then chosen. According to the MM/GBSA analysis, compound 127042987 from the database had the highest binding affinity to the protein with a minimum binding free energy of -77.16 kcal/mole. Compound 127042987 had the most stable RMSD trend and the greatest number of hydrogen bond interactions when these chemical complexes were evaluated further under a 100 ns molecular dynamics simulation. Compound 127042987 displayed the best binding free energy (GBind) of -96.50 kcal/mol, surpassing the native ligand binding energy (-66.17 kcal/mole). Thereafter, an MM/GBSA binding free energy study was conducted to validate the stability of selected chemical complexes. Overall, this study illustrated that compound 127042987 showed preferred binding free energies, suggesting a possible inhibitory mechanism against ZIKV-RdRp. As per this study, it was proposed that compound 127042987 could be used as a therapeutic option to prevent Zika virus infection. These compounds need to be tested in experiments for further validation.

Isoflavonoid and Furanochromone Natural Products as Potential DNA Gyrase Inhibitors: Computational, Spectral, and Antimycobacterial Studies

In pursuit of new antitubercular agents, we here report the antimycobacterial (H₃₇Rv) and DNA gyrase inhibitory potential of daidzein and khellin natural products (NPs). We procured a total of 16 NPs based on their pharmacophoric similarities with known antimycobacterial compounds. The H₃₇Rv strain of M. tuberculosis was found to be susceptible to only two out of the 16 NPs procured; specifically, daidzein and khellin each exhibited an MIC of 25 μg/mL. Moreover, daidzein and khellin inhibited the DNA gyrase enzyme with IC₅₀ values of 0.042 and 0.822 μg/mL, respectively, compared to ciprofloxacin with an IC₅₀ value of 0.018 μg/mL. Daidzein and khellin were found to have lower toxicity toward the vero cell line, with IC₅₀ values of 160.81 and 300.23 μg/mL, respectively. Further, molecular docking study and MD simulation of daidzein indicated that it remained stable inside the cavity of DNA GyrB domain for 100 ns.

Synthesis, Biological Evaluation, and Molecular Modeling Studies of New 1,3,4-Thiadiazole Derivatives as Potent Antimicrobial Agents

In this work, the synthesis, characterization, and biological activities of a new series of 1,3,4-thiadiazole derivatives were investigated. The structures of final compounds were identified using ¹ H-NMR, ¹³ C-NMR, elemental analysis, and HRMS. All the new synthesized compounds were then screened for their antimicrobial activity against four types of pathogenic bacteria and one fungal strain, by application of the MIC assays, using Ampicilin, Gentamycin, Vancomycin, and Fluconazole as standards. Among the compounds, the MIC values of 4 and 8 μg/mL of the compounds 3f and 3g, respectively, are remarkable and indicate that these compounds are good candidates for antifungal activity. The docking experiments were used to identify the binding forms of produced ligands with sterol 14-demethylase to acquire insight into relevant proteins. The MD performed about 100 ns simulations to validate selected compounds' theoretical studies. Finally, using density functional theory (DFT) to predict reactivity, the chemical characteristics and quantum factors of synthesized compounds were computed. These results were then correlated with the experimental data. Furthermore, computational estimation was performed to predict the ADME properties of the most active compound 3f.

Design, synthesis, molecular modeling, DFT, ADME and biological evaluation studies of some new 1,3,4-oxadiazole linked benzimidazoles as anticancer agents and aromatase inhibitors

Breast cancer is the most frequent female cancer and second cause of cancer-related deaths among women around the world. Two thirds of breast cancer patients have hormone-dependent tumors, which is very likely be treated with hormonal therapy. Aromatase is involved in the biosynthesis of estrogen thus a critical target for breast cancer. In this study, in order to identify new aromatase enzyme inhibitors, a series of benzimidazole-1,3,4-oxadiazole derivatives were synthesized and characterized by ¹H NMR, ¹³C NMR, and MS spectra analyses. In the in vitro anticancer assay, all the compounds tested anticancer activities using MTT-based assay against five cancer cell lines (MCF-7, A549, HeLa, C6, and HepG2). Among them, compound 5a exhibited the most potent activity with IC₅₀ values of 5.165 ± 0.211 μM and 5.995 ± 0.264 μM against MCF-7 and HepG2 cell lines. Compound 5a was included in the BrdU test to determine the DNA synthesis inhibition effects for both cell types. Furthermore, compound 5c was also found to be more effective than doxorubicin on the HeLa cell line. The selectivity of anticancer activity was evaluated in NIH3T3 cell line. In vitro, enzymatic inhibition assays of aromatase enzyme were performed for compound 5a acting on the MCF-7 cell line. For compound 5a, in silico molecular docking and dynamics simulations against aromatase enzyme was performed to determine possible protein-ligand interactions and stability. DFT study was performed to evaluate the quantum mechanical and electronic properties of compound 5a. Finally, the theoretical ADME properties of the potential aromatase inhibitor compound 5a were analyzed by calculations.Communicated by Ramaswamy H. Sarma.

QSAR Studies, Molecular Docking, Molecular Dynamics, Synthesis, and Biological Evaluation of Novel Quinolinone-Based Thiosemicarbazones against Mycobacterium tuberculosis

In this study, a series of novel quinolinone-based thiosemicarbazones were designed in silico and their activities tested in vitro against Mycobacterium tuberculosis (M. tuberculosis). Quantitative structure-activity relationship (QSAR) studies were performed using quinolinone and thiosemicarbazide as pharmacophoric nuclei; the best model showed statistical parameters of R² = 0.83; F = 47.96; s = 0.31, and was validated by several different methods. The van der Waals volume, electron density, and electronegativity model results suggested a pivotal role in antituberculosis (anti-TB) activity. Subsequently, from this model a new series of quinolinone-thiosemicarbazone 11a-e was designed and docked against two tuberculosis protein targets: enoyl-acyl carrier protein reductase (InhA) and decaprenylphosphoryl-β-D-ribose-2'-oxidase (DprE1). Molecular dynamics simulation over 200 ns showed a binding energy of -71.3 to -12.7 Kcal/mol, suggesting likely inhibition. In vitro antimycobacterial activity of quinolinone-thiosemicarbazone for 11a-e was evaluated against M. bovis, M. tuberculosis H37Rv, and six different strains of drug-resistant M. tuberculosis. All compounds exhibited good to excellent activity against all the families of M. tuberculosis. Several of the here synthesized compounds were more effective than the standard drugs (isoniazid, oxafloxacin), 11d and 11e being the most active products. The results suggest that these compounds may contribute as lead compounds in the research of new potential antimycobacterial agents.

The Struggle to End a Millennia-Long Pandemic: Novel Candidate and Repurposed Drugs for the Treatment of Tuberculosis

This article provides an encompassing review of the current pipeline of putative and developed treatments for tuberculosis, including multidrug-resistant strains. The review has organized each compound according to its site of activity. To provide context, mention of drugs within current recommended treatment regimens is made, thereafter followed by discussion on recently developed and upcoming molecules at established and novel targets. The review is designed to provide a clinically applicable understanding of the compounds that are deemed most currently relevant, including those already under clinical study and those that have shown promising pre-clinical results. An extensive review of the efficacy and safety data for key contemporary drugs already incorporated into treatment regimens, such as bedaquiline, pretomanid, and linezolid, is provided. The three levels of the bacterial cell wall (mycolic acid, arabinogalactan, and peptidoglycan layers) are highlighted and important compounds designed to target each layer are delineated. Amongst others, the highly optimistic and potent anti-mycobacterial activity of agents such as BTZ-043, PBTZ 169, and OPC-167832 are emphasized. The evolving spectrum of oxazolidinones, such as sutezolid, delpazolid, and TBI-223, all aiming to exceed the efficacy achieved with linezolid yet offer a safer alternative to the potential toxicity, are reviewed. New and exciting prospective agents with novel mechanisms of impact against TB, including 3-aminomethyl benzoxaboroles and telacebec, are underscored. We describe new diaryloquinolines in development, striving to build on the immense success of bedaquiline. Finally, we discuss some of these compounds that have shown encouraging additive or synergistic benefit when used in combination, providing some promise for the future in treating this ancient scourge.

Application of Computational Biology and Artificial Intelligence in Drug Design

Traditional drug design requires a great amount of research time and developmental expense. Booming computational approaches, including computational biology, computer-aided drug design, and artificial intelligence, have the potential to expedite the efficiency of drug discovery by minimizing the time and financial cost. In recent years, computational approaches are being widely used to improve the efficacy and effectiveness of drug discovery and pipeline, leading to the approval of plenty of new drugs for marketing. The present review emphasizes on the applications of these indispensable computational approaches in aiding target identification, lead discovery, and lead optimization. Some challenges of using these approaches for drug design are also discussed. Moreover, we propose a methodology for integrating various computational techniques into new drug discovery and design.

Structural modelling and in silico pharmacology of β-carboline alkaloids as potent 5-HT1A receptor antagonists and reuptake inhibitors

Serotonin (5-HT) antagonists and reuptake inhibitors (SARIs) are atypical antidepressants for managing major depressive disorder. They are oftentimes applied as adjuvants for ameliorating aftereffects of SSRI antidepressants including insomnia and sexual dysfunction. The few available candidates of this class including lorpiprazole and trazodone also display some daunting side effects, making a continuous search for improved alternatives essential. Natural β-carboline alkaloids (NβCs) are interestingly renowned with broad pharmacological spectrum against several neuropsychiatric disorders including depression. However, their potentials as SARIs remain underexplored. In this study, 982 NβCs retrieved from the Ambinter-Greenpharma (Amb) database were virtually screened for potent SARI alternatives using computational and biocheminformatics approaches: homology modelling of 5-HT1A receptor, Glide HTVS, SP and XP molecular docking, molecular dynamics (MD) simulation, ADMET and mutagenicity predictions. The homology receptor was validated as a good representative of human 5HT1A receptor using the RCSB structure validation and quality protocols. From the virtual screening against the 5-HT1A receptor, Amb ligands, Amb18709727 and Amb37857532 showed higher binding affinities by XP scores of -8.725 and -7.976 kcal/mol, and MMGBSA of -87.972 and -107.585 kcal/mol respectively compared to lorpiprazole, a reference SARI with XP score and MMGBSA of -6.512 and -62.788 kcal/mol respectively. They maintained ideal contacts with pharmacologically essential amino acid residues implicated in SARI mechanisms and expressed higher stability and compactness than lorpiprazole throughout the trajectories of 100 ns MD simulation. They also displayed interesting ADME, druggability, low toxicity and mutagenicity profiles, ideal for CNS drug prospects, thus, recommended as putative SARI candidates for further study.

Synthesis, Molecular Modeling Study, and Quantum-Chemical-Based Investigations of Isoindoline-1,3-diones as Antimycobacterial Agents

The condensation of phthalic anhydride afforded structurally modified isoindoline-1,3-dione derivatives with selected amino-containing compounds. The title compounds (2-30) have been characterized by thin-layer chromatography (TLC), infrared spectroscopy, ¹H and ¹³C NMR spectroscopy, and mass spectroscopy. All of the compounds were assessed for their antimycobacterial activity toward the H37Rv strain by a dual read-out assay method. Among the synthesized compounds, compound 27 possessed a significant IC₅₀ of 18 μM, making it the most potent compound of the series. The InhA inhibitory (IC₅₀) activity of compound 27 was 8.65 μM in comparison to Triclosan (1.32 μM). Computational studies like density functional theory (DFT) study, molecular docking, and dynamic simulation studies illustrated the reactivity and stability of the synthesized compounds as InhA inhibitors. A quantum-mechanics-based DFT study was carried out to investigate the molecular and electronic properties, reactivities, and nature of bonding present in the synthesized compounds and theoretical vibrational (IR) and isotropic value (¹H and ¹³C NMR) calculations.

Optimizing cardio, hepato and phospholipidosis toxicity of the Bedaquiline by chemoinformatics and molecular modelling approach

The FDA granted expedited approval for Johnson and Johnson's Bedaquiline to treat pulmonary multidrug resistant tuberculosis on 28 December 2012 which is more common in China, Russian Federation and India. Bedaquiline is the first anti-tubercular drug approved by the FDA in the last 40 years, and it has become a cynosure in the circles of synthetic chemists researching new anti-tubercular drugs. Bedaquiline's highly lipophilic nature raises major concerns like suppression of the hERG gene, hepatotoxicity, and phospholipidosis despite its potential antitubercular profile. To address these toxicity concerns, in the present work, we have employed the structural optimization of Bedaquiline using the ADMETopt web server, which optimizes lead with scaffold hopping and ADMET screening. The ADMETopt web server yielded the 476 structures through optimization of three sites in Bedaquiline. Further, we have validated the optimized structures for their activity by performing molecular docking and molecular dynamics (MD) simulations against the mycobacterial ATP synthase enzyme and density functional theory (DFT) study further provides insight into the reactivity of the compounds. After screening and analysis, compound #449 was observed to be the most promising mycobacterial ATP synthase inhibitor with minimal cardiotoxicity, hepatotoxicity and phospholipidosis.

Structural exploration of selected C6 and C7-substituted coumarin isomers as selective MAO-B inhibitors

Monoamine Oxidase B is considered a successful target for developing antiparkinsonian drugs. Due to the side effects of current MAO-B inhibitors, there's an urgent need for novel potent and highly selective MAO-B inhibitors. A recent study has shown that coumarins tend to be more selective towards MAO-B than MAO-A when connected to a hex-5-ynyloxy chain at position 6 in contrast to their C7-isomers. The present study describes the mode of interaction of the C6 and C7-substituted coumarin isomers characterized by their difference in selectivity towards MAO-B through molecular docking and molecular dynamics simulations in an effort to elucidate the structural components and molecular interactions that may be responsible for MAO-B selectivity. Three isomeric coumarin pairs connected to ether chain at position 6 or 7 were taken from the literature and modelled according to their IUPAC nomenclature. Molecular docking study revealed one π- π stacking interaction with Tyr-326 in common between the selective coumarin C6-isomers. Resulting complexes of one isomeric coumarin pair that displayed the highest selectivity shift towards MAO-B were subject to 100 ns molecular dynamics simulations study to analyze the stability of the docked complexes. Molecular dynamics revealed that the C7-isomer is relatively stable in both MAO isoforms through the simulation duration, whereas the C6-isomer deemed unstable for MAO-A which may be due to the bulky Phe-208 residue in MAO-A. Our results might be applied for further development and optimization of coumarin derivatives into a successful drug against Parkinson's disease.Communicated by Ramaswamy H. Sarma.

In vitro and computational investigations of novel synthetic carboxamide-linked pyridopyrrolopyrimidines with potent activity as SARS-CoV-2-MPro inhibitors

An essential target for COVID-19 is the main protease of SARS-CoV-2 (M^pro). With the objective of targeting this receptor, a novel set of pyrido[1,2-a]pyrrolo[2,3-d]pyrimidines with terminal carboxamide fragments was designed, synthesized, and considered as an initial motif for the creation of effective pan-coronavirus inhibitors. Accordingly, nine derivatives (21–29) have been introduced for in vitro assay to evaluate their antiviral activity and cytotoxicity effect against COVID-19 virus using Vero cells. The obtained data revealed that the majority of these derivatives showed potent cellular anti-COVID-19 activity and prevent viral growth by more than 90% at two different concentrations with weak or even no detectable cytotoxic effect on Vero cells. Extensive molecular docking simulations highlighted proper non-covalent interaction of new compounds within the binding pocket of M^pro as a potential target for their antiviral activity. In vitro assay for all the synthesized derivatives against the viral M^pro target indicated that compounds 25 and 29 have promising inhibitory activity with IC₅₀ values at low micromolar concentrations. The molecular dynamic simulation results predicted the stability of compound 29 in the binding cavity of SARS-CoV-2 M^pro and hence supported the high inhibitory activity shown by the In vitro assay. These results suggested that compounds 25 and 29 merit further investigations as promising drug candidates for the management of SARS-CoV-2.

An essential target for COVID-19 is the main protease of SARS-CoV-2 (M^pro).

Identification of an AXL kinase inhibitor in triple-negative breast cancer by structure-based virtual screening and bioactivity test

Breast cancer is a malignant tumor that occurs in the glandular epithelium of the breast, and more than 15% of the patients are triple-negative breast cancer (TNBC). Therefore, finding new targets and targeted therapeutic drugs for TNBC is urgent. Overexpression of the AXL is associated with motility and invasiveness of the TNBC cells, which is a potential target for breast cancer therapy. A compound Y041-5921 (IC₅₀  = 6.069 μm for AXL kinase and IC₅₀  = 4.1 μm for MDA-MB-231 cell line) was identified through structure-based virtual screening and bioassay test for the first time. The compound Y041-5921 could significantly inhibit the proliferation and invasion of the TNBC cells and the toxicity of Y041-5921 to normal immortalized breast epithelial cells was far lower than that of commonly used clinical chemotherapy drugs. Besides, it also had well inhibitory effect on the proliferation of many other malignant tumor cell lines (the IC₅₀  value are 10.0 m, 7.1 m, 10.3 m, 11.4 m and 5.8 m for U251 cell, COLO cell, PC-9 cell, CAKI-1 cell and MG63 cell, respectively). The interaction mechanism between Y041-5921 and AXL was studied by molecular dynamics (MD) simulations and binding free energy calculation, and the key residues whose energy contribution mainly comes from non-polar solvation interaction (such as Ala565, Lys567, Met598, Leu620, Pro621, Met623, Lys624, Arg676, Asn677 and Met679) were identified. The small molecule inhibitors Y041-5921 targeting AXL reported in this work will lay a foundation and provide a theoretical basis for the development of the TNBC.

Computational identification of 2,4-disubstituted amino-pyrimidines as L858R/T790M-EGFR double mutant inhibitors using pharmacophore mapping, molecular docking, binding free energy calculation, DFT study and molecular dynamic simulation

Pharmacophore modelling studies have been performed for a series of 2,4-disubstituted-pyrimidines derivatives as EGFR L858R/T790M tyrosine kinase inhibitors. The high scoring AARR.15 hypothesis was selected as the best pharmacophore model with the highest survival score of 3.436 having two hydrogen bond acceptors and two aromatic ring features. Pharmacophore-based virtual screening followed by structure-based yielded the six molecules (ZINC17013227, ZINC17013215, ZINC9573324, ZINC9573445, ZINC24023331 and ZINC17013503) from the ZINC database with significant in silico predicted activity and strong binding affinity towords the EGFR L858R/T790M tyrosine kinase. In silico toxicity and cytochrome profiling indicates that all the 06 virtually screened compounds were substrate/inhibitors of the CYP-3A4 metabolizing enzyme and were non-carcinogenic and devoid of Ames mutagenesis. Density functional theory (DFT) and molecular dynamic (MD) simulation further validated the obtained hits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00113-x.

Tomatidine and Patchouli Alcohol as Inhibitors of SARS-CoV-2 Enzymes (3CLpro, PLpro and NSP15) by Molecular Docking and Molecular Dynamics Simulations

Considering the current dramatic and fatal situation due to the high spreading of SARS-CoV-2 infection, there is an urgent unmet medical need to identify novel and effective approaches for prevention and treatment of Coronavirus disease (COVID 19) by re-evaluating and repurposing of known drugs. For this, tomatidine and patchouli alcohol have been selected as potential drugs for combating the virus. The hit compounds were subsequently docked into the active site and molecular docking analyses revealed that both drugs can bind the active site of SARS-CoV-2 3CLpro, PLpro, NSP15, COX-2 and PLA2 targets with a number of important binding interactions. To further validate the interactions of promising compound tomatidine, Molecular dynamics study of 100 ns was carried out towards 3CLpro, NSP15 and COX-2. This indicated that the protein-ligand complex was stable throughout the simulation period, and minimal backbone fluctuations have ensued in the system. Post dynamic MM-GBSA analysis of molecular dynamics data showed promising mean binding free energy 47.4633 ± 9.28, 51.8064 ± 8.91 and 54.8918 ± 7.55 kcal/mol, respectively. Likewise, in silico ADMET studies of the selected ligands showed excellent pharmacokinetic properties with good absorption, bioavailability and devoid of toxicity. Therefore, patchouli alcohol and especially, tomatidine may provide prospect treatment options against SARS-CoV-2 infection by potentially inhibiting virus duplication though more research is guaranteed and secured.

β-Carboline alkaloids induce structural plasticity and inhibition of SARS-CoV-2 nsp3 macrodomain more potently than remdesivir metabolite GS-441524: computational approach

The nsp3 macrodomain is implicated in the viral replication, pathogenesis and host immune responses through the removal of ADP-ribosylation sites during infections of coronaviruses including the SARS-CoV-2. It has ever been modulated by macromolecules including the ADP-ribose until Ni and co-workers recently reported its inhibition and plasticity enhancement unprecedentedly by remdesivir metabolite, GS-441524, creating an opportunity for investigating other biodiverse small molecules such as β-Carboline (βC) alkaloids. In this study, 1497 βC analogues from the HiT2LEAD chemical database were screened, using computational approaches of Glide XP docking, molecular dynamics simulation and pk-CSM ADMET predictions. Selectively, βC ligands, 129, 584, 1303 and 1323 demonstrated higher binding affinities to the receptor, indicated by XP docking scores of -10.72, -10.01, -9.63 and -9.48 kcal/mol respectively than remdesivir and GS-441524 with -4.68 and -9.41 kcal/mol respectively. Consistently, their binding free energies were -36.07, -23.77, -24.07 and -17.76 kcal/mol respectively, while remdesivir and GS-441524 showed -21.22 and -24.20 kcal/mol respectively. Interestingly, the selected βC ligands displayed better stability and flexibility for enhancing the plasticity of the receptor than GS-441524, especially 129 and 1303. Their predicted ADMET parameters favour druggability and low expressions for toxicity. Thus, they are recommended as promising adjuvant/standalone anti-SARS-CoV-2 candidates for further study.Key words: SARS-CoV-2, nsp3 macrodomain, ADP-ribose, β-carboline, bioinformatics, drug design.

The Mur Enzymes Chink in the Armour of Mycobacterium tuberculosis cell wall

TUBERCULOSIS: (TB) transmitted by Mycobacterium tuberculosis (Mtb) is one of the top 10 causes of death globally. Currently, the widespread occurrence of resistance toward Mtb strains is becoming a significant concern to public health. This scenario exaggerated the need for the discovery of novel targets and their inhibitors. Targeting the "Mtb cell wall peptidoglycan synthesis" is an attractive strategy to overcome drug resistance. Mur enzymes (MurA-MurF) play essential roles in the peptidoglycan synthesis by catalyzing the ligation of key amino acid residues to the stem peptide. These enzymes are unique and confined to the eubacteria and are absent in humans, representing potential targets for anti-tubercular drug discovery. Mtb Mur ligases with the same catalytic mechanism share conserved amino acid regions and structural features that can conceivably exploit for the designing of the inhibitors, which can simultaneously target more than one isoforms (MurC-MurF) of the enzyme. In light of these findings in the current review, we have discussed the recent advances in medicinal chemistry of Mtb Mur enzymes (MurA-MurF) and their inhibitors, offering attractive multi-targeted strategies to combat the problem of drug-resistant in M. tuberculosis.

Effect of Aminosilane Coupling Agent-Modified Nano-SiO2 Particles on Thermodynamic Properties of Epoxy Resin Composites

From the perspective of improving the thermodynamic properties of epoxy resin, it has become the focus of research to enhance the operational stability of GIS (Gas Insulated Substation) basin insulators for UHV (Ultra-High Voltage) equipment. In this paper, three aminosilane coupling agents with different chain lengths, (3-Aminopropyl)trimethoxysilane (KH550), Aminoethyl)-γ-aminopropyltrimethoxysilane (KH792) and 3-[2-(2-Aminoethylamino)ethylamino]propyl-trimethoxysilane (TAPS), were used to modify nano-SiO2 and doped into epoxy resin, respectively, using a combination of experimental and molecular dynamics simulations. The experimental results showed that the surface-grafted KH792 model of nano-SiO2 exhibited the most significant improvement in thermal properties compared with the undoped nanoparticle model. The storage modulus increased by 276 MPa and the Tg increased by 61 K. The simulation results also showed that the mechanical properties of the nano-SiO2 surface-grafted KH792 model were about 3 times higher than that of the undoped nanoparticle model, the Tg increased by 36.5 K, and the thermal conductivity increased by 24.5%.