Targeting the protein-protein interface pocket of Aurora-A-TPX2 complex: rational drug design and validation

Unveiling Potential Targeted Therapeutic Opportunities for Co-Overexpressed Targeting Protein for Xklp2 and Aurora-A Kinase in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the most prevalent and leading causes of cancer deaths globally, with limited diagnostic and clinically significant therapeutic targets. Identifying the genes and processes involved in developing and progressing LUAD is crucial for developing effective targeted therapeutics and improving patient outcomes. Therefore, the study aimed to explore the RNA sequencing data of LUAD from The Cancer Genome Atlas (TCGA) and gene expression profile datasets involving GSE10072, GSE31210, and GSE32863 from the Gene Expression Omnibus (GEO) databases. The differential gene expression and the downstream analysis determined clinically significant biomarkers using a network-based approach. These therapeutic targets predominantly enriched the dysregulation of mitotic cell cycle regulation and revealed the co-overexpression of Aurora-A Kinase (AURKA) and Targeting Protein for Xklp2 (TPX2) with high survival risk in LUAD patients. The hydrophobic residues of the AURKA-TPX2 interaction were considered as the target site to block the autophosphorylation of AURKA during the mitotic cell cycle. The tyrosine kinase inhibitor (TKI) dacomitinib demonstrated the strong binding potential to hinder TPX2, shielding the AURKA destabilization. This in silico study lays the foundation for repurposing targeted therapeutic options to impede the Protein-Protein Interactions (PPIs) in LUAD progression and aid in future translational investigations.

Targeting COVID-19 and varicocele by blocking inflammasome: Ligand-based virtual screening

COVID-19 is a new generation of outbreaks that invade not only local emerging region, continental but also the whole globe. Varicocele on the other hand, is a testicular vascular disease that underlies 40 % of male infertility cases. Fortunately, the two diseases can be blocked through targeting one common target, NLRP3 inflammasome. Upon searching for similar drugs that gained FDA-approval in ChEMBL library along with examining their potential blockade of the receptor through docking using CB-DOCK-2, three potential approved drugs can be repurposed, ChEMBL 4297185, ChEMBL 1201749, ChEMBL 1200545 which had binding energy of -9.8 and -9.7 kcal/mol (stronger than the reference inhibitor, -9.3 kcal/mol). Also, ADME profile of the top 3 drugs showed better attributes. Also, the simulated proteins exhibited stable pattern with strong free binding energies. Among the potential inhibitor drugs ChEMBL 4297185 was found to remain inside the binding site of the protein during the 200 ns simulation time. Hence, it is anticipated to have the highest binding and thus inhibition potential against the protein. The suggested drugs, especially ChEMBL 4297185, are potentially repurposable toward treating COVID-19 and varicocele which deserve further experimental validation.

Label-Free Mapping of Multivalent Binding Pathways with Ligand-Receptor-Anchored Nanopores

Understanding single-molecule multivalent ligand-receptor interactions is crucial for comprehending molecular recognition at biological interfaces. However, label-free identifications of these transient interactions during multistep binding processes remains challenging. Herein, we introduce a ligand-receptor-anchored nanopore that allows the protein to maintain structural flexibility and favorable orientations in native states, mapping dynamic multivalent interactions. Using a four-state Markov chain model, we clarify two concentration-dependent binding pathways for the Omicron spike protein (Omicron S) and soluble angiotensin-converting enzyme 2 (sACE2): sequential and concurrent. Real-time kinetic analysis at the single-monomeric subunit level reveals that three S1 monomers of Omicron S exhibit a consistent and robust binding affinity toward sACE2 (-13.1 ± 0.2 kcal/mol). These results highlight the enhanced infectivity of Omicron S compared to other homologous spike proteins (WT S and Delta S). Notably, the preceding binding of sACE2 to Omicron S facilitates the subsequent binding steps, which was previously obscured in bulk measurements. Our single-molecule studies resolve the controversy over the disparity between the measured spike protein binding affinity with sACE2 and the viral infectivity, offering valuable insights for drug design and therapies.

Mechanistic studies on substrate inhibition and substrate activation of ribonuclease A: experimental and in silico investigation

Ribonuclease A (RNase A) is an endonuclease that plays a significant role in antimicrobial activity by the cleavage and hydrolysis of ssRNA. In this study, the interactions between RNase A and cCMP have been investigated, through molecular docking and a 200 ns molecular dynamics simulation. The enzyme kinetic properties were analyzed using saturation curve, Eadie-Hofstee, and Hill plots. The docking results indicate that the cCMP-RNase A complexes are stabilized through hydrophobic interaction, hydrogen bonding, and π-π stacking interaction. The most binding site is observed in the catalytic cleft, especially at residue His12 and His119. Enzyme-ligand docking study indicates that cCMP binds to residues located in the internal cavity of the catalytic site and shows three phases of binding modes. The analysis of MD simulations shows that RMSD, Rg, and RMSF reach equilibrium. The ΔGbinding of the cCMP-RNase A was negative (-619.673 kJ/mol), The comparison between the results pre and post MD simulation showed that ΔGbinding after MD simulation causes to more stable the structure than before simulation. Experimental methods such as saturation, Eadie-Hofstee, and Hill plots confirm theoretical data and show three phases of binding modes as well. Two significant events are demonstrated in the interaction between RNase A and cCMP: substrate activation and substrate inhibition are observed in concentrations below and above 0.8 mM, respectively, for cCMP. Choosing the appropriate concentration of cCMP is very important in investigation of ribonuclease A's catalytic behavour, especially for exploration of the effects of some drugs on biological behaviours related to ribonuclease A.Communicated by Ramaswamy H. Sarma.

In silico analysis of the different variable domain oriented single-chain variable fragment antibody-antigen complexes

Single-chain variable fragment (scFv) antibodies hold great potential as diagnostic tools and therapeutic agents, especially for tumor cells. Since these applications require their production with improved properties, the design strategy of scFvs is crucial for their active, soluble, and high yield expression with high affinity towards their antigens. The order of VL and VH domains is one of the important parameters that affect the expression and binding affinity properties of scFvs. In addition, the optimum order of VL and VH domains could change for each scFv. In the present study, we used computer simulation tools to evaluate the effect of variable domain orientation on structure, stability, interacting residues of scFvs, and binding free energies of scFv-antigen complexes. We selected anti-HER2 scFv, which is specific for human epidermal growth receptor 2 (HER2) overexpressed in breast cancer, and anti-IL-1β scFv against IL-1β which is an important inflammatory biomarker, as model scFvs. Molecular dynamics simulations of the scFv-antigen complexes for 100 ns resulted in stability and compactness for both scFv constructs. Interaction and binding free energies calculated by the Molecular Mechanics-Poisson-Boltzmann Surface Area (MM-PBSA) approach suggested that the relative binding energies of anti-HER2 scFv-VLVH and anti-HER2 scFv-VHVL constructs had similar binding affinity towards HER2, while a relatively more negative binding free energy obtained between anti-IL-1β scFv-VHVL and IL-1β pointed to a higher binding affinity. The in silico approach and the results obtained here could be applied as a guide for future experimental interaction studies for highly specific scFvs used as biotechnological tools.Communicated by Ramaswamy H. Sarma.

Design of multi-epitope chimeric vaccine against Monkeypox virus and SARS-CoV-2: A vaccinomics perspective

