Computational investigation on effect of mutations in PCNA resulting in structural perturbations and inhibition of mismatch repair pathway

Engineered β-glycosidase from Hyperthermophilic Sulfolobus solfataricus with Improved Rd-hydrolyzing Activity for Ginsenoside Compound K Production

Hyperthermophilic Sulfolobus solfataricus β-glycosidase (SS-βGly), with higher stability and activity than mesophilic enzymes, has potential for industrial ginsenosides biotransformation. However, its relatively low ginsenoside Rd-hydrolyzing activity limits the production of pharmaceutically active minor ginsenoside compound K (CK). In this study, first, we used molecular docking to predict the key enzyme residues that may hypothetically interact with ginsenoside Rd. Then, based on sequence alignment and alanine scanning mutagenesis approach, key variant sites were identified that might improve the enzyme catalytic efficiency. The enzyme catalytic efficiency (kcat/Km) and substrate affinity (Km) of the N264D variant enzyme for ginsenoside Rd increased by 60% and decreased by 17.9% compared with WT enzyme, respectively, which may be due to a decrease in the binding free energy (∆G) between the variant enzyme and substrate Rd. In addition, Markov state models (MSM) analysis during the whole 1000-ns MD simulations indicated that altering N264 to D made the variant enzyme achieve a more stable SS-βGly conformational state than the wild-type (WT) enzyme and corresponding Rd complex. Under identical conditions, the relative activities and the CK conversion rates of the N264D enzyme were 1.7 and 1.9 folds higher than those of the WT enzyme. This study identified an excellent hyperthermophilic β-glycosidase candidate for industrial biotransformation of ginsenosides.

In-silico identification of deleterious non-synonymous SNPs of TBX1 gene: Functional and structural impact towards 22q11.2DS

