Structural based study to identify new potential inhibitors for dual specificity tyrosine-phosphorylation- regulated kinase

Evaluation of Pyrrolone-Fused Benzosuberene MK2 Inhibitors as Promising Therapeutic Agents for HNSCC: In Vitro Efficacy, In-Vivo Safety, and Pharmacokinetic Profiling

MAPKAPK2/MK2 is well implicated in the progression of Head and Neck Squamous Cell Carcinoma (HNSCC), and potent MK2-inhibitors are required to suppress its activity. Several MK2-inhibitors have been developed in recent years to combat its effects on cancer. However, inadequate solubility, insufficient cellular permeability, systemic toxicity-mediated side effects, and low bioavailability have severely impeded the advancement of MK2-inhibitors to clinical trials. This void necessitates research to develop less toxic and more bioavailable potent MK2-inhibitors in HNSCC. In the present article, we have evaluated the in-vitro efficacy, in-vivo single-dose acute toxicity, and in-vivo pharmacokinetic profiling of recently developed PfBS (pyrrolone-fused benzosuberene) MK2-inhibitor analogues against HNSCC. The PfBS MK2 inhibitor analogues impeded HPV+ and HPV- HNSCC cell proliferation and two-dimensional migration. Moreover, MK2-inhibitors lowered HNSCC cell clonogenic survival in a dose-dependent manner, significantly enhancing radiation-induced cell death via exerting radio-sensitization effects. Furthermore, γ-H2AX immunostaining revealed that PfBS analogues impaired DNA damage repair in HNSCC cells exposed to gamma radiation. In mice, PfBS MK2 inhibitors at 300 mg/kg were well-tolerated without any lethal effects. Pharmacokinetic studies showed that PfBS analogues exhibited rapid absorption (Tmax), adequate plasma concentration above the micromolar level (C0 or Cmax), limited tissue distribution (Vd), and faster elimination from the body (Cl). Overall, this study summarizes in-vitro efficacy, safety, and pharmacokinetics of developed MK2-inhibitors and opens doors for pharmacodynamics and mechanism of action study of most effective leads in HNSCC.

Structure-based virtual screening of anti-breast cancer compounds from Artemisia absinthium-insights through molecular docking, pharmacokinetics, and molecular dynamic simulations

Breast cancer is the world's second most frequent malignancy, with a significant mortality and morbidity rate. Nowadays, natural breast cancer medicine has piqued attention as disease-curing agent with low side effects. Herein, the leaf powder of Artemisia absinthium was extracted with ethanol, and GC-MS and LC-MS methods were employed to identify the phytocompounds. Using commercial software SeeSAR-9.2 and StarDrop, identified phytocompounds were docked with estrogen and progesterone breast cancer receptors as they promote breast cancer growth to find the binding affinity of the ligands, drugability, and toxicity. Hormone-mediated breast cancer accounts for about 80% of all cases of breast cancer. Cancer cells proliferate when estrogen and progesterone hormones are attached to these receptors. The molecular docking results demonstrated that 3',4',5,7-Tetrahydroxyisoflavanone (THIF) has stronger binding efficacy than standard drugs and other phytocompounds with -28.71 (3 hydrogen bonds) and -24.18 kcal/mol (6 hydrogen bonds) binding energies for estrogen and progesterone receptors, respectively. Pharmacokinetics and toxicity analysis were done to predict the drug-likeness of THIF which results in good drugability and less toxicity. The best fit THIF was subjected to a molecular dynamics simulation analysis by using Gromacs to analyze the conformational changes that occurred during protein-ligand interaction and found that, the structural changes were observed. The results from MD simulation and pharmacokinetic studies suggested that THIF can be expected that in vitro and in vivo research on this compound may lead to the development of a potent anti-breast cancer drug in the future.Communicated by Ramaswamy H. Sarma.

Computational investigation of the impact of potential AT2R polymorphism on small molecule binding

For more than a century, the renin-angiotensin system (RAS) has been acknowledged for playing a crucial part in the physiological control of arterial pressure, as well as sodium and fluid balance. It is now generally acknowledged that one of the receptor of RAS system i.e. angiotensin type 2 receptor (AT2R) functions as a repair system during pathophysiologic circumstances and performs a significant protective role. Efforts have been made previously to design suitable agonist and antagonist molecules to potentially modulate AT2R. One of the agonists and antagonists, named C21 and EMA401, has been studied in a number of pathological conditions. Additionally, a wide panel of single nucleotide polymorphisms (SNPs) has been reported for AT2R, which might potentially affect the efficacy of these molecules. Therefore, computational investigations have been carried out to analyze all the SNPs (1151) reported in NCBI to find potential SNPs affecting the active site of AT2R, as this domain is still unexplored. Structures of these polymorphic forms were modeled, and in silico drug interaction studies with C21 and EMA401 were carried out. The two mutants (rs868939201 and rs1042852794) that significantly affect the binding affinity as that of the wild type were subjected to molecular dynamics simulations. Our analysis of native and mutant AT2R and their complexes with C21 and EMA401 indicated that the occurrence of these mutations affects the conformation of the protein and has affected the binding of these ligand molecules. The study's findings will aid in the development of better, more versatile medications in the near future, and also in vitro and in vivo studies might be planned in accordance with recent findings.Communicated by Ramaswamy H. Sarma.

Design of novel isoxazole derivatives as tubulin inhibitors using computer-aided techniques: QSAR modeling, in silico ADMETox, molecular docking, molecular dynamics, biological efficacy, and retrosynthesis

In the current work, computational methods were used to investigate new isoxazole derivatives that could be used as tubulin inhibitors. The study aims to develop a reliable quantitative structure-activity relationship (QSAR) model, following the criteria set by Golbraikh, Tropsha, and Roy. As a result, seven candidate compounds were developed, all having higher activity than the well-established anticancer agent Cisplatin (Cisp). According to the ADMETox in silico test, the candidates Pr4, Pr5, and P6 can be toxic. As a result, we have chosen to focus our study on compounds Pr1, Pr2, and Pr3. Molecular docking analysis revealed that drug candidate Pr2 exhibits the highest stability within the oxidized quinone reductase 2 (PDB ID: 4zvm), target receptor (ΔG(Pr2) = ΔG(Pr3) = -10.4 < ΔG(Pr1) = -10.0 < ΔG(Cisp) = -7.3 kcal/mol). This finding aligns with the activity predictions made by the QSAR model. Furthermore, molecular dynamics simulations of the Pr2-4zvm complex over 100 ns confirm the ligand's robust stability within the receptor's active site, supporting the results obtained from molecular docking and the QSAR model predictions. The CaverDock software was utilized to identify the tunnels likely to be followed by ligands moving from the active site to the receptor surface. This analysis also helped in determining the biological efficacy of the target compounds. The results indicated that the Pr2 compound is more effective than the others. Finally, the computer-assisted retrosynthesis process of two high confidence sequences was used to synthesize drug candidates.

Novel pyrrolone-fused benzosuberene MK2 inhibitors: synthesis, pharmacophore modelling, molecular docking, and anti-cancer efficacy evaluation in HNSCC cells

