Relation between free energy landscapes of proteins and dynamics

Development of multi-epitope mRNA vaccine against Clostridioides difficile using reverse vaccinology and immunoinformatics approaches

Clostridioides difficile (C. difficile), as the major pathogen of diarrhea in healthcare settings, has become increasingly prevalent within community populations, resulting in significant morbidity and mortality. However, the therapeutic options for Clostridioides difficile infection (CDI) remain limited, and as of now, no authorized vaccine is available to combat this disease. Therefore, the development of a novel vaccine against C. difficile is of paramount importance. In our study, the complete proteome sequences of 118 strains of C. difficile were downloaded and analyzed. We found four antigenic proteins that were highly conserved and can be used for epitope identification. We designed two vaccines, WLcd1 and WLcd2, that contain the ideal T-cell and B-cell epitopes, adjuvants, and the pan HLA DR-binding epitope (PADRE) sequences. The biophysical and chemical assessments of these vaccine candidates indicated that they were suitable for immunogenic applications. Molecular docking analyses revealed that WLcd1 bonded with higher affinity to Toll-like receptors (TLRs) than WLcd2. Furthermore, molecular dynamics (MD) simulations, performed using Gmx_MMPBSA v1.56, confirmed the binding stability of WLcd1 with TLR2 and TLR4. The preliminary findings suggested that this multi-epitope vaccine could be a promising candidate for protection against CDI; however, experimental studies are necessary to confirm these predictions.

Insight of the molecular mechanism of inhibitors located at different allosteric sites regulating the activity of wild type and mutant KRAS (G12)

As the important hub of many cellular signaling networks, KRAS (Kirsten rat sarcoma viral oncogene homologue) has been identified as a tumor biomarker. It is the frequently mutated oncogene in human cancers, and KRAS protein activation caused by mutations, such as G12D, has been found in many human tumors tissues. Although, there are two specific allosteric sites (AS1 and AS2) on the KRAS protein that can be used as the targets for inhibitor development, the difference of regulatory mechanisms between two individual allosteric sites still not be reported. Here, using molecular dynamics simulations combined with molecular mechanics generalized born surface area (MM/GBSA) analysis, we found that both of the inhibitors, located at AS1 and AS2, were able to reduce the binding free energy between wild type, mutant KRAS (G12/D/V/S/C) and GTP remarkably, however the effect of inhibitors on the binding free energy between wild type, mutant KRAS and GDP was limited. In addition, the degree of decrease of binding free energy between KRAS and GTP caused by inhibitors at AS2 was significantly greater than that caused by inhibitors at AS1. Further analysis revealed that both inhibitors at AS1 and AS2 were able to regulate the fluctuation of Switch Ⅰ and Switch Ⅱ to expand the pocket of the orthosteric site (GTP binding site), thereby reducing the binding of KRAS to GTP. Noteworthy there was significant differences in the regulatory preferences on Switch Ⅰ and Switch Ⅱ between two type inhibitor. The inhibitor at AS2 mainly regulated Switch Ⅱ to affect the pocket of the orthosteric site, while the inhibitor at AS1 mainly expand the pocket of the orthosteric site by regulating the fluctuation of Switch Ⅰ. Our study compared the differences between two type inhibitors in regulating the KRAS protein activity and revealed the advantages of the AS2 as the small molecule drug target, aiming to provide theoretical guidance for the research of novel KRAS protein inhibitors.

In-silico assessment of novel peptidomimetics inhibitor targeting STAT3 and STAT4 N-terminal domain dimerization: A comprehensive study using molecular docking, molecular dynamics simulation, and binding free energy analysis

Dysregulation in Janus kinase-Signal Transducer and Activation of Transcription (JAK-STAT) pathway is closely linked to various cancer types. The N-terminal domain (NTD) of STAT proteins, upon dimerization, assumes a multifaceted role with remarkable adaptability in mediating interactions between proteins. Consequently, the strategic targeting of the N-terminal domain of STATs has emerged as a promising tactic for disrupting dimerization and impeding the translocation of STAT proteins. In this study, we have deployed an integrated in-silico methodology to rationally design Peptidomimetic foldamers as inhibitors of the N-terminal domains of STAT3 and STAT4, with the objective of disrupting protein dimerization. Consequently, we have judiciously designed a series of peptidomimetics that encompass β3-amino acids, bearing side chains that mimic the residues within interface II of the dimeric structures of the NTDs. Employing molecular docking techniques; we have assessed the binding affinity of these designed peptidomimetics toward both the NTDs. Furthermore, we have conducted an evaluation of the stability and conformational alterations within the docked complexes over an extensive Molecular Dynamics, subsequently computing the binding free energy utilizing MM/PBSA calculations. Our findings unequivocally demonstrate that the peptidomimetic foldamers we have devised (Peptide-A, Peptide-B, and Peptide-C) exhibit a propensity to bind to and impede the dimerization process of the NTDs of both STAT3 and STAT4. These outcomes serve to underscore the potential of these meticulously designed peptidomimetics as potential candidates meriting further exploration in the realm of cancer prevention and management.

Marine fungal diversity unlocks potent antivirals against monkeypox through methyltransferase inhibition revealed by molecular dynamics and free energy landscape

The escalating threat posed by the Monkeypox virus (MPXV) to global health necessitates the urgent discovery of effective antiviral agents, as there are currently no specific drugs available for its treatment, and existing inhibitors are hindered by toxicity and poor pharmacokinetic profiles. This study aimed to identify potent MPXV inhibitors by screening a diverse library of small molecule compounds derived from marine fungi, focusing on the viral protein VP39, a key methyltransferase involved in viral replication. An extensive virtual screening process identified four promising compounds-CMNPD15724, CMNPD28811, CMNPD30883, and CMNPD18569-alongside a control molecule. Rigorous evaluations, including re-docking, molecular dynamics (MD) simulations, and hydrogen bond analysis, were conducted to assess their inhibitory potential against MPXV VP39. CMNPD15724 and CMNPD30883, in particular, demonstrated a superior binding affinity and stable interactions within the target protein's active site throughout the MD simulations, suggesting a capacity to overcome the limitations associated with sinefungin. The stability of these VP39-compound complexes, corroborated by MD simulations, provided crucial insights into the dynamic behavior of these interactions. Furthermore, Principal Component Analysis (PCA) based free energy landscape assessments offered a detailed understanding of the dynamic conformational changes and energetic profiles underlying these compounds' functional disruption of VP39. These findings establish CMNPD15724, CMNPD28811, CMNPD30883, and CMNPD18569 as promising MPXV inhibitors and highlight marine fungi as a valuable source of novel antiviral agents. These compounds represent potential candidates for further experimental validation, advancing the development of safer and more effective therapeutic options to combat this emerging viral infection.

Virtual screening and molecular dynamics simulation of natural compounds as potential inhibitors of serine/threonine kinase 16 for anticancer drug discovery

Serine/threonine kinase 16 (STK 16) is involved in many facets of cellular regulation; activation of STK 16 plays a crucial role in the migration of cancer cells. Therefore, it is a novel target for the discovery of anticancer agents. Herein, virtual screening and dynamics simulation were used to screen a large library of natural compounds against STK 16 using Schrodinger suit 2021-2 and GROMACS 2021.6. The results predicted five molecules with high binding affinity against the target, with NPC132329 (Arcyriaflavin C) and NPC160898 having higher binding affinity and molecular mechanics generalized born surface area (MM/GBSA), suggesting that it is better than the standard inhibitor. The molecular dymanics (MD) simulation studies showed that the STK 16-NPC132329 complex has the lowest root mean square deviation, and STK 16-NPC160898 was the most stable compared with the standard drug and selective STK 16 inhibitor. The minimal fluctuation was observed in the STK 16-NPC132329 and STK 16-NPC160898 complexes based on the root mean square fluctuation trajectory with NPC132329 and NPC160898 forming 2 and 3 hydrogen bonds respectively with the amino acid residue of the target's binding site. Overall, NPC132329 and NPC160898 are better STK 16 inhibitors than the standard drug and selective inhibitor, which can be further studied to discover novel anticancer drugs.

In vitro and in silico investigation of effects of antimicrobial peptides from Solanaceae plants against rice sheath blight pathogen Rhizoctinia solani