The TBX1 gene plays a critical role in the development of 22q11.2 deletion syndrome (22q11.2DS), a complex genetic disorder associated with various phenotypic manifestations. In this study, we performed in-silico analysis to identify potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) within the TBX1 gene and evaluate their functional and structural impact on 22q11.2DS. A comprehensive analysis pipeline involving multiple computational tools was employed to predict the pathogenicity of nsSNPs. This study assessed protein stability and explored potential alterations in protein-protein interactions. The results revealed the rs751339103(C>A), rs780800634(G>A), rs1936727304(T>C), rs1223320618(G>A), rs1248532217(T>C), rs1294927055 (C>T), rs1331240435 (A>G, rs1601289406 (A>C), rs1936726164 (G>A), and rs911796187(G>A) with a high-risk potential for affecting protein function and stability. These nsSNPs were further analyzed for their impact on post-translational modifications and structural characteristics, indicating their potential disruption of molecular pathways associated with TBX1 and its interacting partners. These findings provide a foundation for further experimental studies and elucidation of potential therapeutic targets and personalized treatment approaches for individuals affected by 22q11.2DS.

Analysis of cholesterol-recognition motifs of the plasma membrane Ca2+-ATPase

The plasma membrane Ca2+-ATPase (PMCA) is crucial for the fine tuning of intracellular calcium levels in eukaryotic cells. In this study, we show the presence of CARC sequences in all human and rat PMCA isoforms and we performed further analysis by molecular dynamics simulations. This analysis focuses on PMCA1, containing three CARC motifs, and PMCA4, with four CARC domains. In PMCA1, two CARC motifs reside within transmembrane domains, while the third is situated at the intracellular interface. The simulations depict more stable RMSD values and lower RMSF fluctuations in the presence of cholesterol, emphasizing its potential stabilizing effect. In PMCA4, a distinct dynamic was found. Notably, the total energy differences between simulations with cholesterol and phospholipids are pronounced in PMCA4 compared to PMCA1. RMSD values for PMCA4 indicate a more energetically favorable conformation in the presence of cholesterol, suggesting a robust interaction between CARCs and this lipid in the membranes. Furthermore, RMSF analysis for CARCs in both PMCA isoforms exhibit lower values in the presence of cholesterol compared to POPC alone. The analysis of H-bond occupancy and total energy values strongly suggests the potential interaction of CARCs with cholesterol. Given the crucial role of PMCAs in physiological calcium regulation and their involvement in diverse pathological processes, this study underscores the significance of CARC motifs and their interaction with cholesterol in elucidating PMCA function. These insights into the energetic preferences associated with CARC-cholesterol interactions offer valuable implications for understanding PMCA function in maintaining calcium homeostasis and addressing potential associated pathologies.

Binding patterns of inhibitors to different pockets of kinesin Eg5

The kinesin-5 family member, Eg5, plays very important role in the mitosis. As a mitotic protein, Eg5 is the target of various mitotic inhibitors. There are two targeting pockets in the motor domain of Eg5, which locates in the α2/L5/α3 region and the α4/α6 region respectively. We investigated the interactions between the different inhibitors and the two binding pockets of Eg5 by using all-atom molecular dynamics method. Combined the conformational analysis with the free-energy calculation, the binding patterns of inhibitors to the two binding pockets are shown. The α2/L5/α3 pocket can be divided into 4 regions. The structures and binding conformations of inhibitors in region 1 and 2 are highly conserved. The shape of α4/α6 pocket is alterable. The space of this pocket in ADP-binding state of Eg5 is larger than that in ADP·Pi-binding state due to the limitation of a hydrogen bond formed in the ADP·Pi-binding state. The results of this investigation provide the structural basis of the inhibitor-Eg5 interaction and offer a reference for the Eg5-targeted drug design.

Screening approach by a combination of computational and in vitro experiments: identification of fluvastatin sodium as a potential SIRT2 inhibitor

BACKGROUND

Sirtuins (SIRTs) are NAD+-dependent deacetylases that play various roles in numerous pathophysiological processes, holding promise as therapeutic targets worthy of further investigation. Among them, the SIRT2 subtype is closely associated with tumorigenesis and malignancies. Dysregulation of SIRT2 activation can regulate the expression levels of related genes in cancer cells, leading to tumor occurrence and metastasis.

METHODS

In this study, we used computer simulations to screen for novel SIRT2 inhibitors from the FDA database, based on which 10 compounds with high docking scores and good interactions were selected for in vitro anti-pancreatic cancer metastasis testing and enzyme binding inhibition experiments. The results showed that fluvastatin sodium may possess inhibitory activity against SIRT2. Subsequently, fluvastatin sodium was subjected to molecular docking experiments with various SIRT isoforms, and the combined results from Western blotting experiments indicated its potential as a SIRT2 inhibitor. Next, molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations were performed, revealing the binding mode of fluvastatin sodium at the SIRT2 active site, further validating the stability and interaction of the ligand-protein complex under physiological conditions.

RESULTS

Overall, this study provides a systematic virtual screening workflow for the discovery of SIRT2 activity inhibitors, identifies the potential inhibitory effect of fluvastatin sodium as a lead compound on SIRT2, and opens up a new direction for developing highly active and selectively targeted SIRT2 inhibitors.

Acute and sub-chronic toxicity study of novel polyherbal formulation in non-alcoholic fatty liver using Wistar rats

Aim: This study assessed the acute and sub-chronic toxicity of a novel polyherbal formulation tablet in Wistar rats Materials & methods: Acute toxicity and sub-chronic toxicity was assessed following OECD (Organisation for the Economic Co-operation and Development) guidelines based on 423 and 408. Results & conclusion: No mortality and toxicity showed in rats during acute toxicity. The LD50 of the extract was at 2000 mg/kg. In sub-chronic study, both sex rats were orally administered at 250, 500,1000 and 2000 mg/kg for 90 days and revealed no significant difference (p < 0.05) in hematological and other parameters compared with the control. Histopathology study did not reveal morphological alteration. The No observed adverse effect level of the tablet was observed until 2000 mg/kg.

Molecular dynamics simulation of the effect of temperature on the conformation of ubiquitin protein

CONTENT

Ubiquitin, a ubiquitous small protein found in all living organisms, is crucial for tagging proteins earmarked for degradation and holds pivotal importance in biomedicine. Protein functionality is intricately linked to its structure. To comprehend the impact of diverse temperatures on ubiquitin protein structure, our study delved into the energy landscape, hydrogen bonding, and overall structural stability of ubiquitin protein at varying temperatures. Through meticulous analysis of root mean square deviation and root mean square fluctuation, we validated the robustness of the simulation conditions employed. Within our simulated system, the bonding energy and electrostatic potential energy exhibited linear augmentation, while the van der Waals energy demonstrated a linear decline. Additionally, our findings highlighted that the α-Helix secondary structure of the ubiquitin protein gradually transitions toward helix destabilization under high-temperature conditions. The secondary structure of ubiquitin protein experiences distinct changes under varying temperatures. The outcomes of our molecular simulations offer a theoretical framework that enhances our comprehension of how temperature impacts the structural stability of ubiquitin protein. These insights contribute not only to a deeper understanding of iniquity's behavior but also hold broader implications in the realm of biomedicine and beyond.

METHODS

All the MD simulations were performed using the GROMACS software with GROMOS96 force field and SPC for water. The ubiquitin protein was put in the center of a cubic box with a length of 8 nm, a setting that allowed > 0.8 nm in the minimal distance between the protein surface and the box wall. To remove the possible coordinate collision of the configurations, in the beginning, the steepest descent method was used until the maximum force between atoms was under 100 kJ/mol·nm with a 0.01 nm step size. Minimization was followed by 30 ps of position-restrained MD simulation. The protein was restrained to its initial position, and the solvent was freely equilibrated. The product phase was obtained with the whole system simulated for 10 ns without any restraint using an integral time step of 1 fs with different temperatures. The cutoff for short-range electronic interaction was set to 1.5 nm. The long-range interactions were treated with a particle-mesh Ewald (PME) method with a grid width of 1.2 nm.

Discovery of a novel SHP2 allosteric inhibitor using virtual screening, FMO calculation, and molecular dynamic simulation

CONTEXT

SHP2 is a non-receptor protein tyrosine phosphatase to remove tyrosine phosphorylation. Functionally, SHP2 is an essential bridge to connect numerous oncogenic cell-signaling cascades including RAS-ERK, PI3K-AKT, JAK-STAT, and PD-1/PD-L1 pathways. This study aims to discover novel and potent SHP2 inhibitors using a hierarchical structure-based virtual screening strategy that combines molecular docking and the fragment molecular orbital method (FMO) for calculating binding affinity (referred to as the Dock-FMO protocol). For the SHP2 target, the FMO method prediction has a high correlation between the binding affinity of the protein-ligand interaction and experimental values (R2 = 0.55), demonstrating a significant advantage over the MM/PBSA (R2 = 0.02) and MM/GBSA (R2 = 0.15) methods. Therefore, we employed Dock-FMO virtual screening of ChemDiv database of ∼2,990,000 compounds to identify a novel SHP2 allosteric inhibitor bearing hydroxyimino acetamide scaffold. Experimental validation demonstrated that the new compound (E)-2-(hydroxyimino)-2-phenyl-N-(piperidin-4-ylmethyl)acetamide (7188-0011) effectively inhibited SHP2 in a dose-dependent manner. Molecular dynamics (MD) simulation analysis revealed the binding stability of compound 7188-0011 and the SHP2 protein, along with the key interacting residues in the allosteric binding site. Overall, our work has identified a novel and promising allosteric inhibitor that targets SHP2, providing a new starting point for further optimization to develop more potent inhibitors.

METHODS

All the molecular docking studies were employed to identify potential leads with Maestro v10.1. The protein-ligand binding affinities of potential leads were further predicted by FMO calculations at MP2/6-31G* level using GAMESS v2020 system. MD simulations were carried out with AmberTools18 by applying the FF14SB force field. MD trajectories were analyzed using VMD v1.9.3. MM/GB(PB)SA binding free energy analysis was carried out with the mmpbsa.py tool of AmberTools18. The docking and MD simulation results were visualized through PyMOL v2.5.0.

Discovery of potential TAAR1 agonist targeting neurological and psychiatric disorders: An in silico approach

Trace amine-associated receptor 1 (TAAR1) is a G-protein-coupled receptor which is primarily expressed in the brain. It is activated by trace amines which play a role in regulating neurotransmitters like dopamine, serotonin and norepinephrine. TAAR1 agonists have potential applications in the treatment of neurological and psychiatric disorders, especially schizophrenia. In this study, we have used a structure-based virtual screening approach to identify potential TAAR1 agonist(s). We have modelled the structure of TAAR1 and predicted the binding pocket. Further, molecular docking of a few well-known antipsychotic drugs was carried out with TAAR1 model, which showed key interactions with the binding pocket. From screening a library of 5 million compounds from the Enamine REAL Database using structure-based virtual screening method, we shortlisted 12 compounds which showed good docking score, glide energy and interactions with the key residues. One lead compound (Z31378290) was finally selected. The lead compound showed promising binding affinity and stable interactions with TAAR1 during molecular dynamics simulations and demonstrated better van der Waals and binding energy than the known agonist, ulotaront. Our findings suggest that the lead compound may serve as a potential TAAR1 agonist, offering a promising avenue for the development of new therapies for neurological and psychiatric disorders.

Gallic Acid-Modified Polyethylenimine-Polypropylene Carbonate-Polyethylenimine Nanoparticles: Synthesis, Characterization, and Anti-Periodontitis Evaluation

The aim of the research was to develop novel gallic acid (GA)-modified amphiphilic nanoparticles of polyethylenimine (PEI)-polypropylene carbonate (PPC)-PEI (PEPE) and comprehensively assess its properties as an antiperiodontitis nanoparticle targeting the Toll-like receptor (TLR). The first step is to evaluate the binding potential of GA to the core trigger receptors TLR2 and TLR4/MD2 for periodontitis using molecular docking techniques. Following this, we conducted NMR, transmission electron microscopy, and dynamic light scattering analyses on the synthesized PEPE nanoparticles. As the final step, we investigated the synthetic results and in vitro antiperiodontitis properties of GA-PEPE nanoparticles. The investigation revealed that GA exhibits potential for targeted binding to TLR2 and the TLR4/MD2 complex. Furthermore, we successfully developed 91.19 nm positively charged PEPE nanoparticles. Spectroscopic analysis indicated the successful synthesis of GA-modified PEPE. Additionally, CCK8 results demonstrated that GA modification significantly reduced the biotoxicity of PEPE. The in vitro antiperiodontitis properties assessment illustrated that 6.25 μM of GA-PEPE nanoparticles significantly reduced the expression of pro-inflammatory factors TNF-α, IL-1β, and IL-6. The GA-PEPE nanoparticles, with their targeted TLR binding capabilities, were found to possess excellent biocompatibility and antiperiodontitis properties. GA-PEPE nanoparticles will provide highly innovative input into the development of anti- periodontitis nanoparticles.

Exosomes derived from human umbilical cord mesenchymal stem cells decrease neuroinflammation and facilitate the restoration of nerve function in rats suffering from intracerebral hemorrhage

Exosomes derived from human umbilical cord mesenchymal stem cells (hUCMSC-ex) have become a hopeful substitute for whole-cell therapy due to their minimal immunogenicity and tumorigenicity. The present study aimed to investigate the hypothesis that hUCMSC-ex can alleviate excessive inflammation resulting from intracerebral hemorrhage (ICH) and facilitate the rehabilitation of the nervous system in rats. In vivo, hemorrhagic stroke was induced by injecting collagenase IV into the striatum of rats using stereotactic techniques. hUCMSC-ex were injected via the tail vein at 6 h after ICH model establishment at a dosage of 200 µg. In vitro, astrocytes were pretreated with hUCMSC-ex and then stimulated with hemin (20 μmol/mL) to establish an ICH cell model. The expression of TLR4/NF-κB signaling pathway proteins and inflammatory factors, including TNF-α, IL-1β, and IL-10, was assessed both in vivo and in vitro to investigate the impact of hUCMSC-ex on inflammation. The neurological function of the ICH rats was evaluated using the corner turn test, forelimb placement test, Longa score, and Bederson score on the 1st, 3rd, and 5th day. Additionally, RT-PCR was employed to examine the mRNA expression of TLR4 following hUCMSC-ex treatment. The findings demonstrated that hUCMSC-ex downregulated the protein expression of TLR4, NF-κB/P65, and p-P65, reduced the levels of pro-inflammatory cytokines TNF-α and IL-1β, and increased the expression of the anti-inflammatory cytokine IL-10. Ultimately, the administration of hUCMSC-ex improved the behavioral performance of the ICH rats. However, the results of PT-PCR indicated that hUCMSC-ex did not affect the expression of TLR4 mRNA induced by ICH, suggesting that hUCMSCs-ex may inhibit TLR4 translation rather than transcription, thereby suppressing the TLR4/NF-κB signaling pathway. We can conclude that hUCMSC-ex mitigates hyperinflammation following ICH by inhibiting the TLR4/NF-κB signaling pathway. This study provides preclinical evidence for the potential future application of hUCMSC-ex in the treatment of cerebral injury.

Structural insights into the pathogenicity of point mutations in human acyl-CoA dehydrogenase homotetramers

Acyl-CoA dehydrogenase deficiency (ACAD) is an inherited and potentially fatal disorder with variable clinical symptoms. The relationship between pathogenicity and deleterious point mutations is investigated here in ACAD structures of short (SCAD) and medium-chain (MCAD) types. Structures and dynamic features of native and mutant forms of enzymes models were compared. A total of 2.88 µs molecular dynamics simulations were performed at four different temperatures. Total energy, RMSD, protein ligand interactions and affinity, RMSF measures, secondary structure changes, and important interactions were studied. Mutations in the three main domains of ACADs are pathogenic, while those located at linker turns are not. Mutations affect mostly tetramer formations, secondary structures, and many contacts and interactions. In R206H (MCAD mutant) which is experimentally known to cause a huge turnover decrease, the lack of a single H-bond between substrate and FAD was observed. Secondary structures showed temperature-dependent changes, and SCAD activity was found to be highly correlated to the enzyme helix 3-10 content. Finally, RMSF patterns pointed to one important loop that maintains the substrate close to the active site and is a cause of substrate wobbling upon mutation. Despite similar structure, function, and cellular location, SCAD and MCAD may have different optimum temperatures that are related to the structure taken at that specific temperature. In conclusion, new insight has been provided on the effect of various SCAD and MCAD pathogenic mutations on the structure and dynamical features of the enzymes.

A discovery of potent kaempferol derivatives as multi-target medicines against diabetes as well as bacterial infections: an in silico approach

Flavonoids demonstrate beneficial effects on human health because flavonoids contain important biological properties. Kaempferol is a flavonol, type of flavonoid found in eatable plants and in plants usually employed in ancient drugs (Moringa oleifera, Tilia spp., fern genus spp. and gingko etc.). Some medicinal studies have shown that the use of foods full of kaempferol decreases the risk of many (cancer, vascular) diseases. All the data of 50 kaempferol derivatives were collected from PubChem database. Through Schrödinger software, 3D-QSAR study was performed for 50 compounds by using method of field base. Conformer of kaempferol derivatives was docked against anti-diabetic, anti-microbial co-crystal structures and protein. To monitor the best anti-diabetic and antibacterial agent, particular kaempferol derivatives were downloaded from PubChem database. Virtual screening by molecular docking provided four lead compounds with four different proteins. These hit compounds were found to be potent inhibitor for diabetic enzymes alpha-amylase and DPP IV and had the potential to suppress DNA gyrase and dihydrofolate reductase synthesis. Molecular dynamic simulation of docked complexes evaluates the value of root mean square fluctuation by iMOD server. Kaempferol 3-O-alpha-L-(2, 3-di-Z-p-coumaroyl) rhamnoside (42) compound used as anti-diabetic and kaempferol 3-O-gentiobioside (3) as antibacterial with good results can be used for drug discovery.Communicated by Ramaswamy H. Sarma.

Genes adaptability and NOL6 protein inhibition studies of fabricated flavan-3-ols lead skeleton intended to treat breast carcinoma

Breast cancer invasive 2.3 million women worldly and second prominent factor of cancer-related mortality. Finding a new site-specific and safe small molecule is a current need in this field. With the aid of deep learning Algorithms, we analyzed the published big database from cancer CBioportal to find the best target protein. Further, Multi-omics analysis such as enrichment analysis, scores of molecular, RNA biological function at a cellular level, and protein domain were obtained and matched to find the better hit molecules. The gene analysis output shows nucleolar protein 6 plays a significant responsibility in breast carcinoma and 354 natural and synthetic lead molecules are docked inside the active site. Docking result gave the output hit molecule falavan-3-ols with a binding score of -5.325 (Kcal/mol) and interaction analysis illustrates, 13 active amino acids favoring the binding interaction with functional groups of the hit molecule compared to the standard molecule Abemacilib (-2.857 (Kcal/mol)). Best docked complex of flavan-3-ols and NOL6 protein subjected to dynamic simulation 100 ns to study the stability. The results proved that π-π stacked, carbon‑hydrogen and electrostatic interactions are stable throughout the 100 ns simulation. The overall results conclude the hit molecule flavan-3-ol will be a safe and potent lead molecule to generate and treat breast carcinoma patients.

Lysine-Specific Demethylase 1 Inhibitors: A Comprehensive Review Utilizing Computer-Aided Drug Design Technologies

Lysine-specific demethylase 1 (LSD1/KDM1A) has emerged as a promising therapeutic target for treating various cancers (such as breast cancer, liver cancer, etc.) and other diseases (blood diseases, cardiovascular diseases, etc.), owing to its observed overexpression, thereby presenting significant opportunities in drug development. Since its discovery in 2004, extensive research has been conducted on LSD1 inhibitors, with notable contributions from computational approaches. This review systematically summarizes LSD1 inhibitors investigated through computer-aided drug design (CADD) technologies since 2010, showcasing a diverse range of chemical scaffolds, including phenelzine derivatives, tranylcypromine (abbreviated as TCP or 2-PCPA) derivatives, nitrogen-containing heterocyclic (pyridine, pyrimidine, azole, thieno[3,2-b]pyrrole, indole, quinoline and benzoxazole) derivatives, natural products (including sanguinarine, phenolic compounds and resveratrol derivatives, flavonoids and other natural products) and others (including thiourea compounds, Fenoldopam and Raloxifene, (4-cyanophenyl)glycine derivatives, propargylamine and benzohydrazide derivatives and inhibitors discovered through AI techniques). Computational techniques, such as virtual screening, molecular docking and 3D-QSAR models, have played a pivotal role in elucidating the interactions between these inhibitors and LSD1. Moreover, the integration of cutting-edge technologies such as artificial intelligence holds promise in facilitating the discovery of novel LSD1 inhibitors. The comprehensive insights presented in this review aim to provide valuable information for advancing further research on LSD1 inhibitors.

Discovery of Palbociclib as a potent c-Myc G4 stabilizer for lung cancer treatment using molecular docking, molecular dynamics simulation, and in vitro activity evaluation

Despite significant progress in lung cancer treatment, this disease remains a prevalent and serious global malignancy, leading to high rates of illness and death. Urgent research is needed to discover new or alternative therapies that can improve clinical outcomes for lung cancer patients. In our study, we successfully demonstrated the effectiveness of Palbociclib, a CDK4/6 inhibitor, in suppressing the growth of lung cancer cells. The IC50 values obtained were 11.00 μM and 11.74 μM for H1299 and A549 cells, respectively. Furthermore, our findings indicate that Palbociclib may possess strong c-Myc G4 stabilizing properties by significantly reducing both protein and mRNA expression levels of c-Myc. Additionally, Palbociclib induces apoptosis and causes cell cycle arrest at the G2/M phase in two cells. Through circular dichroism (CD), molecular docking, and molecular dynamics (MD) simulation, we have provided evidence that Palbociclib enhances the structural stability of c-Myc G4 while exhibiting a high binding affinity to its ligand's binding site on c-Myc G4. These results suggest that Palbociclib holds promise as a novel c-Myc G4 stabilizer for treating cancers associated with abnormal c-Myc activity; further optimization and development are warranted.

Identification, characterization, and CADD analysis of Plasmodium DMAP1 reveals it as a potential molecular target for new anti-malarial discovery

Developing drug resistance in the malaria parasite is a reason for apprehension compelling the scientific community to focus on identifying new molecular targets that can be exploited for developing new anti-malarial compounds. Despite the availability of the Plasmodium genome, many protein-coding genes in Plasmodium are still not characterized or very less information is available about their functions. DMAP1 protein is known to be essential for growth and plays an important role in maintaining genomic integrity and transcriptional repression in vertebrate organisms. In this study, we have identified a homolog of DMAP1 in P. falciparum. Our sequence and structural analysis showed that although PfDMAP1 possesses a conserved SANT domain, parasite protein displays significant structural dissimilarities from human homolog at full-length protein level as well as within its SANT domain. PPIN analysis of PfDMAP1 revealed it to be vital for parasite and virtual High-throughput screening of various pharmacophore libraries using BIOVIA platform-identified compounds that pass ADMET profiling and showed specific binding with PfDMAP1. Based on MD simulations and protein-ligand interaction studies two best hits were identified that could be novel potent inhibitors of PfDMAP1 protein.Communicated by Ramaswamy H. Sarma.

In-silico bioprospecting of secondary metabolites from endophytic Streptomyces spp. against Magnaporthe oryzae, a cereal killer fungus

Rice blast disease, caused by Magnaporthe oryzae, is the most devastating cereal killer worldwide. Note that melanin pigment is an essential factor of M. oryzae virulence, thus fungicides interfering with melanin biosynthesizing enzymes would reduce the pathogenicity. Scytalone dehydratase (SDH) is the key target for commercial fungicides, like carpropamid, due to its role in the dehydration reaction of the fungal melanin pathway. However, a single-point mutation (V75M) in SDH elicits resistance to carpropamid. A lack of effective fungicides against this resistant strain expedited the quest for novel bioactive inhibitors. Currently, bacterial endophytes like Streptomyces have been heralded for synthesizing bioactive metabolites to protect plants from phytopathogens. The literature search led to the identification of 21 Streptomyces spp. symbionts of paddy that can suppress M. oryzae growth. An antiSMASH server was used to explore Streptomyces spp. gene clusters and found 4463 putative metabolites. Besides, 745 unique metabolites were subjected to a series of virtual screening techniques. Ideally, this process identified five potential SDH inhibitors. The docking result highlights that the metabolite pseudopyronine A interacted hydrophobically with both Val75 of SDHWT and Met75 of SDHV75M targets. Moreover, pseudopyronine A has a higher binding free energy with SDHWT (- 89.94 kcal/mol) and SDHV75M (- 71.95 kcal/mol). Interestingly, the pyranones scaffold of pseudopyronine A was reported for antifungal activity against phytopathogens. Dynamic behavior confirms that pseudopyronine A has excellent conformational states with both SDHWT and SDHV75M. Altogether, we hope that this study creates a new avenue for the discovery of novel phytopathogen inhibitors from endophytes.

Supplementary Information

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

Brucine D restrains colorectal cancer tumorigenesis and autophagy by downregulating circ_0068464

Bruceine D (BD) from Brucea javanica (L) exerts an antitumor effect in several human cancers. At present, it has not been reported whether BD inhibits the malignancy of colorectal cancer (CRC) cells. Therefore, investigating the role and regulatory mechanisms of BD in CRC is the main thrust of this study. Effect of BD on CRC cell viability, proliferation, apoptosis, invasion, and autophagy was determined by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazolium bromide, 5-ethynyl-2'-deoxyuridine, flow cytometry, transwell invasion, and western blotting assays. Expression changes of has_circ_0068464 (circ_0068464) were detected using real time quantitative polymerase chain reaction. The molecular mechanisms related to circ_0068464 were predicted through online prediction websites Starbase 2.0, circinteractome, and CircBank and validated using dual-luciferase reporter and RNA pull-down assays. The tumorigenic ability of BD and circ_0068464 on CRC was confirmed by xenograft experiments. The results showed that BD lessened CRC cell proliferation, invasion, autophagy, and prompted cell apoptosis. Circ_0068464 was overexpressed in CRC samples and cells. BD led to a significant reduction in circ_0068464 levels in cells of this carcinoma, but circ_0068464 overexpression partially rescued these effects urged by BD. Also, the combination of BD and circ_0068464 silencing decreased xenograft tumor growth compared to BD alone. Importantly, circ_0068464 could regulate ATG5 expression by functioning as a miR-520h molecular sponge. In conclusion, BD might suppress CRC growth by inhibiting the circ_0068464/miR-520h/ATG5 axis, providing a new perspective for the molecular pathogenesis of CRC and preliminarily indicating that BD may be a promising drug for CRC treatment.

Computational insights into NIMA-related kinase 6: unraveling mutational effects on structure and function

The NEK6 (NIMA-related kinase 6) serine/threonine kinase is a pivotal player in a multitude of cellular processes, including the regulation of the cell cycle and the response to DNA damage. Its significance extends to disease pathogenesis, as changes in NEK6 activity have been linked to the development of cancer. Non-synonymous single nucleotide polymorphisms (nsSNPs) in NEK6 have been linked to cancer as they alter the protein's native structure and function. The association between NEK6 activity and cancer development has prompted researchers to explore the effects of genetic variations within the NEK6 gene. Therefore, we utilized advanced computational tools to analyze 155 high-confidence nsSNPs in the NEK6 gene. From this analysis, 21 nsSNPs were identified as potentially harmful, raising concerns about their impact on NEK6 activity and cancer risk. These 21 mutations were then examined for structural alterations, and eight of nsSNPs (I51M, V76A, I134N, Y152D, R171Q, V186G, L237R, and C285S) were found to destabilize the protein. Among the destabilizing mutations screened, a specific mutation, R171Q, stood out due to its conserved nature. To understand its impact on the protein and conformation, all-atom molecular dynamics simulations (MDS) for 100 ns were performed for both Wildtype NEK6 (WT-NEK6) and R171Q. The simulations revealed that the R171Q variant was unstable and led to significant conformational changes in NEK6. This study provides valuable insights into NEK6 dysfunction caused by single amino acid alterations, offering a novel understanding of the molecular mechanisms underlying NEK6-related cancer progression.

In Vitro Antioxidant and In Vivo Hepatoprotective Properties of Wissadula periplocifolia Extract

Wissadula periplocifolia (L.) Thwaites is a traditional medicinal plant belonging to the family Malvaceae, used in folk medicine for inflamed snake bites and bee stings. The current study was designed to investigate the in vitro antioxidant and in vivo anti-inflammatory and hepatoprotective activities of 80% ethanol extract of W. periplocifolia and its different fractions. The crude ethanolic extract (CEE) was then serially fractionated with petroleum ether fraction (PEF), chloroform fraction (CHF), and aqueous fraction (AQF). The antioxidant activity was assessed using 1,1-diphenyl-2-picrylhydrazyl (DPPH) radical assay, anti-inflammatory activity was determined in the xylene-induced ear edema model, and hepatoprotective activity was measured in the paracetamol-induced hepatic injury model. PEF showed a significant scavenging effect with an IC50 value of 33.5 μg/mL, followed by CEE (IC50 = 42.2 μg/mL), CHF (IC50 = 77 μg/mL), and AQF (IC50 = 80 μg/mL), compared to standard butylated hydroxytoluene (IC50 = 14.8 μg/mL). Both doses of CEE (250 and 500 mg/kg) could reduce ear edema by 41.3 and 50%, respectively, compared to standard diclofenac sodium (76.09%). Moreover, CEE significantly reduces the elevated liver enzymes (ALT, AST, and ALP), compared to control. Nevertheless, it elevated blood protein and reduced the blood bilirubin level (p < 0.01), compared to control. Histopathological studies also indicated significant protection of the liver from paracetamol-induced liver damage. In conclusion, W. periplocifolia could be a good source of antioxidant and hepatoprotective phytochemicals; meanwhile, toxicological and pharmacokinetic studies are recommended.

Nearly complete structure of bacteriophage DT57C reveals architecture of head-to-tail interface and lateral tail fibers

The T5 family of viruses are tailed bacteriophages characterized by a long non-contractile tail. The bacteriophage DT57C is closely related to the paradigmal T5 phage, though it recognizes a different receptor (BtuB) and features highly divergent lateral tail fibers (LTF). Considerable portions of T5-like phages remain structurally uncharacterized. Here, we present the structure of DT57C determined by cryo-EM, and an atomic model of the virus, which was further explored using all-atom molecular dynamics simulations. The structure revealed a unique way of LTF attachment assisted by a dodecameric collar protein LtfC, and an unusual composition of the phage neck constructed of three protein rings. The tape measure protein (TMP) is organized within the tail tube in a three-stranded parallel α-helical coiled coil which makes direct contact with the genomic DNA. The presence of the C-terminal fragment of the TMP that remains within the tail tip suggests that the tail tip complex returns to its original state after DNA ejection. Our results provide a complete atomic structure of a T5-like phage, provide insights into the process of DNA ejection as well as a structural basis for the design of engineered phages and future mechanistic studies.

Decoding the deactivation mechanism of R192W mutation of ZAP-70 using molecular dynamics simulations and binding free energy calculations

CONTEXT

ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation.

METHODS

In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

3D based generative PROTAC linker design with reinforcement learning

Proteolysis targeting chimera (PROTAC), has emerged as an effective modality to selectively degrade disease-related proteins by harnessing the ubiquitin-proteasome system. Due to PROTACs' hetero-bifunctional characteristics, in which a linker joins a warhead binding to a protein of interest (POI), conferring specificity and a E3-ligand binding to an E3 ubiquitin ligase, this could trigger the ubiquitination and transportation of POI to the proteasome, followed by degradation. The rational PROTAC linker design is challenging due to its relatively large molecular weight and the complexity of maintaining the binding mode of warhead and E3-ligand in the binding pockets of counterpart. Conventional linker generation method can only generate linkers in either 1D SMILES or 2D graph, without taking into account the information of ternary structures. Here we propose a novel 3D linker generative model PROTAC-INVENT which can not only generate SMILES of PROTAC but also its 3D putative binding conformation coupled with the target protein and the E3 ligase. The model is trained jointly with the RL approach to bias the generation of PROTAC structures toward pre-defined 2D and 3D based properties. Examples were provided to demonstrate the utility of the model for generating reasonable 3D conformation of PROTACs. On the other hand, our results show that the associated workflow for 3D PROTAC conformation generation can also be used as an efficient docking protocol for PROTACs.

In silico identification of potential protein kinase C alpha inhibitors from phytochemicals from IMPPAT database for anticancer therapeutics: a virtual screening approach

Protein Kinase C alpha (PKCα) is a critical signaling molecule that plays a crucial role in various physiological processes, including cell growth, differentiation, and survival. Over the years, there has been a growing interest in targeting PKCα as a promising drug target for the treatment of various diseases, including cancer. Targeting PKCα can, therefore, serve as a potential strategy to prevent cancer progression and enhance the efficacy of conventional anticancer therapies. We conducted a systematic search for promising compounds for their anticancer potential that target PKCα using natural compounds from the IMPPAT database. The initial compounds were screened through various tests, including analysis of their physical and chemical properties, PAINS filter, ADMET analysis, PASS analysis, and specific interaction analysis. We selected those that showed high binding affinity and specificity to PKCα from the screened compounds, and we further analyzed them using molecular dynamics simulations (MDS) and principal component analysis (PCA). Various systematic parameters from the MDS analyses suggested that the protein-ligand complexes were stabilized throughout the simulation trajectories of 100 nanoseconds (ns). Our findings indicated that compounds Nicandrenone and Withaphysalin D bind to PKCα with high stability and affinity, making them potential candidates for further research in cancer therapeutics innovation in clinical contexts.Communicated by Ramaswamy H. Sarma.

Elucidation of Pharmacological Mechanism Underlying the Anti-Alzheimer's Disease Effects of Evodia rutaecarpa and Discovery of Novel Lead Molecules: An In Silico Study

Alzheimer's disease (AD) is a brain disease with a peculiarity of multiformity and an insidious onset. Multiple-target drugs, especially Chinese traditional medicine, have achieved a measure of success in AD treatment. Evodia rutaecarpa (Juss.) Benth. (Wuzhuyu, WZY, i.e., E. rutaecarpa), a traditional Chinese herb, has been identified as an effective drug to cure migraines. To our surprise, our in silico study showed that rather than migraines, Alzheimer's disease was the primary disease to which the E. rutaecarpa active compounds were targeted. Correspondingly, a behavioral experiment showed that E. rutaecarpa extract could improve impairments in learning and memory in AD model mice. However, the mechanism underlying the way that E. rutaecarpa compounds target AD is still not clear. For this purpose, we employed methods of pharmacology networking and molecular docking to explore this mechanism. We found that E. rutaecarpa showed significant AD-targeting characteristics, and alkaloids of E. rutaecarpa played the main role in binding to the key nodes of AD. Our research detected that E. rutaecarpa affects the pathologic development of AD through the serotonergic synapse signaling pathway (SLC6A4), hormones (PTGS2, ESR1, AR), anti-neuroinflammation (SRC, TNF, NOS3), transcription regulation (NR3C1), and molecular chaperones (HSP90AA1), especially in the key nodes of PTGS2, AR, SLCA64, and SRC. Graveoline, 5-methoxy-N, N-dimethyltryptamine, dehydroevodiamine, and goshuyuamide II in E. rutaecarpa show stronger binding affinities to these key proteins than currently known preclinical and clinical drugs, showing a great potential to be developed as lead molecules for treating AD.

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

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

A lesson for the maestro of the replication fork: Targeting the protein-binding interface of proliferating cell nuclear antigen for anticancer therapy

The proliferating cell nuclear antigen (PCNA) has emerged as a promising candidate for the development of novel cancer therapeutics. PCNA is a nononcogenic mediator of DNA replication that regulates a diverse range of cellular functions and pathways through a comprehensive list of protein-protein interactions. The hydrophobic binding pocket on PCNA offers an opportunity for the development of inhibitors to target various types of cancers and modulate protein-protein interactions. In the present study, we explored the binding modes and affinity of molecule I1 (standard molecule) with the previously suggested dimer interface pocket and the hydrophobic pocket present on the frontal side of the PCNA monomer. We also identified potential lead molecules from the library of in-house synthesized 3-methylenisoindolin-1-one based molecules to inhibit the protein-protein interactions of PCNA. Our results were based on robust computational methods, including molecular docking, conventional, steered, and umbrella sampling molecular dynamics simulations. Our results suggested that the standard inhibitor I1 interacts with the hydrophobic pocket of PCNA with a higher affinity than the previously suggested binding site. Also, the proposed molecules showed better or comparable binding free energies as calculated by the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach and further validated by enhanced umbrella sampling simulations. In vitro and in vivo methods could test the computationally suggested molecules for advancement in the drug discovery pipeline.

Do non-thermal effects exist in microwave heating of glucose aqueous solutions? Evidence from molecular dynamics simulations

The existence of microwave non-thermal effects in food processing is still debated. In this study, molecular dynamics (MD) simulations were performed to investigate the conformation, electrostatic profiles, and intramolecular hydrogen bonds (intra-HB) of glucose in aqueous solution under alternating electric fields of microwaves ranging from 0 to 10⁹ V/m at 2.45 GHz. The results showed a field-induced threshold of 10⁹ V/m. At the threshold, alternating microwaves reoriented the flexible moieties and thus enhanced the intra-HB. The conformational transition among gg, gt, and tg conformers at 10⁹ V/m possibly resulted from the uneven electrostatic potential and the increased intra-HB. In practice, the maximum electric field of microwaves is several times weaker than the threshold, verifying the absence of microwave non-thermal effects for glucose molecules in food processing. This study provides a novel strategy to evaluate the potential non-thermal effects of microwaves in food processing and the related underlying food safety issues.

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.

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

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

Synthesis, Molecular Docking and Molecular Dynamics Simulation of 2- Thioxothiazolidin-4-One Derivatives against Gp41

INTRODUCTION

Gp41 and its conserved hydrophobic groove on the N-terminal heptad repeat region are attractive targets in the design of HIV-1 entry inhibitors. Linearly extended molecules have shown potent anti-HIV-1 activity for their effective interactions with the gp41 binding pocket. Rhodanine ring attached to substituted pyrrole or furan rings has been proved a preferred moiety to be inserted inside the molecular structure of the gp41 inhibitors.

OBJECTIVES

Based on the previous findings we are going to describe some rhodanine derivatives in which a substituted imidazole ring is introduced in place of the pyrrole or furan rings. The compounds' flexibility is increased by inserting methylene groups inside the main scaffold.

METHODS

Molecular docking and molecular dynamics simulations approaches were exploited to investigate the chemical interactions and the stability of the designed ligands-gp41 complex. All compounds were synthesized and their chemical structures were elucidated by 1HNMR, 13CNMR, FTIR and Mass spectroscopy. Biological activities of the compounds against HIV-1 and HIV-2 and their cellular toxicities against the T-lymphocyte (MT-4) cell line were determined.

RESULTS

All the designed compounds showed proper and stable chemical interactions with gp41 according to the in silico studies. The results of the biological tests proved none of the compounds active against HIV-1 replication in cell cultures.

CONCLUSION

Since all the studied compounds were potently toxic for the host cell; it was therefore not possible to assess their anti-HIV activities.

Molecular docking and simulation studies of natural compounds of Vitex negundo L. against papain-like protease (PLpro) of SARS CoV-2 (coronavirus) to conquer the pandemic situation in the world

The severe acute respiratory syndrome coronavirus-2 (SARS CoV-2) is β-coronavirus that is responsible for the pandemic coronavirus disease 2019 (COVID-19) all over the world. The rapid spread of the novel SARS CoV-2 worldwide is raising a significant global public health issue with nearly 61.86 million people infected and 1.4 million deaths. To date, no specific drugs are available for the treatment of COVID-19. The inhibition of proteases essential for the proteolytic treatment of viral polyproteins is a conventional therapeutic strategy for conquering viral infections. In the study, molecular docking approach was used to screen potential drug compounds among the phytochemicals of Vitex negundo L. against COVID-19 infection. Molecular docking analysis showed that oleanolic acid forms a stable complex and other phyto-compounds ursolic acid, 3β-acetoxyolean-12-en-27-oic acid and isovitexin of V. negundo natural compounds form a less-stable complex. When compared with the control the synergistic interaction of these compounds shows inhibitory activity against papain-like protease (PL^pro) of SARS CoV-2 (COVID-19). The molecular dynamics (MD) simulation (50 ns) were performed on the complexes of PL^pro and the phyto-compounds viz. oleanolic acid, ursolic acid, 3β-acetoxyolean-12-en-27-oic acid and isovitexin followed by the binding free energy calculations using MM-GBSA and these molecules have stable interactions with PL^pro protein binding site. The MD simulation study provides more insight into the functional properties of the protein-ligand complex and suggests that these molecules can be considered as a potential drug molecule against COVID-19. In this pandemic situation, these herbal compounds provide a rich resource to produce new antivirals against COVID-19.Communicated by Ramaswamy H. Sarma.

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

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

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

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

Himalayan bioactive molecules as potential entry inhibitors for the human immunodeficiency virus

The human immunodeficiency virus interacts with the cluster of differentiation 4 receptors and one of the two chemokine receptors (CCR5 and CXCR4) to gain entry in human cells. Both the co-receptors are essential for viral entry, replication, and are considered critical targets for antiviral drugs. In this study, bioactive molecules from different Himalayan plants were screened considering their potential to bind with the CCR5 and CXCR4 co-receptors. We utilized computational and thermodynamic parameters to validate the binding of the selected biomolecules to the active site of the co-receptors. The molecules Butyl 2-ethylhexyl phthalate and Dactylorhin-A showed a higher binding affinity with CCR5 co-receptor than the standard antagonist Maraviroc. Moreover, Pseudohypericin, Amarogentin, and Dactylorhin-E exhibited stronger interactions with CXCR4 than the co-crystallized inhibitor Isothiourea-1 t. Hence, we suggest that these molecules could be developed as potential inhibitors of the CCR5 and CXCR4 co-receptors. However, this require further in-vitro and in-vivo validation.

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

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

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

Sensing the interactions between carbohydrate-binding agents and N-linked glycans of SARS-CoV-2 spike glycoprotein using molecular docking and simulation studies

A recent surge in finding new candidate vaccines and potential antivirals to tackle atypical pneumonia triggered by the novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) needs new and unexplored approaches in solving this global pandemic. The homotrimeric transmembrane spike (S) glycoprotein of coronaviruses which facilitates virus entry into the host cells is covered with N-linked glycans having oligomannose and complex sugars. These glycans provide a unique opportunity for their targeting via carbohydrate-binding agents (CBAs) which have shown their antiviral potential against coronaviruses and enveloped viruses. However, CBA-ligand interaction is not fully explored in developing novel carbohydrate-binding-based antivirals due to associated unfavorable responses with CBAs. CBAs possess unique carbohydrate-binding specificity, therefore, CBAs like mannose-specific plant lectins/lectin-like mimic Pradimicin-A (PRM-A) can be used for targeting N-linked glycans of S glycoproteins. Here, we report studies on the binding and stability of lectins (NPA, UDA, GRFT, CV-N and wild-type and mutant BanLec) and PRM-A with the S glycoprotein glycans via docking and MD simulation. MM/GBSA calculations were also performed for docked complexes. Interestingly, stable BanLec mutant (H84T) also showed similar docking affinity and interactions as compared to wild-type BanLec, thus, confirming that uncoupling the mitogenic activity did not alter the lectin binding activity of BanLec. The stability of the docked complexes, i.e. PRM-A and lectins with SARS-CoV-2 S glycoprotein showed favorable intermolecular hydrogen-bond formation during the 100 ns MD simulation. Taking these together, our predicted in silico results will be helpful in the design and development of novel CBA-based antivirals for the SARS-CoV-2 neutralization.Communicated by Ramaswamy H. Sarma.

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

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

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

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

