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

Computer-aided drug design approaches for the identification of potent inhibitors targeting elongation factor G of Mycobacterium tuberculosis

Elongation factor G (EF-G) is essential for protein synthesis in Mycobacterium tuberculosis (Mtb), positioning it as a promising target for anti-tubercular drug development. This study employs Structure-Based Drug Design (SBDD) to identify potential small molecule inhibitors that specifically target EF-G. Initially, binding hotspots on EF-G were pinpointed, and the binding modes of various compounds were analyzed. Through protein-protein interaction studies, several promising candidates were validated. Virtual screening and molecular docking techniques were utilized to evaluate the binding affinities and interactions of 20 candidate molecules with Mtb EF-G. Additionally, toxicity profiles of these compounds were assessed using predictive models, which indicated non-carcinogenic properties. To further refine the selection process, Support Vector Machine (SVM) and Random Forest models were applied to predict cell wall permeability. Notably, Asinex (8853) and Asinex (102619) emerged as top candidates, boasting high probability scores for effective permeability. Molecular docking and molecular dynamics (MD) simulations revealed that Asinex (8853), Asinex (102619), and Otava (79226) exhibited strong binding affinities and favorable conformations within the active site of Mtb EF-G. These findings suggest that these compounds have significant potential as inhibitors, warranting further investigation into their efficacy as novel anti-tubercular agents. Overall, this study emphasizes the value of Structure-Based Drug Design in identifying promising therapeutic candidates against tuberculosis by targeting essential bacterial factors like EF-G.

Pharmacoinformatics-based prediction of Checkpoint kinase-1 inhibitors from Momordica charantia Linn. for cancer

Checkpoint kinase 1 (Chk-1), a serine/threonine kinase family protein, is an emerging target in cancer research owing to its crucial role in cell cycle arrest. Therefore, we aimed to predict potential Chk-1 inhibitors from Momordica charantia Linn., using high-throughput molecular docking. We used a graph theoretical network approach to determine the target protein, Chk-1. Among 86 compounds identified from M. charantia L., five molecules such as α-spinasterol (-9.7 kcal × mol-1), stigmasterol (-9.6 kcal × mol-1), stigmasta-7,22,25-trienol (-9.5 kcal × mol-1), campesterol (-9.5 kcal × mol-1), and stigmasta-7,25-dien-3beta-ol (-9.5 kcal × mol-1) and standard drug CCT245737 (-8.3 kcal × mol-1) displayed highest binding affinity with Chk-1. Besides, pharmacokinetic studies have demonstrated the non-toxic and drug-like properties of these compounds. Furthermore, molecular dynamics (MD) simulation studies confirmed the strong intermolecular interactions and stability of the compounds with Chk-1. The estimation of binding free-energy derived from molecular docking was fully recognized by the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) produced from the MD simulation paths. Altogether, these five compounds may serve as effective inhibitors of Chk-1, thereby could be used to develop new medications for cancer treatment.

Pyrrolo[2,3-d]pyrimidines as potential kinase inhibitors in cancer drug discovery: A critical review

Pyrrolo[2,3-d]pyrimidine-based kinase inhibitors have emerged as an important class of targeted therapeutics to combat various types of cancer. The distinctive structural feature of pyrrolopyrimidine ring system offers an adaptable platform for designing potent inhibitors of various kinases, crucial in regulating cellular processes. The deazapurine framework inherent to pyrrolopyrimidines bears a conspicuous resemblance to adenine, the natural ligand ATP. The structural mimicry enhances their appeal as potent inhibitors of key kinases. This review reconnoitres the intricate process of designing and developing pyrrolopyrimidine based derivatives, accentuating their structural diversity and the strategic modifications employed to enhance selectivity, potency, and pharmacokinetic properties. The discussion delves into medicinal chemistry strategies, highlighting successful examples that have been progressed to clinical evaluation. Furthermore, the review highlights the promise of pyrrolopyrimidine scaffolds in revolutionizing targeted cancer therapy and provides a pioneering perspective on future directions.

An update on small molecule compounds targeting synthetic lethality for cancer therapy