The current era we experience is full with pandemic infectious agents that no longer threatens the major local source but the whole globe. Almost the most emerging infectious agents are severe acute respiratory syndrome coronavirus-2 (SARS CoV-2), followed by monkeypox virus (MPXV). Since no approved antiviral drugs nor licensed active vaccines are yet available, we aimed to utilize immunoinformatics approach to design chimeric vaccine against the two mentioned viruses. This is the first study to deal with design divalent vaccine against SARS-CoV-2 and MPXV. ORF8, E and M proteins from Omicron SARS-CoV-2 and gp182 from MPXV were used as the protein precursor from which multi-epitopes (inducing B-cell, helper T cells, cytotoxic T cells and interferon-ɣ) chimeric vaccine was contrived. The structure of the vaccine construct was predicted, validated, and docked to toll-like receptor-2 (TLR-2). Moreover, its sequence was also used to examine the immune simulation profile and was then inserted into the pET-28a plasmid for in silico cloning. The vaccine construct was probable antigen (0.543) and safe (non-allergen) with strong binding energy to TLR-2 (-1169.8 kcal/mol) and found to have significant immune simulation profile. In conclusion, the designed chimeric vaccine was potent and safe against SARS-CoV-2 and MPXV, which deserves further consideration.

Structure-based virtual screening and molecular dynamics approaches to identify new inhibitors of Staphylococcus aureus sortase A

Staphylococcus aureus is a prevalent Gram-positive bacteria leading cause of a wide range of human pathologies. Moreover, antibiotic résistance of pathogenesis bacteria is one of the worldwide health problems. In Gram-positive bacteria, the enzyme of SrtA, is responsible for the anchoring of surface-exposed proteins to the cell wall peptidoglycan. Because of its critical role in Gram-positive bacterial pathogenesis, SrtA is an attractive target for anti-virulence during drug development. To date, some SrtA inhibitors have been discovered most of them being derived from flavonoid compounds, like Myricetin. In order to provide potential hit molecules against SrtA for clinical use, we obtained a total of 293 compounds by performing in silico shape-based screening of compound libraries against Myristin as a reference structure. Employing molecular docking and scoring functions, the top 3 compounds Apigenin, Efloxate, and Compound 8261032 were screened by comparing their docking scores with Myricetin. Furthermore, MD simulations and MM-PBSA binding energy calculation studies revealed that only Compound 8261032 strongly binds to the catalytic core of the SrtA enzyme than Myricetin, and stable behavior was consistently observed in the docking complex. Compound 8261032 showed a good number of hydrogen bonds with SrtA and higher MM-PBSA binding energy when compared to all three molecules. Also, it makes strength interactions with Arg139 and His62, which are critical for SrtA biological activity. This study showed that the development of this inhibitor could be a fundamental strategy against resistant bacteria, but further studies in vitro are needed to confirm this claim.Communicated by Ramaswamy H. Sarma.

Alantolactone Induced Apoptosis and DNA Damage of Cervical Cancer through ATM/CHK2 Signaling Pathway

Traditional Chinese Medicine, known for its minimal side effects and significant clinical efficacy, has attracted considerable interest for its potential in cancer therapy. In particular, Inula helenium L. has demonstrated effectiveness in inhibiting a variety of cancers. This study focuses on alantolactone (ALT), a prominent compound from Inula helenium L., recognized for its anti-cancer capabilities across multiple cancer types. The primary objective of this study is to examine the influence of ALT on the proliferation, apoptosis, cell cycle, and tumor growth of cervical cancer (CC) cells, along with its associated signaling pathways. To determine protein expression alterations, Western blot analysis was conducted. Furthermore, an in vivo model was created by subcutaneously injecting HeLa cells into nude mice to assess the impact of ALT on cervical cancer. Our research thoroughly investigates the anti-tumor potential of ALT in the context of CC. ALT was found to inhibit cell proliferation and induce apoptosis in SiHa and HeLa cell lines, particularly targeting ataxia-telangiectasia mutated (ATM) proteins associated with DNA damage. The suppression of DNA damage and apoptosis induction when ATM was inhibited underscores the crucial role of the ATM/cell cycle checkpoint kinase 2 (CHK2) axis in ALT's anti-tumor effects. In vivo studies with a xenograft mouse model further validated ALT's effectiveness in reducing CC tumor growth and promoting apoptosis. This study offers new insights into how ALT combats CC, highlighting its promise as an effective anti-cervical cancer agent and providing hope for improved treatment outcomes for CC patients.

Molecular Drug Simulation and Experimental Validation of the CD36 Receptor Competitively Binding to Long-Chain Fatty Acids by 7-Ketocholesteryl-9-carboxynonanoate

Long-chain fatty acids (LCFAs) are one of the main energy-supplying substances in the body. LCFAs with different lengths and saturations may have contrasting biological effects that exacerbate or alleviate progress against a variety of systemic disorders of lipid metabolism in organisms. Nonalcoholic fatty liver disease is characterized by chronic inflammation and steatosis, mainly caused by the ectopic accumulation of lipids in the liver, especially LCFAs. CD36 is a scavenger receptor that recognizes and mediates the transmembrane absorption of LCFAs and is expressed in a variety of cells throughout the body. In previous studies, our group found that 7-ketocholesteryl-9-carboxynonanoate (oxLig-1) has the biological effect of targeting CD36 to inhibit oxidized low-density lipoprotein lipotoxicity-induced lipid metabolism disorder; it has an ω-carboxyl physiologically active center and is structurally similar to LCFAs. However, the biological mechanism of oxLig-1 binding to CD36 and competing for binding to different types of LCFAs is still not clear. In this study, molecular docking and molecular dynamics simulation were utilized to simulate and analyze the binding activity between oxLig-1 and different types of LCFAs to CD36 and confirmed by the enzyme-linked immunosorbent assay (ELISA) method. Absorption, distribution, metabolism, excretion, and toxicity (ADMET) platform was applied to predict the drug-forming properties of oxLig-1, and HepG2 cells model of oleic acid and nonalcoholic fatty liver disease (NAFLD) model mice were validated to verify the biological protection of oxLig-1 on lipid lowering. The results showed that there was a co-binding site of LCFAs and oxLig-1 on CD36, and the binding driving forces were mainly hydrogen bonding and hydrophobic interactions. The binding abilities of polyunsaturated LCFAs, oxLig-1, monounsaturated LCFAs, and saturated LCFAs to CD36 showed a decreasing trend in this order. There was a similar decreasing trend in the stability of the molecular dynamics simulation. ELISA results similarly confirmed that the binding activity of oxLig-1 to CD36 was significantly higher than that of typical monounsaturated and saturated LCFAs. ADMET prediction results indicated that oxLig-1 had a good drug-forming property. HepG2 cells model of oleic acid and NAFLD model mice study results demonstrated the favorable lipid-lowering biological effects of oxLig-1. Therefore, oxLig-1 may have a protective effect by targeting CD36 to inhibit the excessive influx and deposition of lipotoxicity monounsaturated LCFAs and saturated LCFAs in hepatocytes.

Active components and molecular mechanisms of Sagacious Confucius' Pillow Elixir to treat cognitive impairment based on systems pharmacology

BACKGROUND

Sagacious Confucius' Pillow Elixir (SCPE) is a common clinical prescription to treat cognitive impairment (CI) in East Asia.

OBJECTIVE

To predict the active components of SCPE, identify the associated signaling pathway, and explore the molecular mechanism using systems pharmacology and an animal study.

METHODS