Conformational behavior of coat protein in plants and association with coat protein-mediated resistance against TMV

Tobacco mosaic virus (TMV) coat protein (CP) self assembles in viral RNA deprived transgenic plants to form aggregates based on the physical conditions of the environment. Transgenic plants in which these aggregates are developed show resistance toward infection by TMV referred to as CP-MR. This phenomenon has been extensively used to protect transgenic plants against viral diseases. The mutants T42W and E50Q CP confer enhanced CP-MR as compared to the WT CP. The aggregates, when examined, show the presence of helical discs in the case of WT CP; on the other hand, mutants show the presence of highly stable non-helical long rods. These aggregates interfere with the accumulation of MP as well as with the disassembly of TMV in plant cells. Here, we explored an atomic level insight to the process of CP-MR through MD simulations. The subunit-subunit interactions were assessed with the help of MM-PBSA calculations. Moreover, classification of secondary structure elements of the protein also provided unambiguous information about the conformational changes occurring in the two chains, which indicated toward increased flexibility of the mutant protein and seconded the other results of simulations. Our finding indicates the essential structural changes caused by the mutation in CP subunits, which are critically responsible for CP-MR and provides an in silico insight into the effects of these transitions over CP-MR. These results could further be utilized to design TMV-CP-based small peptides that would be able to provide appropriate protection against TMV infection.