Targeting cancer-specific vulnerabilities through synthetic lethality (SL) is an emerging paradigm in precision oncology. A SL strategy based on PARP inhibitors has demonstrated clinical efficacy. Advances in DNA damage response (DDR) uncover novel SL gene pairs. Beyond BRCA-PARP, emerging SL targets like ATR, ATM, DNA-PK, CHK1, WEE1, CDK12, RAD51, and RAD52 show clinical promise. Selective and bioavailable small molecule inhibitors have been developed to induce SL, but optimization for potency, specificity, and drug-like properties remains challenging. This article illuminated recent progress in the field of medicinal chemistry centered on the rational design of agents capable of eliciting SL specifically in neoplastic cells. It is envisioned that innovative strategies harnessing SL for small molecule design may unlock novel prospects for targeted cancer therapeutics going forward.

PIM Kinases as Potential Biomarkers and Therapeutic Targets in Inflammatory Arthritides

The Proviral Integration site for the Moloney murine leukemia virus (PIM)-1 kinase and its family members (PIM-2 and PIM-3) regulate several cellular functions including survival, proliferation, and apoptosis. Recent studies showed their involvement in the pathogenesis of rheumatoid arthritis RA, while no studies are available on psoriatic arthritis (PsA) and axial spondyloarthritis (axSpA). The main objective of this study is to assess the expression of PIM kinases in inflammatory arthritides, their correlation with proinflammatory cytokines, and their variation after treatment with biologic disease-modifying anti-rheumatic drugs or JAK inhibitors. We evaluated PIM-1, -2, and -3 expression at the gene and protein level, respectively, in the peripheral blood mononuclear cells and serum of patients with RA, PsA, axSpA, and healthy individuals (CTR). All the samples showed expression of PIM-1, -2, and -3 kinases both at the gene and protein level. PIM-1 was the most expressed protein, PIM-3 the least. PIM kinase levels differed between controls and disease groups, with reduced PIM-1 protein and increased PIM-3 protein in all disease samples compared to controls. No difference was found in the expression of these molecules between the three different pathologies. PIM levels were not modified after 6 months of therapy. In conclusion, our preliminary data suggest a deregulation of the PIM pathway in inflammatory arthritides. In-depth studies on the role of PIM kinases in this field are warranted.

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.

Molecular modeling and dynamics studies of the synthetic small molecule agonists with GPR17 and P2Y1 receptor

The human Guanine Protein coupled membrane Receptor 17 (hGPR17), an orphan receptor that activates uracil nucleotides and cysteinyl leukotrienes is considered as a crucial target for the neurodegenerative diseases. Yet, the detailed molecular interaction of potential synthetic ligands of GPR17 needs to be characterized. Here, we have studied a comparative analysis on the interaction specificity of GPR17-ligands with hGPR17 and human purinergic G protein-coupled receptor (hP2Y1) receptors. Previously, we have simulated the interaction stability of synthetic ligands such as T0510.3657, AC1MLNKK, and MDL29951 with hGPR17 and hP2Y1 receptor in the lipid environment. In the present work, we have comparatively studied the protein-ligand interaction of hGPR17-T0510.3657 and P2Y1-MRS2500. Sequence analysis and structural superimposition of hGPR17 and hP2Y1 receptor revealed the similarities in the structural arrangement with the local backbone root mean square deviation (RMSD) value of 1.16 Å and global backbone RMSD value of 5.30 Å. The comparative receptor-ligand interaction analysis between hGPR17 and hP2Y1 receptor exposed the distinct binding sites in terms of geometrical properties. Further, the molecular docking of T0510.3657 with the hP2Y1 receptor have shown non-specific interaction. The experimental validation also revealed that Gi-coupled activation of GPR17 by specific ligands leads to the adenylyl cyclase inhibition, while there is no inhibition upon hP2Y1 activation. Overall, the above findings suggest that T0510.3657-GPR17 binding specificity could be further explored for the treatment of numerous neuronal diseases. Communicated by Ramaswamy H. Sarma.

Design, synthesis, and molecular dynamics simulation studies of quinoline derivatives as protease inhibitors against SARS-CoV-2

