Molecular Dynamics: New Frontier in Personalized Medicine

Identification of potential inhibitors for drug resistance in acute lymphoblastic leukemia through differentially expressed gene analysis and in silico screening

Acute lymphoblastic leukemia (ALL) is a disease of lymphocyte origin predominantly diagnosed in children. While its 5-year survival rate is high, resistance to chemotherapy drugs is still an obstacle. Our aim is to determine differentially expressed genes (DEGs) related to Asparaginase, Daunorubicin, Prednisolone, and Vincristine resistance and identify potential inhibitors via docking. Three datasets were accessed from the Gene Expression Omnibus database; GSE635, GSE19143, and GSE22529. The microarray data was analyzed using R4.2.0 and Bioconductor packages, and pathway and protein-protein interaction analysis were performed. We identified 1294 upregulated DEGs, with 12 genes consistently upregulated in all four resistant groups. KEGG analysis revealed an association with the PI3K-Akt pathway. Among DEGs, 33 hub genes including MDM2 and USP7 were pinpointed. Within common genes, CLDN9 and HS3ST3A1 were subjected to molecular docking against 3556 molecules. Following ADMET analysis, three drugs emerged as potential inhibitors: Flunarizine, Talniflumate, and Eltrombopag. Molecular dynamics analysis for HS3ST3A1 indicated all candidates had the potential to overcome drug resistance, Eltrombopag displaying particularly promising results. This study promotes a further understanding of drug resistance in ALL, introducing novel genes for consideration in diagnostic screening. It also presents potential inhibitor candidates to tackle drug resistance through repurposing.

Computational exploration of natural compounds targeting Staphylococcus aureus: inhibiting AgrA promoter binding for antimicrobial intervention

Staphylococcus aureus is a highly virulent nosocomial pathogen that poses a significant threat to individuals exposed to healthcare settings. Due to its sophisticated machinery for producing virulence factors, S. aureus can cause severe and potentially fatal infections in humans. This study focuses on the response regulator AgrA, which plays a crucial role in regulating the production of virulence factors in S. aureus. The objective is to identify natural compounds that can inhibit the binding of AgrA to its promoter site, thus inhibiting the expression of virulence genes. To achieve this, a pharmacophore model was generated using known drugs and applied to screen the ZINC natural product database. The resulting compounds were subjected to molecular docking-based virtual screening against the C-terminal DNA binding domain of AgrA. Three compounds, namely ZINC000077269178, ZINC000051012304, and ZINC000004266026, were shortlisted based on their strong affinity for key residues involved in DNA binding and transcription initiation. Subsequently, the unbound and ligand-bound complexes were subjected to a 200 ns molecular dynamics simulation to assess their conformational stability. Various analyses, including RMSD, RMSF, Rg, SASA, Principal Component Analysis, and Gibbs free energy landscape, were conducted on the simulation trajectory. The RMSD profile indicated similar fluctuations in both bound and unbound structures, while the Rg profile demonstrated the compactness of the protein without any unfolding during the simulation. Furthermore, Principal component analysis revealed that ligand binding reduced the overall atomic motion of the protein whereas free energy landscape suggested the energy variations obtained in complexes.Communicated by Ramaswamy H. Sarma.

Synergistic Effect Mechanism of Binary Sweet Taste Compounds

In our previous study, phloridzin, sucrose, l-alanine, and dulcitol presented synergistic effects in Camellia nanchuanica black tea (NCBT). This study aims to verify the synergistic effects of the aforementioned sweet taste compounds and the mechanism involved. By conducting σ-τ plot analysis, phloridzin at the recognition threshold concentration (phl) exhibited synergistic effects with different concentrations of sucrose (Lsuc-6suc). Various concentrations of sucrose, phloridzin, and their combinations were selected to investigate the impact on sweet taste receptor cells. The results revealed that sucrose/phloridzin significantly increased the calcium signal compared to phloridzin and sucrose alone, attributed to the greater stability of the sucrose/phloridzin combination when binding to Taste 1 Receptor Member 3 (TAS1R3; one subunit of sweet taste receptor proteins). Ultimately, the sweet taste signal of sucrose/phloridzin was transmitted to the brain, triggering the activation of more brain regions associated with sweet taste perception (right insular, postcentral, and amygdala).

Mutations in LIFR rewire the JAK/STAT signaling pathway: A study unveiling mechanistic details of Stüve-Wiedemann syndrome

Stüve-Wiedemann syndrome (SWS), a rare autosomal recessive disorder, characterized by diminutive size, curvature of the elongated bones, bent fingers, episodes of heightened body temperature, respiratory distress or periods of breath-holding, and challenges with feeding, especially causes fatality in infants. SWS is an outcome of potential missense mutations in the leukemia inhibitory factor receptor gene reflected as numerous amino acid mutations at protein level. Employing in silico tools and techniques like mutational screening with Pred_MutHTP, I-Mutant2.0, PANTHER.db, PolyPhen, to classify mutations as deleterious/destabilizing, in conjunction with experimental data analysis, P136A and S279P emerged as 'effect'-causing mutations. Pre-existing knowledge suggests, SWS progression is effectuated conformationally altered and dysfunctional LIFR, unable to bind to LIF and further form the LIF/LIFR/gp130 signalling complex. To gain functional insights into the effect of the said mutations on the wild type protein, an all-atom, explicit, solvent molecular dynamics simulation was performed following docking approaches. Consequently, referring to the RMSD, RMSF, protein dynamic network analysis, energy landscape plots and domain motion analysis, it was revealed that unbound LIFR_WT was more prone to LIF binding as usual whereas the mutants exhibited considerable domain closure to inhibit LIF binding. We conducted binding affinity analysis via MM/GBSA and dissociation constant estimation after LIFR-LIF docking and found the WT_complex to be more stable and compact as a whole when compared to the flexible mutant complexes thus being associated with SWS. Our study offers a route for understanding molecular level implications upon LIFR mutations which opens an avenue for therapeutic interventions.

Postharvest ripening of newly harvested corn: Weakened interactions between starch and protein

The effects of postharvest ripening of corn on the mechanisms of starch and protein interactions were investigated using molecular dynamics and several chemical substances. Sodium dodecyl sulfate (SDS) treatment all significantly affected the starch content, molecular weight of proteins, relative crystallinity, pasting characteristics and dynamic viscoelasticity in samples before and after postharvest ripening. In the corn that had not undergone postharvest ripening, there were also significant electrostatic interactions and hydrogen bonds between starch and protein. In addition, molecular dynamics had demonstrated that the forces between starch and protein in corn were mainly hydrophobic interactions, electrostatic interaction, and hydrogen bonds. Compared with zein, corn glutelin was more tightly bound to starch. The binding energy of starch to both proteins was reduced in after postharvest-ripened corn. This study laid a rationale for investigating the change mechanism of corn postharvest ripening quality and improving processing property and edible quality of corn.

Synthesis of 2-aminothiazole sulfonamides as potent biological agents: Synthesis, structural investigations and docking studies

A simplified synthetic approach involving sulfonylation followed by amino group alkylation produced new 2-aminothiazole derivatives. UV/Vis, infrared, and NMR spectroscopies confirmed their structures. Compounds 36, 22, 34, and 35 showed strong inhibition against Jack bean and Bacillus Pasteurii urease, with IC50 values from 14.06 to 20.21 μM/mL. Compounds 20, 26, 21, 29, 30, 31, and 32 exhibited potent inhibitory effects against α-glucosidase and α-amylase, with IC50 values between 20.34 and 37.20 μM/mL. Compounds 33, 26, and 27 demonstrated potent DPPH scavenging, with IC50 values around 34.4-39.2 μM/mL. FMO analysis showed compounds 21, 22, 24, and 25 having parallel aromatic ring systems due to π cloud interactions, while compounds 32 and 38 had distinct electronic density distributions. Compound 22 had HOMO and LUMO energy gaps of 5.805 eV, with bromo and fluoro substitutions in compounds 21 and 24 slightly increasing the gaps to 6.089 eV and 6.078 eV, respectively. Nitro groups in compounds 25 and 32 reduced the gaps to 0.384 eV and 1.187 eV. All compounds demonstrated high gastrointestinal absorption, non-permeability to the blood-brain barrier, and optimal skin permeation (Log Kp between -5.83 and -6.54 cm/s). Compounds 22, 24, and 38 had promising QED scores of 0.719, 0.707, and 0.860, respectively, with synthetic accessibility scores from 2.057 to 2.517. ADMET predictions indicated minimal toxicity, cardiovascular safety, and significant inhibitory potential for CYP enzymes. Strong in silico binding affinities (binding energies -5.75 to -7.63 kcal/mol) and metabolic stability suggest these derivatives are promising candidates for further drug development.