Systems pharmacology and Python programming language-based molecular docking were used to select and analyze the active components and targets. Senescence-accelerated prone 8 mice were used as a CI model. The molecular mechanism was evaluated using the water maze test, neuropathological observation, cerebrospinal fluid microdialysis, and Western blotting.

RESULTS

Thirty active components were revealed by screening relevant databases and performing topological analysis. Additionally, 376 differentially expressed genes for CI were identified. Pathway enrichment analysis, protein-protein interaction (PPI) network analysis and molecular docking indicated that SCPE played a crucial role in modulating the PI3K/Akt/mTOR signaling pathway, and 23 SCPE components interacted with it. In the CI model, SCPE improved cognitive function, increased the levels of the neurotransmitter 5-hydroxytryptamine (5-HT) and metabolite 5-hydroxyindole acetic acid (5-HIAA), ameliorated pathological damage and regulated the PI3K/AKT/mTOR signaling pathway. SCPE increased the LC3-II/LC3-I, p-PI3K p85/PI3K p85, p-AKT/AKT, and p-mTOR/mTOR protein expression ratios and inhibited P62 expression in the hippocampal tissue of the CI model.

CONCLUSION

Our study revealed that 23 active SCPE components improve CI by increasing the levels of the neurotransmitter 5-HT and metabolite 5-HIAA, suppressing pathological injury and regulating the PI3K/Akt/mTOR signaling pathway to improve cognitive function.

Design of a Novel and Potent Multi-Epitope Chimeric Vaccine against Human Papillomavirus (HPV): An Immunoinformatics Approach

In the current era, our experience is full of pandemic infectious agents that no longer threaten the major local source but the whole globe. One such infectious agent is HPV, a sexually transmitted disease that can cause various clinical disorders, including benign lesions and cervical cancer. Since available vaccines still need further improvements in order to enhance efficacy, our goal was to design a chimeric vaccine against HPV using an immunoinformatics approach. For designing the vaccine, the sequence of HPV was retrieved, and then phylogenetic analysis was performed. Several CTL epitopes, HTL epitopes, and LBL epitopes were all predicted using bioinformatics tools. After checking the antigenicity, allergenicity, and toxicity scores, the best epitopes were selected for vaccine construction, and then physicochemical and immunological properties were analyzed. Subsequently, vaccine 3D structure prediction, refinement, and validation were performed. Molecular docking and dynamics simulation techniques were used to explore the interactions between the Toll-like receptor 2 and the modeled vaccine construct. To ensure the vaccine protein was expressed at a higher level, the construct was computationally cloned into the pET28a (+) plasmid. The molecular docking and simulation results showed that our designed vaccine is stable, of immunogenic quality, and has considerable solubility. Through in silico immune simulation, it was predicted that the designed polypeptide vaccine construct would trigger both humoral and cellular immune responses. The developed vaccine showed significant affinity for the TLR2 receptor molecule. However, additional laboratory research is required to evaluate its safety and efficacy.

Evaluation of Pathogenicity and Structural Alterations for the Mutations Identified in the Conserved Region of the C-Terminal Kinase Domain of Human-Ribosomal S6 Kinase 1

Human-ribosomal s6 kinase 1 (h-RSK1) is an effector kinase of the Ras/MAPK signaling pathway, which is involved in the regulation of the cell cycle, proliferation, and survival. RSKs comprise two functionally distinct kinase domains at the N-terminal (NTKD) and C-terminal (CTKD) separated by a linker region. The mutations in RSK1 may have the potential to provide an extra benefit to the cancer cell to proliferate, migrate, and survive. The present study focuses on evaluating the structural basis for the missense mutations identified at the C-terminal kinase domain of human-RSK1. A total of 139 mutations reported on RSK1 were retrieved from cBioPortal, where 62 were located at the CTKD region. Furthermore, 10 missense mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, Arg726Gln, His533Asn, Pro613Leu, Ser720Cys, Arg725Gln, and Ser732Phe were predicted to be deleterious using in silico tools. To our observation, these mutations are located in the evolutionarily conserved region of RSK1 and shown to alter the inter- and intramolecular interactions and also the conformational stability of RSK1-CTKD. The molecular dynamics (MD) simulation study further revealed that the five mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, and Arg726Gln showed maximum structural alterations in RSK1-CTKD. Thus, based on the in silico and MD simulation analysis, it can be concluded that the reported mutations may serve as potential candidates for further functional studies.

Structure restoration and aggregate inhibition of V30M mutant transthyretin protein by potential quinoline molecules

Transthyretin (TTR) is a tetrameric protein found in human plasma and cerebrospinal fluid that functions as a transporter of thyroxine (T4) and retinol. A mutation resulting in the substitution of valine to methionine at position 30 (V30M) is the most common mutation that destabilizes the tetramer structure of TTR protein resulting in a fatal neuropathy known as TTR amyloidosis. The V30M TTR-induced neuropathy can be inhibited through stabilization of the TTR tetramer by the binding of small molecules. We accessed the potential of in-house synthesized quinoline molecules to stabilize the V30M TTR structure and analyzed the impact of protein-ligand interactions through molecular docking, molecular dynamics (MD) simulations, steered MD, and umbrella sampling simulations. This study revealed that the binding of quinoline molecules reverted back the structural changes including the residual flexibility, changes in secondary structural elements, and also restored the alterations in the electrostatic surface potential induced by the V30M mutation. Further, the top-most 4G and 4R molecules were compared with an FDA-approved drug (Tafamidis) and a reference quinoline molecule 14C. Here, we intend to suggest that the quinoline molecules could revert the structural changes, cease tetramer dissociation, prevent abnormal oligomerization and therefore could be developed as an effective therapeutics against TTR amyloidosis.

miR-212-5p Suppresses Glioma Development via Targeting SUMO2

Recently, increasing studies have suggested that miRNAs play a significant role in the occurrence and development of glioma. More researches are needed to explore the role of miRNAs in glioma, which will help to find new therapeutic targets. miR-212-5p has been reported to be involved in the progression in many cancers. However, whether miR-212-5p has a regulative effect on glioma remains un clear. In this study, we aimed to explore the effect of miR-212-5p on glioma development and its mechanism. Here, we demonstrated that miR-212-5p was lowly expressed in glioma cell. miR-212-5p suppressed the glioma cell proliferation, inhibited the migratory and invasive capabilities and promoted apoptosis in glioma cells. Besides, miR-212-5p also inhibited tumor growth in vivo. We found small ubiquitin-like modifier 2 (SUMO2) was the target of miR-212-5p, and miR-212-5p suppressed SUMO2 expression to regulate the proliferation, migration, and apoptosis of glioma cells. These findings indicated that miR-212-5p may be a possible therapeutic target for the treatment for glioma.

Identification of Activated Cdc42-Associated Kinase Inhibitors as Potential Anticancer Agents Using Pharmacoinformatic Approaches

BACKGROUND

Activated Cdc42-associated kinase (ACK1) is essential for numerous cellular functions, such as growth, proliferation, and migration. ACK1 signaling occurs through multiple receptor tyrosine kinases; therefore, its inhibition can provide effective antiproliferative effects against multiple human cancers. A number of ACK1-specific inhibitors were designed and discovered in the previous decade, but none have reached the clinic. Potent and selective ACK1 inhibitors are urgently needed.

METHODS