A new series of quinoline derivatives has been designed and synthesized as probable protease inhibitors (PIs) against severe acute respiratory syndrome coronavirus 2. In silico studies using DS v20.1.0.19295 software have shown that these compounds behaved as PIs while interacting at the allosteric site of target Mpro enzyme (6LU7). The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with His41, His164, Glu166, Tyr54, Asp187, and showed π-interaction with His41, the highly conserved amino acids in the target protein. Toxicity Prediction by Komputer Assisted Technology results confirmed that the compounds were found to be less toxic than the reference drug. Further, molecular dynamics simulations were performed on compound 5 and remdesivir with protease enzyme. Analysis of conformational stability, residue flexibility, compactness, hydrogen bonding, solvent accessible surface area (SASA), and binding free energy revealed comparable stability of protease:5 complex to the protease: remdesivir complex. The result of hydrogen bonding showed a large number of intermolecular hydrogen bonds formed between protein residues (Glu166 and Gln189) and ligand 5, indicating strong interaction, which validated the docking result. Further, compactness analysis, SASA and interactions like hydrogen-bonding demonstrated inhibitory properties of compound 5 similar to the existing reference drug. Thus, the designed compound 5 might act as a potential inhibitor against the protease enzyme.Communicated by Ramaswamy H. SarmaHighlightsQuinoline derivatives have been designed as protease inhibitors against SARS-CoV-2.The compounds were docked at the allosteric site of SARS-CoV-2-Mpro enzyme (PDB ID: 6LU7) to study the stability of protein-ligand complex.Docking studies indicated the stable ligand-protein complexes for all designed compounds.The Toxicity Prediction by Komputer Assisted Technology protocol in DS v20.1.0.19295 software was used to evaluate the toxicity of the designed quinoline derivatives.Molecular dynamics studies indicated the formation of stable ligand-Mpro complexes.

Imaging of prostate cancer: optimizing affinity to prostate specific membrane antigen by spacer modifications in a tumor spheroid model

Early diagnosis of prostate cancer (PCa) is crucial for staging, treatment and management of patients. Prostate specific membrane antigen (PSMA), highly over-expressed on PCa cells, is an excellent target for selective imaging of PCa. In recent years, various scaffolds have been explored as potential carriers to target diagnostic and therapeutic agents to PSMA⁺ tumour cells. Numerous fluorescent or radioisotope probes linked via a peptide linker have been developed that selectively binds to PCa cells. However, there are very few reports that examine the effects of chemical modifications in the peptide linker of an imaging probe on its affinity to PSMA protein. This report systematically investigates the impact of hydrophobic aromatic moieties in the peptide linker on PSMA affinity and in vitro performance. For this, a series of fluorescent bioconjugates 12-17 with different aromatic spacers were designed, synthesized, and their interactions within the PSMA pocket were first analysed in silico. Cell uptake studies were then performed for 12-17 in PSMA⁺ cell lines and 3D tumour models in vitro. Binding affinity values of 12-17 were found to be in the range of 36 to 157.9 nM, and 12 with three aromatic groups in the spacer exhibit highest affinity (K_D = 36 nM) compared to 17 which is devoid of aromatic groups. These studies suggest that aromatic groups in the spacer region can significantly affect deep tissue imaging of fluorescent bioconjugates. Bioconjugate 12 can be a promising diagnostic tool, and conjugation to near-infrared agents would further its applications in deep-tissue imaging and surgery. 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.

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.

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.

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.

Structure-based study to identify alkaloids as promising cytochrome P450 (CYP1A1) inhibitors: An in silico approach using virtual screening, molecular dynamic simulations, and binding free energy calculation

Carcinogens present in smokeless tobacco (SLT) like tobacco-specific nitrosamines can be metabolized by the cytochrome P450 (CYP450) enzyme. Functionally, the CYP450 enzyme resides in a heme pigment to perform the catalytic activity. The CYP1A1 is one of the main extrahepatic CYP450 enzymes known to detoxify toxic substances and activate carcinogens. The CYP1A1 inhibition by potential inhibitors reduce the chance of oral cancer. The current study aimed to explore more about the inhibitor binding site and identification of lead alkaloids, that could work as putative inhibitors against target CYP1A1. In respect, we have performed docking studies, virtual screening of alkaloids, and natural product libraries against CYP1A1 followed by molecular dynamic simulations and binding free energy calculations. Docking studies of tobacco-specific nitrosamine (TSNA) products and their similar carcinogen analogs revealed that the heme group is bound to the floor of the bowl-shaped cavity whereas carcinogens are bound to the roof of the rounded shape cavity. Furthermore, virtual screening and binding free energy calculations revealed Tomatidine as a putative inhibitor against CYP1A1. On the basis of altogether outcomes of the current study, we have concluded that the addition of lead-hit alkaloid Tomatidine and others in SLT products may be working as a supplement that could be able to reduce the expression of human CYP1A1 and suppresses carcinogenic by-products formations.