Structural modifications and kinetic effects of KRAS interactions with HRAS and NRAS: an in silico comparative analysis of KRAS mutants

The RAS genes which code for KRAS, HRAS, and NRAS are three of the most frequently mutated oncogenes responsible for cancer deaths. Tumorigenesis is one of the most significant outcomes of deregulation of RAS GTPases. Although the structures have been extensively studied, there is still more to be discovered about the actual binding conformations of the three isoforms, especially when mutated, to design an inhibitory drug. Recent studies have identified important interactions between the three isoforms that affect the oncogenic strength of the others when they are mutated. In this study, we utilize molecular dynamics simulations to examine the modifications of the structural property, mechanism, and kinetic energy of KRAS when interacting individually and with HRAS and NRAS. Notably, we found that WT-KRAS' orientation when bound to WT-HRAS vs. WT-NRAS is rotated 180°, with mutants demonstrating a similar binding pattern. The binding sites of the isoforms with KRAS share similarities with those involved in the GDP/GTP active site and site of KRAS dimerization. Thus, the isoform interaction can serve as an inhibitory method of KRAS actions. This study advances the understanding of inhibiting RAS-driven cancers through a novel isoform interaction approach only recently discovered, which has been proven to be an effective alternate therapeutic approach. We developed a blueprint of the interaction which would be beneficial in the development of KRAS mutant-specific and pan-KRAS mutant inhibitory drugs that mimic the isoform interactions. Our results support the direct interaction inhibition mechanism of mutant KRAS when bound to WT-HRAS and WT-NRAS by the isoforms' hypervariable region binding to the G-domain of KRAS. Furthermore, our results support the approach of reducing the effects of oncogenic KRAS by altering the concentration of the isoforms or a drug alternative based on the overall structural and kinetic stability, as well as the binding strength of the mutant-isoform complexes.

Optimizing expression, purification, structural and functional assessments of a novel dimeric incretin (GLP-1cpGLP-1)

Glucagon-like peptide-1 (GLP-1) is an incretin hormone that reduces postprandial glycemic excursions by enhancing insulin secretion. In this study, a new dimeric GLP-1 analogue (GLP-1cpGLP-1) was designed by inserting human insulin C-peptide (CP) in the middle of a dimer of [Gly8] GLP-1 (7-36). Then, the dimeric incretin (GLP-1cpGLP-1) was ligated to human αB-crystallin (αB-Cry) to create a hybrid protein, abbreviated as αB-GLP-1cpGLP-1. The constructed gene was well expressed in the bacterial host system. After specific chemical release from the hybrid protein, the dimeric incretin was purified by size exclusion chromatography (SEC). Finally, the RP-HPLC analysis indicated a purity of >99 % for the dimeric incretin. The secondary structure assessments by various spectroscopic methods, and in silico analysis suggested that the dimeric incretin has α-helical rich structure. The dynamic light scattering (DLS) analysis indicates that our dimeric incretin forms large oligomeric structures. This incretin analogue significantly reduced blood glucose levels in both healthy and diabetic mice while effectively triggering insulin release. The size exclusion HPLC also indicates the interaction of the new incretin analogue with human serum albumin, the main carrier protein in the bloodstream. Consistent with the results obtained from the biological activity assessments, this significant interaction indicates its potential as a viable therapeutic agent with a long-lasting effect. The results of our research represent a significant breakthrough in the successful design of an active incretin dimer capable of effectively controlling blood sugar levels and inducing insulin secretion in the realm of diabetes treatment.

Unraveling the binding behavior of the antifouling biocide 4,5-dichloro-2-n-octyl-4-isothiazolin-3-one with human serum albumin: Multi-spectroscopic, atomic force microscope, computational simulation, and esterase activity

As a marine antifouling biocide, 4,5-dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT) exhibited high toxicity to marine organisms. This study investigated the interaction between DCOIT and human serum albumin (HSA) using several spectroscopic techniques combined with computer prediction methods. The UV-vis absorption spectra, Stern-Volmer constant (KSV) and fluorescence resonance energy transfer (FRET) results indicated that DCOIT caused static quenching of HSA fluorescence. The ΔG°, ΔH° and ΔS° values were -31.03 ± 0.17 kJ·mol-1, -133.54 ± 0.88 kJ·mol-1 and -348.46 ± 2.86 J.mol-1·K-1, respectively, suggesting that van der Waals forces and hydrogen bonds governed the spontaneous formation of the complex. Synchronous fluorescence and circular dichroism (CD) spectroscopy observed the burial of Trp residues within HSA and the unfolding of HSA secondary structure induced by DCOIT. Three-dimensional (3D) fluorescence and Atomic Force Microscopy (AFM) further detected DCOIT-induced loosening of HSA peptide chain structure. Site displacement experiments indicated that DCOIT binding at site I of HSA. Computational predictions indicated that hydrophobic interactions were also essential in the complex. The increased RMSD, Rg, SASA, and RMSF confirmed that DCOIT weakened the stability and compactness of HSA, rendering residues more flexible. Lastly, esterase activity assays demonstrated that DCOIT inhibited esterase activity and interfered with the human detoxification process.

Unveiling potential PET degrading eukaryotes through in silico bioprospecting of PETases

This study addresses the environmental problem of PET plastic through in silico bioprospecting for the identification and experimental validation of novel PET degrading eukaryotes through the in silico bioprospectingI of PETases, employing a methodology that combines Hidden Markov Models (HMMs), clustering techniques, molecular docking, and dynamic simulations. A total of 424 putative PETase sequences were identified from 219 eukaryotic organisms, highlighting six sequences with low affinity energies. The Aspergillus luchuensis sequence showed the lowest Gibbs free energy and exhibited stability at different temperatures in molecular dynamics assays. Experimental validation, through a plate clearance assay and HPLC, confirmed PETase activity in three wild-type fungal strains, with A. luchuensis showing the highest efficiency. The results obtained demonstrate the effectiveness of combining computational and experimental approaches as proof of concept to discover and validate eukaryotes with PET-degrading capabilities opening new perspectives for the sustainable management of this type of waste and contributing to its environmental mitigation.

Potent anti-angiogenic component in Kaempferia galanga L. and its mechanism of action

ETHNOPHARMACOLOGICAL RELEVANCE

Traditionally, the roots of Kaempferia galanga has been used to treat high blood pressure, chest pain, headache, toothache, rheumatism, indigestion, cough, inflammation and cancer in Asia. Nevertheless, most of its pharmacological studies were focused on ethanolic extracts and volatile oils. The exact active chemical constituents and their underlying mechanisms are still poorly understood, especially towards its anti-cancer treatment. Inhibition of angiogenesis is an important atrategy to inhibit tumor growth. It has been reported that the low polar component of the plant possessed anti-angiogenic activity. Yet, the potent compound which is responsible for the effect and its molecular mechanism has not been reported.

AIM OF THE STUDY

To determine the potent anti-angiogenic component in K.galanga and its mechanism of action.

MATERIAL AND METHODS

The low polar components of the plant were concentrated using the methods of supercritical fluid extraction (SFE), subcritical extraction (SCE) and steam distillation (SD). The anti-angiogenic activity of the three extracts was evaluated using a zebrafish model. The content of the active compound in those extracts was determined with HPLC analysis. The in-vitro and in-vivo activity of the isolated compound was evaluated using human umbilical vein endothelial cells (HUVECs) model, the aortic ring assay and the matrigel plug assay, respectively. Its molecular mechanism was further studied by the western blotting assay and computer-docking experiments. Besides, its cytotoxicity on cancer and normal cell lines was evaluated using the cell-counting kit.

RESULTS

HPLC results showed that trans-ethyl p-methoxycinnamate (TEM) was the major component of the extracts. The extract of SFE showed the best effect as it has the highest content of TEM. TEM could inhibit vascular endothelial growth factor (VEGF)-induced viability, migration, invasion and tube formation in human umbilical vein endothelial cells (HUVECs) in vitro. Moreover, it inhibited VEGF-induced sprout formation ex vivo and vessel formation in vivo. Mechanistic study showed that it could suppress tyrosine kinase activity of the receptor of VEGF (VEGFR2) and alter its downstream signaling pathways. In addition, the molecular docking showed that the binding of TEM and VEGFR2 is stable, which mainly attributed to the non-covalent binding interaction. Beside, TEM possessed little toxicity to both cancer and normal cells.