Rhizoctonia solani, the causative agent of sheath blight disease in rice, poses a significant threat to agricultural productivity. Traditional management approaches involving chemical fungicides have been effective but come with detrimental consequences for the ecosystem. This study aimed to investigate sustainable alternatives in the form of antifungal peptides derived from Solanaceous plant species as potential agents against R. solani. Peptide extracts were obtained using an optimized antimicrobial peptide (AMP) extraction method and desalted using the solid-phase extraction technique. The antifungal potential of peptide-rich extracts from Solanum tuberosum and Capsicum annum was assessed through in vitro tests employing the agar well diffusion method. Furthermore, peptide-protein docking analysis was performed on HPEPDOCK and HDOCK server; and molecular dynamics simulations (MDS) of 100 ns period were performed using the Gromacs 2020.4. The results demonstrated significant inhibition zones for both extracts at concentrations of 100 mg/mL. Additionally, the extracts of Solanum tuberosum and Capsicum annum had minimum inhibitory concentrations of 50 mg/mL and 25 mg/mL, respectively with minimum fungicidal concentrations of 25 mg/mL. Insights into the potential mechanisms of key peptides inhibiting R. solani targets were gleaned from in-silico studies. Notably, certain AMPs exhibited favorable free energy of binding against pathogenicity-related targets, including histone demethylase, sortin nexin, and squalene synthase, in protein-peptide docking simulations. Extended molecular dynamics simulations lasting 100 ns and MM-PBSA calculations were performed on select protein-peptide complexes. AMP10 displayed the most favorable binding free energy against all target proteins, with AMP3, AMP12b, AMP6, and AMP15 also exhibiting promising results against specific targets of R. solani. These findings underscore the potential of peptide extracts from S. tuberosum and C. annum as effective antifungal agents against rice sheath blight caused by R. solani.

Impact of Phosphorylation at Various Sites on the Active Pocket of Human Ferrochelatase: Insights from Molecular Dynamics Simulations

Ferrochelatase (FECH) is the terminal enzyme in human heme biosynthesis, catalyzing the insertion of ferrous iron into protoporphyrin IX (PPIX) to form protoheme IX (Heme). Phosphorylation increases the activity of FECH, and it has been confirmed that the activity of FECH phosphorylated at T116 increases. However, it remains unclear whether the T116 site and other potential phosphorylation modification sites collaboratively regulate the activity of FECH. In this study, we identified a new phosphorylation site, T218, and explored the allosteric effects of unphosphorylated (UP), PT116, PT218, and PT116 + PT218 states on FECH in the presence and absence of substrates (PPIX and Heme) using molecular dynamics (MD) simulations. Binding free energies were evaluated with the MM/PBSA method. Our findings indicate that the PT116 + PT218 state exhibits the lowest binding free energy with PPIX, suggesting the strongest binding affinity. Additionally, this state showed a higher binding free energy with Heme compared to UP, which facilitates Heme release. Moreover, employing multiple analysis methods, including free energy landscape (FEL), principal component analysis (PCA), dynamic cross-correlation matrix (DCCM), and hydrogen bond interaction analysis, we demonstrated that phosphorylation significantly affects the dynamic behavior and binding patterns of substrates to FECH. Insights from this study provide valuable theoretical guidance for treating conditions related to disrupted heme metabolism, such as various porphyrias and iron-related disorders.

A computational odyssey: uncovering classical β-lactamase inhibitors in dry fruits

In the antibacterial arsenal, β-lactams have held a prominent position, but increasing resistance due to unauthorized use and genetic factors requires new strategies. Combining β-lactamase inhibitors with broad-spectrum β-lactams proves effective in combating this resistance. ESBL producers demand new inhibitors, leading to the exploration of plant-derived secondary metabolites for potent β-lactam antibiotics or alternative inhibitors. Using virtual screening, molecular docking, ADMET analysis, and molecular dynamic simulation, this study actively analyzed the inhibitory activity of figs, cashews, walnuts, and peanuts against SHV-1, NDM-1, KPC-2, and OXA-48 β-lactamases. Using AutoDock Vina, the docking affinities of various compounds for target enzymes were initially screened, revealing 12 bioactive compounds with higher affinities for the target enzymes compared to Avibactam and Tazobactam. Top-scoring metabolites, including Oleanolic acid, Protocatechuic acid, and Tannin, were subjected to MD simulation studies to further analyze the stability of the docked complexes using WebGro. The simulation coordinates, in terms of RMSD, RMSF, SASA, Rg, and hydrogen bonds formed, showed that these phytocompounds are stable enough to retain in the active sites at various orientations. The PCA and FEL analysis also showed the stability of the dynamic motion of Cα residues of phytochemical-bound enzymes. The pharmacokinetic analysis of the top phytochemicals was performed to analyze their bioavailability and toxicity. This study provides new insights into the therapeutic potential of phytochemicals of selected dry fruits and contributes to future experimental studies to identify βL inhibitors from plants.Communicated by Ramaswamy H. Sarma.

Structural Dynamics Analysis of USP14 Activation by AKT-Mediated Phosphorylation

Ubiquitin-specific protease 14 (USP14), one of the three major proteasome-associated deubiquitinating enzymes (DUBs), is known to be activated by the AKT-mediated phosphorylation at Ser432. Thereby, AKT can regulate global protein degradation by controlling the ubiquitin-proteasome system (UPS). However, the exact molecular mechanism of USP14 activation by AKT phosphorylation at the atomic level remains unknown. By performing the molecular dynamics (MD) simulation of the USP14 catalytic domain at three different states (inactive, active, and USP14-ubiquitin complex), we characterized the change in structural dynamics by phosphorylation. We observed that the Ser432 phosphorylation induced substantial conformational changes of USP14 in the blocking loop (BL) region to fold it from an open loop into a β-sheet, which is critical for USP14 activation. Furthermore, phosphorylation also increased the frequency of critical hydrogen bonding and salt bridge interactions between USP14 and ubiquitin, which is essential for DUB activity. Structural dynamics insights from this study pinpoint the important local conformational landscape of USP14 by the phosphorylation event, which would be critical for understanding USP14-mediated proteasome regulation and designing future therapeutics.

Investigating the role of functional mutations in leucine binding to Sestrin2 in aging and age-associated degenerative pathologies using structural and molecular simulation approaches

Leucine is the native known ligand of Sestrin2 (Sesn2) and its interaction with Sesn2 is particularly important, as it influences the activity of mTOR in aging and its associated pathologies. It is important to find out how leucine interacts with Sesn2 and how mutations in the binding pocket of leucine affect the binding of leucine. Therefore, this study was committed to investigating the impact of non-synonymous mutations by incorporating a broad spectrum of simulation techniques, from molecular dynamics to free energy calculations. Our study was designed to model the atomic-scale interactions between leucine and mutant forms of Sesn2. Our results demonstrated that the interaction paradigm for the mutants has been altered thus showing a significant decline in the hydrogen bonding network. Moreover, these mutations compromised the dynamic stability by altering the conformational flexibility, sampling time, and leucine-induced structural constraints that consequently caused variation in the binding and structural stability. Molecular dynamics-based flexibility analysis revealed that the regions 217-339 and 371-380 demonstrated a higher fluctuation. Noteworthy, these regions correspond to a linker (217-339) and a loop (371-380) that cover the leucine binding cavity that is critical for the 'latch' mechanism in the N-terminal, which is essential for leucine binding. Further validation of reduced binding and modified internal motions caused by the mutants was obtained through binding free energy calculations, principal components analysis (PCA), and free energy landscape (FEL) analysis. By unraveling the molecular intricacies of Sesn2-leucine interactions and their mutations, we hope to pave the way for innovative strategies to combat the inevitable tide of aging and its associated diseases.Communicated by Ramaswamy H. Sarma.

Deciphering the Impact of Rare Missense Variants in EGFR-TKI-Resistant Non-Small-Cell Lung Cancer through Whole Exome Sequencing: A Computational Approach

Targeted therapy revolutionizes the treatment of non-small-cell lung cancer (NSCLC), harboring molecular change. Epidermal growth factor receptor(EGFR) mutations play a crucial role in the development of NSCLC, serving as a pivotal factor in its pathogenesis. We elucidated the mechanisms of resistance and potential therapeutic strategies in NSCLC resistant to the EGFR-tyrosine kinase inhibitor (EGFR-TKI). This is achieved by identifying rare missense variants through whole exome sequencing (WES). The goal is to enhance our understanding, identify biomarkers, and lay the groundwork for targeted interventions, thereby offering hope for an improved NSCLC treatment landscape. We conducted WES analysis on 16 NSCLC samples with EGFR-TKI-resistant NSCLC obtained from SRA-NCBI (PRJEB50602) to reveal genomic profiles within the EGFR-TKI. Our findings showed that 48% of the variants were missense, and after filtering with the Ensembl variant effect predictor, 53 rare missense variants in 23 genes were identified as highly deleterious. Further examination using pathogenic tools like PredictSNP revealed 12 deleterious rare missense variants in 7 genes: ZNF717, PSPH, ESRRA, SEMA3G, PTPN7, CAVIN4, and MYBBP1A. Molecular dynamics simulation (MDS) suggested that the L385P variant alters the structural flexibility of ESRRA, potentially leading to unfolding of ERRα proteins. This could impact their function and alter ERRα expression. These insights from MDS enhance our understanding of the structural and dynamic consequences of the L385P ESRRA variant and provide valuable implications for subsequent therapeutic considerations and targeted interventions.

Antifungal plant flavonoids identified in silico with potential to control rice blast disease caused by Magnaporthe oryzae