Integrating Network Pharmacology and Experimental Validation to Explore the Key Mechanism of Gubitong Recipe in the Treatment of Osteoarthritis

Background

Gubitong Recipe (GBT) is a prescription based on the Traditional Chinese Medicine (TCM) theory of tonifying the kidney yang and strengthening the bone. A previous multicentral randomized clinical trial has shown that GBT can effectively relieve joint pain and improve quality of life with a high safety in treating osteoarthritis (OA). This study is aimed at elucidating the active compounds, potential targets, and mechanisms of GBT for treating OA.

Method

The network pharmacology method was used to predict the key active compounds, targets, and mechanisms of GBT in treating OA. An OA rat model was established with Hulth surgery, and the pathological changes of articular cartilage were observed to evaluate the effects of GBT. Chondrocytes were stimulated with LPS to establish in vitro models, and key targets and mechanisms predicted by network pharmacology were verified via qRT-PCR, ELISA, western blot, and immunofluorescence. The Contribution Index Model and molecular docking were used to determine the key active compounds of GBT and the major nodes affecting predicted pathways.

Result

A total of 500 compounds were acquired from related databases, where 87 active compounds and their 254 corresponding targets were identified. 2979 OA-related genes were collected from three databases, 150 of which were GBT-regulating OA genes. The compound-target network weight analysis and PPI results showed that IL-6 and PGE2 are key targets of GBT in treating OA. KEGG results showed that PI3K/AKT, Toll-like receptor, NFκB, TNF, and HIF-1 are the key signaling pathways. An in vivo experiment showed that GBT could effectively suppress cartilage degradation of OA rats. In vitro experiments demonstrated that GBT can inhibit the key targets of KEGG-related pathways. Molecular-docking results suggested that luteolin, licochalcone A, and β-carotene were key targets of GBT, and the mechanisms may be associated with the NFκB signaling pathway. Blockage experiments showed that the NFκB pathway is the key pathway of GBT in treating OA.

Conclusion

This study verified that GBT can effectively protect articular cartilage through multitarget and multipathway, and its inhibitory effect on the NFκB pathway is the most key mechanism in treating OA.

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.

In silico analyses of acetylcholinesterase (AChE) and its genetic variants in interaction with the anti-Alzheimer drug Rivastigmine

Alzheimer's disease (AD) is the leading cause of dementia worldwide. Despite causing great social and economic impact, there is currently no cure for AD. The most effective therapy to manage AD symptoms is based on acetylcholinesterase inhibitors (AChEi), from which rivastigmine presented numerous benefits. However, mutations in AChE, which affect approximately 5% of the population, can modify protein structure and function, changing the individual response to Alzheimer's treatment. In this study, we performed computer simulations of AChE wild type and variants R34Q, P135A, V333E, and H353N, identified by one or more genome-wide association studies, to evaluate their effects on protein structure and interaction with rivastigmine. The functional effects of AChE variants were predicted using eight machine learning algorithms, while the evolutionary conservation of AChE residues was analyzed using the ConSurf server. Autodock4.2.6 was used to predict the binding modes for the hAChE-rivastigmine complex, which is still unknown. Molecular dynamics (MD) simulations were performed in triplicates for the AChE wild type and mutants using the GROMACS packages. Among the analyzed variants, P135A was classified as deleterious by all the functional prediction algorithms, in addition to occurring at highly conserved positions, which may have harmful consequences on protein function. The molecular docking results suggested that rivastigmine interacts with hAChE at the upper active-site gorge, which was further confirmed by MD simulations. Our MD findings also suggested that the complex hAChE-rivastigmine remains stable over time. The essential dynamics revealed flexibility alterations at the active-site gorge upon mutations P135A, V333E, and H353N, which may lead to strong and nonintuitive consequences to hAChE binding. Nonetheless, similar binding affinities were registered in the MMPBSA analysis for the hAChE wild type and variants when complexed to rivastigmine. Finally, our findings indicated that the rivastigmine binding to hAChE is an energetically favorable process mainly driven by negatively charged amino acids.

Molecular Dynamics Simulation and Essential Dynamics of Deleterious Proline 12 Alanine Single-Nucleotide Polymorphism in PPARγ2 Associated with Type 2 Diabetes, Cardiovascular Disease, and Nonalcoholic Fatty Liver Disease