Head and neck Squamous Cell Carcinoma (HNSCC) is a growing concern worldwide and MAPKAPK2/MK2 (Mitogen-Activated Protein Kinase Activated Protein Kinase 2) is crucially involved in HNSCC progression. Increased disease burden and lacuna of targeted therapies require novel and safe pharmacological inhibitors to suppress the well-explored molecular targets in HNSCC. Here, we used dibromo-substituted benzosuberene synthesized from the mixture of α, β, γ-himachalenes and utilized as a precursor for the synthesis of Pyrrolone-fused benzosuberenes (PfBS) as MK2 inhibitors through aminocarbonylation approach in a single-pot reaction. The devised protocol provides a broad substrate scope, facile recovery, recyclability of Polystyrene-supported palladium (Pd@PS) nanoparticle catalyst, and fewer synthesis steps. In-silico molecular docking, pharmacophore modeling, and ADMET revealed MK2-inhibitory potential and drug-likeliness of PfBS analogues. Surface plasmon resonance (SPR) analysis revealed effective high binding affinity (KD) and kinetics of PfBS analogues with MK2. Additionally, the SPR-mediated in-solution inhibition assay established the MK2-inhibition properties of PfBS analogues through abrogation of MK2-Hsp27 interaction. Further, in-vitro studies validate the findings in HNSCC cells. PfBS analogues exhibited significant anti-proliferative effects on CAL 27 tongue squamous carcinoma cells and were found safe on IEC-6 intestinal epithelial cells. Moreover, immunofluorescence analysis and western-blot assays potentiated, that selected analogues inhibited the inflammatory cytokine TNF-α induced activation of MK2 on cellular and molecular levels in HNSCC cells. In conclusion, this study presents novel MK2-inhibitors and opens the avenue for further pre-clinical and clinical efficacy evaluation of developed PfBS analogues in the treatment of HNSCC.Communicated by Ramaswamy H. Sarma.

Effect of co-pigments on anthocyanins of Rhododendron arboreum and insights into interaction mechanism

The impact of intermolecular copigmentation between five phenolic acids, two flavonoid and three amino acids with R. arboreum anthocyanins (ANS) and its isolated cyanidin-3-O-monoglycosides were investigated through experimental and theoretical approach. On addition of different copigments, phenolic acid induced strong hyperchromic (0.26-0.55 nm) and bathochromic shift (6.6-14.2 nm). The color intensity and stability of ANS with, storage at 4 °C & 25 °C, sunlight, oxidation and heat were evaluated by chromaticity, anthocyanin content, kinetic and structural simulation analysis. The strongest copigmentation reaction was observed with narningin (NA) and also showed high thermostability and highest half-life i.e. 3.39 h-1.24 h at 90-160 °C. The cyanidin-3-O-monoglycosides were analysed for their copigmentation effect and observations revealed that NA displayed best copigmentation effect to cyanidin-3-O-arabinoside (B) followed by cyanidin-3-O-galactoside (A), and cyanidin-3-O-rhamnoside (C). Additionally, structural simulation and steered molecular dynamics insights NA is the most favourable co-pigment involving π-π stacking and H-bonding.

Structure-based virtual screening for identification of potential CDC20 inhibitors and their therapeutic evaluation in breast cancer

Cell division cycle 20 (CDC20), a critical partner of anaphase promoting complex (APC/C), is indispensably required for metaphase-to-anaphase transition. CDC20 overexpression in TNBC breast cancer patients has been found to be correlated with poor prognosis, hence, we aimed to target CDC20 for TNBC therapeutics. In silico molecular docking of large-scale chemical libraries (phytochemicals/synthetic drugs) against CDC20 protein structure identified five synthetic drugs and four phytochemicals as potential hits interacting with CDC20 active site. The molecular selection was done based on docking scores, binding interactions, binding energies and MM/GBSA scores. Further, we analysed ADME profiles for all the hits and identified lidocaine, an aminoamide anaesthetic group of synthetic drug, with high drug-likeness properties. We explored the anti-tumorigenic effects of lidocaine on MDA-MB-231 TNBC breast cancer cells, which resulted in increased growth inhibition in dose-dependent manner. The molecular mechanism behind the cell viability defect mediated by lidocaine was found to be induction of G2/M cell cycle arrest and cellular apoptosis. Notably, lidocaine treatment of TNBC cells also resulted in downregulation of CDC20 gene expression. Thus, this study identifies lidocaine as a potential anti-neoplastic agent for TNBC cells emphasizing CDC20 as a suitable therapeutic target for breast cancer.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03554-7.

β-Sitosterol, a phytocompound from Parthenium hysterophorus, reveals anti-diabetic properties through α-Amylase inhibition: an in-silico and in-vitro analysis

The study aims to identify and validate a potential α-Amylase inhibitor from the leaf extract of the Parthenium hysterophorus. Molecular docking and dynamics analyses were performed to test the anti-diabetic efficacy of the compound by focusing on α-Amylase inhibition. The molecular docking study using AutoDock Vina (PyRx) and SeeSAR tools identified β-Sitosterol as an effective α-Amylase inhibitory compound. Among the analysed fifteen phytochemicals, β-Sitosterol demonstrated the most appreciable binding energy (-9.0 Kcal/mol) and is comparatively higher than the binding energy of the standard α-Amylase inhibitor, the Acarbose (-7.6 Kcal/mol). The significance of the interaction between β-Sitosterol and α-Amylase was further investigated using Molecular Dynamics Simulation (MDS) for 100 ns via GROMACS. The data reveals that the compound could exhibit the highest stability with α-Amylase regarding RMSD, RMSF, SASA and Potential Energy analysis. The key residue of α-Amylase (Asp -197) shows a significantly low fluctuation of 0.7 Å while interacting with β-Sitosterol. The data obtained from MDS results strongly suggested the potential inhibitory impact of β-Sitosterol on α-Amylase. In addition, the proposed phytochemical was purified from the leaf extracts of P.hysterophorus using the silica gel column chromatography and identified by GC-MS analysis. The purified β-Sitosterol demonstrated a significant 42.30% in-vitro α-Amylase enzyme inhibition property under 400 µg/ml concentration and thus supported the in-silico predictions. Further in-vivo investigations are necessary to analyse the efficiency of β-Sitosterol on α-Amylase inhibition to help the anti-diabetic potential of the phytocompound.Communicated by Ramaswamy H. Sarma.

Designing of kinase hinge binders: A medicinal chemistry perspective

Protein kinases are key regulators of cellular signaling and play a critical role in oncogenesis. Inhibitors of protein kinases are pursued by both industry and academia as a promising target for cancer therapy. Within the protein kinases, the ATP site has produced more than 40 FDA-approved drugs. The ATP site is broadly composed of a hinge region, gatekeeper residues, DFG-loop, ribose pocket, and other hydrophobic regions. The hinge region in the ATP site can be used for designing potent inhibitors. In this review, we discuss some representative studies that will highlight the interactions of heterocyclic compounds with hinge regions of different kinases like BRAF kinase, EGRF kinase, MAP kinase, and Mps1 kinase.

Molecular docking analysis reveals differential binding affinities of multiple classes of selective inhibitors towards cancer-associated KRAS mutants

KRAS is the most frequently mutated oncogene in solid cancers, and inhibitors that specifically target the KRAS-G12C mutant were recently approved for clinical use. The limited availability of experimental data pertaining to the sensitivity of KRAS-non-G12C mutants towards RAS inhibitors made it difficult to predict the response of KRAS-mutated cancers towards RAS-targeted therapies. The current study aims at evaluating sensitivity profiles of KRAS-non-G12C mutations towards clinically approved sotorasib and adagrasib, and experimental RAS inhibitors based on binding energies derived through molecular docking analysis. Computationally predicted sensitivities of KRAS mutants conformed with the available but limited experimental data, thus validating the usefulness of molecular docking approach in predicting clinical response towards RAS inhibitor treatment. Our results indicate differential sensitivity of KRAS mutants towards both clinical and experimental therapeutics; while certain mutants exhibited broad cross-resistance to most inhibitors, some mutants showed resistance towards specific inhibitors. These results thus suggest the potential of emergence of more resistance mutations in future towards RAS-targeted therapy and points to an urgent need to develop novel classes of inhibitors that are able to overcome both primary and secondary drug resistance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03407-9.

Novel benzylidene benzofuranone analogues as potential anticancer agents: design, synthesis and in vitro evaluation based on CDK2 inhibition assays