In the present investigation, the pharmacophore model (PM) was rationally built utilizing two distinct inhibitors coupled with ACK1 crystal structures. The generated PM was utilized to screen the drug-like database generated from the four chemical databases. The binding mode of pharmacophore-mapped compounds was predicted using a molecular docking (MD) study. The selected hit-protein complexes from MD were studied under all-atom molecular dynamics simulations (MDS) for 500 ns. The obtained trajectories were ranked using binding free energy calculations (ΔG kJ/mol) and Gibb's free energy landscape.

RESULTS

Our results indicate that the three hit compounds displayed higher binding affinity toward ACK1 when compared with the known multi-kinase inhibitor dasatinib. The inter-molecular interactions of Hit1 and Hit3 reveal that compounds form desirable hydrogen bond interactions with gatekeeper T205, hinge region A208, and DFG motif D270. As a result, we anticipate that the proposed scaffolds might help in the design of promising selective ACK1 inhibitors.

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.

Computational investigation of natural compounds as potential main protease (Mpro) inhibitors for SARS-CoV-2 virus

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is significantly impacting human lives, overburdening the healthcare system and weakening global economies. Plant-derived natural compounds are being largely tested for their efficacy against COVID-19 targets to combat SARS-CoV-2 infection. The SARS-CoV-2 Main protease (Mpro) is considered an appealing target because of its role in replication in host cells. We curated a set of 7809 natural compounds by combining the collections of five databases viz Dr Duke's Phytochemical and Ethnobotanical database, IMPPAT, PhytoHub, AromaDb and Zinc. We applied a rigorous computational approach to identify lead molecules from our curated compound set using docking, dynamic simulations, the free energy of binding and DFT calculations. Theaflavin and ginkgetin have emerged as better molecules with a similar inhibition profile in both SARS-CoV-2 and Omicron variants.

Immunoinformatic-based design of immune-boosting multiepitope subunit vaccines against monkeypox virus and validation through molecular dynamics and immune simulation

Monkeypox virus is the causative agent of monkeypox disease, belonging to an orthopoxvirus genus, with a disease pattern similar to that of smallpox. The number of monkeypox cases have robustly increased recently in several countries around the world, potentially causing an international threat. Therefore, serious measures are indispensable to be taken to mitigate the spread of the disease and hence, under these circumstances, vaccination is the best choice to neutralize the monkeypox virus. In the current study, we used immunoinformatic approaches to target the L1R, B5R, and A33R proteins of the monkeypox virus to screen for immunogenic cytotoxic T-lymphocyte (CTL), helper T-lymphocyte (HTL), and B-cell epitopes to construct multiepitope subunit vaccines. Various online tools predicted the best epitope from immunogenic targets (L1R, B5R, and A33R) of monkeypox virus. The predicted epitopes were joined together by different linkers and subjected to 3D structure prediction. Molecular dynamics simulation analysis confirmed the proper folding of the modeled proteins. The strong binding of the constructed vaccines with human TLR-2 was verified by the molecular docking and determination of dissociation constant values. The GC content and codon adaptation index (CAI) values confirmed the high expression of the constructed vaccines in the pET-28a (+) expression vector. The immune response simulation data delineated that the injected vaccines robustly activated the immune system, triggering the production of high titers of IgG and IgM antibodies. In conclusion, this study provided a solid base of concept to develop dynamic and effective vaccines that contain several monkeypox virus-derived highly antigenic and nonallergenic peptides to control the current pandemic of monkeypox virus.

In vitro and in silico assessment of new beta amino ketones with antiplasmodial activity

BACKGROUND

Based on the current need for new drugs against malaria, our study evaluated eight beta amino ketones in silico and in vitro for potential antimalarial activity.

METHODS

Using the Brazilian Malaria Molecular Targets (BraMMT) and OCTOPUS® software programs, the pattern of interactions of beta-amino ketones was described against different proteins of P. falciparum and screened to evaluate their physicochemical properties. The in vitro antiplasmodial activities of the compounds were evaluated using a SYBR Green-based assay. In parallel, in vitro cytotoxic data were obtained using the MTT assay.

RESULTS

Among the eight compounds, compound 1 was the most active and selective against P. falciparum (IC50 = 0.98 µM; SI > 60). Six targets were identified in BraMMT that interact with compounds exhibiting a stronger binding energy than the crystallographic ligand: P. falciparum triophosphate phosphoglycolate complex (1LYX), P. falciparum reductase (2OK8), PfPK7 (2PML), P. falciparum glutaredoxin (4N0Z), PfATP6, and PfHT.

CONCLUSIONS

The physicochemical properties of compound 1 were compatible with the set of criteria established by the Lipinski rule and demonstrated its potential as a drug prototype for antiplasmodial activity.

Prospecting for Cressa cretica to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2

The severe acute respiratory syndrome COVID-19 declared as a global pandemic by the World Health Organization has become the present wellbeing worry to the whole world. There is an emergent need to search for possible medications. Cressa cretica is reported to show antitubercular, antibacterial and expectorant property. In this research, we aim to prospect the COVID-19 main protease crystal structure (Mpro; PDB ID: 6LU7) and the active chemical constituents from Cressa cretica in order to understand the structural basis of their interactions. We examined the binding potential of active constituents of Cressa cretica plant to immensely conserved protein Mpro of SARS-CoV-2 followed by exploration of the vast conformational space of protein-ligand complexes by molecular dynamics (MD) simulations. The results suggest the effectiveness of 3,5-Dicaffeoylquinic acid and Quercetin against standard drug Remdesivir. The active chemical constituents exhibited good docking scores, and interacts with binding site residues of Mpro by forming hydrogen bond and hydrophobic interactions. 3,5-Dicaffeoylquinic acid showed the best affinity towards Mpro receptor which is one of the target enzymes required by SARS CoV-2 virus for replication suggesting it to be a novel research molecule. The potential of the active chemical constituents from Cressa cretica against the SARS-CoV-2 virus has best been highlighted through this study. Therefore, these chemical entities can be further scrutinized and provides direction for further consideration for in-vivo and in-vitro validations for the treatment of covid-19. Communicated by Ramaswamy H. Sarma.

COVID-19: Pathophysiology, transmission, and drug development for therapeutic treatment and vaccination strategies

COVID-19, a dreaded and highly contagious pandemic, is flagrantly known for its rapid prevalence across the world. Till date, none of the treatments are distinctly accessible for this life-threatening disease. Under the prevailing conditions of medical emergency, one creative strategy for the identification of novel and potential antiviral agents gaining momentum in research institutions and progressively being leveraged by pharmaceutical companies is target-based drug repositioning/repurposing. A continuous monitoring and recording of results offer an anticipation that this strategy may help to reveal new medications for viral infections. This review recapitulates the neoteric illation of COVID-19, its genomic dispensation, molecular evolution via phylogenetic assessment, drug targets, the most frequently worldwide used repurposed drugs and their therapeutic applications, and a recent update on vaccine management strategies. The available data from solidarity trials exposed that the treatment with several known drugs, viz. lopinavir-ritonavir, chloroquine, hydroxychloroquine, etc had displayed various antagonistic effects along with no impactful result in diminution of mortality rate. The drugs like remdesivir, favipiravir, and ribavirin proved to be quite safer therapeutic options for treatment against COVID-19. Similarly, dexamethasone, convalescent plasma therapy and oral administration of 2DG are expected to reduce the mortality rate of COVID-19 patients.

Comparative study of the unbinding process of some HTLV-1 protease inhibitors using unbiased molecular dynamics simulations