CONCLUSION

TEM is the major anti-angiogenic component present in K. galanga and its anti-angiogenic property rather than toxicity provides scientific basis for the traditional use of K. galanga in cancer treatment.

Identification of phytocompounds as potent inhibitors of sodium/glucose cotransporter-2 leading to diabetes treatment

Type-II diabetes, a major metabolic disorder has threatened the very existence of a healthy life since long ago. Commercially available antidiabetic drugs are known for several adverse effects. The present study attempted to identify potential phytocompounds as inhibitors of sodium/glucose cotransporter-2 (SGLT2), a major protein that helps in glucose re-absorption from renal tubules. A total of 28 phytocompounds were collected based on the literature survey. 3D co-ordinates of phytocompounds were collected from PubChem database. Molecular docking was carried out with SGLT2 protein and the best 3 docking complexes were subjected to molecular dynamics simulation for 100 ns. Free energy changes were also analyzed using MM/PBSA analysis. Phytocompounds were also analyzed for their drug-likeness and ADMET properties. Docking study observed a strong binding affinity of phytocompounds (> -7.0 kcal/mol). More than 10 phytocompounds showed better binding affinity compared to reference drugs. Further analysis of three best docking complexes when analyzed by MD simulation showed better stability and compactness of the complexes compared to reference drug, empagliflozin. MM/PBSA analysis also revealed that van der Waals force and electrostatic energy are the major binding energy involved in the complex formation. Like docking energy, free energy analysis also observed stronger binding energies (ΔGGAS) in SGLT2-phytocompound complexes compared to empagliflozin complex. All the phytocompounds showed drug-likeness and considerable ADMET properties. The study, therefore, suggests that Trifolirhizin-6'-monoacetate, Aspalathin, and Quercetin-3-glucoside could be a possible inhibitor of SGLT2 protein. However, further studies need to be carried out to reveal the exact mode of activity.Communicated by Ramaswamy H. Sarma.

Prediction of the active compounds and mechanism of Biochanin A in the treatment of Legg-Calvé-Perthes disease based on network pharmacology and molecular docking

BACKGROUND

Legg-Calvé-Perthes disease is a special self-limited disease in pediatric orthopedics with a high disability rate and a long-term course, and there is still no clear and effective therapeutic drug in clinic. This study aimed to investigate the potential efficacy of biochanin A, a kind of oxygen-methylated isoflavone compound, in treating Perthes disease based on network pharmacology, molecular docking and in vitro experiments.

METHODS

IL-6 was used to stimulate human umbilical vein endothelial cells to construct endothelial cell dysfunction model. We demonstrated whether biochanin A could alleviate endothelial dysfunction through CCK8 assay, immunofluorescence. Targets of biochanin A from pharmMappeer, SWISS, and TargetNet databases were screened. Targets of endothelial dysfunction were obtained from Genecards and OMIM databases. Protein-protein interaction, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomics analyses were used to analyze the potential target and the key pathway of the anti-endothelial dysfunction activity of biochanin A. To validate the potential target-drug interactions, molecular docking and molecular dynamics simulations were performed and the result was proved by western blot.

RESULTS

It was found that biochanin A can promote the expression of ZO-1, reduce the expression of ICAM-1, which means improving endothelial dysfunction. A total of 585 targets of biochanin A from pharmMappeer, SWISS, and TargetNet databases were screened. A total of 10,832 targets of endothelial dysfunction were obtained from Genecards and OMIM databases. A total of 527 overlapping targets of endothelial dysfunction and biochanin A were obtained. AKT1, TNF-α, VCAM1, ICAM1, and NOS3 might be the key targets of the anti-endothelial dysfunction activity of biochanin A, and the key pathways might be PI3K-Akt and TNF signaling pathways. Molecular docking results indicated that the AKT1 and TNF-α had the highest affinity binding with biochanin A.

CONCLUSION

This study indicates that biochanin A can target AKT1 and TNF-α to alleviate endothelial dysfunction induced by IL-6 in Perthes disease, which provides a theoretical basis for the treatment of Perthes disease by using biochanin A.

An in silico analysis of the interaction of marine sponge-derived bioactive compounds with type 2 diabetes mellitus targets DPP-4 and PTP1B

Type 2 diabetes is a medical condition involving elevated blood glucose levels resulting from impaired or improper insulin utilization. As the number of type 2 diabetes cases increases each year, there is an urgent need to develop novel drugs having new targets and/or complementing existing therapeutic protocols. In this regard, marine sponge-derived compounds hold great potential due to their potent biological activity and structural diversity. In this study, a small library of 50 marine sponge-derived compounds were examined for their activity towards type 2 diabetes targets, namely dipeptidyl peptidase-4 (DPP-4) and protein tyrosine phosphatase 1B (PTP1B). The compounds were first subjected to molecular docking on protein models based on their respective co-crystal structures to assess binding free energies (BFE) and conformations. Clustering analysis yielded BFE that ranged from 24.54 kcal/mol to -9.97 kcal/mol for DPP-4, and from -4.98 kcal/mol to -8.67 kcal/mol for PTP1B. Interaction analysis on the top ten compounds with the most negative BFE towards each protein target showed similar intermolecular interactions and key interacting residues as in the previously solved co-crystal structure. These compounds were subjected to absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling to characterize drug-likeness and combining the results from these analyses, (S)-6'-debromohamacanthin B was identified as a potential multi-target inhibitor of DPP-4 and PTP1B, having favorable protein interaction, no Lipinski violations, good gastrointestinal (GI) tract absorption, blood-brain barrier (BBB) penetration, and no predicted toxicity. Finally, the interaction of (S)-6'-debromohamacanthin B with the two proteins was validated using molecular dynamics simulations over 100 ns through RMSD, radius of gyration, PCA, and molecular mechanics Poisson-Boltzmann surface area (MMPBSA) confirming favorable interactions with the respective proteins.Communicated by Ramaswamy H. Sarma.

Molecular Design of Novel Inhibitor by Targeting IL-6Rα using Combined Pharmacophore and Experimentally Verified Plant Products with Scaffold-Hopping Techniques: A Dual Therapeutic Strategy for COVID-19 and Cancer

The IL-6/IL-6R/gp130 complex serves as a significant indicator of cytokine release syndrome in COVID-19 and chronic inflammation, increasing the risk of cancer. Therefore, we identified IL-6Rα as a potential target to block gp130 interaction. Notably, there has been no reception of approval for an orally available drug to serve this purpose, to date. In this study, we targeted IL-6Rα to inhibit IL-6Rα/gp130 interaction. The selection of the lead candidate L821 involved the amalgamation of three drug discovery approaches. This library was screened employing tertiary structure-based pharmacophore models followed by molecular docking models, scaffold-hopping, MM/PBSA as well as MM/GBSA analysis, and assessments of pKi and ADMET properties. After evaluating the binding interactions with key amino acids, 15 potential ligands were chosen, with the top ligand undergoing further investigation by means of molecular dynamics simulations. Considering the stability of the complexes, the strong interactions observed between ligand and residues of IL-6Rα/gp130, and the favorable binding free energy calculations, L821 emerged as the prime candidate for inhibiting IL-6Rα. Notably, L821 exhibited a docking-based binding affinity of -9.5 kcal/mol. Our study presents L821 as a promising inhibitor for future in vitro analysis, potentially combatting SARS-CoV-2-related cytokine storms and serving as an oncogenic drug therapy.

Polyphenolic flavonoid compounds act as the inhibitory potential of aggregation process: Implications for the prevention and therapeutics against FALS-associated D101G SOD1 mutant

Aggregation of proteins is a biological phenomenon caused by misfolded proteins. Human superoxide dismutase (hSOD1) misfolding and aggregation underlie the neurological illness amyotrophic lateral sclerosis (ALS). The most significant contributing factor to ALS is genetic point mutations in SOD1. particularly, D101G mutant is the most harmful because it significantly reduces the life expectancy of patients. Subsequently, the use of natural polyphenolic flavonoids is strongly recommended to reduce the amyloidogenic behavior of protopathic proteins. In this study, using computational parameters such as protein-ligand interaction and molecular dynamics (MD) simulation analyses, we are trying to identify a pharmacodynamically promising flavonoid compound that can effectively inhibit the pathogenic behavior of the D101G mutant. Epigallocatechin-gallate (EGCG), Hesperidin, Isorhamnetin, and Diosmetin were identified as potential leads in a preliminary screening of flavonoids to anti-amyloid action. The results of MD showed that the binding of flavonoids to D101G mutant caused changes in stability, hydrophobicity of protein, and flexibility, as well as significantly led to the restoration of lost hydrogen bonds. Secondary structure analysis showed that protein destabilization and the increased propensity of β-sheet caused by the mutation were restored to the wild-type state upon binding of flavonoids. Besides, to differentiate aggregation, we elucidated alterations in the free energy landscape (FEL) and dynamic cross-correlation matrix (DCCM) of WT-SOD1 and mutant (unbound /bound) states. Among flavonoids, Epigallocatechin-gallate and Hesperidin had the most therapeutic efficacy against the D101G mutant. Therefore, Epigallocatechin-gallate and Hesperidin promise considerable therapeutic potential to develop highly effective inhibitors in reducing fatal and irreversible ALS.