The classical anticancer agents do not have their efficacy on inhibiting the G2 phase of the cell cycle. There are a very few reports available on drugs that work at G2 phase. Flavopiridol is one such drug candidate. In the current study, we sought to make analogues of flavopiridol. Still, the conditions used during their synthesis were unfavourable for the formation of flavopiridol and led to the generation of benzofuranones. In the present work, a new series of benzylidene benzofuranones were designed, synthesized and evaluated for their antioxidant, anti-colorectal cancer activity. Molecular docking, MMGBSA and molecular dynamics studies were conducted to assess their binding affinity at the active site of CDK2. Based on the cytotoxicity exhibited by test compounds, the compound NISOA4 (from isopropyl series) was further selected for mechanistic anticancer studies on HCT 116 cell lines. The compound selected was evaluated by comet assay, DNA fragmentation assay, cell cycle analysis, apoptosis detection by annexin FITC, semi-quantitative RTPCR based gene expression studies and FRET assay on the target CDK2/Cyclin A. Compound NISOA4 exhibited marked olive moments in comet assay studies. The apoptotic DNA fragmentation for compound NISOA4 demonstrated a marked change in the DNA fragmentation. The compound exhibited cell cycle arrest at G2/M phase at both the test concentrations. Apoptosis induction was observed at both the test concentrations and the compound was found to be a potent proapoptotic agent. It exhibited marked inhibition for the CDK2 gene expression and did not show any effect on CyclinA gene expression. However, the compound NISOA4 along with other analogues showed appreciable inhibition for the CDK2/Cyclin A target enzyme.

Exploring the Role of Chemical Reactions in the Selectivity of Tyrosine Kinase Inhibitors

A variety of diseases are associated with tyrosine kinase enzymes that activate many proteins via signal transduction cascades. The similar ATP-binding pockets of these tyrosine kinases make it extremely difficult to design selective covalent inhibitors. The present study explores the contribution of the chemical reaction steps to the selectivity of the commercialized inhibitor acalabrutinib over the Bruton's tyrosine kinase (BTK) and the interleukin-2-inducible T-cell kinase (ITK). Ab initio and empirical valence bond (EVB) simulations of the two kinases indicate that the most favorable reaction path involves a water-assisted mechanism of the 2-butynamide reactive group of acalabrutinib. BTK reacts with acalabrutinib with a substantially lower barrier than ITK, according to our calculated free-energy profile and kinetic simulations. Such a difference is due to the microenvironment of the active site, as further supported by a sequence-based analysis of specificity determinants for several commercialized inhibitors. Our study involves a new approach of simulating directly the IC50 and inactivation efficiency k_eff, instead of using the standard formulas. This new strategy is particularly important in studies of covalent inhibitors with a very exothermic bonding step. Overall, our results demonstrate the importance of understanding the chemical reaction steps in designing selective covalent inhibitors for tyrosine kinases.

Theaflavin 3-gallate inhibits the main protease (Mpro) of SARS-CoV-2 and reduces its count in vitro

The main protease (M^pro) of SARS-CoV-2 has been recognized as an attractive drug target because of its central role in viral replication. Our previous preliminary molecular docking studies showed that theaflavin 3-gallate (a natural bioactive molecule derived from theaflavin and found in high abundance in black tea) exhibited better docking scores than repurposed drugs (Atazanavir, Darunavir, Lopinavir). In this study, conventional and steered MD-simulations analyses revealed stronger interactions of theaflavin 3-gallate with the active site residues of M^pro than theaflavin and a standard molecule GC373 (a known inhibitor of M^pro and novel broad-spectrum anti-viral agent). Theaflavin 3-gallate inhibited M^pro protein of SARS-CoV-2 with an IC₅₀ value of 18.48 ± 1.29 μM. Treatment of SARS-CoV-2 (Indian/a3i clade/2020 isolate) with 200 μM of theaflavin 3-gallate in vitro using Vero cells and quantifying viral transcripts demonstrated reduction of viral count by 75% (viral particles reduced from Log10^6.7 to Log10^6.1). Overall, our findings suggest that theaflavin 3-gallate effectively targets the M^pro thus limiting the replication of the SARS-CoV-2 virus in vitro.

Natural Compounds or Their Derivatives against Breast Cancer: A Computational Study

BACKGROUND

Breast cancer is one of the most common types of cancer diagnosed and the second leading cause of death among women. Breast cancer susceptibility proteins of type 1 and 2 are human tumor suppressor genes. Genetic variations/mutations in these two genes lead to overexpression of human breast tumor suppressor genes (e.g., BRCA1, BRCA2), which triggers uncontrolled duplication of cells in humans. In addition, multidrug resistance protein 1 (MDR1), an important cell membrane protein that pumps many foreign substances from cells, is also responsible for developing resistance to cancer chemotherapy. Aim of the Study. The aim of this study was to analyze some natural compounds or their derivatives as part of the development of strong inhibitors for breast cancer. Methodology. Molecular docking studies were performed using compounds known in the literature to be effective against BRCA1 and BRCA2 and MDR1, with positive control being 5-fluorouracil, an antineoplastic drug as a positive control.

RESULTS

The binding affinity of the compounds was analyzed, and it was observed that they had a better binding affinity for the target proteins than the standard drug 5-fluorouracil. Among the compounds analyzed, α-hederin, andrographolide, apigenin, asiatic acid, auricular acid, sinularin, curcumin, citrinin, hispolon, nerol, phytol, retinol palmitate, and sclareol showed the best binding affinity energy to the BRCA1, BRCA2, and MDR1 proteins, respectively.

CONCLUSIONS

α-Hederin, andrographolide, apigenin, asiatic acid, auricular acid, hispolon, sclareol, curcumin, citrinin, and sinularin or their derivatives can be a good source of anticancer agents in breast cancer.

Virtual repurposing of ursodeoxycholate and chenodeoxycholate as lead candidates against SARS-Cov2-Envelope protein: A molecular dynamics investigation

Drug repurposing is an apt choice to combat the currently prevailing global threat of COVID-19, caused by SARS-Cov2in absence of any specific medication/vaccine. The present work employs state of art computational methods like homology modelling, molecular docking and molecular dynamics simulations to evaluate the potential of two widely used surfactant drugs namely chenodeoxycholate(CDC) and ursodeoxycholate (UDC), to bind to the envelope protein of SARS-Cov2(SARS-Cov2-E).The monomeric unit of SARS-Cov2-E was modelled from a close homologue (>90% sequence identity) and a pentameric assembly was modelled using symmetric docking, followed by energy minimization in a DPPC membrane environment. The minimized structure was used to generate best scoring SARS-Cov2-E-CDC/UDC complexes through blind docking. These complexes were subjected to 230 ns molecular dynamics simulations in triplicates in a DPPC membrane environment. Comparative analyses of structural properties and molecular interaction profiles from the MD trajectories revealed that, both CDC and UDC could stably bind to SARS-Cov2-E through H-bonds, water-bridges and hydrophobic contacts with the transmembrane-channelresidues.T30 was observed to be a key residue for CDC/UDC binding. CDC/UDC binding affected the H-bonding pattern between adjacent monomeric chains, slackening the compact transmembrane region of SARS-Cov2-E. Additionally, the polar functional groups of CDC/UDC facilitated entry of a large number of water molecules into the channel. These observations suggest CDC/UDC as potential candidates to hinder the survival of SARS-Cov2 by disrupting the structure of SARS-Cov2-E and facilitating the entry of solvents/polar inhibitors inside the viral cell.Communicated by Ramaswamy H. Sarma.

Prediction of potential inhibitors against SARS-CoV-2 endoribonuclease: RNA immunity sensing