Background. Peroxisome proliferator-activated receptor-γ (PPARγ) gene is located at 3p25 position. PPARγ functions as the master regulator of glucose homeostasis and lipoprotein metabolism, and recent studies have reported that it is involved in various metabolic diseases such as diabetes mellitus, hyperlipidemia, coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). PPARγ1 and PPARγ2 are necessary for the development of adipose tissue and insulin sensitivity regulation. But PPARγ2 is the isoform that was controlled in response to nutrient intake and obesity. Methodology. In this study, we used computational techniques to show the impact of Pro12Ala polymorphism on PPARγ2. The 3-D structure of PPARγ2 was modeled using I-TASSER server. The modeled structure was validated with the ZLab server, and the mutation was created with SPDB viewer. Stability prediction tools were used. Molecular dynamics simulation (MDS) was used to understand the structural and functional behavior of the wild type and mutant. Essential dynamics was also applied. Results and Conclusions. Stability prediction tools were showed that this mutation has a destabilizing effect on the PPARγ2 structure. The RMSD, RMSF, Rg, SASA, and DSSP were in line with H-bond results that showed less flexibility in the mutant structure. Essential dynamics was used to verify MDS results. Pro12Ala polymorphism leads to rigidity of the PPARγ2 protein and might disturb the conformational changes and interactions of PPARγ2 and results in type 2 diabetes mellitus (T2DM), CAD, and NAFLD. This study can help scientists to develop a drug therapy against these diseases.

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.

Interaction analyses of hTAAR1 and mTAAR1 with antagonist EPPTB

Trace amine-associated receptor 1 (TAAR1) plays a critical role in regulating monoaminergic activity. EPPTB is the only known selective potent antagonist of the mouse (m) TAAR1 presently, while it was shown to be weak at antagonizing human (h) TAAR1. The lack of high-resolution structure of TAAR1 hinders the understanding of the differences in the interaction modes between EPPTB and m/hTARR1. The purpose of this study is to probe these interaction modes using homology modeling, molecular docking, molecular dynamics (MD) simulations, and molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations. Eight populated conformers of hTAAR1-EPPTB complex were observed during the MD simulations and could be used in structure-based virtual screening in future. The MM-GBSA binding energy of hTAAR1-EPPTB complex (-96.5 kcal/mol) is larger than that of mTAAR1-EPPTB complex (-106.7 kcal/mol), which is consistent with the experimental finding that EPPTB has weaker binding affinity to hTAAR1. The several residues in binding site of hTAAR1 (F154^4.56, T194^5.42 and I290^7.39) are different from these of mTAAR1 (Y153^4.56, A193^5.42 and Y287^7.39), which might contribute to the binding affinity difference. Our docking analysis on another hTAAR1 antagonist Compound 3 has found that 1). this compound binds in different pockets of our mTAAR1 and hTAAR1 homology models with a slightly stronger binding affinity to hTAAR1; 2). both antagonists bind to a very similar pocket of hTAAR1.

Identification of 11β-HSD1 inhibitors through enhanced sampling methods

Aminoarylbenzosuberene (AAB) molecules were chosen for in silico analysis to develop effective and more competent 11β-hydroxysteroid dehydrogenase (11β-HSD1) protein inhibitors. The AAB4 molecule was shown to have stronger interactions and binding affinity than standard inhibitors (co-crystallized molecules). These results were based on conventional, steered and enhanced umbrella sampling simulations.

The Computer Simulation of Therapy with the NMDA Antagonist in Excitotoxic Neurodegeneration in an Alzheimer's Disease-like Pathology