The HTLV-1 protease is one of the major antiviral targets to overwhelm this virus. Several research groups have developed protease inhibitors, but none has been successful. In this regard, developing new HTLV-1 protease inhibitors to fix the defects in previous inhibitors may overcome the lack of curative treatment for this oncovirus. Thus, we decided to study the unbinding pathways of the most potent (compound 10, PDB ID 4YDF, Ki = 15 nM) and one of the weakest (compound 9, PDB ID 4YDG, Ki = 7900 nM) protease inhibitors, which are very structurally similar. We conducted 12 successful short and long simulations (totaling 14.8 μs) to unbind the compounds from two monoprotonated (mp) forms of protease using the Supervised Molecular Dynamics (SuMD) without applying any biasing force. The results revealed that Asp32 or Asp32′ in the two forms of mp state similarly exert powerful effects on maintaining both potent and weak inhibitors in the binding pocket of HTLV-1 protease. In the potent inhibitor’s unbinding process, His66′ was a great supporter that was absent in the weak inhibitor’s unbinding pathway. In contrast, in the weak inhibitor’s unbinding process, Trp98/Trp98′ by pi-pi stacking interactions were unfavorable for the stability of the inhibitor in the binding site. In our opinion, these results will assist in designing more potent and effective inhibitors for the HTLV-1 protease.

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.

Antiviral Effect of Polyphenolic Substances in Geranium wilfordii Maxim against HSV-2 Infection Using in vitro and in silico Approaches

Background

Herpes simplex virus type 2 (HSV-2) infestation was the most widespread STD (sexually transmitted diseases) among humans and was the leading cause of infectious recurrent genital herpes. Existing therapies against HSV-2 did incompletely restrain the comeback of activated HSV-2 infestation. Geranium wilfordii Maxim had long been used as traditional Chinese medicine for treating the diseases owing to its anti-inflammatory and antiviral effects. Herein, the study was designed to investigate the antiviral activity of G.wilfordii and its potential effect in regulating the host's immune response.

Methods

To identify the stage of infection at which the compounds inhibited HSV-2, we performed virucidal, therapeutic, and prophylactic assays. The antiviral efficacy was evaluated by the analysis of viral components HSV-2 gD and VP16. The antiviral activities of these compounds were also evaluated by phenotypic analysis, such as cell proliferation and apoptosis. Molecular docking studies on candidate compounds were done to indicate binding interactions between the compounds and adopted compound targets.

Results

Quercetin, corilagin, and geraniin inhibited the replication of HSV-2, with geraniin showing greater TI. The obtained IC₅₀ value of quercetin was 204.7 μM and TI (IC₅₀/EC₅₀) was 5.1, whereas the obtained IC₅₀ value of corilagin was 118.0 μg/ml and TI was 4.05. Geraniin exhibited prominent antiviral activity with an IC₅₀ of 212.4 μM and an EC₅₀ of 18.37 μM, resulting in a therapeutic index (TI) of 11.56. Geraniin showed important in vitro virucidal activity through blocking viral attachment. Compared with the virus group, the apoptosis rates in quercetin-, corilagin-, and geraniin-treated groups were significantly decreased (p < 0.001).The expressions at the transcription genes of virus own replication key factors (including HSV-2 gD and VP16) and cytokines (including TBK1) of infected cells treated with quercetin, corilagin, and geraniin were inhibited. The in silico approaches demonstrated a high number of potential strong intermolecular interactions as hydrogen bonds between geraniin, corilagin, and the activity site of HSV-2 gD. Molecular docking studies demonstrated the effects of corilagin by targeting TBK1.

Conclusions

Together, these results highlighted the importance of G.wilfordii treatment in HSV-2 infection and underscored its therapeutic potential. However, additional in vitro and in vivo research was required to validate our findings.

Emodin Ameliorates Intestinal Dysfunction by Maintaining Intestinal Barrier Integrity and Modulating the Microbiota in Septic Mice

Sepsis-induced inflammatory response leads to intestinal damage and secondary bacterial translocation, causing systemic infections and eventually death. Emodin is a natural anthraquinone derivative in many plants with promising bioactivities. However, the effects and mechanisms of emodin on sepsis-induced intestinal dysfunctions have not been well clarified yet. We found that emodin treatment suppressed the inflammatory response in the intestines of septic mice. Intestinal barrier function was also improved by emodin through enhancing ZO-1 and occludin expression, which prevented the secondary translocation of Escherichia coli. By proteome microarray investigation, JNK2 was identified as a direct target of emodin. In vitro study also showed that emodin inhibited LPS-induced inflammatory response in intestinal epithelial cells. Nuclear factors including NF-κB and AP-1 were further identified as downstream effectors of JNK2. Bioinformatic analysis based on 16s rRNA gene sequencing illustrated that emodin treatment significantly increased the alpha- and beta-diversity of gut microbiota in septic mice. Moreover, data according to functional prediction showed that emodin decreased the abundance of potential pathogenic bacteria in gut. Our findings have shown that emodin treatment prevented inflammatory induced barrier dysfunction and decreased the potential pathogenicity of lumen bacteria, reducing the hazard of lumen bacterial translocation during sepsis.

Computer aided identification of potential SARS CoV-2 main protease inhibitors from diterpenoids and biflavonoids of Torreya nucifera leaves

SARS CoV-2 is the causative agent of the pandemic disease COVID-19. There is an urgent need for effective drugs or vaccines which can effectively combat this outbreak. The main protease (Mpro), a key component for the SARS CoV-2 replication, is considered to be one of the important drug targets for developing anti-COVID-19 drugs. This SARS CoV-2 Mpro/cysteine protease has high sequence similarity with the same protease from SARS CoV-1. Previously, it has been shown experimentally that eight diterpenoids and four biflavonoids derived from the leaf of Torreya nucifera show inhibitory effect on the cleavage/catalytic activity of the SARS CoV-1 Mpro. But whether these phytochemicals exhibit any inhibitory effect on SARS CoV-2 Mpro is unclear. To understand this fact, here, we have adopted various in-silico approaches. Diterpenoids and biflavonoids those qualified pharmacological test (hinokiol, amentoflavone, bilobetin and ginkgetin) and two well-known Mpro inhibitors (N3 and lopinavir) were subjected for molecular docking studies. Only three biflavonoids (amentoflavone, bilobetin and ginkgetin) were selected by comparing their binding affinities with N3 and lopinavir. They interacted with two most important catalytic residues of Mpro (His41 and Cys145). Molecular dynamics studies further revealed that these three Mpro-biflavonoid complexes are highly stable and share a similar degree of compactness. Besides, these complexes experience less conformational fluctuations and more expansion than Mpro-N3 and/or Mpro-lopinavir complex. MM-GBSA and H-bond analysis further corroborated these findings. Altogether, our study suggested that these three biflavonoids could possibly inhibit the proteolytic/catalytic activity of SARS CoV-2 Mpro and might be useful for COVID-19 treatment.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.

Structural analysis, virtual screening and molecular simulation to identify potential inhibitors targeting 2'-O-ribose methyltransferase of SARS-CoV-2 coronavirus