The World Health Organization has classified the COVID-19 outbreak a pandemic which is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) and declared it a global health emergency. Repurposing drugs with minimum side effects are one approach to quickly respond in attempt to prevent the spread of COVID-19. SARS-CoV-2 encodes several RNA processing enzymes that are unusual and unique for single-stranded RNA viruses, including Nsp15, a hexameric endoribonuclease that discriminatory cleaves immediately 3' of uridines. The structure of SARS-CoV-2 Nsp15 is reported to be homologous to that of the Nsp15 endoribonucleases of SARS-CoV and MERS-CoV, but it exhibits differences that may contribute to the greater virulence of SARS-CoV-2. This study aimed to identify drugs that targeted SARS-COV-2 Nsp15 using a molecular docking-based virtual screening of a library containing 10,000 approved and experimental drugs. The molecular docking results revealed 19 medications that demonstrated a good ability to inhibit Nsp15. Among all the candidated 19 drugs only five FDA approved drugs were used for further investigation by molecular dynamics simulation, the stability of Nsp15-ligand system was evaluated by calculating the RMSD, RMSF, radius of gyration and hydrogen bond profile. Furthermore, MM-PBSA method was employed to validate the binding affinity. According to the obtained results of MD, the complex of Olaparib was showed more stability and lower binding free energy than the control inhibitor during MD simulation time. Finally, we suggest that Olaparib is a potential drug for treating patients infected with SARS-CoV-2 and provide insight into the host immune response to viral RNA.Communicated by Ramaswamy H. Sarma.

Identifying New Ligands for JNK3 by Fluorescence Thermal Shift Assays and Native Mass Spectrometry

The c-Jun N-terminal kinases (JNKs) are evolutionary highly conserved serine/threonine kinases. Numerous findings suggest that JNK3 is involved in the pathogenesis of neurodegenerative diseases, so the inhibition of JNK3 may be a potential therapeutic intervention. The identification of novel compounds with promising pharmacological properties still represents a challenge. Fluorescence thermal shift screening of a chemically diversified lead-like scaffold library of 2024 pure compounds led to the initial identification of seven JNK3 binding hits, which were classified into four scaffold groups according to their chemical structures. Native mass spectrometry validated the interaction of 4 out of the 7 hits with JNK3. Binding geometries and interactions of the top 2 hits were evaluated by docking into a JNK3 crystal structure. Hit 5 had a K _d of 21 μM with JNK3 suggested scaffold 5-(phenylamino)-1H-1,2,3-triazole-4-carboxamide as a novel and selective JNK3 binder.

Computational studies on the cholinesterase, beta-secretase 1 (BACE1) and monoamine oxidase (MAO) inhibitory activities of endophytes-derived compounds: towards discovery of novel neurotherapeutics

Cholinesterases, beta-secretase 1 (BACE1) and monoamine oxidase (MAO) are significant in the etiology of neurodegenerative diseases. Inhibition of these enzymes is therefore a major strategy for the development of neurotherapeutics. Even though, this strategy has birthed some approved synthetic drugs, they are characterized by adverse effects. It is therefore, imperative to explore promising alternatives. Consequently, we assessed the inhibitory activities of some endophytes-derived compounds against selected targets towards discovery of novel neurotherapeutics. Standard inhibitors and 83 endophytes-derived compounds were docked against acetylcholinesterase (AChE), butyrylcholinesterase (BChE), BACE 1 and MAO using AutodockVina while the molecular interactions between the selected targets and the compounds with notable binding affinity were viewed through Discovery Studio Visualizer. Druglikeness and Absorption-Distribution-Metabolism-Excretion-Toxicity (ADMET) and blood brain barrier (BBB) properties of the top 4 compounds were evaluated using the Swiss online ADME web tool and OSIRIS server; ligands-enzymes complex stability was assessed through molecular dynamics (MD) simulation. From the 83 compounds, asperflavin, ascomfurans C, camptothecine and corynesidone A exhibited remarkable inhibitory activity against all the four target enzymes compared to the respective standard inhibitors. However, only corynesidone A could transverse the BBB and predicted to be safe. MD simulation of the unbound and complexed enzymes with corynesidone A showed that the complexes were stable throughout the simulation time. Given the exceptional inhibitory activity of endophytes-derived corynesidone A against the four selected targets, its ability to permeate the BBB, excellent drugability properties as well as its stability when complexed with the enzymes, it is a good candidate for further studies towards development of new neurotherapeutics.Communicated by Ramaswamy H. Sarma.

A computational approach to drug repurposing against SARS-CoV-2 RNA dependent RNA polymerase (RdRp)

The spread of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) caused a worldwide outbreak of coronavirus disease 19 (COVID-19), which rapidly evolved as a global concern. The efforts of the scientific community are pointed towards the identification of promptly available therapeutic options. RNA-dependent RNA polymerase (RdRp) is a promising target for developing small molecules to contrast SARS-CoV-2 replication. Modern computational tools can boost identification and repurposing of known drugs targeting RdRp. We here report the results regarding the screening of a database containing more than 8800 molecules, including approved, experimental, nutraceutical, illicit, withdrawn and investigational compounds. The molecules were docked against the cryo-electron microscopy structure of SARS-CoV-2 RdRp, optimized by means of molecular dynamics (MD) simulations. The adopted three-stage ensemble docking study underline that compounds formerly developed as kinase inhibitors may interact with RdRp.Communicated by Ramaswamy H. Sarma.

Identification of polypharmacological anticancerous molecules against Aurora kinase family of proteins

The Human Aurora Kinase (AURK) protein family is the key player of cell cycle events including spindle assembly, kinetochore formation, chromosomal segregation, centrosome separation, microtubule dynamics, and cytokinesis. Their aberrant expression has been extensively linked with chromosomal instability in addition to derangement of multiple tumor suppressors and oncoprotein regulated pathways. Therefore, the AURK family of kinases is a promising target for the treatment of various types of cancer. Over the past few decades, several potential inhibitors of AURK proteins have been identified and have reached various phases of clinical trials. But very few molecules have currently crossed the safety criteria due to their various toxic side effects. In the present study, we have adopted a computational polypharmacological strategy and identified four novel molecules that can target all three AURKs. These molecules were further investigated for their binding stabilities at the ATP binding pocket using molecular dynamics based simulation studies. The molecules selected adopting a multipronged computational approach can be considered as potential AURKs inhibitors for cancer therapeutics.

Computational studies to identify the common type-I and type-II inhibitors against the CDK8 enzyme

In this study, multicomplex-based pharmacophore modeling was conducted on the structural proteome of the two states of CDK8 protein, that is, DMG-in and out. Three pharmacophores having six, five, and four features were selected as the representative models to conduct the virtual screening process using the prepared drug-like natural product database. The screened candidates were subjected to molecular docking studies on DMG-in (5XS2) and out (4F6U) conformation of the CDK8 protein. Subsequently, the common four docked candidates of 5XS2 and 4F6U were selected to perform the molecular dynamics simulation studies. Apart from one of the complexes of DMG-in (5XS2-UNPD163102), all other complexes displayed stable dynamic behavior. The interaction and stability studies of the docked complexes were compared with the references selected from the two conformations (DMG-in and out) of the protein. The current work leads to the identification of three common DMG-in and out hits with diverse scaffolds which can be employed as the initial leads for the design of the novel CDK8 inhibitors.

Elucidating the Glucokinase Activating Potentials of Naturally Occurring Prenylated Flavonoids: An Explicit Computational Approach

Glucokinase activators are considered as new therapeutic arsenals that bind to the allosteric activator sites of glucokinase enzymes, thereby maximizing its catalytic rate and increasing its affinity to glucose. This study was designed to identify potent glucokinase activators from prenylated flavonoids isolated from medicinal plants using molecular docking, molecular dynamics simulation, density functional theory, and ADMET analysis. Virtual screening was carried out on glucokinase enzymes using 221 naturally occurring prenylated flavonoids, followed by molecular dynamics simulation (100 ns), density functional theory (B3LYP model), and ADMET (admeSar 2 online server) studies. The result obtained from the virtual screening with the glucokinase revealed arcommunol B (−10.1 kcal/mol), kuwanon S (−9.6 kcal/mol), manuifolin H (−9.5 kcal/mol), and kuwanon F (−9.4 kcal/mol) as the top-ranked molecules. Additionally, the molecular dynamics simulation and MM/GBSA calculations showed that the hit molecules were stable at the active site of the glucokinase enzyme. Furthermore, the DFT and ADMET studies revealed the hit molecules as potential glucokinase activators and drug-like candidates. Our findings suggested further evaluation of the top-ranked prenylated flavonoids for their in vitro and in vivo glucokinase activating potentials.