Rice blast disease, caused by the fungus Magnaporthe oryzae, poses a severe threat to rice production, particularly in Asia where rice is a staple food. Concerns over fungicide resistance and environmental impact have sparked interest in exploring natural fungicides as potential alternatives. This study aimed to identify highly potent natural fungicides against M. oryzae to combat rice blast disease, using advanced molecular dynamics techniques. Four key proteins (CATALASE PEROXIDASES 2, HYBRID PKS-NRPS SYNTHETASE TAS1, MANGANESE LIPOXYGENASE, and PRE-MRNA-SPLICING FACTOR CEF1) involved in M. oryzae's infection process were identified. A list of 30 plant metabolites with documented antifungal properties was compiled for evaluation as potential fungicides. Molecular docking studies revealed that 2-Coumaroylquinic acid, Myricetin, Rosmarinic Acid, and Quercetin exhibited superior binding affinities compared to reference fungicides (Azoxystrobin and Tricyclazole). High throughput molecular dynamics simulations were performed, analyzing parameters like RMSD, RMSF, Rg, SASA, hydrogen bonds, contact analysis, Gibbs free energy, and cluster analysis. The results revealed stable interactions between the selected metabolites and the target proteins, involving important hydrogen bonds and contacts. The SwissADME server analysis indicated that the metabolites possess fungicide properties, making them effective and safe fungicides with low toxicity to the environment and living beings. Additionally, bioactivity assays confirmed their biological activity as nuclear receptor ligands and enzyme inhibitors. Overall, this study offers valuable insights into potential natural fungicides for combating rice blast disease, with 2-Coumaroylquinic acid, Myricetin, Rosmarinic Acid, and Quercetin standing out as promising and environmentally friendly alternatives to conventional fungicides. These findings have significant implications for developing crop protection strategies and enhancing global food security, particularly in rice-dependent regions.

Characterization of Plant-Derived Natural Inhibitors of Dipeptidyl Peptidase-4 as Potential Antidiabetic Agents: A Computational Study

Diabetes, characterized by elevated blood sugar levels, poses significant health and economic risks, correlating with complications like cardiovascular disease, kidney failure, and blindness. Dipeptidyl peptidase-4 (DPP-4), also referred to as T-cell activation antigen CD26 (EC 3.4.14.5.), plays a crucial role in glucose metabolism and immune function. Inhibiting DPP-4 was anticipated as a potential new therapy for diabetes. Therefore, identification of plant-based natural inhibitors of DPP-4 would help in eradicating diabetes worldwide. Here, for the identification of the potential natural inhibitors of DPP-4, we developed a phytochemicals library consisting of over 6000 phytochemicals detected in 81 medicinal plants that exhibited anti-diabetic potency. The library has been docked against the target proteins, where isorhamnetin, Benzyl 5-Amino-5-deoxy-2,3-O-isopropyl-alpha-D-mannofuranoside (DTXSID90724586), and 5-Oxo-7-[4-(trifluoromethyl) phenyl]-4H,6H,7H-[1,2]thiazolo[4,5-b]pyridine 3-carboxylic acid (CHEMBL3446108) showed binding affinities of -8.5, -8.3, and -8.3 kcal/mol, respectively. These compounds exhibiting strong interactions with DPP-4 active sites (Glu205, Glu206, Tyr547, Trp629, Ser630, Tyr662, His740) were identified. ADME/T and bioactivity predictions affirmed their pharmacological safety. Density functional theory calculations assessed stability and reactivity, while molecular dynamics simulations demonstrated persistent stability. Analyzing parameters like RMSD, RG, RMSF, SASA, H-bonds, MM-PBSA, and FEL confirmed stable protein-ligand compound formation. Principal component analysis provided structural variation insights. Our findings suggest that those compounds might be possible candidates for developing novel inhibitors targeting DPP-4 for treating diabetes.

Identification of β-cycloidal-derived mono-carbonyl curcumin analogs as potential interleukin-6 inhibitor to treat wound healing through QSAR, molecular docking, MD simulation, MM-GBSA calculation

Interleukin-6 (IL-6) is a cytokine that involved in the different phases of wound healing. It is responsible for promoting inflammation, regulating tissue repair scar formation, stimulating the production of extracellular matrix components and recruiting immune cells to the wound site. Therefore, suppressing IL-6 is beneficial for wound healing. However, no small molecules are currently available in the market against the IL-6. As a result, this research gap motivates us to find a potential inhibitor. This study aimed to investigate the wound healing potential of novel β-cycloidal-derived mono-carbonyl curcumin analogs reported in the literature through screening a series of computational studies. The calculated pIC50 value of 18 compounds (below 10) showed that all compounds may have potential therapeutic efficacy. Molecular docking studies revealed that compound C12 (-45.6044 kcal/mol) bound most strongly in the active site of IL-6 compared to the FDA-approved drug clindamycin (-42.3223). The Molecular Dynamic (MD) simulation displayed that lead compound C12 had the highest stability in the active site of IL-6 compared to the reference drug clindamycin. Furthermore, MMGBSA results indicated that C12 (-20.28 kcal/mol) had the highest binding energy compared to clindamycin (-8.36 kcal/mol). The ADMET analysis predicted that C12 are favourable for drug candidates. This study recommended compound C12 as a lead IL-6 inhibitor for future testing and development as therapeutics for wound healing.Communicated by Ramaswamy H. Sarma.

De novo antioxidant peptide design via machine learning and DFT studies

Antioxidant peptides (AOPs) are highly valued in food and pharmaceutical industries due to their significant role in human function. This study introduces a novel approach to identifying robust AOPs using a deep generative model based on sequence representation. Through filtration with a deep-learning classification model and subsequent clustering via the Butina cluster algorithm, twelve peptides (GP1-GP12) with potential antioxidant capacity were predicted. Density functional theory (DFT) calculations guided the selection of six peptides for synthesis and biological experiments. Molecular orbital representations revealed that the HOMO for these peptides is primarily localized on the indole segment, underscoring its pivotal role in antioxidant activity. All six synthesized peptides exhibited antioxidant activity in the DPPH assay, while the hydroxyl radical test showed suboptimal results. A hemolysis assay confirmed the non-hemolytic nature of the generated peptides. Additionally, an in silico investigation explored the potential inhibitory interaction between the peptides and the Keap1 protein. Analysis revealed that ligands GP3, GP4, and GP12 induced significant structural changes in proteins, affecting their stability and flexibility. These findings highlight the capability of machine learning approaches in generating novel antioxidant peptides.

In silico development of a novel anti-mutation, multi-epitope mRNA vaccine against MPXV variants of emerging lineage and sub-lineages by using immunoinformatics approaches

Over the past year, an unexpected surge in human monkeypox (hMPX) cases has been observed. This outbreak differs from previous ones, displaying distinct epidemiological characteristics and transmission patterns, believed to be influenced by a newly emerging monkeypox virus (MPXV) lineage. Notably, this emerging MPXV lineage has exhibited several non-synonymous mutations, some of which are linked to immunomodulatory activities and antigenic characteristics that aid in host detection. However, specific treatments or vaccines for human monkeypox are currently lacking. Hence, we aim to develop a multi-epitope mRNA vaccine by using immunoinformatics approaches against the MPXV, particularly its emerging variants. Six proteins (A29L, A35R, B6R, M1R, H3L, and E8L) were chosen for epitope and mutation site identification. Seventeen top-performing epitopes and eight epitopes containing mutation sites were selected and combined with adjuvants, the PADRE sequence, and linkers for vaccine development. The molecular and physical properties of the designed vaccine (WLmpx) were favorable. Immunological characteristics of WLmpx were assessed through molecular docking, molecular dynamics (MD) simulations, and immune simulations. Finally, the vaccine sequence was utilized to formulate an mRNA-based vaccine. The informatics-based predicted results indicated that the designed vaccine exhibits significant potential in eliciting high-level humoral and cellular immune responses, but further validation through in vivo and vitro studies is warranted.Communicated by Ramaswamy H. Sarma.

Repurposing of SARS-CoV-2 compounds against Marburg Virus using MD simulation, mm/GBSA, PCA analysis, and free energy landscape

The significant mortality rate associated with Marburg virus infection made it the greatest hazard among infectious diseases. Drug repurposing using in silico methods has been crucial in identifying potential compounds that could prevent viral replication by targeting the virus's primary proteins. This study aimed at repurposing the drugs of SARS-CoV-2 for identifying potential candidates against the matrix protein VP40 of the Marburg virus. Virtual screening was performed where the control compound, Nilotinib, showed a binding score of -9.99 kcal/mol. Based on binding scores, hit compounds 9549298, 11960895, 44545852, 51039094, and 89670174 were selected that had a lower binding score than the control. Subsequent molecular dynamics (MD) simulation revealed that compound 9549298 consistently formed a hydrogen bond with the residue Gln290. This was observed both in molecular docking and MD simulation poses, indicating a strong and significant interaction with the protein. 11960895 had the most stable and consistent RMSD pattern exhibited in 100 ns simulation, while 9549298 had the most identical RMSD plot compared to the control molecule. MM/PBSA analysis showed that the binding free energy (ΔG) of 9549298 and 11960895 was lower than the control, with -30.84 and -38.86 kcal/mol, respectively. It was observed by the PCA (principal component analysis) and FEL (free energy landscape) analysis that compounds 9549298 and 11960895 had lesser conformational variation. Overall, this study proposed 9549298 and 11960895 as potential binders of VP40 MARV that can cause its inhibition, however it inherently lacks experimental validation. Furthermore, the study proposes in-vitro experiments as the next step to validate these computational findings, offering a practical approach to further explore these compounds' potential as antiviral agents.Communicated by Ramaswamy H. Sarma.