Computational analysis of the functional and structural impact of the most deleterious missense mutations in the human Protein C

Protein C (PC) is a vitamin K-dependent factor that plays a crucial role in controlling anticoagulant processes and acts as a cytoprotective agent to promote cell survival. Several mutations in human PC are associated with decreased protein production or altered protein structure, resulting in PC deficiency. In this study, we conducted a comprehensive analysis of nonsynonymous single nucleotide polymorphisms in human PC to prioritize and confirm the most high-risk mutations predicted to cause disease. Of the 340 missense mutations obtained from the NCBI database, only 26 were classified as high-risk mutations using various bioinformatic tools. Among these, we identified that 12 mutations reduced the stability of protein, and thereby had the greatest potential to disturb protein structure and function. Molecular dynamics simulations revealed moderate alterations in the structural stability, flexibility, and secondary structural organization of the serine protease domain of human PC for five missense mutations (L305R, W342C, G403R, V420E, and W444C) when compared to the native structure that could maybe influence its interaction with other molecules. Protein-protein interaction analyses demonstrated that the occurrence of these five mutations can affect the regular interaction between PC and activated factor V. Therefore, our findings assume that these mutants can be used in the identification and development of therapeutics for diseases associated with PC dysfunction, although assessment the effect of these mutations need to be proofed in in-vitro.

Fighting against amyotrophic lateral sclerosis (ALS) with flavonoids: a computational approach to inhibit superoxide dismutase (SOD1) mutant aggregation

Protein aggregation is a biological process that occurs when proteins misfold. Misfolding and aggregation of human superoxide dismutase (hSOD1) cause a neurodegenerative disease called amyotrophic lateral sclerosis (ALS). Among the mutations occurring, targeting the E21K mutation could be a good choice to understand the pathological mechanism of SOD1 in ALS, whereof it significantly reduces life hopefulness in patients. Naturally occurring polyphenolic flavonoids have been suggested as a way to alleviate the amyloidogenic behavior of proteins. In this study, computational tools were used to identify promising flavonoid compounds that effectively inhibit the pathogenic behavior of the E21K mutant. Initial screening identified Pelargonidin, Curcumin, and Silybin as promising leads. Molecular dynamics (MD) simulations showed that the binding of flavonoids to the mutated SOD1 caused changes in the protein stability, hydrophobicity, flexibility, and restoration of lost hydrogen bonds. Secondary structure analysis indicated that the protein destabilization and the increased propensity of β-sheet caused by the mutation were restored to the wild-type state upon binding of flavonoids. Free energy landscape (FEL) analysis was also used to differentiate aggregation, and results showed that Silybin followed by Pelargonidin had the most therapeutic efficacy against the E21K mutant SOD1. Therefore, these flavonoids hold great potential as highly effective inhibitors in mitigating ALS's fatal and insuperable effects.Communicated by Ramaswamy H. Sarma.

In silico engineering and simulation of RNA interferences nanoplatforms for osteoporosis treating and bone healing promoting

Osteoporosis is a bone condition characterized by reduced bone mineral density (BMD), poor bone microarchitecture/mineralization, and/or diminished bone strength. This asymptomatic disorder typically goes untreated until it presents as a low-trauma fracture of the hip, spine, proximal humerus, pelvis, and/or wrist, requiring surgery. Utilizing RNA interference (RNAi) may be accomplished in a number of ways, one of which is by the use of very tiny RNA molecules called microRNAs (miRNAs) and small interfering RNAs (siRNAs). Several kinds of antagomirs and siRNAs are now being developed to prevent the detrimental effects of miRNAs. The goal of this study is to find new antagonists for miRNAs and siRNAs that target multiple genes in order to reduce osteoporosis and promote bone repair. Also, choosing the optimum nanocarriers to deliver these RNAis appropriately to the body could lighten up the research road. In this context, we employed gene ontology analysis to search across multiple datasets. Following data analysis, a systems biology approach was used to process it. A molecular dynamics (MD) simulation was used to explore the possibility of incorporating the suggested siRNAs and miRNA antagonists into polymeric bioresponsive nanocarriers for delivery purposes. Among the three nanocarriers tested [polyethylene glycol (PEG), polyethylenimine (PEI), and PEG-PEI copolymer], MD simulations show that the integration of PEG-PEI with has-mIR-146a-5p is the most stable (total energy = -372.84 kJ/mol, Gyration radius = 2.1084 nm), whereas PEI is an appropriate delivery carrier for has-mIR-7155. The findings of the systems biology and MD simulations indicate that the proposed RNAis might be given through bioresponsive nanocarriers to accelerate bone repair and osteoporosis treatment.

Current stage and future development of Belgrade collisional and radiative databases/datasets of importance for molecular dynamics

Atomic and molecular (A&M) databases that contain information about species, their identities and radiative/collisional processes are essential and helpful tools that are utilized in many fields of physics, chemistry, and chem/phys-informatics. Errors or inconsistencies in the datasets are a serious issue since they can lead to inaccurate predictions and generate problems with the modeling. This demonstrates that data curation efforts around A&M databases are still indispensable and that in the curation process studious attention is required. Therefore, we herein present research activities around Belgrade "nodes" - datasets of collision/radiative cross-sections and rates needed for spectroscopy analysis in various A&M, optical and plasma physics fields. Methodologies of our research and both present and future aspects of the applications are explained. We explored the possibility to extend our nodes towards building a new database on Judd-Ofelt parameters by using machine learning in order to predict optical properties of luminescence materials. In addition, we hope that public availability of our datasets and their graphical representations will also motivate others to investigate the potential of these data.

Deciphering the Molecular Mechanism of Escin against Neuropathic Pain: A Network Pharmacology Study

Background

Escin is the main active component in Aesculus hippocastanum. It has been demonstrated that escin has anti-inflammatory properties. This study combined the methods of network pharmacology, molecular docking, and molecular dynamics to explore the molecular mechanism of escin against neuropathic pain (NP).

Methods

The Swiss Target Prediction and the Pharm Mapper database were employed for predicting the targets of escin. Also, the candidate targets of NP were gathered via the databases including Therapeutic Targets, DisGeNet, GeneCards, DrugBank, and OMIM. Subsequently, the network of protein-protein interaction was screened for the key targets by the software Cytoscape 3.8.0. Then, the intersection of these targets was analysed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Additionally, we further investigated the ligand-target interactions by molecular docking and molecular dynamics simulations.

Results

In total, 94 escin targets were predicted by network pharmacology. Among them, SRC, MMP9, PTGS2, and MAPK1 were the core candidate targets. Subsequently, the analysis of GO and KEGG enrichment revealed that escin affected NP by regulating protein kinase C, MAP kinase, TRP channels, the T-cell receptors signaling pathway, and the TNF signaling pathway. The results of molecular docking and molecular dynamics simulation confirmed that escin not only had a strong binding activity with the four core target proteins but also stably combined in 50 ns.

Conclusions

Our study revealed that escin acts on the core targets SRC, MMP9, PTGS2, MAPK1, and associated enrichment pathways to alleviate neuronal inflammation and regulate the immune response, thus exerting anti-NP efficacy. This study provided innovative ideas and methods for the promising treatment of escin in relieving NP.

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.

Association of IL-17F rs2397084 (E126G), rs11465553 (V155I) and rs763780 (H161R) variants with rheumatoid arthritis and their effects on the stability of protein