Target-based drug repurposing against Candida albicans-A computational modeling, docking, and molecular dynamic simulations study

The emergence of multidrug-resistant strains of Candida albicans has become a global threat mostly due to co-infection with immune-compromised patients leading to invasive candidiasis. The life-threatening form of the disease can be managed quickly and effectively by drug repurposing. Thus, the study used in silico approaches to evaluate Food and Drug Administration (FDA) approved drugs against three drug targets-TRR1, TOM40, and YHB1. The tertiary structures of three drug targets were modeled, refined, and evaluated for their structural integrity based on PROCHECK, ERRAT, and PROSA. High-throughput virtual screening of FDA-approved drugs (8815), interaction analysis, and energy profiles had revealed that DB01102 (Arbutamine), DB01611 (Hydroxychloroquine), and DB09319 (Carindacillin) exhibited better binding affinity with TRR1, TOM40, and YHB1, respectively. Notably, the molecular dynamic simulation explored that Gln45, Thr119, and Asp288 of TRR1; Thr107 and Ser121 of TOM40; Arg193, Glu213, and Ser228 of YHB1 are crucial residues for stable drug-target interaction. Additionally, it also prioritized Arbutamine-TRR1 as the best drug-target complex based on MM-PBSA (-52.72 kcal/mol), RMSD (2.43 Å), and radius of gyration (-21.49 Å) analysis. In-depth, PCA results supported the findings of molecular dynamic simulations. Interestingly, the conserved region (>70%) among the TRR1 sequences from pathogenic Candida species indicated the effectiveness of Arbutamine against multiple species of Candida as well. Thus, the study dispenses new insight and enriches the understanding of developing an advanced technique to consider potential antifungals against C. albicans. Nonetheless, a detailed experimental validation is needed to investigate the efficacy of Arbutamin against life-threatening candidiasis.

Pathogenic genetic variants from highly connected cancer susceptibility genes confer the loss of structural stability

Genetic polymorphisms in DNA damage repair and tumor suppressor genes have been associated with increasing the risk of several types of cancer. Analyses of putative functional single nucleotide polymorphisms (SNP) in such genes can greatly improve human health by guiding choice of therapeutics. In this study, we selected nine genes responsible for various cancer types for gene enrichment analysis and found that BRCA1, ATM, and TP53 were more enriched in connectivity. Therefore, we used different computational algorithms to classify the nonsynonymous SNPs which are deleterious to the structure and/or function of these three proteins. The present study showed that the major pathogenic variants (V1687G and V1736G of BRCA1, I2865T and V2906A of ATM, V216G and L194H of TP53) might have a greater impact on the destabilization of the proteins. To stabilize the high-risk SNPs, we performed mutation site-specific molecular docking analysis and validated using molecular dynamics (MD) simulation and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) studies. Additionally, SNPs of untranslated regions of these genes affecting miRNA binding were characterized. Hence, this study will assist in developing precision medicines for cancer types related to these polymorphisms.

In Silico Evaluation of Iranian Medicinal Plant Phytoconstituents as Inhibitors against Main Protease and the Receptor-Binding Domain of SARS-CoV-2

The novel coronavirus disease 2019 (COVID-19) is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which initially appeared in Wuhan, China, in December 2019. Elderly individuals and those with comorbid conditions may be more vulnerable to this disease. Consequently, several research laboratories continue to focus on developing drugs to treat this infection because this disease has developed into a global pandemic with an extremely limited number of specific treatments available. Natural herbal remedies have long been used to treat illnesses in a variety of cultures. Modern medicine has achieved success due to the effectiveness of traditional medicines, which are derived from medicinal plants. The objective of this study was to determine whether components of natural origin from Iranian medicinal plants have an antiviral effect that can prevent humans from this coronavirus infection using the most reliable molecular docking method; in our case, we focused on the main protease (Mpro) and a receptor-binding domain (RBD). The results of molecular docking showed that among 169 molecules of natural origin from common Iranian medicinal plants, 20 molecules (chelidimerine, rutin, fumariline, catechin gallate, adlumidine, astragalin, somniferine, etc.) can be proposed as inhibitors against this coronavirus based on the binding free energy and type of interactions between these molecules and the studied proteins. Moreover, a molecular dynamics simulation study revealed that the chelidimerine-Mpro and somniferine-RBD complexes were stable for up to 50 ns below 0.5 nm. Our results provide valuable insights into this mechanism, which sheds light on future structure-based designs of high-potency inhibitors for SARS-CoV-2.

Ebsulfur and Ebselen as highly potent scaffolds for the development of potential SARS-CoV-2 antivirals

The emerging COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has raised a global catastrophe. To date, there is no specific antiviral drug available to combat this virus, except the vaccine. In this study, the main protease (Mpro) required for SARS-CoV-2 viral replication was expressed and purified. Thirty-six compounds were tested as inhibitors of SARS-CoV-2 Mpro by fluorescence resonance energy transfer (FRET) technique. The half-maximal inhibitory concentration (IC50) values of Ebselen and Ebsulfur analogs were obtained to be in the range of 0.074-0.91 μM. Notably, the molecules containing furane substituent displayed higher inhibition against Mpro, followed by Ebselen 1i (IC50 = 0.074 μM) and Ebsulfur 2k (IC50 = 0.11 μM). The action mechanism of 1i and 2k were characterized by enzyme kinetics, pre-incubation and jump dilution assays, as well as fluorescent labeling experiments, which suggested that both compounds covalently and irreversibly bind to Mpro, while molecular docking suggested that 2k formed an SS bond with the Cys145 at the enzymatic active site. This study provides two very potent scaffolds Ebsulfur and Ebselen for the development of covalent inhibitors of Mpro to combat COVID-19.

Structure-based investigation of MARK4 inhibitory potential of Naringenin for therapeutic management of cancer and neurodegenerative diseases

MAP/microtubule affinity-regulating kinase 4 (MARK4) is a member of serine/threonine kinase family and considered an attractive drug target for many diseases. Screening of Indian Medicinal Plants, Phytochemistry, and Therapeutics (IMPPAT) using virtual high-throughput screening coupled with enzyme assay suggested that Naringenin (NAG) could be a potent inhibitor of MARK4. Structure-based molecular docking analysis showed that NAG binds to the critical residues found in the active site pocket of MARK4. Furthermore, molecular dynamics (MD) simulation studies for 100 ns have delineated the binding mechanism of NAG to MARK4. Results of MD simulation suggested that binding of NAG further stabilizes the structure of MARK4 by forming a stable complex. In addition, no significant conformational change in the MARK4 structure was observed. Fluorescence binding and isothermal titration calorimetric measurements revealed an excellent binding affinity of NAG to MARK4 with a binding constant (K) = 0.13 × 10⁶ M^-1 obtained from fluorescence binding studies. Further, enzyme inhibition studies showed that NAG has an admirable IC₅₀ value of 4.11 µM for MARK4. Together, these findings suggest that NAG could be an effective MARK4 inhibitor that can potentially be used to treat cancer and neurodegenerative diseases.

A computational approach for rational discovery of inhibitors for non-structural protein 1 of SARS-CoV-2

Background

Non-structural protein 1 (Nsp1), a virulence agent of SARS-CoV-2, has emerged as an important target for drug discovery. Nsp1 shuts down the host gene function by associating with the 40S ribosomal subunit.

Methods

Molecular interactions, drug-likeness, physiochemical property predictions, and robust molecular dynamics (MD) simulations were employed to discover novel Nsp1 inhibitors. In this study, we evaluated a series of molecules based on the plant (Cedrus deodara) derived α,β,γ-Himachalenes scaffolds.

Results