In Silico Structural Modeling and Binding Site Analysis of Cerebroside Sulfotransferase (CST): A Therapeutic Target for Developing Substrate Reduction Therapy for Metachromatic Leukodystrophy

Cerebroside sulfotransferase (CST) is emerging as an important therapeutic target to develop substrate reduction therapy (SRT) for metachromatic leukodystrophy (MLD), a rare neurodegenerative lysosomal storage disorder. MLD develops with progressive impairment and destruction of the myelin sheath as a result of accumulation of sulfatide around the nerve cells in the absence of its recycling mechanism with deficiency of arylsulfatase A (ARSA). Sulfatide is the product of the catalytic action of cerebroside sulfotransferase (CST), which needs to be regulated under pathophysiological conditions by inhibitor development. To carry out in silico-based preliminary drug screening or for designing new drug candidates, a high-quality three-dimensional (3D) structure is needed in the absence of an experimentally derived three-dimensional crystal structure. In this study, a 3D model of the protein was developed using a primary sequence with the SWISS-MODEL server by applying the top four GMEQ score-based templates belonging to the sulfotransferase family as a reference. The 3D model of CST highlights the features of the protein responsible for its catalytic action. The CST model comprises five β-strands, which are flanked by ten α-helices from both sides as well as form the upside cover of the catalytic pocket of CST. CST has two catalytic regions: PAPS (-sulfo donor) binding and galactosylceramide (-sulfo acceptor) binding. The catalytic action of CST was proposed via molecular docking and molecular dynamic (MD) simulation with PAPS, galactosylceramide (GC), PAPS-galactosylceramide, and PAP. The stability of the model and its catalytic action were confirmed using molecular dynamic simulation-based trajectory analysis. CST response against the inhibition potential of the experimentally reported competitive inhibitor of CST was confirmed via molecular docking and molecular dynamics simulation, which suggested the suitability of the CST model for future drug discovery to strengthen substrate reduction therapy for MLD.

An amalgamated molecular dynamic and Gaussian based 3D-QSAR study for the design of 2,4-thiazolidinediones as potential PTP1B inhibitors

Overexpression of protein tyrosine phosphatase 1B (PTP1B) is the major cause of various diseases such as diabetes, obesity, and cancer. PTP1B has been identified as a negative regulator of the insulin signaling cascade, thereby causing diabetes. Numerous anti-diabetic medications based on thiazolidinedione have been successfully developed; however, 2,4-thiazolidinedione (2,4-TZD) scaffolds have been reported as potential PTP1B inhibitors for the manifestation of type 2 diabetes mellitus involving insulin resistance. In the present study, we have employed amalgamated approach involving MD-simulation studies (100 ns) as well as Gaussian field-based 3D-QSAR to develop a pharmacophoric model of 2,4-TZD as potent PTP1B inhibitors. MD simulation studies of the most potent compound in the PTP1B (PDB Id: 2QBS) binding pocket revealed that compound 43 was stable in the binding pocket and demonstrated excellent binding efficacy within the active site pocket. MM/GBSA results revealed that compound 43, bearing C-5 arylidine substitution, strongly bound to the target as compared to rosiglitazone with ΔGMM/GBSA difference of -11.13 kcal/mol. PCA, Rg, RMSF, RMSD, and SASA were analyzed from the complex's trajectories to anticipate the simulation outcome. We have suggested a series of 2,4-TZD as possible PTP1B inhibitors based on the results of MD simulation and 3D-QSAR studies.

Marine Brown Algae-Derived Compounds as Potential Inhibitors of Japanese Encephalitis Virus RNA-Dependent RNA Polymerase

The Japanese encephalitis virus (JEV) is a mosquito-borne flavivirus that primarily affects people in Asia and seriously threatens public health. Considering the rising occurrence rates and lack of targeted antiviral treatments, it is essential to comprehend and tackle obstacles related to JEV in order to lessen its influence on world health. This investigation explores compounds derived from marine brown algae (Phaeophyceae) as potential inhibitors of JEV RNA-dependent RNA polymerase (RdRp), a critical enzyme in the virus's replication cycle. Employing the computational virtual screen approach, four compounds, i.e., CMNPD16749, CMNPD2606, CMNPD27817, and CMNPD23662, with favorable binding energies ranging from -15.7 Kcal/mol to -13.9 kcal/mol were identified. Subsequently, through molecular docking analysis, the interactions responsible for the binding stability between the target protein and hit molecules compared to the reference molecule Galidesvir were studied. Further, through extensive molecular dynamic (MD) simulation studies at 200 ns, it was confirmed that each docked complex showed acceptable dynamic stability compared to the reference molecule. These findings were further validated using MM/PBSA free binding energy calculations, PCA analysis and free energy landscape construction. These computational findings suggested that the brown algae-derived compounds may act as an antiviral drug against JEV infection and lay a crucial foundation for future experimental studies against JEV.

Identification of amentoflavone as a potent SARS-CoV-2 Mpro inhibitor: a combination of computational studies and in vitro biological evaluation