(1) Background: The use of uncompetitive N-methyl-D-aspartate (NMDA) receptor antagonists results in neuroprotective benefits in patients with moderate to severe Alzheimer’s disease. In this study, we demonstrated mathematical and computer modelling of the excitotoxicity phenomenon and performed virtual memantine therapy. (2) Methods: A computer simulation environment of the N-methyl-D-aspartate receptor combining biological mechanisms of channel activation by means of excessive extracellular glutamic acid concentration in three models of excitotoxicity severity. The simulation model is based on sliding register tables, where each table is associated with corresponding synaptic inputs. Modelling of the increase in extracellular glutamate concentration, through over-stimulation of NMDA receptors and exacerbation of excitotoxicity, is performed by gradually increasing the parameters of phenomenological events by the power function. Pathological models were virtually treated with 3−30 µM doses of memantine compared to controls. (3) Results: The virtual therapy results of memantine at doses of 3−30 µM in the pathological models of excitotoxicity severity show statistically significant neuroprotective benefits in AD patients with moderate severity, 1.25 (95% CI, 1.18−1.32) vs. 1.76 (95% CI, 1.71−1.80) vs. 1.53 (95% CI, 1.48−1.59), (p < 0.001), to severe, 1.32 (95% CI, 1.12−1.53) vs. 1.77 (95% CI, 1.72−1.82) vs. 1.73 (95% CI, 1.68−1.79), (p < 0.001), in the area of effects on memory. A statistically significant benefit of memantine was demonstrated for all neuronal parameters in pathological models. In the mild severity model, a statistically significant increase in frequency was obtained relative to virtual memantine treatment with a dose of 3 µM, which was 23.5 Hz (95% CI, 15.5−28.4) vs. 38.8 Hz (95% CI, 34.0−43.6), (p < 0.0001). In the intermediate excitotoxicity severity model, a statistically significant increase in frequency was obtained relative to virtual memantine therapy with a 3 µM dose of 26.0 Hz (95% CI, 15.7−36.2) vs. 39.0 Hz (95% CI, 34.2−43.8) and a 10 µM dose of 26.0 Hz (95% CI, 15.7−36.2) vs. 30.9 Hz (95% CI, 26.4−35.4), (p < 0.0001). A statistically significant increase in frequency was obtained in the advanced excitotoxicity severity model as in the medium. (4) Conclusions: The NMDA antagonist memantine causes neuroprotective benefits in patients with moderate to severe AD. One of the most important benefits of memantine is the improvement of cognitive function and beneficial effects on memory. On the other hand, memantine provides only symptomatic and temporary support for AD patients. Memantine is prescribed in the US and Europe if a patient has moderate to severe AD. Memantine has also been approved for mild to moderate AD patients. However, its very modest effect provides motivation for further research into new drugs in AD. We are the first to present a mathematical model of the NMDA receptor that allows the simulation of excitotoxicity and virtual memantine therapy.

Isolation and acetylcholinesterase inhibitory activity of asterric acid derivatives produced by Talaromyces aurantiacus FL15, an endophytic fungus from Huperzia serrata

Alzheimer's disease (AD) is a neurodegenerative disease and the fourth leading cause of death after cardiovascular disease, tumors, and stroke. Acetylcholinesterase (AChE) inhibitors, which are based on cholinergic damage, remain the mainstream drugs to alleviate AD-related symptoms. This study aimed to explore novel AChE inhibitors produced by the endophytic fungus FL15 from Huperzia serrata. The fungus was identified as Talaromyces aurantiacus FL15 according to its morphological characteristics and ITS, 18S rDNA, and 28S rDNA sequence analysis. Subsequently, seven natural metabolites were isolated from strain FL15, and identified as asterric acid (1), methyl asterrate (2), ethyl asterrate (3), emodin (4), physcion (5), chrysophanol (6), and sulochrin (7). Compounds 1-3, which possess a diphenyl ether structure, exhibited highly selective and moderate AChE inhibitory activities with IC₅₀ values of 66.7, 23.3, and 20.1 μM, respectively. The molecular docking analysis showed that compounds 1-3 interacted with the active catalytic site and peripheral anionic site of AChE, and the esterification substitution groups at position 8 of asterric acid may contribute to its bioactivity. The asterric acid derivatives showed highly selective and moderate AChE inhibitory activities, probably via interaction with the peripheral anionic site and catalytic site of AChE. To the best of our knowledge, this study was the first report of the AChE inhibitory activity of asterric acid derivatives, which opens new perspectives for the design of more effective derivatives that could serve as a drug carrier for new chemotherapeutic agents to treat AD.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03125-2.

High throughput virtual screening and molecular dynamics simulation analysis of phytomolecules against BfmR of Acinetobacter baumannii: anti-virulent drug development campaign