The results obtained from estimated affinity and ligand efficiency suggested that BCH10, BCH15, BCH16, and BCH17 could act as potential inhibitors of Nsp1. Moreover, MD simulations comprising various MD driven time-dependent analyses and thermodynamic free energy calculations also suggested stable protein-ligand complexes and strong interactions with the binding site. Furthermore, the selected molecules passed drug likeliness parameters and the physiochemical property analysis showed acceptable bioactivity scores.

Conclusion

The structural parameters of dynamic simulations revealed that the reported molecules could act as lead compounds against SARS-CoV-2 Nsp1 protein.

Xanthone glucoside 2-β-D-glucopyranosyl-1,3,6,7-tetrahydroxy-9H-xanthen-9-one binds to the ATP-binding pocket of glycogen synthase kinase 3β and inhibits its activity: implications in prostate cancer and associated cardiovascular disease risk

Glycogen synthase kinase 3 (GSK3) is a serine/threonine kinase which in the presence of ATP in its ATP-binding pocket transfers a phosphate to a primed substrate. GSK3β is an isoform of GSK3 which has been projected as a potent therapeutic target in human diseases including cancers and metabolic syndrome. Incidentally, cardiovascular disease is a common cause of non-cancer related deaths in prostate cancer (PCa) patients, mainly due to the effects of androgen-deprivation therapy (ADT), a mainstay for PCa treatment. Several small molecular inhibitors of GSK3 are either ATP-competitive (bind to the ATP-binding pocket), or non-ATP-competitive inhibitors (binding to the substrate-binding site of the enzyme). In this study, 2-β-D-glucopyranosyl-1,3,6,7-tetrahydroxy-9H-xanthen-9-one (βDGT), a natural xanthonoid present in many plant species, is reported to bind to the ATP-binding pocket of GSK3β and inhibit its activity, as demonstrated by the molecular docking and molecular dynamics simulation analysis and experimental validation in vitro. A comparison of the binding affinities with five known ATP-competitive inhibitors of GSK3β suggested similarity in binding site residues in the ATP-binding pocket of the enzyme. The optimum inhibitory concentration of the xanthonoid as determined by the luminescent kinase assay was 200 µM. The study envisages the use of βDGT as a natural ATP-competitive inhibitor of GSK3β and implicates its use in PCa patients on ADT, a cardiovascular disease risk, and other pathological conditions where GSK3 inhibition may be clinically important. HighlightsGSK3β is a multifaceted kinase known for its role in cancers, cardiovascular, and other diseases.In this study, βDGT, a xanthonoid, is reported to bind to the ATP-binding pocket of GSK3β.A comparison of βDGT binding with 5 known ATP-competitive inhibitors of GSK3β suggested the involvement of residues at the ATP binding site.The binding site analysis suggested an ATP-competitive mechanism of enzyme inhibition.Study envisages the use of βDGT as a natural ATP-competitive inhibitor of GSK3β and implicates its use in prostate cancer patients on androgen-deprivation therapy, a cardiovascular disease risk, and other pathological conditions.Communicated by Ramaswamy H. Sarma.

The molecular interaction of human anti-apoptotic proteins and in silico ADMET, drug-likeness and toxicity computation of N-cyclohexylmethacrylamide

Cancer is an uncontrolled growth of normal cells and apoptosis has an important role in cancer progression and cancer treatment. Antiapoptotic proteins are overexpressed in several tumors including breast, brain, lung cancer cells. The protein-ligand interaction has a critical role in drug designing. The present study aims to evaluate the interaction of synthesized N-cyclohexylmethacrylamide (NCMA) with proteins using in silico molecular docking and toxicity analyses. The NCMA monomer was synthesized and characterized by our team, previously. Kinetics stability, binding affinities and toxic potential of protein-NCMA complex were examined with the aid of molecular simulation. The toxicity results of this study indicate that NCMA is a sample with low toxic potential. According to the docking results, NCMA may be a drug active substance with chemical modifications and toxicity results support this situation. The drug-likeness and ADMET parameters were screened properties of NCMA.

In silico drug designing for COVID-19: an approach of high-throughput virtual screening, molecular, and essential dynamics simulations

Severe acute respiratory syndrome-coronavirus2 (SARS-CoV2), a new coronavirus has emerged in Wuhan city of China, December 2019 causing pneumonia named Coronavirus disease-19 (COVID-19), which has spread to the entire world. By January 2021, number of confirmed cumulative cases crossed ∼104 million worldwide. Till date, no effective treatment or drug is available for this virus. Availability of X-ray structures of SARS-CoV2 main protease (M^pro) provides the potential opportunity for structure-based drug designing. Here, we have made an attempt to do computational drug design by targeting main protease of SARS-CoV2. High-throughput virtual screening of million molecules and natural compounds databases were performed followed by docking. After that, the protein-ligand complexes were optimized and rescoring of binding energies were accomplished through molecular dynamics simulation and Molecular mechanics Poisson Boltzmann surface area approaches, respectively. In addition, conformational effect of various ligands on protein was also examined through essential dynamics simulation. Three compounds namely ZINC14732869, ZINC19774413, and ZINC19774479 were finally filtered that displayed better binding affinities than N3 (known) inhibitor and formed conformationally stable complexes. Hence, the current study features the potential novel inhibitors against main protease of SARS-CoV2 which might provide an effective therapeutic strategy against COVID-19.Communicated by Ramaswamy H. Sarma.

Investigation of Cyc1 protein structure stability after H53I mutation using computational approaches to improve redox potential

Acidithiobacillus ferrooxidans (Af) is an acidophilic bacterium that grows in rigid surroundings and gets its own energy from the oxidation of Fe²⁺ to Fe³⁺. These bacteria are involved in the bioleaching process. Cyc₁ is a periplasmic protein with a crucial role in electron transportation in the respiratory chain. His53 of the Cyc₁ protein, involved in electron transfer to CoxB, was selected for mutation and bioinformatics studies. His53 was substituted by Ile using PyMol software. Molecular dynamics simulations were performed for wild and mutant types of Cyc₁ protein. The conformational changes of mutated protein were studied by analyzing RMSD, RMSF, SASA, Rg, H Bond, and DSSP. The results of the RMSF analysis indicated an increase in the flexibility of the ligand in the mutant. Finally, active site instability leads to an increase in the value of E₀ at the mutation point and improving electron transfer. On the other, His53 in Cyc₁ is interconnected to Glu126 in CoxB through the water molecule (W76) and hydrogen bonding. In the H53I mutation, there was a decrease in the distance between H2O 2030, 2033, and isoleucine 53, and subsequently, the distance to the water molecule 76 between the two proteins was reduced and strengthens the hydrogen bond between Cyc₁ and CoxB, finally improves electron transfer and the bioleaching process.

Effects of the association of the αv β8 lower legs on integrin ligand binding

Many integrins transmit signals through global conformational changes. However, it is unclear whether integrin α_v β₈ adopts a similar mechanism during integrin activation and signaling on the cell surface. Here, we showed that disulfide-bonded mutants, which prevented integrin α_v β₈ lower leg dissociation, bound ligands with similar level as the wild-type protein, suggesting that α_v β₈ ligand binding did not require lower leg disassociation. We further showed that the N-glycosylation mutant at the interface between the β I and hybrid domains did not affect ligand binding, suggesting that the α_v β₈ open headpiece was not present on the cell surface. We proposed that α_v β₈ integrin may adopt only one state, that is, the extended conformation with a closed headpiece. Our results showed that two lower legs retained heterodimeric interfaces, and this association might be important for stabilizing integrin in the extended conformation. Therefore, α_v β₈ may not transmit bidirectional signals across the plasma membrane but instead may serve as an anchoring site with high affinity and high accessibility for extracellular ligands.

In silico drug repositioning based on the integration of chemical, genomic and pharmacological spaces

BACKGROUND