Interleukin-17F (IL-17F), considered a pro-inflammatory cytokine, has been shown to contribute to skeletal tissue degradation and hence chronic inflammation in rheumatoid arthritis (RA). In this study we utilized bioinformatics tools to analyze the effect of three exonic SNPs (rs2397084, rs11465553, and rs763780) on the structure and function of the IL-17F gene, and evaluated their association with RA in Pakistani patients. The predicted deleterious and damaging effects of identified genetic variants were assessed through the utilization of multiple bioinformatics tools including PROVEAN, SNP&GO, SIFT, and PolyPhen2. Structural and functional effects of these variants on protein structures were evaluated through the use of additional tools such as I-Mutant, MutPred, and ConSurf. Three-dimensional (3D) models of both the wild-type and mutant proteins were constructed through the utilization of I-TASSER software, with subsequent structural comparisons between the models conducted through the use of the TM-align score. A total of 500 individuals, 250 cases and 250 controls, were genotyped through Tri-ARMS-PCR method and the resultant data was statistically analyzed using various inheritance models. Our bioinformatics analysis showed significant structural differences for wild type and mutant protein (TM-scores and RMSD values were 0.85934 and 2.34 for rs2397084 (E126G), 0.87388 and 2.49 for rs11465553 (V155I), and 0.86572 and 0.86572 for rs763780 (H161R) with decrease stability for the later. Overall, these tools enabled us to predict that these variants are crucial in causing disease phenotypes. We further tested each of these single nucleotide variants for their association with RA. Our analysis revealed a strong positive association between the genetic variant rs763780 and the risk of developing rheumatoid arthritis (RA) at both the genotypic and allelic levels. The genotypic association was statistically significant[χ2 = 111.8; P value <0.0001], as was the allelic level [OR 3.444 (2.539-4.672); P value 0.0008]. These findings suggest that the presence of this genetic variant may increase the susceptibility to RA. Similarly, we observed a significant distribution of the genetic variant rs11465553 at the genotypic level [χ2 = 25.24; P value = 0.0001]. However, this variant did not show a significant association with RA at the allelic level [OR = 1.194 (0.930-1.531); P value = 0.183]. However, the distribution of variant rs2397084 was more or less random across our sample with no significant association either at genotypic and or allelic level. Put together, our association study and in silico prediction of decreasing of IL17-F protein stabilty confirmed that two SNPs, rs11465553 and rs763780 are crucial to the suscetibility of and showed that these RA in Pakistani patients.

In silico prediction of mutation sites for anthranilate synthase from Serratia marcesens to deregulate tryptophan feedback inhibition

Allosteric feedback inhibition of the committed step in amino acid biosynthetic pathways is a major concern for production of amino acids at industrial scale. Anthranilate synthase (AS) catalyzes the first reaction of tryptophan biosynthetic pathway found in microorganisms and is feedback inhibited by its own product i.e. tryptophan. Here, we identified new mutant sites in AS using computational mutagenesis approach. MD simulations (20 ns) followed by MMPBSA and per residue decomposition energy analysis identified seven amino acid residues with best binding affinity for tryptophan. All 19 mutant structures were generated for each identified amino acid residue followed by simulation to evaluate effect of mutation on protein stability. Later, molecular docking studies were employed to generate mutant-tryptophan complex and structures with binding energies (kcal/mol) much higher than wild-type AS were selected. Finally, two mutants i.e., S37W and S37H were identified on the basis of positive binding scores and loss of tryptophan binding inside pocket. Further, MD simulations run for 200 ns were performed over these mutant-tryptophan complexes followed by RMSD, RMSF, radius of gyration , solvent accessible surface area , intra-protein hydrogen bond numbers, principal component analysis, free energy landscape (FEL) and secondary structure analysis to rationale effect of mutations on stability of protein. Cross correlation analysis of mutant site amino acids (S37W) with key residues of catalytic site (G325, T326, H395 and G482) was done to evaluate the effect of mutations on catalytic site conformation. Current computational mutagenesis approach predicted two mutants S37W and S37H with proposed deregulated feedback inhibition by tryptophan and retained catalytic activity.Communicated by Ramaswamy H. Sarma.

Antibacterial Potency of an Active Compound from Sansevieria trifasciata Prain: An Integrated In Vitro and In Silico Study

Sansevieria trifasciata Prain holds great potential as a valuable asset in pharmaceutical development. In this study, our focus is to explore and assess the antibacterial activity of various components derived from this plant, including extracts, fractions, subfractions, and isolates, explicitly targeting two common bacteria: Escherichia coli and Streptococcus aureus. The isolated compound, identified as a derivative pyridone alkaloid (5-methyl-11-(2-oxopyridin-1(2H)-yl)undecaneperoxoicacid), demonstrates notable antibacterial effects. The extracts, fractions, subfractions, and isolates reveal significant bacterial growth reductions (p < 0.05). The minimum inhibitory concentration (MIC) values for Escherichia coli were 1.95 ppm, 3.9 ppm, 15.62 ppm, and 7.81 ppm, respectively, while the MIC values for Streptococcus aureus were 1.95 ppm, 1.95 ppm, 15.62 ppm, and 7.81 ppm, respectively. Computational analysis showed the isolates' interaction with key residues on the active site of β-ketoacyl-ACP synthase from Escherichia coli and TyrRS from Streptococcus aureus. The findings indicate that the isolates exhibit a strong affinity for specific residues, including His333, Cys163, and Phe392 in β-ketoacyl-ACP synthase, as well as Arg88, His117, Glu160, and Gln213 in TyrRS. Comparative energy calculations using MMPBSA demonstrate the isolates' favorable binding energy (-104,101 kJ/mol for β-ketoacyl-ACP synthase and -81,060 kJ/mol for TyrRS) compared to ciprofloxacin. The elucidated antibacterial activity and molecular interactions of the isolates present valuable knowledge for future in vitro studies, facilitating the development of novel antibacterial agents targeting diverse bacterial strains.

Concatenation of molecular docking and dynamics simulation of human papillomavirus type 16 E7 oncoprotein targeted ligands: In quest of cervical cancer's treatment

The Human papillomaviruses type 16 E7 oncoprotein is a 98-amino-acid, 11-kilodalton acidic oncoprotein with three conserved portions. Due to its interaction with the pRb-E2F complex, CKII, CKI (mostly p21), and even HDAC1, it possesses strong transformative and carcinogenic qualities that inhibit normal differentiation and cell cycle regulation. Here, we target the E7 oncoprotein using two prior research active compounds: asarinin and thiazolo[3,2-a]benzimidazole-3(2H)-one,2-(2-fluorobenzylideno)-7,8-dimethyl (thiazolo), and valproic acid as a control. We are performing molecular docking followed by molecular dynamic analysis. By acting as competitive inhibitors in the binding site, it was hypothesized that both drugs would inhibit E7-mediated pRb degradation and E7-mediated p21 degradation, resulting in decreased cell cycle progression, immortalization, and proliferation. In addition, we expect that the direct inhibitory action of valproic acid in E7 will target the CKII-mediated phosphorylation pathway necessary for destabilizing p130 and pRb. According to the results of the dynamic simulation, stable interactions exist between every compound. Despite the instability of E7 protein, stability results indicate that both natural chemicals are preferable, with thiazolo outperforming valproic acid.

Probing the mechanism of interaction between capsaicin and myofibrillar proteins through multispectral, molecular docking, and molecular dynamics simulation methods

The interaction between myofibrillar proteins (MPs) and capsaicin (CAP) was investigated using multispectral, molecular docking, and molecular dynamics simulation methods. The resulting complex increased the hydrophobicity of the tryptophan and tyrosine microenvironment as revealed by fluorescence spectral analysis. The fluorescence burst mechanism study indicated that the fluorescence burst of CAP on the MPs was a static one (Kq = 1.386 × 10¹² m^-1s^-1) and that CAP could bind with MPs well (Ka = 3.31 × 10⁴ L/mol, n = 1.09). The analysis of circular dichroism demonstrated that the interaction between CAP and MPs caused a decrease in the α-helical structure of MPs. The complexes formed exhibited lower particle size and higher absolute ζ potential. Furthermore, hydrogen bonding, van der Waals forces, and hydrophobic interactions were found to be the primary factors facilitating the interaction between CAP and MPs, as suggested by molecular docking models and molecular dynamics simulations.

DSF inactivator RpfB homologous FadD upregulated in Bradyrhizobium japonicum under iron limiting conditions

Phytopathogenic bacteria Xanthomonas campestris pv. campestris (Xcc) causes black rot and other plant diseases. Xcc senses diffusible signal factor (DSF) as a quorum-sensing (QS) signal that mediates mainly iron uptake and virulence. RpfB deactivates DSF in this DSF-QS circuit. We examined differential gene expression profiles of Bradyrhizobium japonicum under low versus high iron conditions and found that fadD and irr were upregulated under low iron (log2 fold change 0.825 and 1.716, respectively). In addition to having similar protein folding patterns and functional domain similarities, FadD shared 58% sequence similarity with RpfB of Xcc. The RpfB-DSF and FadD-DSF complexes had SWISSDock molecular docking scores of - 8.88 kcal/mol and - 9.85 kcal/mol, respectively, and the 100 ns molecular dynamics simulation results were in accord with the docking results. However, significant differences were found between the binding energies of FadD-DSF and RpfB-DSF, indicating possible FadD-dependent DSF turnover. The protein-protein interaction network showed that FadD connected indirectly with ABC transporter permease (ABCtp), which was also upregulated (log2 fold change 5.485). We speculate that the low iron condition may be a mimetic environmental stimulus for fadD upregulation in B. japonicum to deactivate DSF, inhibit iron uptake and virulence of DSF-producing neighbors. This finding provides a new option of using B. japonicum or a genetically improved B. japonicum as a potential biocontrol agent against Xcc, with the added benefit of plant growth-promoting properties.

Ligand-Based Virtual Screening, Molecular Docking, and Molecular Dynamic Simulations of New β-Estrogen Receptor Activators with Potential for Pharmacological Obesity Treatment

Obesity is a pandemic and a serious health problem in developed and undeveloped countries. Activation of estrogen receptor beta (ERβ) has been shown to promote weight loss without modifying caloric intake, making it an attractive target for developing new drugs against obesity. This work aimed to predict new small molecules as potential ERβ activators. A ligand-based virtual screening of the ZINC15, PubChem, and Molport databases by substructure and similarity was carried out using the three-dimensional organization of known ligands as a reference. A molecular docking screening of FDA-approved drugs was also conducted as a repositioning strategy. Finally, selected compounds were evaluated by molecular dynamic simulations. Compounds 1 (-24.27 ± 0.34 kcal/mol), 2 (-23.33 ± 0.3 kcal/mol), and 6 (-29.55 ± 0.51 kcal/mol) showed the best stability on the active site in complex with ERβ with an RMSD < 3.3 Å. RMSF analysis showed that these compounds do not affect the fluctuation of the Cα of ERβ nor the compactness according to the radius of gyration. Finally, an in silico evaluation of ADMET showed they are safe molecules. These results suggest that new ERβ ligands could be promising molecules for obesity control.

Multi-Responsive Eight-State Bis(acridinium-Zn(II) porphyrin) Receptor

A multi-responsive receptor consisting of two (acridinium-Zn(II) porphyrin) conjugates has been designed. The binding constant between this receptor and a ditopic guest has been modulated (i) upon addition of nucleophiles converting acridinium moieties into the non-aromatic acridane derivatives and (ii) upon oxidation of the porphyrin units. A total of eight states has been probed for this receptor resulting from the cascade of the recognition and responsive events. Moreover, the acridinium/acridane conversion leads to a significant change of the photophysical properties, switching from electron to energy transfer processes. Interestingly, for the bis(acridinium-Zn(II) porphyrin) receptor, charge-transfer luminescence in the near-infrared has been observed.

A computational study of structural analysis of Class I human glucose-6-phosphate dehydrogenase (G6PD) variants: Elaborating the correlation to chronic non-spherocytic hemolytic anemia (CNSHA)

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common human enzyme defect that affects more than 500 million people worldwide. Individuals affected with G6PD deficiency may occasionally suffer mild-to-severe chronic hemolytic anemia. Chronic non-spherocytic hemolytic anemia (CNSHA) is a potential result of the Class I G6PD variants. This comparative computational study attempted to correct the defect in variants structure by docking the AG1 molecule to selected Class I G6PD variants [G6PD^Nashville (Arg393His), G6PD^Alhambra (Val394Leu), and G6PD^Durham (Lys238Arg)] at the dimer interface and structural NADP⁺ binding site. It was followed by an analysis of the enzyme conformations before and after binding to the AG1 molecule using the molecular dynamics simulation (MDS) approach, while the severity of CNSHA was determined via root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), hydrogen bonds, salt bridges, radius of gyration (Rg), solvent accessible surface area analysis (SASA), and principal component analysis (PCA). The results revealed that G6PD^Nashville (Arg393His) and G6PD^Durham (Lys238Arg) had lost the direct contact with structural NADP⁺ and salt bridges at Glu419 - Arg427 and Glu206 - Lys407 were disrupted in all selected variants. Furthermore, the AG1 molecule re-stabilized the enzyme structure by restoring the missing interactions. Bioinformatics approaches were also used to conduct a detailed structural analysis of the G6PD enzyme at a molecular level to understand the implications of these variants toward enzyme function. Our findings suggest that despite the lack of treatment for G6PDD to date, AG1 remains a novel molecule that promotes activation in a variety of G6PD variants.