SARS-CoV-2, an emerging coronavirus, has spread rapidly around the world, resulting in over ten million cases and more than half a million deaths as of July 1, 2020. Effective treatments and vaccines for SARS-CoV-2 infection do not currently exist. Previous studies demonstrated that nonstructural protein 16 (nsp16) of coronavirus is an S-adenosyl methionine (SAM)-dependent 2'-O-methyltransferase (2'-O-MTase) that has an important role in viral replication and prevents recognition by the host innate immune system. In the present study, we employed structural analysis, virtual screening, and molecular simulation approaches to identify clinically investigated and approved drugs which can act as promising inhibitors against nsp16 2'-O-MTase of SARS-CoV-2. Comparative analysis of primary amino acid sequences and crystal structures of seven human CoVs defined the key residues for nsp16 2-O'-MTase functions. Virtual screening and docking analysis ranked the potential inhibitors of nsp16 from more than 4,500 clinically investigated and approved drugs. Furthermore, molecular dynamics simulations were carried out on eight top candidates, including Hesperidin, Rimegepant, Gs-9667, and Sonedenoson, to calculate various structural parameters and understand the dynamic behavior of the drug-protein complexes. Our studies provided the foundation to further test and repurpose these candidate drugs experimentally and/or clinically for COVID-19 treatment.Communicated by Ramaswamy H. Sarma.

Excavating phytochemicals from plants possessing antiviral activities for identifying SARS-CoV hemagglutinin-esterase inhibitors by diligent computational workflow

Coronaviruses (CoVs) belong to a group of RNA viruses that cause diseases in vertebrates including. Newer and deadlier than SARS CoV-2 are sought to appear in future for which the scientific community must be prepared with the strategies for their control. Spike protein (S-protein) of all the CoVs require angiotensin-converting enzyme2 (ACE2), while CoVs also require hemagglutinin-acetylesterase (HE) glycoprotein receptor to simultaneously interact with O-acetylated sialic acids on host cells, both these interactions enable viral particle to enter host cell leading to its infection. Target inhibition of viral S-protein and HE glycoprotein receptor can lead to a development of therapy against the SARS CoV-2. The proposition is to recognize molecules from the bundle of phytochemicals of medicinal plants known to possess antiviral potentials as a lead that could interact and mask the active site of, HE glycoprotein which would ideally bind to O-acetylated sialic acids on human host cells. Such molecules can be addressed as 'HE glycoprotein blockers'. A library of 110 phytochemicals from Withania somnifera, Asparagus racemosus, Zinziber officinalis, Allium sativum, Curcuma longa and Adhatoda vasica was constructed and was used under present study. In silico analysis was employed with plant-derived phytochemicals. The molecular docking, molecular dynamics simulations over the scale of 1000 ns (1 μs) and ADMET prediction revealed that the Withania somnifera (ashwagandha) and Asparagus racemosus (shatavari) plants possessed various steroidal saponins and alkaloids which could potentially inhibit the COVID-19 virus and even other CoVs targeted HE glycoprotein receptor.

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.

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.

Using mathematical modeling to estimate time-independent cancer chemotherapy efficacy parameters

One of the primary cancer treatment modalities is chemotherapy. Unfortunately, traditional anti-cancer drugs are often not selective and cause damage to healthy cells, leading to serious side effects for patients. For this reason more targeted therapeutics and drug delivery methods are being developed. The effectiveness of new treatments is initially determined via in vitro cell viability assays, which determine the IC 50  of the drug. However, these assays are known to result in estimates of IC 50  that depend on the measurement time, possibly resulting in over- or under-estimation of the IC 50 . Here, we test the possibility of using cell growth curves and fitting of mathematical models to determine the IC 50  as well as the maximum efficacy of a drug ( ε max ). We measured cell growth of MCF-7 and HeLa cells in the presence of different concentrations of doxorubicin and fit the data with a logistic growth model that incorporates the effect of the drug. This method leads to measurement time-independent estimates of IC 50  and ε max , but we find that ε max  is not identifiable. Further refinement of this methodology is needed to produce uniquely identifiable parameter estimates.

Alantolactone: A Natural Plant Extract as a Potential Therapeutic Agent for Cancer

Alantolactone (ALT) is a natural compound extracted from Chinese traditional medicine Inula helenium L. with therapeutic potential in the treatment of various diseases. Recently, in vitro and in vivo studies have indicated cytotoxic effects of ALT on various cancers, including liver cancer, colorectal cancer, breast cancer, etc. The inhibitory effects of ALT depend on several cancer-associated signaling pathways and abnormal regulatory factors in cancer cells. Moreover, emerging studies have reported several promising strategies to enhance the oral bioavailability of ALT, such as combining ALT with other herbs and using ALT-entrapped nanostructured carriers. In this review, studies on the anti-tumor roles of ALT are mainly summarized, and the underlying molecular mechanisms of ALT exerting anticancer effects on cells investigated in animal-based studies are also discussed.

Molecular simulation probes the potency of resveratrol in regulating the toxic aggregation of mutant V30M TTR fibrils in Transthyretin mediated amyloidosis

Transthyretin (TTR) mediated amyloidosis is a highly ruinous illness that affects various organs by aggravating the deposition of misfolded or mutated TTR protein aggregates in tissues. Hence, hindering the formation of TTR amyloid aggregates could be a key strategy in finding an effective cure towards the aggravating disorder. In this analysis, we examined the subversive nature of point mutation, V30M, in TTR that promotes amyloidogenicity using discrete molecular dynamics (DMD) simulations. Besides, we probed the association of naturally occurring polyphenols: EGCG (a proven anti TTR aggregation agent as positive control), resveratrol and curcumin in mitigating the pathogenic repercussions of mutant TTR. Results from the computational studies endorsed that the resveratrol constitutes a restorative potential to subjugate TTR mediated amyloidosis, besides EGCG. Hence, this study could be a reminiscent aspect in understanding the inhibitory role of key polyphenols against the mutant TTR aggregates, which could be an aid towards structure-based drug design in the upcoming research era on familial amyloid disorders.

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.

Immunoinformatics and Immunogenetics-Based Design of Immunogenic Peptides Vaccine against the Emerging Tick-Borne Encephalitis Virus (TBEV) and Its Validation through In Silico Cloning and Immune Simulation

Tick-borne encephalitis virus (TBEV), belonging to the Flaviviridae family, is transmitted to humans via infected tick bites, leading to serious neurological complications and, in some cases, death. The available vaccines against the TBEV are reported to have low immunogenicity and are associated with adverse effects like swelling, redness and fever. Moreover, these vaccines are whole-organism-based, carry a risk of reactivation and potential for significant mortality. Consequently, to design a potential antigenic and non-allergenic multi-epitope subunit vaccine against the TBEV, we used an immunoinformatic approach to screen the Tick-borne virus proteome for highly antigenic CTL, HTL and B cell epitopes. The proper folding of the constructed vaccine was validated by a molecular dynamic simulation. Additionally, the molecular docking and binding free energy (−87.50 kcal/mol) further confirmed the strong binding affinity of the constructed vaccine with TLR-4. The vaccine exhibited a CAI value of 0.93 and a GC content of 49%, showing a high expression capability in E coli. Moreover, the analysis of immune simulation demonstrated robust immune responses against the injected vaccine and clearance of the antigen with time. In conclusion, our vaccine candidate shows promise for both in vitro and in vivo analyses due to its high immunogenicity, non-allergenicity and stable interaction with the human TLR-4 receptor.

Drug repurposing studies targeting SARS-CoV-2: an ensemble docking approach on drug target 3C-like protease (3CLpro)