Acinetobacter baumannii is a notorious multidrug resistant bacterium responsible for several hospital acquired infections assisted by its capacity to develop biofilms. A. baumannii BfmR (RstA), a response regulator from the BfmR/S two-component signal transduction system, is the major controller of A. baumannii biofilm development and formation. As a result, BfmR represents a novel target for anti-biofilm treatment against A. baumannii. The discovery of the high-resolution crystal structure of BfmR provides a good chance for computational screening of its probable inhibitors. Therefore, in this study we aim to search new, less toxic, and natural BfmR inhibitors from 8450 phytomolecules available in the Indian Medicinal Plants, Phytochemistry and Therapeutic (IMPPAT) database by analyzing molecular docking against BfmR (PDB ID: 6BR7). Out of these 8450 phytomolecules 6742 molecules were successfully docked with BfmR with the docking score range -6.305 kcal/mol to +5.120 kcal/mol. Structure based-molecular docking (SB-MD) and ADMET (absorption, distribution, metabolism, excretion, & toxicity) profile examination revealed that Norepinephrine, Australine, Calystegine B3, 7,7 A-Diepialexine, and Alpha-Methylnoradrenaline phytocompounds strongly binds to the active site residues of BfmR. Furthermore, molecular dynamics simulation (MDS) studies for 100 ns and the binding free energy (MM/GBSA) analysis elucidated the binding mechanism of Calystegine B3, 7,7 A-Diepialexine, and Alpha-Methylnoradrenaline to BfmR. In summary, these phytocompounds seems to have the promising molecules against BfmR, and thus necessitates further verification by both in vitro and in vivo experiments. HighlightsBfmR plays a key role in biofilm development and exopolysaccharide (EPS) synthesis in A. baumannii.Computational approach to search for promising BfmR inhibitors from IMPAAT database.The lead phytomolecules such as Calystegine B3, 7,7 A-Diepialexine, and Alpha-Methylnoradrenaline displayed significant binding with BfmR active site.The outcome of BfmR binding phytomolecules has broadened the scope of hit molecules validation.Communicated by Ramaswamy H. Sarma.

Network Pharmacology Prediction and Molecular Docking-Based Strategy to Discover the Potential Pharmacological Mechanism of Wen-Yu-Jin against Pulmonary Fibrosis in a Mouse Model

BACKGROUND

Pulmonary fibrosis (PF) is a devastating lung disease, resulting in gas exchange dysfunction until death. The two drugs approved by the FDA, pirfenidone and nintedanib, have obvious side effects. Wen-yu-jin (WYJ), one of the commonly used herbs in China, can treat respiratory diseases. The potential effects and the underlying mechanism of WYJ against PF are unclear.

PURPOSE

Employing network pharmacology, molecular docking, and in vivo and in vitro experiments to explore the potential effects and underlying mechanisms of WYJ in the treatment of PF.

METHODS

Ultra-high pressure liquid chromatography combined with linear ion trap-orbital tandem mass spectrometry (UHPLC-LTQ-orbital trap) was used to identify compounds of WYJ. We got PF-related targets and WYJ compounds-related targets from public databases and further completed critical targets exploration, network construction, and pathway analysis by network pharmacology. Molecular docking predicted binding activity of WYJ compounds and critical targets. Based on the above results, in vivo and in vitro experiments validated the potential effects and mechanisms of WYJ against PF.

RESULTS

23 major compositions of WYJ were identified based on UHPLC-LTQ-Orbitrap. According to the results of network pharmacology, STAT3, SRC, IL6, MAPK1, AKT1, EGFR, MAPK8, MAPK14, and IL1B are critical therapeutic targets. Molecular docking results showed that most of the compounds have good binding activities with critical targets. The results of in vivo and in vitro experiments showed that WYJ alleviated the process of fibrosis by targeting MAPK and STAT3 pathways.

CONCLUSION

Network pharmacology, molecular docking, and in vivo and in vitro experiments showed the potential effects and mechanisms of WYJ against PF, which provides a theoretical basis for the treatment of WYJ with PF.

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.

Hexaconazole exposure ravages biosynthesis pathway of steroid hormones: revealed by molecular dynamics and interaction

Widespread application of hexaconazole for agriculture purpose poses a threat to human health by disrupting normal endocrine homeostasis. To avoid adverse health effects on human, it is crucial to identify the effects of hexaconazole on key enzymes responsible for steroidal hormone synthesis. In view of this, present study was conducted to investigate the interaction mechanisms of hexaconazole with key enzymes in comparison with their food drug administration (FDA) approved inhibitor by molecular docking and molecular dynamics simulations. Results indicate that hexaconazole contacts with the active site of the key enzymes required for steroidal hormonal synthesis. Results pertaining to root-mean-square deviation, root-mean-square calculation, radius of gyration, hydrogen bonding and solvent accessible surface area exhibited that the interaction pattern and stability of interaction of hexaconazole was similar to enzyme specific inhibitor. In addition, ligand and enzyme complex interaction energy of hexaconazole was almost similar to key enzyme and FDA-approved enzyme specific inhibitor complex. This study offers a molecular level of understanding of hexaconazole with different enzymes required for steroidal hormonal synthesis. Findings of the study clearly suggest that hexaconazole has efficacy to stably interact with various enzyme required to progress the pathway of hormonal synthesis. If incessant exposure of hexaconazole occurs during agricultural work it may lead to ravage hormonal synthesis or potent endocrine disruption. The result of binding energy and complex interaction energy is depicted in the graphical abstract.

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.