Bi-terminal fusion of intrinsically-disordered mussel foot protein fragments boosts mechanical strength for protein fibers

Microbially-synthesized protein-based materials are attractive replacements for petroleum-derived synthetic polymers. However, the high molecular weight, high repetitiveness, and highly-biased amino acid composition of high-performance protein-based materials have restricted their production and widespread use. Here we present a general strategy for enhancing both strength and toughness of low-molecular-weight protein-based materials by fusing intrinsically-disordered mussel foot protein fragments to their termini, thereby promoting end-to-end protein-protein interactions. We demonstrate that fibers of a ~60 kDa bi-terminally fused amyloid-silk protein exhibit ultimate tensile strength up to 481 ± 31 MPa and toughness of 179 ± 39 MJ*m^-3, while achieving a high titer of 8.0 ± 0.70 g/L by bioreactor production. We show that bi-terminal fusion of Mfp5 fragments significantly enhances the alignment of β-nanocrystals, and intermolecular interactions are promoted by cation-π and π-π interactions between terminal fragments. Our approach highlights the advantage of self-interacting intrinsically-disordered proteins in enhancing material mechanical properties and can be applied to a wide range of protein-based materials.

In Silico Analyses of a Promising Drug Candidate for the Treatment of Amyotrophic Lateral Sclerosis Targeting Superoxide Dismutase I Protein

Amyotrophic lateral sclerosis (ALS) is the most prevalent motor neuron disorder in adults, which is associated with a highly disabling condition. To date, ALS remains incurable, and the only drugs approved by the FDA for its treatment confer a limited survival benefit. Recently, SOD1 binding ligand 1 (SBL-1) was shown to inhibit in vitro the oxidation of a critical residue for SOD1 aggregation, which is a central event in ALS-related neurodegeneration. In this work, we investigated the interactions between SOD1 wild-type and its most frequent variants, i.e., A4V (NP_000445.1:p.Ala5Val) and D90A (NP_000445.1:p.Asp91Val), with SBL-1 using molecular dynamics (MD) simulations. The pharmacokinetics and toxicological profile of SBL-1 were also characterized in silico. The MD results suggest that the complex SOD1-SBL-1 remains relatively stable and interacts within a close distance during the simulations. This analysis also suggests that the mechanism of action proposed by SBL-1 and its binding affinity to SOD1 may be preserved upon mutations A4V and D90A. The pharmacokinetics and toxicological assessments suggest that SBL-1 has drug-likeness characteristics with low toxicity. Our findings, therefore, suggested that SBL-1 may be a promising strategy to treat ALS based on an unprecedented mechanism, including for patients with these frequent mutations.

A comprehensive insight into the effects of caffeic acid (CA) on pepsin: Multi-spectroscopy and MD simulations methods

The interaction between caffeic acid (CA) and pepsin was investigated using multi-spectroscopy approaches and molecular dynamic simulations (MDS). The effects of CA on the structure, stability, and activity of pepsin were studied. Fluorescence emission spectra and UV-vis absorption peaks all represented the static quenching mechanism of pepsin by CA. Moreover, the fluorescence spectra displayed that the interaction of CA exposed the tryptophan chromophores of pepsin to a more hydrophilic micro-environment. Consistent with the simulation results, thermodynamic parameters revealed that CA was bound to pepsin with a high binding affinity. The Van der Waals force and Hydrogen bond interaction were the dominant driving forces during the binding process. The circular dichroism (CD) spectroscopy analysis showed that the CA binding to pepsin decreased the contents of α-Helix and Random Coil but increased the content of β-sheet in the pepsin structure. Accordingly, MD simulations confirmed all the experimental results. As a result, CA is considered an inhibitor with adverse effects on pepsin activity.

Molecular docking and dynamic approach to screen the drug candidate against the Imipenem-resistant CarO porin in Acinetobacter baumannii