Virtual screening, ADMET prediction and dynamics simulation of potential compounds targeting the main protease of SARS-CoV-2

The coronavirus disease-2019 caused by a novel SARS CoV-2 virus has emerged as a global threat. Still, no drugs are available for its treatment. The main protease is the most conserved structure responsible for the posttranslational processing of non-structural polyproteins of this virus. Therefore, it can be the potential target for drug discovery against SARS CoV-2. Twenty-one thousand two hundred and seven chemical compounds used for sequential virtual screening studies including coronavirus screening compounds (Life Chemical database) and antiviral compounds (Asinex database). The Schrodinger suite 2019 employed for high throughput screening, molecular docking and MM-GBSA through the Glide module. Subsequently, 23 compounds were selected in the phase first selection criteria for re-docking with AutoDock and iDock followed by ADMET prediction. The drug-likeness predicted through Lipinski’s rule of five, Veber’s rule and Muegge’s rule. Finally, three ligands were selected for molecular dynamics simulation studies over 150 ns against the main protease of the SARS CoV-2. They showed promising docking scores on Glide, iDock and AutoDock Vina algorithms (ligand F2679-0163: −10.75, −10.29 and −9.2; ligand F6355-0442: −9.38, −8.61 and −7.6; ligand 8250: −9.795, −7.94 and −7.5), respectively. The RMSD parameter remained stable at 2.5 Å for all the three ligands for 150 ns. The high RMSF fluctuations, RoG of around 22 Å and the binding free energy were favorable in each case. The hydrogen bond interactions of 8250, F6355-0442 and F2679-0163 were six, five and three, respectively. These compounds can be further explored for in vitro experimental validation against SARS-CoV-2.