Genomic and in silico protein structural analyses provide insights into marine polysaccharide-degrading enzymes in the sponge-derived Pseudoalteromonas sp. PA2MD11

Active heterotrophic metabolism is a critical metabolic role performed by sponge-associated microorganisms, but little is known about their capacity to metabolize marine polysaccharides (MPs). Here, we investigated the genome of the sponge-derived Pseudoalteromonas sp. strain PA2MD11 focusing on its macroalgal carbohydrate-degrading potential. Carbohydrate-active enzymes (CAZymes) for the depolymerization of agar and alginate were found in PA2MD11's genome, including glycoside hydrolases (GHs) and polysaccharide lyases (PLs) belonging to families GH16, GH50 and GH117, and PL6 and PL17, respectively. A gene potentially encoding a sulfatase was also identified, which may play a role in the strain's ability to consume carrageenans. The complete metabolism of agar and alginate by PA2MD11 could also be predicted and was consistent with the results obtained in physiological assays. The polysaccharide utilization locus (PUL) potentially involved in the metabolism of agarose contained mobile genetic elements from other marine Gammaproteobacteria and its unusual larger size might be due to gene duplication events. Homology modelling and structural protein analyses of the agarases, alginate lyases and sulfatase depicted clear conservation of catalytic machinery and protein folding together with suitable industrially-relevant features. Pseudoalteromonas sp. PA2MD11 is therefore a source of potential MP-degrading biocatalysts for biorefinery applications and in the preparation of pharmacologically-active oligosaccharides.

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.

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.

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.

Targeting SARS-CoV-2 nonstructural protein 15 endoribonuclease: an in silico perspective

The newly emerged human coronavirus, SARS-CoV-2, had begun to spread last year and sparked worldwide. In this study, molecular docking is utilized to test some previously approved drugs against the SARS-CoV-2 nonstructural protein 15 (Nsp15). We screened 23 drugs, from which three (saquinavir, valrubicin and aprepitant) show a paramount predicted binding affinity (-9.1, -9.6 and -9.2 kcal/mol, respectively) against SARS-CoV-2 Nsp15. Moreover, saquinavir and aprepitant make nonbonded interactions with Leu201 in the active site cavity of Nsp15, while the drug valrubicin interacts with Arg199 and Leu201. This binding pattern may be effective against the targeted protein, leading to Nsp15 blockage and virus abolition. Additionally, the pharmacological properties of the screened drugs are known since they have been approved against different viruses.

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.

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.

Synthesis of magnetic nanoparticles functionalized with histidine and nickel to immobilize His-tagged enzymes using β-galactosidase as a model

The aim of this study was to synthesize iron magnetic nanoparticles functionalized with histidine and nickel (Fe₃O₄-His-Ni) to be used as support materials for oriented immobilization of His-tagged recombinant enzymes of high molecular weight, using β-galactosidase as a model. The texture, morphology, magnetism, thermal stability, pH and temperature reaction conditions, and the kinetic parameters of the biocatalyst obtained were assessed. In addition, the operational stability of the biocatalyst in the lactose hydrolysis of cheese whey and skim milk by batch processes was also assessed. The load of 600 U_enzyme/g_support showed the highest recovered activity value (~50%). After the immobilization process, the recombinant β-galactosidase (HisGal) showed increased substrate affinity and greater thermal stability (~50×) compared to the free enzyme. The immobilized β-galactosidase was employed in batch processes for lactose hydrolysis of skim milk and cheese whey, resulting in hydrolysis rates higher than 50% after 15 cycles of reuse. The support used was obtained in the present study without modifying chemical agents. The support easily recovered from the reaction medium due to its magnetic characteristics. The iron nanoparticles functionalized with histidine and nickel were efficient in the oriented immobilization of the recombinant β-galactosidase, showing its potential application in other high-molecular-weight enzymes.