Drug repositioning refers to the identification of new indications for existing drugs. Drug-based inference methods for drug repositioning apply some unique features of drugs for new indication prediction. Complementary information is provided by these different features. It is therefore necessary to integrate these features for more accurate in silico drug repositioning.

RESULTS

In this study, we collect 3 different types of drug features (i.e., chemical, genomic and pharmacological spaces) from public databases. Similarities between drugs are separately calculated based on each of the features. We further develop a fusion method to combine the 3 similarity measurements. We test the inference abilities of the 4 similarity datasets in drug repositioning under the guilt-by-association principle. Leave-one-out cross-validations show the integrated similarity measurement IntegratedSim receives the best prediction performance, with the highest AUC value of 0.8451 and the highest AUPR value of 0.2201. Case studies demonstrate IntegratedSim produces the largest numbers of confirmed predictions in most cases. Moreover, we compare our integration method with 3 other similarity-fusion methods using the datasets in our study. Cross-validation results suggest our method improves the prediction accuracy in terms of AUC and AUPR values.

CONCLUSIONS

Our study suggests that the 3 drug features used in our manuscript are valuable information for drug repositioning. The comparative results indicate that integration of the 3 drug features would improve drug-disease association prediction. Our study provides a strategy for the fusion of different drug features for in silico drug repositioning.

Could Dermaseptin Analogue be a Competitive Inhibitor for ACE2 Towards Binding with Viral Spike Protein Causing COVID19?: Computational Investigation

Initial phase of COVID-19 infection is associated with the binding of viral spike protein S1 receptor binding domain (RBD) with the host cell surface receptor, ACE2. Peptide inhibitors typically interact with spike proteins in order to block its interaction with ACE2, and this knowledge would promote the use of such peptides as therapeutic scaffolds. The present study examined the competitive inhibitor activity of a broad spectrum antimicrobial peptide, Dermaseptin-S4 (S4) and its analogues. Three structural S4 analogues viz., S4 (K₄), S4 (K₂₀) and S4 (K₄K₂₀) were modelled by substituting charged lysine for non-polar residues in S4 and subsequently, docked with S1. Further, the comparative analysis of inter-residue contacts and non-covalent intermolecular interactions among S1–S4 (K₄), S1–S4 (K₄K₂₀) and S1–ACE2 complexes were carried out to explore their mode of binding with S1. Interestingly, S1–S4 (K₄) established more inter-molecular interactions compared to S4 (K₄K₂₀) and S1–ACE2. In order to substantiate this study, the normal mode analysis (NMA) was conducted to show how the structural stability of the flexible loop region in S1 is affected by atomic displacements in unbound S1 and docked complexes. Markedly, the strong interactions consistently maintained by S1–S4 (K₄) complex revealed their conformational transition over the harmonic motion period. Moreover, S1–S4 (K₄) peptide complex showed a higher energy deformation profile compared to S1–S4 (K₄K₂₀), where the higher energy deformation suggests the rigidity of the docked complex and thus it’s harder deformability, which is also substantiated by molecular dynamics simulation. In conclusion, S1–S4 (K₄) complex has definitely exhibited a functionally significant dynamics compared to S1–ACE2 complex; this peptide inhibitor, S4 (K₄) will need to be considered as the best therapeutic scaffold to block SARS-CoV-2 infection.

A computational perspective on the dynamic behaviour of recurrent drug resistance mutations in the pncA gene from Mycobacterium tuberculosis

Tuberculosis is still one of the top 10 causes of death worldwide, particularly with the emergence of multidrug-resistant tuberculosis. As the most effective first-line anti-tuberculosis drug, pyrazinamide also develops resistance due to the mutation in the pncA gene. Among these mutations, seven mutations at positions F94L, F94S, K96N, K96R, G97C, G97D, and G97S are classified as high-level resistance mutations. However, the resistance mechanism of Mtb to PZA (pyrazinamide) caused by these mutations is still unclear. In this work, we combined molecular dynamics simulation, molecular mechanics/generalized-Born surface area calculation, principal component analysis, and free energy landscape analysis to explore the resistance mechanism of Mtb to PZA due to F94L, F94S, K96N, K96R, G97C, G97D, and G97S mutations, as well as compare interaction changes in wild-type and mutant PZA-bound complexes. The results of molecular mechanics/generalized-Born surface area calculations indicated that the binding free energy of PZA with seven mutants decreased. In mutant systems, the most significant interactions with His137 and Cys138 were lost. Besides, PCA and FEL confirmed significant variations in the protein dynamics during the simulation specifically by altering the Fe2+ binding and its destabilization. Furthermore, mutants also flipped the β-sheet 2, which also affects the binding of Fe2+ and PZA. In the G97D mutant, reported as most lethal, mutation causes the binding pocket to change considerably, so that the position of PZA has a large movement in the binding pocket. In this study, the resistance mechanism of PZA at the atomic level is proposed. The proposed drug-resistance mechanism will provide valuable guidance for the design of anti-tuberculosis drugs.

Mushroom-derived bioactive compounds potentially serve as the inhibitors of SARS-CoV-2 main protease: An in silico approach

Background and aim

Coronavirus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has now become the world pandemic. There is a race to develop suitable drugs and vaccines for the disease. The anti-HIV protease drugs are currently repurposed for the potential treatment of COVID-19. The drugs were primarily screened against the SARS-CoV-2 main protease. With an urgent need for safe and effective drugs to treat the virus, we have explored natural products isolated from edible and medicinal mushrooms that have been reported to possess anti-HIV protease.

Experimental procedures

We have examined 36 compounds for their potential to be SARS-CoV-2 main protease inhibitors using molecular docking study. Moreover, drug-likeness properties including absorption, distribution, metabolism, excretion and toxicity were evaluated by in silico ADMET analysis.

Results

Our AutoDock study showed that 25 of 36 candidate compounds have the potential to inhibit the main viral protease based on their binding affinity against the enzyme's active site when compared to the standard drugs. Interestingly, ADMET analysis and toxicity prediction revealed that 6 out of 25 compounds are the best drug-like property candidates, including colossolactone VIII, colossolactone E, colossolactone G, ergosterol, heliantriol F and velutin.

Conclusion

Our study highlights the potential of existing mushroom-derived natural compounds for further investigation and possibly can be used to fight against SARS-CoV-2 infection.

Taxonomy classification by evise

Disease, Infectious Disease, Respiratory System Disease, Covid-19, Traditional Medicine, Traditional Herbal Medicine, Phamaceutical Analysis.

Protein–ligand free energies of binding from full-protein DFT calculations: convergence and choice of exchange–correlation functional

Quantum mechanical binding free energies based on thousands of full-protein DFT calculations are tractable, reproducible and converge well.

Enhancing the catalytic activity of a GH5 processive endoglucanase from Bacillus subtilis BS-5 by site-directed mutagenesis

Processive endoglucanases possess both endo- and exoglucanase activity, making them attractive discovery and engineering targets. Here, a processive endoglucanase EG5C-1 from Bacillus subtilis was employed as the starting point for enzyme engineering. Referring to the complex structure information of EG5C-1 and cellohexaose, the amino acid residues in the active site architecture were identified and subjected to alanine scanning mutagenesis. The residues were chosen for a saturation mutagenesis since their variants showed similar activities to EG5C-1. Variants D70Q and S235W showed increased activity towards the substrates CMC and Avicel, an increase was further enhanced in D70Q/S235W double mutant, which displayed a 2.1- and 1.7-fold improvement in the hydrolytic activity towards CMC and Avicel, respectively. In addition, kinetic measurements showed that double mutant had higher substrate affinity (K_m) and a significantly higher catalytic efficiency (k_cat/K_m). The binding isotherms of wild-type EG5C-1 and double mutant D70Q/S235W suggested that the binding capability of EG5C-1 for the insoluble substrate was weaker than that of D70Q/S235W. Molecular dynamics simulations suggested that the collaborative substitutions of D70Q and S235W altered the hydrogen bonding network within the active site architecture and introduced new hydrogen bonds between the enzyme and cellohexaose, thus enhancing both substrate affinity and catalytic efficiency.