The multidrug-resistant Acinetobacter baumannii is an emerging nosocomial pathogen in the healthcare sector. Intrinsic resistance in A. baumannii is a significant problem framing a perfect treatment regimen. Also, this organism showed more resistance towards the carbapenem antibiotics, especially for imipenem and meropenem. The development of carbapenem-resistant Acinetobacter baumannii is mainly due to the alteration or loss of the porin region in the outer membrane. The most well-known porin in Acinetobacter baumannii is CarO (carbapenem-associated outer membrane protein). The CarO protein, which functions as a porin channel for carbapenem inflow, may contribute to carbapenem resistance. The current study identifies a potent drug candidate with a better binding affinity to the carbapenem-resistant outer membrane protein. We investigated the specificity of carbapenems such as imipenem, meropenem, ertapenem, biapenem, doripenem, and fluoroquinolone drugs such as sitafloxacin against the imipenem-resistant CarO protein was demonstrated using the computational approaches molecular docking and dynamic simulation for 50 ns. As a result, the high to low enzyme-ligand complex's binding affinity exhibited a greater binding affinity for ertapenem -7.76 kcal·mol^-1 and sitafloxacin -7.75 kcal·mol^-1 than biapenem, doripenem, meropenem, and imipenem. The molecular dynamic simulation and the MMPBSA analysis depicted ertapenem -55.431±25.908 kJ/mol and sitafloxacin -47.154 ± 11.052 kJ/mol with better binding affinity and more stability against the imipenem resistant CarO protein when it compared to other antibiotics.

Mechanistic Insight into the Inhibition of Choline Acetyltransferase by Proton Pump Inhibitors

Various pharmacoepidemiological investigational studies have indicated that Proton Pump Inhibitors (PPIs) may increase the likelihood of developing Alzheimer's disease (AD) and non-AD related dementias. Previously, we have reported the inhibition of the acetylcholine biosynthesizing enzyme choline acetyltransferase (ChAT) by PPIs, for which omeprazole, lansoprazole, and pantoprazole exhibited IC₅₀ values of 0.1, 1.5, and 5.3 μM, respectively. In this study we utilize a battery of computational tools to perceive a mechanistic insight into the molecular interaction of PPIs with the ChAT binding pocket that may further help in designing novel ChAT ligands. Various in-silico tools make it possible for us to elucidate the binding interaction, conformational stability, and dynamics of the protein-ligand complexes within a 200 ns time frame. Further, the binding free energies for the PPI-ChAT complexes were explored. The results suggest that the PPIs exhibit equal or higher binding affinity toward the ChAT catalytic tunnel and are stable throughout the simulated time and that the pyridine ring of the PPIs interacts primarily with the catalytic residue His₃₂₄. A free energy landscape analysis showed that the folding process was linear, and the residue interaction network analysis can provide insight into the roles of various amino acid residues in stabilization of the PPIs in the ChAT binding pocket. As a major factor for the onset of Alzheimer's disease is linked to cholinergic dysfunction, our previous and the present findings give clear insight into the PPI interaction with ChAT. The scaffold can be further simplified to develop novel ChAT ligands, which can also be used as ChAT tracer probes for the diagnosis of cholinergic dysfunction and to initiate timely therapeutic interventions to prevent or delay the progression of AD.

Quantum computing led innovation for achieving a more sustainable Covid-19 healthcare industry

Involvement of multiple stakeholders in healthcare industry, even the simple healthcare problems become complex due to classical approach to treatment. In the Covid-19 era where quick and accurate solutions in healthcare are needed along with quick collaboration of stakeholders such as patients, insurance agents, healthcare providers and medicine supplier etc., a classical computing approach is not enough. Therefore, this study aims to identify the role of quantum computing in disrupting the healthcare sector with the lens of organizational information processing theory (OIPT), creating a more sustainable (less strained) healthcare system. A semi-structured interview approach is adopted to gauge the expectations of professionals from healthcare industry regarding quantum computing. A structured approach of coding, using open, axial and selective approach is adopted to map the themes under quantum computing for healthcare industry. The findings indicate the potential applications of quantum computing for pharmaceutical, hospital, health insurance organizations along with patients to have precise and quick solutions to the problems, where greater accuracy and speed can be achieved. Existing research focuses on the technological background of quantum computing, whereas this study makes an effort to mark the beginning of quantum computing research with respect to organizational management theory.

In silico drug screen reveals potential competitive MTHFR inhibitors for clinical repurposing

MTHFR (Methylenetetrahydrofolate reductase) is a pivotal enzyme involved in one-carbon metabolism, which is critical for the proliferation of cancer cells. In line with this, published literature showed that MTHFR knockdown caused impaired growth of multiple types of cancer cells. Moreover, higher MTHFR expression levels were linked to shorter overall survival in hepatocellular carcinoma, adrenocortical carcinoma, and low-grade glioma, bringing the need to design MTHFR inhibitors as a possible treatment option. No competitive inhibitors of MTHFR have been reported as of today. This study aimed to identify potential competitive MTHFR inhibitor candidates using an in silico drug screen. A total of 30470 molecules containing biogenic compounds, FDA-approved drugs, and those in clinical trials were screened against the catalytic pocket of MTHFR in the presence and absence of cofactors. Binding energy and ADMET analysis revealed that Vilanterol (β2-adrenergic agonist), Selexipag (prostacyclin receptor agonist), and Ramipril Diketopiperazine (ACE inhibitor) are potential competitive inhibitors of MTHFR. Molecular dynamics analyses and MM-PBSA calculations with these compounds particularly revealed the amino acids between 285-290 for ligand binding and highlighted Vilanterol as the strongest candidate for MTHFR inhibition. Our results could guide the development of novel MTHFR inhibitor compounds, which could be inspired by the drugs brought into the spotlight here. More importantly, these potential candidates could be quhickly tested as a repurposing strategy in pre-clinical and clinical studies of the cancers mentioned above.Communicated by Ramaswamy H. Sarma.

Identification of Diosmin and Flavin Adenine Dinucleotide as Repurposing Treatments for Monkeypox Virus: A Computational Study

The World Health Organization declared monkeypox a global public health emergency on 23 July 2022. This disease was caused by the monkeypox virus (MPXV), which was first identified in 1958 in Denmark. The MPXV is a member of the Poxviridae family, the Chordopoxvirinae subfamily, and the genus Orthopoxvirus, which share high similarities with the vaccinia virus (the virus used to produce the smallpox vaccine). For the initial stage of infection, the MPXV needs to attach to the human cell surface glycosaminoglycan (GAG) adhesion molecules using its E8 protein. However, up until now, neither a structure for the MPXV E8 protein nor a specific cure for the MPXV exists. This study aimed to search for small molecules that inhibit the MPXV E8 protein, using computational approaches. In this study, a high-quality three-dimensional structure of the MPXV E8 protein was retrieved by homology modeling using the AlphaFold deep learning server. Subsequent molecular docking and molecular dynamics simulations (MDs) for a cumulative duration of 2.1 microseconds revealed that ZINC003977803 (Diosmin) and ZINC008215434 (Flavin adenine dinucleotide-FAD) could be potential inhibitors against the E8 protein with the MM/GBSA binding free energies of −38.19 ± 9.69 and −35.59 ± 7.65 kcal·mol⁻¹, respectively.

Control of eriocitrin release from pH-sensitive gelatin-based microgels to inhibit α-glucosidase: an experimental and computational study

α-Glucosidase is among the intestinal epithelial enzymes that produce absorbable glucose in the final stage of glycan catabolism. It leads to an increase in blood glucose levels as a result of high glucose uptake in diabetic patients. However, inhibition of this essential biochemical process can be a useful therapeutic approach to diabetes mellitus (DM). Eriocitrin (ER) is an abundant "flavanone glycoside" in citrus fruits with rich antioxidant properties whose effects on α-Glu inhibition in the small intestine remain to be determined. Herein, pH-sensitive microgels (MGs) were designed based on cross-linked methacrylate with acrylamide (AM) and acrylic acid (AAc) (molar ratio 70 : 30 of AAc : AM) as a controlled release system for sustained delivery of ER into the small intestine. The presence of amide and acrylate in MGs and the mechanical resistance were determined using FT-IR spectroscopy, rheology, and viscoelastometry. In vitro experiments showed that MGs could protect ER against diffusion in the gastric location and adjust its release in the intestinal milieu. The intestinal α-Glu activity was inhibited by ER (IC₅₀ value of 12.50 ± 0.73 μM) in an uncompetitive dose-dependent manner. The presence of ER altered the structure of α-Glu and reduced the hydrophobic pockets of the enzyme. Molecular docking analysis along with molecular dynamics simulation displayed that ER-α-Glu formation is directed by hydrogen binding with Asp69, Asp215, Glu411, Asp307, and Tyr347 residues. Moreover, in vivo assessment showed that rat blood glucose concentration decreased after ER administration compared with the control group. The results highlight that ER-loaded-MGs can be considered as a useful releasing strategy in treating DM via α-Glu inhibition.