Communicated by Ramaswamy H. Sarma

Microsecond Simulation of the Proteoglycan-like Region of Carbonic Anhydrase IX and Design of Chemical Inhibitors Targeting pH Homeostasis in Cancer Cells

Carbonic anhydrase IX (CAIX) is a membrane-bound enzyme associated with tumor hypoxia and found to be over expressed in various tumor conditions. Targeting CAIX catalytic activity is proven to be efficient modality in modulating pH homeostasis in cancer cells. Proteoglycan-like (PG) region is unique to CAIX and is proposed to serve as an antenna enhancing the export of protons in conjunction with facilitated efflux of lactate ions via monocarboxylate transporters. Moreover, the PG region is also reported to contribute to the assembly and maturation of focal adhesion links during cellular attachment and dispersion on solid supports. Thus, drug targeting of this region shall efficiently modulate pH homeostasis and cell adhesion in cancer cells. As the PG region is intrinsically disordered, the complete crystal structure is not elucidated. Hence, in this study, we intend to sample the conformational landscape of the PG region at microsecond scale simulation in order to sample the most probable conformations that shall be utilized for structure-based drug design. In addition, the sampled conformations were subjected to high-throughput virtual screening against NCI and Maybridge datasets to identify potential hits based on consensus scoring and validation by molecular dynamics simulation. Further, the identified hits were experimentally validated for efficacy by in vitro and direct enzymatic assays. The results reveal 5-(2-aminoethyl)-1,2,3-benzenetriol to be the most promising hit as it showed significant CAIX inhibition at all levels of in silico and experimental validation.