Antimicrobial, cytotoxicity, molecular modeling and DNA cleavage/binding studies of zinc-naproxen complex: switching DNA binding mode of naproxen by coordination to zinc ion

The intercalation DNA binding mode of the naproxen, a non-steroidal anti-inflammatory drug, has been reported previously. In this study, calf thymus deoxyribonucleic acid (CT-DNA) binding of zinc-naproxen complex, [Zn(naproxen)₂(MeOH)₂], at physiological pH has been investigated by multi-spectroscopic techniques and molecular docking. Zinc-naproxen complex displays significant binding property to the CT-DNA (K_b = 0.2 × 10⁵ L.mol^-1). All of the experimental results; relative increasing in viscosity of CT-DNA and fluorimetric studies using ethidium bromide (EB) and Hoechst 33258 probes, are indicative of groove binding mode of zinc-naproxen complex to CT-DNA. These results show that the coordination of naproxen to zinc metal switches the mode of binding from intercalation to groove. The molecular modeling also shows that the complex binds to the AT-rich region of minor groove of DNA. Structural and topography changes of DNA in interaction with the complex by atomic force microscopy (AFM) indicated that CT-DNA becomes swollen after interaction. The pUC18 plasmid DNA cleavage ability of zinc-naproxen complex by gel electrophoresis experiments revealed that zinc-naproxen complex cleaved supercoiled pUC18 plasmid DNA to nicked DNA. The cytotoxicity of the zinc complex performed by MTT method on HT29 and MCF7 cancer cell lines and on HEK 293 normal cell lines indicates that zinc complex has no cytotoxic effect on both HT29 and MCF7 cell lines but has better cytotoxicity effect on HEK 293 cell lines compared to cisplatin standard drug. The antimicrobial activity of the complex against Staphylococcus aureus and Escherichia coli bacteria revealed the high antimicrobial activity of the complex.Communicated by Ramaswamy H. Sarma.

RAB5A is associated with genes involved in exosome secretion: Integration of bioinformatics analysis and experimental validation

Exosomes, as cell-cell communicators with an endosomal origin, are involved in the progression of various diseases. RAB5A, a member of the small Rab GTPases family, which is well known as a key regulator of cellular endocytosis, is expected to be involved in exosome secretion. Here, we found the impact of RAB5A on exosome secretion from human hepatocellular carcinoma cell line using a rapid yet reliable bioinformatics approach followed by experimental analysis. Initially, RAB5A and exosome secretion-related genes were gathered from bioinformatics tools, namely, CTD, COREMINE, and GeneMANIA; and published papers. Protein-protein interaction (PPI) was then constructed by the Search Tool for Retrieval of Interacting Genes (STRING) database. Among them, several genes with different combined scores were validated by the real-time quantitative polymerase chain reaction (RT-qPCR) in stable RAB5A knockdown cells. Thereafter, to validate the bioinformatics results functionally, the impact of RAB5A knockdown on exosome secretion was evaluated. Bioinformatics analysis showed that RAB5A interacts with 37 genes involved in exosome secretion regulatory pathways. Validation by RT-qPCR confirmed the association of RAB5A with candidate interacted genes and interestingly showed that even medium to low combined scores of the STRING database could be experimentally valid. Moreover, the functional analysis demonstrated that the stable silencing of RAB5A could experimentally decrease exosome secretion. In conclusion, we suggest RAB5A as a regulator of exosome secretion based on our bioinformatics approach and experimental analysis. Also, we propose the usage of PPI-derived from the STRING database regardless of their combined scores in advanced bioinformatics analysis.

Transcriptomic analysis of glioblastoma multiforme providing new insights into GPR17 signaling communication

Glioblastoma Multiforme (GBM) is one of the most aggressive malignant tumors in the central nervous system, which arises due to the failure or crosstalk in the signaling networks. GPR17, an orphan G protein-coupled receptor is anticipated to be associated with the biology of the GBM disease progression. In the present study, we have identified the differential expressions of around 170 genes along with GPR17 through the RNA-Seq analysis of 169 GBM samples. Coordinated expression patterns of all other gene products with this receptor were analysed using gene ontology and protein-protein interaction data. Several crucial signaling components and genes that play a significant role in tumor progression have been identified among which GPR17 was found to be significantly interacting with about 30 different pathways. High-throughput molecular docking of GPR17 by homology-based model against differentially expressed proteins, showed effective recognition and binding of PX, SH3, and Ig-like domains besides Gi protein. Pathways of PI3, Src, Ptdn, Ras, cytoplasmic tyrosine kinases, phospholipases, nexins and other proteins possessing these structural domains are identified as critical signaling components of the complex GBM signaling network. Our findings also provide a mechanistic insight of GPR17-T0510-3657 interaction, which potentially regulates the interaction of PX domain and helical mPTS recognition domain-containing proteins. Overall, our results demonstrate that GPR17 mediated signaling networks could be used as a therapeutic target for GBM.Communicated by Ramaswamy H. Sarma.

Discovery and in silico evaluation of aminoarylbenzosuberene molecules as novel checkpoint kinase 1 inhibitor determinants

Checkpoint kinase 1 (CHK1) is an essential kinase with a critical function in cell cycle arrest. Several potent inhibitors targeting CHK1 have been published, but most of them have failed in clinical trials. Acknowledging the emerging consequence of CHK1 inhibitors in medication of cancer, there is a demand for widening the chemical range of CHK1 inhibitors. In this research, we considered a set of in-house plant based semi-synthetic aminoarylbenzosuberene molecules as potential CHK1 inhibitors. Based on a combined computational research that consolidates molecular docking and binding free energy computations we recognized the crucial determinants for their receptor binding. The drug likeness of these molecules were also scrutinized based on their toxicity and bioavailibilty profile. The computational strategy indicates that the Bch10 could be regarded as a potential CHK1 inhibitor in comparison with top five co-crystallize molecules. Bch10 signifies a promising outlet for the development of potent inhibitors for CHK1.

Identification of polyphenols from Broussonetia papyrifera as SARS CoV-2 main protease inhibitors using in silico docking and molecular dynamics simulation approaches

The current COVID-19 pandemic is caused by SARS CoV-2. To date, ∼463,000 people died worldwide due to this disease. Several attempts have been taken in search of effective drugs to control the spread of SARS CoV-2 infection. The main protease (Mpro) from SARS CoV-2 plays a vital role in viral replication and thus serves as an important drug target. This Mpro shares a high degree of sequence similarity (>96%) with the same protease from SARS CoV-1 and MERS. It was already reported that Broussonetia papyrifera polyphenols efficiently inhibit the catalytic activity of SARS CoV-1 and MERS Mpro. But whether these polyphenols exhibit any inhibitory effect on SARS CoV-2 Mpro is far from clear. To understand this fact, here we have adopted computational approaches. Polyphenols having proper drug-likeness properties and two repurposed drugs (lopinavir and darunavir; having binding affinity −7.3 to −7.4 kcal/mol) were docked against SARS CoV-2 Mpro to study their binding properties. Only six polyphenols (broussochalcone A, papyriflavonol A, 3'-(3-methylbut-2-enyl)-3',4',7-trihydroxyflavane, broussoflavan A, kazinol F and kazinol J) had interaction with both the catalytic residues (His41 and Cys145) of Mpro and exhibited good binding affinity (−7.6 to −8.2 kcal/mol). Molecular dynamic simulations (100 ns) revealed that all Mpro-polyphenol complexes are more stable, conformationally less fluctuated; slightly less compact and marginally expanded than Mpro-darunavir/lopinavir complex. Even the number of intermolecular H-bond and MM-GBSA analysis suggested that these six polyphenols are more potent Mpro inhibitors than the two repurposed drugs (lopinavir and darunavir) and may serve as promising anti-COVID-19 drugs.

Communicated by Ramaswamy H. Sarma