Mass Spectrometry-Based Chemical Proteomics for Drug Target Discoveries

Chemical proteomics, emerging rapidly in recent years, has become a main approach to identifying interactions between the small molecules and proteins in the cells on a proteome scale and mapping the signaling and/or metabolic pathways activated and regulated by these interactions. The methods of chemical proteomics allow not only identifying proteins targeted by drugs, characterizing their toxicity and discovering possible off-target proteins, but also elucidation of the fundamental mechanisms of cell functioning under conditions of drug exposure or due to the changes in physiological state of the organism itself. Solving these problems is essential for both basic research in biology and clinical practice, including approaches to early diagnosis of various forms of serious diseases or prediction of the effectiveness of therapeutic treatment. At the same time, recent developments in high-resolution mass spectrometry have provided the technology for searching the drug targets across the whole cell proteomes. This review provides a concise description of the main objectives and problems of mass spectrometry-based chemical proteomics, the methods and approaches to their solution, and examples of implementation of these methods in biomedical research.

Phospholipid peroxidation inhibits autophagy via stimulating the delipidation of oxidized LC3-PE

Phospholipid peroxidation of polyunsaturated fatty acids at the bis-allylic position drives ferroptosis. Here we identify a novel role for phospholipid peroxidation in the inhibition of autophagy. Using in vitro and in vivo models, we report that phospholipid peroxidation induced by glutathione peroxidase-4 inhibition and arachidonate 15-lipoxygenase overexpression leads to overload of peroxidized phospholipids and culminate in inhibition of autophagy. Functional and lipidomics analysis further demonstrated that inhibition of autophagy was associated with an increase of peroxidized phosphatidylethanolamine (PE) conjugated LC3. We further demonstrate that autophagy inhibition occurred due to preferential cleavage of peroxidized LC3-PE by ATG4B to yield delipidated LC3. Mouse models of phospholipid peroxidation and autophagy additionally supported a role for peroxidized PE in autophagy inhibition. Our results agree with the recognized role of endoplasmic reticulum as the primary source for autophagosomal membranes. In summary, our studies demonstrated that phospholipid peroxidation inhibited autophagy via stimulating the ATG4B-mediated delipidation of peroxidized LC3-PE.

Probing the polyphenolic flavonoid, morin as a highly efficacious inhibitor against amyloid(A4V) mutant SOD1 in fatal amyotrophic lateral sclerosis

Deposition of misfolded protein aggregates in key areas of human brain is the quintessential trait of various pertinent neurodegenerative disorders including amyotrophic lateral sclerosis (ALS). Genetic point mutations in Cu/Zn superoxide dismutase (SOD1) are found to be the most important contributing factor behind familial ALS. Especially, single nucleotide polymorphism (SNP) A4V is the most nocuous since it substantially decreases life expectancy of patients. Besides, the use of naturally occurring polyphenolic flavonoids is profoundly being advocated for palliating amyloidogenic behavior of proteopathic proteins. In the present analysis, through proficient computational tools, we have attempted to ascertain a pharmacodynamically promising flavonoid compound that effectively curbs the pathogenic behavior of A4V SOD1 mutant. Initial screening of flavonoids that exhibit potency against amyloids identified morin, myricetin and epigallocatechin gallate as promising leads. Further, with the help of feasible and yet adept protein-ligand interaction studies and stalwart molecular simulation analyses, we were able to observe that aforementioned flavonoids were able to considerably divert mutant A4V SOD1 from its distinct pathogenic behavior. Among which, morin showed the most curative potential against A4V SOD1. Therefore, morin holds a great therapeutic potential in contriving highly efficacious inhibitors in mitigating fatal and insuperable ALS.

Computational Study on the Inhibitory Effect of Natural Compounds against the SARS-CoV-2 Proteins

COVID-19 is more virulent and challenging to human life. In India, the Ministry of AYUSH recommended some strategies through Siddha, homeopathy, and other methods to effectively manage COVID-19 (Guidelines for AYUSH Clinical Studies in COVID-19, 2020). Kabasura Kudineer and homeopathy medicines are in use for the prevention and treatment of COVID-19 infection; however, the mechanism of action is less explored. This study aims to understand the antagonist activity of natural compounds found in Kabasura Kudineer and homeopathy medicines against the SARS-CoV-2 using computational methods. Potential compounds were screened against NSP-12, NSP-13, NSP-14, NSP-15, main protease, and spike proteins. Structure-based virtual screening results shows that, out of 14,682 Kabasura Kudineer compounds, the 250395, 129677029, 44259583, 44259584, and 88583189 compounds and, out of 3,112 homeopathy compounds, the 3802778, 320361, 5315832, 14590080, and 74029795 compounds have good scoring function against the SARS-CoV-2 structural and nonstructural proteins. As a result of docking, homeopathy compounds have a docking score ranging from -5.636 to 13.631 kcal/mol, while Kabasura Kudineer compounds have a docking score varying from -8.290 to -13.759 kcal/mol. It has been found that the selected compounds bind well to the active site of SARS-CoV-2 proteins and form hydrogen bonds. The molecular dynamics simulation study shows that the selected compounds have maintained stable conformation in the simulation period and interact with the target. This study supports the antagonist activity of natural compounds from Kabasura Kudineer and homeopathy against SARS-CoV-2's structural and nonstructural proteins.

In Silico Drug Repurposing Approach: Investigation of Mycobacterium tuberculosis FadD32 Targeted by FDA-Approved Drugs

Background: Despite the enormous efforts made towards combating tuberculosis (TB), the disease remains a major global threat. Hence, new drugs with novel mechanisms against TB are urgently needed. Fatty acid degradation protein D32 (FadD32) has been identified as a promising drug target against TB, the protein is required for the biosynthesis of mycolic acids, hence, essential for the growth and multiplication of the mycobacterium. However, the FadD32 mechanism upon the binding of FDA-approved drugs is not well established. Herein, we applied virtual screening (VS), molecular docking, and molecular dynamic (MD) simulation to identify potential FDA-approved drugs against FadD32. Methodology/Results: VS technique was found promising to identify four FDA-approved drugs (accolate, sorafenib, mefloquine, and loperamide) with higher molecular docking scores, ranging from -8.0 to -10.0 kcal/mol. Post-MD analysis showed that the accolate hit displayed the highest total binding energy of -45.13 kcal/mol. Results also showed that the accolate hit formed more interactions with FadD32 active site residues and all active site residues displayed an increase in total binding contribution. RMSD, RMSF, Rg, and DCCM analysis further supported that the presence of accolate exhibited more structural stability, lower bimolecular flexibility, and more compactness into the FadD32 protein. Conclusions: Our study revealed accolate as the best potential drug against FadD32, hence a prospective anti-TB drug in TB therapy. In addition, we believe that the approach presented in the current study will serve as a cornerstone to identifying new potential inhibitors against a wide range of biological targets.

Discovery of small molecular inhibitors for interleukin-33/ST2 protein-protein interaction: a virtual screening, molecular dynamics simulations and binding free energy calculations

The interleukin-1 receptor like ST2 has emerged as a potential drug discovery target since it was identified as the receptor of the novel cytokine IL-33, which is involved in many inflammatory and autoimmune diseases. For the treatment of such IL-33-related disorders, efforts have been made to discover molecules that can inhibit the protein–protein interactions (PPIs) between IL-33 and ST2, but to date no drug has been approved. Although several anti-ST2 antibodies have entered clinical trials, the exploration of small molecular inhibitors is highly sought-after because of its advantages in terms of oral bioavailability and manufacturing cost. The aim of this study was to discover ST2 receptor inhibitors based on its PPIs with IL-33 in crystal structure (PDB ID: 4KC3) using virtual screening tools with pharmacophore modeling and molecular docking. From an enormous chemical space ZINC, a potential series of compounds has been discovered with stronger binding affinities than the control compound from a previous study. Among them, four compounds strongly interacted with the key residues of the receptor and had a binding free energy <  − 20 kcal/mol. By intensive calculations using data from molecular dynamics simulations, ZINC59514725 was identified as the most potential candidate for ST2 receptor inhibitor in this study.