The COVID-19 pandemic has been responsible for several deaths worldwide. The causative agent behind this disease is the Severe Acute Respiratory Syndrome - novel Coronavirus 2 (SARS-CoV-2). SARS-CoV-2 belongs to the category of RNA viruses. The main protease, responsible for the cleavage of the viral polyprotein is considered as one of the hot targets for treating COVID-19. Earlier reports suggest the use of HIV anti-viral drugs for targeting the main protease of SARS-CoV, which caused SARS in the year 2002-2003. Hence, drug repurposing approach may prove to be useful in targeting the main protease of SARS-CoV-2. The high-resolution crystal structure of the main protease of SARS-CoV-2 (PDB ID: 6LU7) was used as the target. The Food and Drug Administration approved and SWEETLEAD database of drug molecules were screened. The apo form of the main protease was simulated for a cumulative of 150 ns and 10 μs open-source simulation data was used, to obtain conformations for ensemble docking. The representative structures for docking were selected using RMSD-based clustering and Markov State Modeling analysis. This ensemble docking approach for the main protease helped in exploring the conformational variation in the drug-binding site of the main protease leading to the efficient binding of more relevant drug molecules. The drugs obtained as top hits from the ensemble docking possessed anti-bacterial and anti-viral properties. This in silico ensemble docking approach would support the identification of potential candidates for repurposing against COVID-19.Communicated by Ramaswamy H. Sarma.

A network pharmacology analysis on drug-like compounds from Ganoderma lucidum for alleviation of atherosclerosis