Small-molecule inhibitors of SARS-CoV-2 Mpro that block the active site pocket of the viral main protease have been considered potential therapeutics for the development of drugs against SARS-CoV-2. Here, we report the identification of amentoflavone (a biflavonoid) through docking-based virtual screening of a library comprised of 231 compounds consisting of flavonoids and isoflavonoids. The docking results were further substantiated through extensive analysis of the data obtained from all-atom 150 ns MD simulation. End-state effective free energy calculations using MM-PBSA calculations further suggested that (Ra)-amentoflavone (C3'-C8''-atropisomer) may show a greater binding affinity towards the Mpro than (Sa)-amentoflavone. In vitro cytotoxicity assay established that amentoflavone showed a high CC50 value indicating much lower toxicity. Further, potent inhibition of the Mpro by amentoflavone was established by studying the effect on HEK293T cells treated with SARS-CoV-2 Mpro expressing plasmid.Communicated by Ramaswamy H. Sarma.

Computational insights into overcoming resistance mechanisms in targeted therapies for advanced breast cancer: focus on EGFR and HER2 co-inhibition

In the present study, the formation of a heterodimer involving both epidermal growth factor receptor (EGFR) and human epidermal growth factor receptor 2 (HER2) has been explored as a potential therapeutic mechanism to inhibit the progression of breast cancer. Virtual screening using molecular docking resulted in the three hit compounds (ZINC08382411, ZINC08382438, and ZINC08382292) with minimum binding scores and commonly binding to both receptors. Further, MD simulation analysis of these complexes illustrated the high stability of these compounds with EGFR and HER2. RMSD showed that ZINC08382411 displayed the most stable RMSD of 2 - 3 Å when bound to both receptors, suggesting to have strong compatibility with the active site of the receptor. Hydrogen bond analysis showed that ZINC08382411 forms the maximum number of H-bonds (2 to 3) in both EGFR and HER2 bound complexes, with the highest occupancy of 62% and 79%, respectively. Binding free energy calculation showed that ZINC08382411 possesses maximum affinity towards both the receptors with ΔGbind = -129.628 and -164.063 kJ/mol, respectively. This approach recognizes the significance of EGFR and HER2 in breast cancer development and aims to disrupt their collaborative signaling, which is known to promote the antagonistic behavior of cancer cells. By focusing on this EGFR/HER2 heterodimer, the study offers a promising avenue for identifying a potential candidate (ZINC08382411) that may inhibit breast cancer cell growth and potentially improve patient outcomes. The study's findings may contribute to the ongoing efforts to advance breast cancer treatment strategies.Communicated by Ramaswamy H. Sarma.

Investigating the effects of four medicinal plants against dengue virus through QSAR modeling and molecular dynamics studies

The Dengue virus (DENV) has been increasingly recognized as a prevalent viral pathogen responsible for global transmission of infection. It has been established that DENV's NS5 methyltransferase (MTase) controls viral replication. As a result, NS5 MTase is considered a potentially useful drug target for DENV. In this study, the two phases of virtual screening were conducted using the ML-based QSAR model and molecular docking to identify potential compounds against NS5 of DENV. Four medicinal plants [Aloe vera, Cannabis sativa (Hemp), Ocimum sanctum (Holy Basil; Tulsi), and Zingiber officinale (Ginger)] that showed anti-viral properties were selected for sourcing the phytochemicals and screening them against NS5. Additionally, re-docking at higher exhaustiveness and interaction analysis were performed which resulted in the identification of the top four hits (135398658, 5281675, 119394, and 969516) which showed comparable results with the control Sinefungin (SFG). Post molecular dynamics simulation, 135398658 showed the lowest RMSD (0.4-0.5 nm) and the maximum number of hydrogen bonds (eight hydrogen bonds) after the control while 5281675 and 969516 showed comparable hydrogen bonds to the control. These compounds showed direct interactions with the catalytic site residues GLU111 and ASP131, in addition to this these compounds showed stable complex formation as depicted by principal component analysis and free energy landscape. 135398658 showed lower total binding free energy (ΔGTotal = -36.56 kcal/mol) than the control, while 5281675 had comparable values to the control (ΔGTotal = -34.1 kcal/mol). Overall, the purpose of this study was to identify phytochemicals that inhibit NS5 function, that could be further tested experimentally to treat dengue virus (DENV).Communicated by Ramaswamy H. Sarma.

Structural and Dynamical Basis of VP35-RBD Inhibition by Marine Fungi Compounds to Combat Marburg Virus Infection

The Marburg virus (MBV), a deadly pathogen, poses a serious threat to world health due to the lack of effective treatments, calling for an immediate search for targeted and efficient treatments. In this study, we focused on compounds originating from marine fungi in order to identify possible inhibitory compounds against the Marburg virus (MBV) VP35-RNA binding domain (VP35-RBD) using a computational approach. We started with a virtual screening procedure using the Lipinski filter as a guide. Based on their docking scores, 42 potential candidates were found. Four of these compounds-CMNPD17596, CMNPD22144, CMNPD25994, and CMNPD17598-as well as myricetin, the control compound, were chosen for re-docking analysis. Re-docking revealed that these particular compounds had a higher affinity for MBV VP35-RBD in comparison to the control. Analyzing the chemical interactions revealed unique binding properties for every compound, identified by a range of Pi-cation interactions and hydrogen bond types. We were able to learn more about the dynamic behaviors and stability of the protein-ligand complexes through a 200-nanosecond molecular dynamics simulation, as demonstrated by the compounds' consistent RMSD and RMSF values. The multidimensional nature of the data was clarified by the application of principal component analysis, which suggested stable conformations in the complexes with little modification. Further insight into the energy profiles and stability states of these complexes was also obtained by an examination of the free energy landscape. Our findings underscore the effectiveness of computational strategies in identifying and analyzing potential inhibitors for MBV VP35-RBD, offering promising paths for further experimental investigations and possible therapeutic development against the MBV.

Inhibition of SARS-CoV-2 NSP-15 by Uridine-5'-Monophosphate Analogues Using QSAR Modelling, Molecular Dynamics Simulations, and Free Energy Landscape

SARS-CoV-2 is accountable for severe social and economic disruption around the world causing COVID-19. Non-structural protein-15 (NSP15) possesses a domain that is vital to the viral life cycle and is known as uridylate-specific endoribonuclease (EndoU). This domain binds to the uridine 5'-monophosphate (U5P) so that the protein may carry out its native activity. It is considered a vital drug target to inhibit the growth of the virus. Thus, in this current study, ML-based QSAR and virtual screening of U5P analogues targeting Nsp15 were performed to identify potential molecules against SARS-CoV-2. Screening of 816 unique U5P analogues using ML-based QSAR identified 397 compounds ranked on their predicted bioactivity (pIC50). Further, molecular docking and hydrogen bond interaction analysis resulted in the selection of the top three compounds (53309102, 57398422, and 76314921). Molecular dynamics simulation of the most promising compounds showed that two molecules 53309102 and 57398422 acted as potential binders of Nsp15. The compound was able to inhibit nsp15 activity as it was successfully bound to the active site of the nsp15 protein. This was achieved by the formation of relevant contacts with enzymatically critical amino acid residues (His235, His250, and Lys290). Principal component analysis and free energy landscape studies showed stable complex formation while MM/GBSA calculation showed lower binding energies for 53309102 (ΔGTOTAL = -29.4 kcal/mol) and 57398422 (ΔGTOTAL = -39.4 kcal/mol) compared to the control U5P (ΔGTOTAL = -18.8 kcal/mol). This study aimed to identify analogues of U5P inhibiting the NSP15 function that potentially could be used for treating COVID-19.

Investigating the bispecific lead compounds against methicillin-resistant Staphylococcus aureus SarA and CrtM using machine learning and molecular dynamics approach

Methicillin-resistant Staphylococcus aureus (MRSA) is a notorious pathogen that has emerged as a serious global health concern over the past few decades. Staphylococcal accessory regulator A (SarA) and 4,4'-diapophytoene synthase (CrtM) play a crucial role in biofilm formation and staphyloxanthin biosynthesis. Thus, the present study used a machine learning-based QSAR model to screen 1261 plant-derived natural organic compounds in order to identify a medication candidate with both biofilm and virulence inhibitory potential. Additionally, the in-silico molecular docking analysis has demonstrated significant binding efficacy of the identified hit compound, that is 85137543, with SarA and CrtM when compared to the control compound, hesperidin. Post-MD simulation analysis of the complexes depicted strong binding of 85137543 to both SarA and CrtM. Moreover, 85137543 showed hydrogen bonding with the key residues of both proteins during docking (ALA138 of SarA and ALA134 of CrtM) and post-MD simulation (LYS273 of CrtM and ASN212 of SarA). The RMSD of 85137543 was stable and consistent when bound to both CrtM and SarA with RMSDs of 1.3 and 1 nm, respectively. In addition, principal component analysis and the free energy landscape showed stable complex formation with both proteins. Low binding free energy (ΔGTotal) was observed by 85137543 for SarA (-47.92 kcal/mol) and CrtM (-36.43 kcal/mol), which showed strong binding. Overall, this study identified 85137543 as a potential inhibitor of both SarA and CrtM in MRSA.Communicated by Ramaswamy H. Sarma.

The influence of the hydrogen-bond network on the structure and dynamics of the RAPRKKG heptapeptide and its mutants

The structural behaviour of the RAPRKKG heptapeptide after individual or multiple mutations was inspected through molecular dynamics simulation. The nature of the mutations provided information on the flexibility of the heptapeptide and on how water molecules establish hydrogen bonds with it. The structural behaviour of the wild-type and the mutated structures were measured through the analysis of protein‒protein and protein‒solvent hydrogen bonds. The conformational behaviours of the different structures were analysed through free energy landscape analysis. The flexibility characteristics of the mutants seem to depend on the reorganization of water molecules and their static or dynamic behaviour around amino acid side chains.

Sappanin-type homoisoflavonoids from Scilla nervosa inhibits acetylcholinesterase enzyme: a combined in silico and in vitro approach

Alzheimer's disease is among the major health challenges that currently attract the attention of health care givers and drug discovery and development experts worldwide. This study investigated the acetylcholinesterase inhibitory activity of sappanin-type homisoflavonoids isolated from the inter-bulb surface of Scilla nervosa. Molecular docking, molecular dynamics simulation, ADMET and in vitro studies were performed to identify the hit molecules, understand their binding mode and interaction, druggability and establish their inhibitory potentials against acetylcholinesterase enzyme. The hit compounds 5, 2, 1 and 4 were identified as the hit molecules through the molecular docking. The molecular dynamics simulation and MM-PBSA analysis showed the hit homoisoflavonoids established stability and good binding affinity against the acetylcholinesterase enzyme. Also, 5 elicited the best inhibitory activity followed by 2, 1 and 4 in the in vitro experiment. Furthermore, the selected homoisoflavonoids exhibit interesting drug likeness and pharmacokinetic properties as drug candidate. The results suggest further investigations towards the development of the phytochemicals as possible acetylcholinesterase inhibitors.Communicated by Ramaswamy H. Sarma.

Synergistic effect of conformational changes in phosphoglycerate kinase 1 product release

In the glycolysis pathway, phosphoglycerate kinase 1 (PGK1) transfers one phosphoryl-group from 1,3-diphosphoglycerate (1,3BPG) to ADP to product 3-phosphoglycerate (3PG) and ATP. The catalytic process is accompanied with the conversion between the open conformation and the closed conformation of PGK1. However, the dynamic collaboration mechanism between the PGK1 conformation transition and the products releasing process remains poorly understood. Here using molecular dynamics simulations combined with molecular mechanics generalized born surface area (MM/GBSA) analysis, we demonstrated that PGK1 in the closed conformation first releases the product ATP to reach a semi-open conformation, and releases the product 3PG to achieve the full open conformation, which could accept new substrates ADP and 1,3BPG for the next cycle. It is noteworthy that the phosphorylation of PGK1 at T243 causes the loop region (residues L248-E260) flip outside the protein, and the phosphorylation of Y324 leads PGK1 become looser. Both modifications cause the exposure of the ADP/ATP binding site, which was beneficial for the substrates/products binding/releasing of PGK1. In addition, the other post translational modifications (PTMs) were also able to regulate the ligands binding/releasing with different effects. Our results revealed the dynamic cooperative molecular mechanism of PGK1 conformational transition with products releasing, as well as the influence of PTMs, which would contribute to the understanding of PGK1 substrates/products conversion process and the development of small molecule drugs targeting PGK1.Communicated by Ramaswamy H. Sarma.

Protein Kinase Structure and Dynamics: Role of the αC-β4 Loop

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, Local Spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

Unraveling the mechanisms of Cefoxitin resistance in methicillin-resistant Staphylococcus aureus (MRSA): structural and molecular simulation-based insights

Methicillin-resistant Staphylococcus aureus (MRSA) severely affects human health, including the skin glands, nasal cavity, wound infections, bone infections, and pneumonia. Among the most effective MRSA drugs, Cefoxitin also develops resistance due to mutations in the mecA gene. Four mutations at positions E229K, E239R, G246K, and E447K are classified as high-level resistance mutations. However, the resistance mechanism of MRSA towards Cefoxitin caused by these mutations is still unclear, as there is less information available regarding the structural and functional effects of the mutations against Cefoxitin. Therefore, our present study was designed to explore the mechanisms of binding interactions between wild-type and mutated PBP2a against Cefoxitin using molecular docking and MD simulations. Subsequently, we identified that the mutant form of PBP2a affects the activity of Cefoxitin. Interestingly, the binding of Cefoxitin with G246K and E239R mutants demonstrates unstable behavior compared to E447K-Cefoxitin and E229K-Cefoxitin. In this study, we propose the resistance mechanism of Cefoxitin at the atomic level. The proposed drug-resistance mechanism will provide valuable guidance for the design of MRSA drugs. This research might provide a new framework for designing new agents against the mutated form of PBP2a.Communicated by Ramaswamy H. Sarma.

Identification of N-(3-(methyl(3-(orotic amido)propyl)amino)propyl) oleanolamide as a novel topoisomerase I catalytic inhibitor by rational design, molecular dynamics simulation, and biological evaluation

DNA topoisomerase I (TOP1) catalytic inhibitors are a promising class of antitumor agents. Oleanolic acid derivatives are potential TOP1 catalytic inhibitors. However, their inhibitory activity still needs to be enhanced, and the stability and hotspot residue sites of their interaction with TOP1 remain to be elucidated. Herein, a novel oleanolic acid derivative, OA4 (N-(3-(methyl(3-(orotic amido)propyl)amino)propyl)oleanolamide), was identified by rational design. Subsequently, molecular dynamics simulations were performed to explore the stability and conformational dynamics of the TOP1-OA4 complex. The molecular mechanics/generalized Born surface area method calculated the binding free energy and predicted Arg488, Ile535, and His632 to be hotspot residues. Biological experiments verified that OA4 is a nonintercalative TOP1 catalytic inhibitor. OA4 exhibits better proliferation inhibitory activity against tumor cells than normal cells. Furthermore, OA4 can induce apoptosis and effectively suppress the proliferation and migration of cancer cells. This work provides new insights for the development of novel TOP1 catalytic inhibitors.

An immunoinformatics and extended molecular dynamics approach for designing a polyvalent vaccine against multiple strains of Human T-lymphotropic virus (HTLV)

Human T-lymphotropic virus (HTLV), a group of retroviruses belonging to the oncovirus family, has long been associated with various inflammatory and immunosuppressive disorders. At present, there is no approved vaccine capable of effectively combating all the highly pathogenic strains of HTLV that makes this group of viruses a potential threat to human health. To combat the devastating impact of any potential future outbreak caused by this virus group, our study employed a reverse vaccinology approach to design a novel polyvalent vaccine targeting the highly virulent subtypes of HTLV. Moreover, we comprehensively analyzed the molecular interactions between the designed vaccine and corresponding Toll-like receptors (TLRs), providing valuable insights for future research on preventing and managing HTLV-related diseases and any possible outbreaks. The vaccine was designed by focusing on the envelope glycoprotein gp62, a crucial protein involved in the infectious process and immune mechanisms of HTLV inside the human body. Epitope mapping identified T cell and B cell epitopes with low binding energies, ensuring their immunogenicity and safety. Linkers and adjuvants were incorporated to enhance the vaccine's stability, antigenicity, and immunogenicity. Initially, two vaccine constructs were formulated, and among them, vaccine construct-2 exhibited superior solubility and structural stability. Molecular docking analyses also revealed strong binding affinity between the vaccine construct-2 and both targeted TLR2 and TLR4. Molecular dynamics simulations demonstrated enhanced stability, compactness, and consistent hydrogen bonding within TLR-vaccine complexes, suggesting a strong binding affinity. The stability of the complexes was further corroborated by contact, free energy, structure, and MM-PBSA analyses. Consequently, our research proposes a vaccine targeting multiple HTLV subtypes, offering valuable insights into the molecular interactions between the vaccine and TLRs. These findings should contribute to developing effective preventive and treatment approaches against HTLV-related diseases and preventing possible outbreaks. However, future research should focus on in-depth validation through experimental studies to confirm the interactions identified in silico and to evaluate the vaccine's efficacy in relevant animal models and, eventually, in clinical trials.

Subtractive proteomics analysis to uncover the potent drug targets for distinctive drug design of Candidaauris

Candida auris is a serious health concern of the current world that possesses a serious global health threat and is emerging at a high rate. Available antifungal drugs are failing to combat this pathogen as they are growing resistant to those drugs and some strains have already shown resistance to all three available antifungal drugs in the market. Hence, finding alternative therapies is essential for saving lives from this enemy. To make the development of new treatments easier, we conducted some in silico study of this pathogen to discover possible targets for drug design and also recommended some possible metabolites to test in vivo circumstances. The complete proteome of the representative strain was retrieved, and the duplicate, non-essential, human homologous, non-metabolic, and druggable proteins were then eliminated. As a result, out of a total of 5441 C. auris proteins, we were able to isolate three proteins (XP 028890156.1, XP 028891672.1, and XP 028891858.1) that are crucial for the pathogen's survival as well as host-non-homolog, metabolic, and unrelated proteins to the human microbiome. Their subcellular locations and interactions with a large number of proteins (10 proteins) further point to them being good candidates for therapeutic targets. Following in silico docking of 29 putative antifungals of plant origin against the three proteins we chose, Caledonixanthone E, Viniferin, Glaucine, and Jatrorrhizine were discovered to be the most effective means of inhibiting those proteins since they displayed higher binding affinities (ranging from -28.97 kcal/mol to -51.99 kcal/mol) than the control fluconazole (which ranged between -28.84 kcal/mol and -41.15 kcal/mol). According to the results of MD simulations and MM-PBSA calculations, Viniferin and Caledonixanthone E are the most effective ligands for the proteins XP 028890156.1, XP 028891672.1, and XP 028891858.1. Furthermore, they were predicted to be safe and also showed proper ADME properties.

Molecular dynamics simulations corroborate recombinant expression studies carried out on three αIIb β-propeller mutations reported in Indian Glanzmann thrombasthenia patients

Mutations in the αIIb β-propeller domain have long been known to disrupt heterodimerization and intracellular trafficking of αIIbβ3 complexes leading to diminished surface expression and/or function, resulting in Glanzmann thrombasthenia. Our previous study on three β-propeller mutations, namely G128S, S287L, and G357S, showed variable defects in protein transport correlated with the patient's clinical phenotypes. Pulse-chase experiments revealed differences in αIIbβ3 complex maturation among the three mutations. Hence, the current study aims to correlate conformational changes caused by each one of them. Evolutionary conservation analysis, stability analysis, and molecular dynamics simulations of the three mutant structures were carried out. Stability analysis revealed that, while G128S and G357S mutations destabilized the β-propeller structure, S287L retained the stability. Wild-type and mutant β-propeller structures, when subjected to molecular dynamics simulations, confirmed that G128S and G357S were both destabilizing in nature when compared with the wild-type and S287L based on several parameters studied, like RMSD, RMSF, Rg, FEL, PCA, secondary structure, and hydrogen bonds. In our previous study, we demonstrated that mutant S287L αIIbβ3 complexes were more stable than the wild-type αIIbβ3 complexes, as evidenced in pulse-chase experiments. These findings corroborate variable intracellular fates of mutant αIIbβ3 complexes as a result of these β-propeller mutations.

Enhanced stability of a disaggregated Aβ fibril on removal of ligand inhibits refibrillation: An all atom Molecular Dynamics simulation study

The extraneuronally deposited senile plaques, composed of neurotoxic aggregates of Aβ fibril, define Alzheimer's disease (AD). Natural compounds have been tested for their destabilization potential on Aβ fibril, thereby curing AD. However, the resultant destabilized Aβ fibril, needs to be checked for its irreversibility to the native organized state after removal of the ligand. Herein, we assessed the stability of a destabilized fibril after the ligand (ellagic acid represented as REF) is removed from the complex. The study has been conducted via Molecular Dynamics (MD) simulation of 1 μs for both Aβ-Water (control) and Aβ-REF″ (test or REF removed) system. The increased value of RMSD, R_g, SASA, lower β-sheet content and reduced number of H-bonds explains enhanced destabilization observed in Aβ-REF″ system. The increased inter-chain distance demonstrates breaking of the residual contacts, testifying the drift of terminal chains from the pentamer. The increased SASA along with the ∆G_ps(polar solvation energy) accounts for the reduced interaction amongst residues, and more with solvent molecules, governing irreversibility to native state. The higher Gibb's free energy of the misaligned structure of Aβ-REF″ ensures irreversibility to the organized structure due to its inability to cross such high energy barrier. The observed stability of the disaggregated structure, despite ligand elimination, establishes the effectiveness of the destabilization technique as a promising therapeutic approach towards treating AD.

Artificial intelligence assisted identification of potential tau aggregation inhibitors: ligand- and structure-based virtual screening, in silico ADME, and molecular dynamics study

Alzheimer's disease (AD) is a severe, growing, multifactorial disorder affecting millions of people worldwide characterized by cognitive decline and neurodegeneration. The accumulation of tau protein into paired helical filaments is one of the major pathological hallmarks of AD and has gained the interest of researchers as a potential drug target to treat AD. Lately, Artificial Intelligence (AI) has revolutionized the drug discovery process by speeding it up and reducing the overall cost. As a part of our continuous effort to identify potential tau aggregation inhibitors, and leveraging the power of AI, in this study, we used a fully automated AI-assisted ligand-based virtual screening tool, PyRMD to screen a library of 12 million compounds from the ZINC database to identify potential tau aggregation inhibitors. The preliminary hits from virtual screening were filtered for similar compounds and pan-assay interference compounds (the compounds containing reactive functional groups which can interfere with the assays) using RDKit. Further, the selected compounds were prioritized based on their molecular docking score with the binding pocket of tau where the binding pockets were identified using replica exchange molecular dynamics simulation. Thirty-three compounds showing good docking scores for all the tau clusters were selected and were further subjected to in silico pharmacokinetic prediction. Finally, top 10 compounds were selected for molecular dynamics simulation and MMPBSA binding free energy calculations resulting in the identification of UNK_175, UNK_1027, UNK_1172, UNK_1173, UNK_1237, UNK_1518, and UNK_2181 as potential tau aggregation inhibitors.

Development of multi-epitope vaccines against the monkeypox virus based on envelope proteins using immunoinformatics approaches

Background

Since May 2022, cases of monkeypox, a zoonotic disease caused by the monkeypox virus (MPXV), have been increasingly reported worldwide. There are, however, no proven therapies or vaccines available for monkeypox. In this study, several multi-epitope vaccines were designed against the MPXV using immunoinformatics approaches.

Methods

Three target proteins, A35R and B6R, enveloped virion (EV) form-derived antigens, and H3L, expressed on the mature virion (MV) form, were selected for epitope identification. The shortlisted epitopes were fused with appropriate adjuvants and linkers to vaccine candidates. The biophysical andbiochemical features of vaccine candidates were evaluated. The Molecular docking and molecular dynamics(MD) simulation were run to understand the binding mode and binding stability between the vaccines and Toll-like receptors (TLRs) and major histocompatibility complexes (MHCs). The immunogenicity of the designed vaccines was evaluated via immune simulation.

Results

Five vaccine constructs (MPXV-1-5) were formed. After the evaluation of various immunological and physicochemical parameters, MPXV-2 and MPXV-5 were selected for further analysis. The results of molecular docking showed that the MPXV-2 and MPXV-5 had a stronger affinity to TLRs (TLR2 and TLR4) and MHC (HLA-A*02:01 and HLA-DRB1*02:01) molecules, and the analyses of molecular dynamics (MD) simulation have further confirmed the strong binding stability of MPXV-2 and MPXV-5 with TLRs and MHC molecules. The results of the immune simulation indicated that both MPXV-2 and MPXV-5 could effectively induce robust protective immune responses in the human body.

Conclusion

The MPXV-2 and MPXV-5 have good efficacy against the MPXV in theory, but further studies are required to validate their safety and efficacy.

Effect of tryptophan mutation on the structure of LOV1 domain of phototropin1 protein of Ostreococcus tauri: A combined molecular dynamics simulation and biophysical approach

BACKGROUND

Light, oxygen and voltage (LOV) proteins detect blue light by formation of a covalent 'photoadduct' between the flavin chromophore and the neighboring conserved cysteine residue. LOV proteins devoid of this conserved photoactive cysteine are unable to form this 'photoadduct' upon light illumination, but they can still elicit functional response via the formation of neutral flavin radical. Recently, tryptophan residue has been shown to be the primary electron donors to the flavin excited state.

METHODS

Photoactive cysteine (Cys42) and tryptophan (Trp68) residues in the LOV1 domain of phototropin1 of Ostreococcus tauri (OtLOV1) was mutated to alanine and threonine respectively. Effect of these mutations have been studied using molecular dynamics simulation and spectroscopic techniques.

RESULTS

Molecular dynamics simulation indicated that W68T did not affect the structure of OtLOV1 protein, but C42A leads to some structural changes. An increase in the fluorescence lifetime and quantum yield values was observed for the Trp68 mutant.

CONCLUSIONS

An increase in the fluorescence lifetime and quantum yield of Trp68 mutant compared to the wild type protein suggests that Trp68 residue participates in quenching of the flavin excited state followed by photoexcitation.

GENERAL SIGNIFICANCE

Enhanced photo-physical properties of Trp68 OtLOV1 mutant might enable its use for the optogenetic and microscopic applications.

Revealing Free Energy Landscape From MD Data via Conditional Angle Partition Tree

Deciphering the free energy landscape of biomolecular structure space is crucial for understanding many complex molecular processes, such as protein-protein interaction, RNA folding, and protein folding. A major source of current dynamic structure data is Molecular Dynamics (MD) simulations. Several methods have been proposed to investigate the free energy landscape from MD data, but all of them rely on the assumption that kinetic similarity is associated with global geometric similarity, which may lead to unsatisfactory results. In this paper, we proposed a new method called Conditional Angle Partition Tree to reveal the hierarchical free energy landscape by correlating local geometric similarity with kinetic similarity. Its application on the benchmark alanine dipeptide MD data showed a much better performance than existing methods in exploring and understanding the free energy landscape. We also applied it to the MD data of Villin HP35. Our results are more reasonable on various aspects than those from other methods and very informative on the hierarchical structure of its energy landscape.

Structure-activity correlations for peptaibols obtained from clade Longibrachiatum of Trichoderma: A combined experimental and computational approach

Integrated disease management and plant protection have been discussed with much fervor in the past decade due to the rising environmental concerns of using industrially produced pesticides. Members of the genus Trichoderma are a subject of considerable research today due to their several properties as biocontrol agents. In our study, the peptaibol production of Trichoderma longibrachiatum SZMC 1775, T. longibrachiatum f. bissettii SZMC 12546, T. reesei SZMC 22616, T. reesei SZMC 22614, T. saturnisporum SZMC 22606 and T. effusum SZMC 22611 were investigated to elucidate structure-activity relationships (SARs) between the properties of peptaibols and their 3D structures. The effects of peptaibol mixtures obtained from every Trichoderma strain were examined against nine commonly known bacteria. The lowest minimum inhibitory concentrations (MIC, mg ml-1) were exerted by T. longibrachiatum f. bissettii SZMC 12546 against Gram-positive bacteria, which was also able to inhibit the plant pathogenic Gram-negative Rhizobium radiobacter. Accelerated molecular dynamics (aMD) simulations were performed in aqueous solvent to explore the folding dynamics of 12 selected peptaibol sequences. The most characteristic difference between the peptaibols from group A and B relies in the 'Gly-Leu-Aib-Pro' and 'Gly-Aib-Aib-Pro' motifs ('Aib' stands for α-aminoisobutyric acid), which imparted a significant effect on the folding dynamics in water and might be correlated with their expressed bioactivity. In our aMD simulation experiments, Group A peptaibols showed more restricted folding dynamics with well-folded helical conformations as the most stable representative structures. This structural stability and dynamics may contribute to their bioactivity against the selected bacterial species.

Gold Nanorods Inhibit Tumor Metastasis by Regulating MMP-9 Activity: Implications for Radiotherapy

Dysregulation of matrix metalloproteinase (MMP) is strongly implicated in tumor invasion and metastasis. Nanomaterials can interact with proteins and have impacts on protein activity, which provides a potential strategy for inhibiting tumor invasion and metastasis. However, the regulation of MMP activity by nanomaterials has not been fully determined. Herein, we have found that gold nanorods (Au NRs) are able to induce the change of the secondary structure of MMP-9 and thereby inhibit their activity. Interestingly, the inhibition of MMP-9 activity is highly dependent on the aspect ratio of Au NRs, and an aspect ratio of 3.3 shows the maximum inhibition efficiency. Molecular dynamics simulations combined with mathematical statistics algorithm reveal the binding behaviors and interaction modes of MMP-9 with Au NRs in atomic details and disclose the mechanism of aspect ratio-dependent inhibition effect of Au NRs on MMP-9 activity. Au NRs with an aspect ratio of 3.3 successfully suppress the X-ray-activated invasion and metastasis of tumor by inhibiting MMP-9 activity. Our findings provide important guidance for the modulation of MMP-9 activity by tuning key parameters of nanomaterials and demonstrate that gold nanorods could be developed as potential MMP inhibitors.

Deciphering the structural and functional impact of missense mutations in Egr1-DNA interacting interface: an integrative computational approach

Early growth response-1 (Egr1) is a zinc-finger transcription factor that plays a critical role in controlling cell growth, proliferation, differentiation, angiogenesis, and apoptosis. Egr1 is induced by many growth factors, cytokines, and stress signals and is also known to be involved in several pathological conditions like cancer, neurological and ocular disorders. The DNA binding domain of Egr1 is a highly conserved Cys2His2 (C2H2) zinc finger (ZNF) domain which specifically binds to GC-rich consensus sequence GcG (G/T) GGGCG and activates transcription. As the C2H2 domain specifically recognizes its DNA target, the mutations spanning this region shall perturb DNA recognition and may hinder transcription of target genes. Therefore, in this study, the missense mutations occurring specifically at the DNA binding domain (DBD) of Egr1 were probed by computational approaches involving in silico screening of pathogenic and functional mutants coupled with extensive molecular dynamics simulations, to determine the mutants that affect its structural stability and interactions with DNA. From the pathogenicity analysis of 38 missense mutations spanning Egr1-DBD, 17 were predicted as pathogenic, and 7 amongst these were found to have functional impact on Egr1. On combined analysis of molecular dynamics simulation, Residue interaction analysis and Egr1-DNA interaction analysis results, the mutants R371C and R375C showed least impact, whilst, H382R tend to increase the structural stability, whereas R360H, H390R, E393V, and H414Y conferred greater impact by altering the structural stability and DNA interactions. Hence, this study exposes the prospects of considering these 4 deleterious mutations for clinical significance, but needs further experimental validation.HighlightsEgr1's DNA binding domain is a highly conserved Cys2His2 (C2H2) zinc finger domain that specifically recognizes its DNA target.Mutations spanning in the DNA binding domain shall perturb DNA recognition and may hinder transcription.Among the missense mutations, mutants R360H, H390R, E393V, and H414Y were inferred to have a greater impact on Egr1 by altering the structural stability and DNA interactions.

Immunoinformatics Approach to Design Novel Subunit Vaccine against the Epstein-Barr Virus

Epstein-Barr virus (EBV) is a lymphotropic virus responsible for numerous epithelial and lymphoid cell malignancies, including gastric carcinoma, Hodgkin's lymphoma, nasopharyngeal carcinoma, and Burkitt lymphoma. Hundreds of thousands of people worldwide get infected with this virus, and in most cases, this viral infection leads to cancer. Although researchers are trying to develop potential vaccines and drug therapeutics, there is still no effective vaccine to combat this virus. In this study, the immunoinformatics approach was utilized to develop a potential multiepitope subunit vaccine against the two most common subtypes of EBV, targeting three of their virulent envelope glycoproteins. Eleven cytotoxic T lymphocyte (CTL) epitopes, 11 helper T lymphocyte (HTL) epitopes, and 10 B-cell lymphocyte (BCL) epitopes were predicted to be antigenic, nonallergenic, nontoxic, and fully conserved among the two subtypes, and nonhuman homologs were used for constructing the vaccine after much analysis. Later, further validation experiments, including molecular docking with different immune receptors (e.g., Toll-like receptors [TLRs]), molecular dynamics simulation analyses (including root means square deviation [RMSD], root mean square fluctuation [RMSF], radius of gyration [Rg], principal-component analysis [PCA], dynamic cross-correlation [DCC], definition of the secondary structure of proteins [DSSP], and Molecular Mechanics Poisson-Boltzmann Surface Area [MM-PBSA]), and immune simulation analyses generated promising results, ensuring the safe and stable response of the vaccine with specific immune receptors after potential administration within the human body. The vaccine's high binding affinity with TLRs was revealed in the docking study, and a very stable interaction throughout the simulation proved the potential high efficacy of the proposed vaccine. Further, in silico cloning was also conducted to design an efficient mass production strategy for future bulk industrial vaccine production. IMPORTANCE Epstein-Barr virus (EBV) vaccines have been developing for over 30 years, but polyphyletic and therapeutic vaccines have failed to get licensed. Our vaccine surpasses the limitations of many such vaccines and remains very promising, which is crucial because the infection rate is higher than most viral infections, affecting a whopping 90% of the adult population. One of the major identifications covers a holistic analysis of populations worldwide, giving us crucial information about its effectiveness for everyone's unique immunological system. We targeted three glycoproteins that enhance the virulence of the virus to design an epitope-based polyvalent vaccine against two different strains of EBV, type 1 and 2. Our methodology in this study is nonconventional yet swift to show effective results while designing vaccines.

Molecular basis of the anchoring and stabilization of human islet amyloid polypeptide in lipid hydroperoxidized bilayers

The molecular structure of membrane lipids is formed by mono- or polyunsaturations on their aliphatic tails that make them susceptible to oxidation, facilitating the incorporation of hydroperoxide (R-OOH) functional groups. Such groups promote changes in both composition and complexity of the membrane significantly modifying its physicochemical properties. Human Langerhans islets amyloid polypeptide (hIAPP) is the main component of amyloid deposits found in the pancreas of patients with type-2 diabetes (T2D). hIAPP in the presence of membranes with oxidized lipid species accelerates the formation of amyloid fibrils or the formation of intermediate oligomeric structures. However, the molecular bases at the initial stage of the anchoring and stabilization of the hIAPP in a hydroperoxidized membrane are not yet well understood. To shed some light on this matter, in this contribution, three bilayer models were modeled: neutral (POPC), anionic (POPS), and oxidized (POPC_OOH), and full atom Molecular Dynamics (MD) simulations were performed. Our results show that the POPC_OOH bilayer increases the helicity in hIAPP when compared to POPC or POPS bilayer. The modification in the secondary structure covers the residues of the so-called amyloidogenic core of the hIAPP. Overall, the hydroperoxidation of the neutral lipids modifies both the anchoring and the stabilization of the peptide hIAPP by reducing the random conformations of the peptide and increasing of hydrogen bond population with the hydroperoxidized lipids.

Insight into the Nucleoside Transport and Inhibition of Human ENT1

The human equilibrative nucleoside transporter 1 (hENT1) is an effective controller of adenosine signaling by regulating its extracellular and intracellular concentration, and has become a solid drug target of clinical used adenosine reuptake inhibitors (AdoRIs). Currently, the mechanisms of adenosine transport and inhibition for hENT1 remain unclear, which greatly limits the in-depth understanding of its inner workings as well as the development of novel inhibitors. In this work, the dynamic details of hENT1 underlie adenosine transport and the inhibition mechanism of the non-nucleoside AdoRIs dilazep both were investigated by comparative long-time unbiased molecular dynamics simulations. The calculation results show that the conformational transitions of hENT1 from the outward open to metastable occluded state are mainly driven by TM1, TM2, TM7 and TM9. One of the trimethoxyphenyl rings in dilazep serves as the adenosyl moiety of the endogenous adenosine substrate to competitively occupy the orthosteric site of hENT1. Due to extensive and various VDW interactions with N30, M33, M84, P308 and F334, the other trimethoxyphenyl ring is stuck in the opportunistic site near the extracellular side preventing the complete occlusion of thin gate simultaneously. Obviously, dilazep shows significant inhibitory activity by disrupting the local induce-fit action in substrate binding cavity and blocking the transport cycle of whole protein. This study not only reveals the nucleoside transport mechanism by hENT1 at atomic level, but also provides structural guidance for the subsequent design of novel non-nucleoside AdoRIs with enhanced pharmacologic properties.

•

The transitions of hENT1 from the outward open to metastable occluded state are mainly driven by TM1, TM2, TM7 and TM9.

•

The induce-fit action by adenosine recognition precedes. inward contraction of the extracellular side.

•

Dilazep exerts its special hENT1 inhibitory function through competitive binding and allosteric regulation.

•

A gating strategy of extracellular loop is revealed to ensure adenosine is firmly located in the transport cavity.

Molecular rationale for the impairment of the MexAB-OprM efflux pump by a single mutation in MexA

Efflux pumps of the Resistance-Nodulation-cell Division (RND) superfamily contribute to intrinsic and acquired resistance in Gram-negative pathogens by expelling chemically unrelated antibiotics with high efficiency. They are tripartite systems constituted by an inner-membrane-anchored transporter, an outer membrane factor protein, and a membrane fusion protein. Multimerization of the membrane fusion protein is an essential prerequisite for full functionality of these efflux pumps. In this work, we employed complementary computational techniques to investigate the stability of a dimeric unit of MexA (the membrane fusion protein of the MexAB-OprM RND efflux pump of Pseudomonas aeruginosa), and to provide a molecular rationale for the effect of the G72S substitution, which affects MexAB-OprM functionality by impairing the assembly of MexA. Our findings indicate that: i) dimers of this protein are stable in multiple µs-long molecular dynamics simulations; ii) the mutation drastically alters the conformational equilibrium of MexA, favouring a collapsed conformation that is unlikely to form dimers or higher order assemblies. Unveiling the mechanistic aspects underlying large conformational distortions induced by minor sequence changes is informative to efforts at interfering with the activity of this elusive bacterial weapon. In this respect, our work further confirms how molecular simulations can give important contribution and useful insights to characterize the mechanism of highly complex biological systems.