Ganoderma lucidum (GL) is known as a potent alleviator against chronic inflammatory disease like atherosclerosis (AS), but its mechanisms against AS have not been unveiled. This research aimed to identify the key compounds(s) and mechanism(s) of GL against AS through network pharmacology. The compounds from GL were identified by gas chromatography-mass spectrum (GC-MS), and SwissADME screened their physicochemical properties. Then, the target(s) associated with the screened compound(s) or AS related targets were identified by public databases, and we selected the overlapping targets using a Venn diagram. The networks between overlapping targets and compounds were visualized, constructed, and analyzed by RStudio. Finally, we performed a molecular docking test (MDT) to explore key target(s), compound(s), on AutoDockVina. A total of 35 compounds in GL were detected via GC-MS, and 34 compounds (accepted by Lipinski's rule) were selected as drug-like compounds (DLCs). A total of 34 compounds were connected to the number of 785 targets, and DisGeNET and Online Mendelian Inheritance in Man (OMIM) identified 2,606 AS-related targets. The final 98 overlapping targets were extracted between the compounds-targets and AS-related targets. On Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, the number of 27 signaling pathways were sorted out, and a hub signaling pathway (MAPK signaling pathway), a core gene (PRKCA), and a key compound (Benzamide, 4-acetyl-N-[2,6-dimethylphenyl]) were selected among the 27 signaling pathways via MDT. Overall, we found that the identified 3 DLCs from GL have potent anti-inflammatory efficacy, improving AS by inactivating the MAPK signaling pathway. PRACTICAL APPLICATIONS: Ganoderma lucidum (GL) has been used as a medicinal or edible mushroom for chronic inflammatory patients: diabetes mellitus and dyslipidemia, especially atherosclerosis (AS). Until now, the majority of mushroom research has been implemented regarding β-glucan derivatives with very hydrophilic physicochemical properties. It implies that β-glucan or its derivatives have poor bioavailability. Hence, we have involved GC-MS in identifying lipophilic compounds from GL, which filtered them in silico to sort drug-like compounds (DLCs). Then, we retrieved targets associated with the DLCs, and identified a key signaling pathway, key targets, and key compounds against AS. In this paper, we utilized bioinformatics and network pharmacology theory to understand the uncovered pharmacological mechanism of GL on AS. To sum things up, our analysis elucidates the relationships between signaling pathways, targets, and compounds in GL. Ultimately, this work provides biochemical evidence to identify the therapeutic effect of GL on AS, and a scientific basis for deciphering the key mechanism on DLCs of GL against AS.

Integrated bioinformatics analysis of differentially expressed genes and immune cell infiltration characteristics in Esophageal Squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a life-threatening thoracic tumor with a poor prognosis. The role of molecular alterations and the immune microenvironment in ESCC development has not been fully elucidated. The present study aimed to elucidate key candidate genes and immune cell infiltration characteristics in ESCC by integrated bioinformatics analysis. Nine gene expression datasets from the Gene Expression Omnibus (GEO) database were analysed to identify robust differentially expressed genes (DEGs) using the robust rank aggregation (RRA) algorithm. Functional enrichment analyses showed that the 152 robust DEGs are involved in multiple processes in the tumor microenvironment (TME). Immune cell infiltration analysis based on the 9 normalized GEO microarray datasets was conducted with the CIBERSORT algorithm. The changes in macrophages between ESCC and normal tissues were particularly obvious. In ESCC tissues, M0 and M1 macrophages were increased dramatically, while M2 macrophages were decreased. A robust DEG-based protein–protein interaction (PPI) network was used for hub gene selection with the CytoHubba plugin in Cytoscape. Nine hub genes (CDA, CXCL1, IGFBP3, MMP3, MMP11, PLAU, SERPINE1, SPP1 and VCAN) had high diagnostic efficiency for ESCC according to receiver operating characteristic (ROC) curve analysis. The expression of all hub genes except MMP3 and PLAU was significantly related to macrophage infiltration. Univariate and multivariate regression analyses showed that a 7-gene signature constructed from the robust DEGs was useful for predicting ESCC prognosis. Our results might facilitate the exploration of potential targeted TME therapies and prognostic evaluation in ESCC.

Promising inhibitors of main protease of novel corona virus to prevent the spread of COVID-19 using docking and molecular dynamics simulation

Coronavirus disease-2019 (COVID-19) is a global health emergency and the matter of serious concern, which has been declared a pandemic by WHO. Till date, no potential medicine/ drug is available to cure the infected persons from SARS-CoV-2. This deadly virus is named as novel 2019-nCoV coronavirus and caused coronavirus disease, that is, COVID-19. The first case of SARS-CoV-2 infection in human was confirmed in the Wuhan city of the China. COVID-19 is an infectious disease and spread from man to man as well as surface to man . In the present work, in silico approach was followed to find potential molecule to control this infection. Authors have screened more than one million molecules available in the ZINC database and taken the best two compounds based on binding energy score. These lead molecules were further studied through docking against the main protease of SARS-CoV-2. Then, molecular dynamics simulations of the main protease with and without screened compounds were performed at room temperature to determine the thermodynamic parameters to understand the inhibition. Further, molecular dynamics simulations at different temperatures were performed to understand the efficiency of the inhibition of the main protease in the presence of the screened compounds. Change in energy for the formation of the complexes between the main protease of novel coronavirus and ZINC20601870 as well ZINC00793735 at room temperature was determined on applying MM-GBSA calculations. Docking and molecular dynamics simulations showed their antiviral potential and may inhibit viral replication experimentally. Communicated by Ramaswamy H. Sarma.

Polymorphisms at site 469 of B-RAF protein associated with skin melanoma may be correlated with dabrafenib resistance: An in silico study

Melanoma is a type of skin cancer. Numerous genes and their proteins are strongly associated with melanoma susceptibility. This study aims to use an in silico method to identify genetic variants in the melanoma susceptibility gene. The COSMIC database was queried for genes and cross-referenced with three environment-gene interaction databases (EGP, SeattleSNPs and CTD) to identify shared genes. The majority of approved skin melanoma drugs were found to act on the protein serine/threonine-protein kinase (B-RAF) encoded by the BRAF gene, which was also present in all three referenced databases. Comprehensive computational analysis was performed to predict deleterious genetic variants associated with skin melanoma, and the nsSNPs G469V and G469E were prioritized based on their predicted deleterious effects. Molecular dynamic simulation analysis of the B-RAF protein mutants G469V and G469E reveals that variations in the amino acid conformation at the drug binding site result in inconsistency in drug interaction. Additionally, this analysis showed that the G469V and G469E mutants have lower binding energy for dabrafenib than the wild type. The population with the highest frequency of each deleterious and pathogenic variant has been determined. The study's findings would support the development of more effective treatment strategies for skin melanoma. Communicated by Ramaswamy H. Sarma.

Universal two-point interaction of mediator KIX with 9aaTAD activation domains

The nine-amino-acid activation domain (9aaTAD) is defined by a short amino acid pattern including two hydrophobic regions (positions p3-4 and p6-7). The KIX domain of mediator transcription CBP interacts with the 9aaTAD domains of transcription factors MLL, E2A, NF-kB, and p53. In this study, we analyzed the 9aaTADs-KIX interactions by nuclear magnetic resonance. The positions of three KIX helixes α1-α2-α3 are influenced by sterically-associated hydrophobic I611, L628, and I660 residues that are exposed to solvent. The positions of two rigid KIX helixes α1 and α2 generate conditions for structural folding in the flexible KIX-L12-G2 regions localized between them. The three KIX I611, L628, and I660 residues interact with two 9aaTAD hydrophobic residues in positions p3 and p4 and together build a hydrophobic core of five residues (5R). Numerous residues in 9aaTAD position p3 and p4 could provide this interaction. Following binding of the 9aaTAD to KIX, the hydrophobic I611, L628, and I660 residues are no longer exposed to solvent and their position changes inside the hydrophobic core together with position of KIX α1-α2-α3 helixes. The new positions of the KIX helixes α1 and α2 allow the KIX-L12-G2 enhanced formation. The second hydrophobic region of the 9aaTAD (positions p6 and p7) provides strong binding with the KIX-L12-G2 region. Similarly, multiple residues in 9aaTAD position p6 and p7 could provide this interaction. In conclusion, both 9aaTAD regions p3, p4 and p6, p7 provide co-operative and highly universal binding to mediator KIX. The hydrophobic core 5R formation allows new positions of the rigid KIX α-helixes and enables the enhanced formation of the KIX-L12-G2 region. This contributes to free energy and is the key for the KIX-9aaTAD binding. Therefore, the 9aaTAD-KIX interactions do not operate under the rigid key-and-lock mechanism what explains the 9aaTAD natural variability.

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.

Identifying structural-functional analogue of GRL0617, the only well-established inhibitor for papain-like protease (PLpro) of SARS-CoV2 from the pool of fungal metabolites using docking and molecular dynamics simulation

Abstract

The non-structural protein (nsp)-3 of SARS-CoV2 coronavirus is sought to be an essential target protein which is also named as papain-like protease (PLpro). This protease cleaves the viral polyprotein, but importantly in human host it also removes ubiquitin-like interferon-stimulated gene 15 protein (ISG15) from interferon responsive factor 3 (IRF3) protein which ultimately downregulates the production of type I interferon leading to weakening of immune response. GRL0617 is the most potent known inhibitor for PLpro that was initially developed for SARS outbreak of 2003. The PLpro of SARS-CoV and CoV2 share 83% sequence identity but interestingly have several identical conserved amino acids that suggests GRL0617 to be an effective inhibitor for PLpro of SARS-CoV2. GRL0617 is a naphthalene-based molecule and interacts with Tyr268 of SARS-CoV2-PLpro (and Tyr269 of SARS-CoV-PLpro). To identify PLpro inhibitors, we prepared a library of secondary metabolites from fungi with aromatic nature and docked them with PLpro of SARS-CoV and SARS-CoV2. We found six hits which interacts with Tyr268 of SARS-CoV2-PLpro (and Tyr269 of SARS-CoV-PLpro). More surprisingly the top hit, Fonsecin, has naphthalene moiety in its structure, which recruits Tyr268 of SARS-CoV2-PLpro (and Tyr269 of SARS-CoV-PLpro) and has binding energy at par with control (GRL0617). Molecular dynamics (MD) simulation showed Fonsecin to interact with Tyr268 of SARS-CoV2-PLpro more efficiently than control (GRL0617) and interacting with a greater number of amino acids in the binding cleft of PLpro.

Graphic abstract

A computational study on active constituents of Habb-ul-aas and Tabasheer as inhibitors of SARS-CoV-2 main protease

A respiratory pandemic known as coronavirus disease-19 (COVID-19) has created havoc since it emerged from Wuhan, China. COVID-19 is caused by a newly emerged SARS coronavirus (SARS-CoV) with increased pathogenicity named severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Due to the lack of understanding of the mechanism of pathogenesis, an effective therapeutic option is unavailable. Epidemics described in Unani ancient literature include nazla-e-wabai and humma-e-wabai, and most of the symptoms of COVID-19 resemble nazla-e-wabai. Hence, in light of Unani literature, the treatment of COVID-19 can be managed with the composites prescribed in Unani medicine for nazla-e-wabai. In this study, a structure-based drug design approach was carried out to check the effectiveness of the pharmacologically active constituents of the Unani composites prescribed to treat nazla-e-wabai against SARS-CoV-2. We performed molecular docking of the active constituents of these composites against the main protease (M^pro), a potential drug target in SARS-CoV-2. Using detailed molecular docking analysis, Habb-ul-aas and Tabasheer were identified as potential inhibitors of SARS-CoV-2 M^pro. The active constituents of both these composites bind to the substrate-binding pocket of SARS-CoV-2 M^pro, forming interactions with key residues of the binding pocket. Molecular dynamics (MD) simulation suggested the binding of active constituents of Habb-ul-aas with SARS-CoV-2 M^pro with a strong affinity as compared to the constituents of Tabasheer. Thus, this study sheds light on the use of these Unani composites in COVID-19 therapeutics.Communicated by Ramaswamy H. Sarma.

Repurposing of anticancer phytochemicals for identifying potential fusion inhibitor for SARS-CoV-2 using molecular docking and molecular dynamics (MD) simulations

The viral particle, SARS-CoV-2 is responsible for causing the epidemic of Coronavirus disease 2019 (COVID-19). To combat this situation, numerous strategies are being thought for either creating its antidote, vaccine, or agents that can prevent its infection. For enabling research on these strategies, several target proteins are identified where, Spike (S) protein is of great potential. S-protein interacts with human angiotensin-converting-enzyme-2 (ACE2) for entering the cell. S-protein is a large protein and a portion of it designated as a receptor-binding domain (RBD) is the key region that interacts with ACE2, following to which the viral membrane fuses with the alveolar membrane to enter the human cell. The hypothesis is to identify molecules from the pool of anticancer phytochemicals as a lead possessing the ability to interact and mask the amino acids of RBD, making them unavailable to form associations with ACE2. Such a molecule is termed as 'fusion inhibitor'. We hypothesized to identify fusion inhibitors from the NPACT library of anticancer phytochemicals. For this, all the molecules from the NPACT were screened using molecular docking, the five top hits (Theaflavin, Ginkgetin, Ursolic acid, Silymarin and Spirosolane) were analyzed for essential Pharmacophore features and their ADMET profiles were studied following to which the best two hits were further analyzed for their interaction with RBD using Molecular Dynamics (MD) simulation. Binding free energy calculations were performed using MM/GBSA, proving these phytochemicals containing anticancer properties to serve as fusion inhibitors.Communicated by Ramaswamy H. Sarma.