Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors

Phytochemical on-line screening and in silico study of Helianthemum confertum: antioxidant activity, DFT, MD simulation, ADME/T analysis, and xanthine oxidase binding

Seven components from the methanol extract of the aerial part of the endemic species Helianthemum confertum were isolated and identified for the first time. Investigating this species and its separated components chemical make-up and radical scavenging capacity, was the main goal. Using an online HPLC-ABTS˙+ test, ORAC, and TEAC assays, the free radical scavenging capacity of the ethyl acetate extract was assessed. The fractionation of these extracts by CC, TLC, and reverse-phase HPLC was guided by the collected data, which was corroborated by TEAC and ORAC assays. Molecular docking studies, DFT at the B3LYP level, and an examination of the ADME/T predictions of all compounds helped to further clarify the phytochemicals' antioxidant potential. Isolation and identification of all components were confirmed through spectroscopy, which revealed a mixture (50-50%) of para-hydroxybenzoic acid 1 and methyl gallate 2, protocatechuic acid 3, astragalin 4, trans-tiliroside 5, cis-tiliroside 6, contaminated by trans-tiliroside and 3-oxo-α-ionol-β-d-glucopyranoside 7, as well as two new compounds for the genus Helianthemum (2 and 7). With a focus on compounds 1, 2, 3, and 4, the results clearly showed that the extract and the compounds tested from this species had a high antioxidant capacity. Within the xanthine oxidase enzyme's pocket, all of the components tested showed strong and stable binding. In light of these findings, the xanthine oxidase/methyl gallate 2 complex was simulated using the Desmond module of the Schrodinger suite molecular dynamics (MD) for 100 ns. Substantially stable receptor-ligand complexes were observed following 1 ns of MD simulation.

Computational method for designing vaccines applied to virus-like particles (VLPs) as epitope carriers

Nonenveloped virus-like particles (VLPs) are self-assembled oligomeric structures composed of one or more proteins that originate from diverse viruses. Because these VLPs have similar antigenicity to the parental virus, they are successfully used as vaccines against cognate virus infection. Furthermore, after foreign antigenic sequences are inserted in their protein components (chimVLPs), some VLPs are also amenable to producing vaccines against pathogens other than the virus it originates from (these VLPs are named platform or epitope carrier). Designing chimVLP vaccines is challenging because the immunogenic response must be oriented against a given antigen without altering stimulant properties inherent to the VLP. An important step in this process is choosing the location of the sequence modifications because this must be performed without compromising the assembly and stability of the original VLP. Currently, many immunogenic data and computational tools can help guide the design of chimVLPs, thus reducing experimental costs and work. In this study, we analyze the structure of a novel VLP that originate from an insect virus and describe the putative regions of its three structural proteins amenable to insertion. For this purpose, we employed molecular dynamics (MD) simulations to assess chimVLP stability by comparing mutated and wild-type (WT) VLP protein trajectories. We applied this procedure to design a chimVLP that can serve as a prophylactic vaccine against the SARS-CoV-2 virus. The methodology described in this work is generally applicable for VLP-based vaccine development.

Investigating Pyruvate Dehydrogenase Kinase 3 Inhibitory Potential of Myricetin Using Integrated Computational and Spectroscopic Approaches

Protein kinases are involved in various diseases and currently represent potential targets for drug discovery. These kinases play major roles in regulating the cellular machinery and control growth, homeostasis, and cell signaling. Dysregulation of kinase expression is associated with various disorders such as cancer and neurodegeneration. Pyruvate dehydrogenase kinase 3 (PDK3) is implicated in cancer therapeutics as a potential drug target. In this current study, a molecular docking exhibited a strong binding affinity of myricetin to PDK3. Further, a 100 ns all-atom molecular dynamics (MD) simulation study provided insights into the structural dynamics and stability of the PDK3-myricetin complex, revealing the formation of a stable complex with minimal structural alterations upon ligand binding. Additionally, the actual affinity was ascertained by fluorescence binding studies, and myricetin showed appreciable binding affinity to PDK3. Further, the kinase inhibition assay suggested significant inhibition of PDK3 by myricetin, revealing an excellent inhibitory potential with an IC50 value of 3.3 μM. In conclusion, this study establishes myricetin as a potent PDK3 inhibitor that can be implicated in therapeutic targeting cancer and PDK3-associated diseases. In addition, this study underscores the efficacy of myricetin as a potential lead to drug discovery and provides valuable insights into the inhibition mechanism, enabling advancements in cancer therapeutics.

A theoretical screening of phytochemical constituents from Millettia brandisiana as inhibitors against acetylcholinesterase

Alzheimer's disease is one of the causes associated with the early stages of dementia. Nowadays, the main treatment available is to inhibit the actions of the acetylcholinesterase (AChE) enzyme, which has been identified as responsible for the disease. In this study, computational methods were used to examine the structure and therapeutic ability of chemical compounds extracted from Millettia brandisiana natural products against AChE. This plant is commonly known as a traditional medicine in Vietnam and Thailand for the treatment of several diseases. Furthermore, machine learning helped us narrow down the choice of 85 substances for further studies by molecular docking and molecular dynamics simulations to gain deeper insights into the interactions between inhibitors and disease proteins. Of the five top-choice substances, γ-dimethylallyloxy-5,7,2,5-tetramethoxyisoflavone emerges as a promising substance due to its large free binding energy to AChE and the high thermodynamic stability of the resulting complex.

Understanding Genetic Risks: Computational Exploration of Human β-Synuclein nsSNPs and their Potential Impact on Structural Alteration

Synucleins are pivotal in neurodegenerative conditions. Beta-synuclein (β-synuclein) is part of the synuclein protein family alongside alpha-synuclein (α-synuclein) and gamma-synuclein (γ-synuclein). These proteins, found mainly in brain tissue and cancers, are soluble and unstructured. β-synuclein shares significant similarity with α-synuclein, especially in their N-terminus, with a 90% match. However, their aggregation tendencies differ significantly. While α-synuclein aggregation is believed to be counteracted by β-synuclein, which occurs in conditions like Parkinson's disease, β-synuclein may counteract α-synuclein's toxic effects on the nervous system, offering potential treatment for neurodegenerative diseases. Under normal circumstances, β-synuclein may guard against disease by interacting with α-synuclein. Yet, in pathological environments with heightened levels or toxic substances, it might contribute to disease. Our research aims to explore potential harmful mutations in the β-synuclein using computational tools to predict their destabilizing impact on protein structure. Consensus analysis revealed rs1207608813 (A63P), rs1340051870 (S72F), and rs1581178262 (G36C) as deleterious. These findings highlight the intricate relationship between nsSNPs and protein function, shedding light on their potential implications in disease pathways. Understanding the structural consequences of nsSNPs is crucial for elucidating their role in pathogenesis and developing targeted therapeutic interventions. Our results offer a robust computational framework for identifying neurodegenerative disorder-related mutations from SNP datasets, potentially reducing the costs associated with experimental characterization.

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.

Structure-guided engineered urethanase from Candida parapsilosis with pH and ethanol tolerance to efficiently degrade ethyl carbamate in Chinese rice wine

Urethane hydrolase can degrade the carcinogen ethyl carbamate (EC) in fermented food, but its stability and activity limit its application. In this study, a mutant G246A and a double mutant N194V/G246A with improved cpUH activity and stability of Candida parapsilosis were obtained by site-directed mutagenesis. The catalytic efficiency (Kcat/Km) of mutant G246A and double mutant N194V/G246A are 1.95 times and 1.88 times higher than that of WT, respectively. In addition, compared with WT, the thermal stability and pH stability of mutant G246A and double mutant N194V/G246A were enhanced. The ability of mutant G246A and double mutant N194V/G246A to degrade EC in rice wine was also stronger than that of WT. The mutation increased the stability of the enzyme, as evidenced by decreased root mean square deviation (RMSD) and increased hydrogen bonds between the enzyme and substrate by molecular dynamics simulation and molecular docking analysis. The molecule modification of new cpUH promotes the industrial process of EC degradation.

Discovery of novel BRD4-BD2 inhibitors via in silico approaches: QSAR techniques, molecular docking, and molecular dynamics simulations

Bromodomain-containing protein 4(BRD4) plays an important role in the occurrence and development of various malignant tumors, which has attracted the attention of scientific research institutions and pharmaceutical companies. The structural modification of most currently available BRD4 inhibitors is relatively simple, but the drug effectiveness is limited. Research has found that the inhibition of BD1 may promote the differentiation of oligodendrocyte progenitor cell; however, the inhibition of BD2 will not cause this outcome. Therefore, newly potential drugs which target BRD4-BD2 need further research. Herein, we initially built QSAR models out of 49 compounds using HQSAR, CoMFA, CoMSIA, and Topomer CoMFA technology. All of the models have shown suitable reliabilities (q2 = 0.778, 0.533, 0.640, 0.702, respectively) and predictive abilities (r2pred = 0.716, 0.6289, 0.6153, 0.7968, respectively) for BRD4-BD2 inhibitors. On the basis of QSAR results and the search of the R-group in the topomer search module, we designed 20 new compounds with high activity that showed appropriate docking score and suitable ADMET. Docking studies and MD simulation were carried out to reveal the amino acid residues (Asn351, Cys347, Tyr350, Pro293, and Asp299) at the active site of BRD4-BD2. Free energy calculations and free energy landscapes verified the stable binding results and indicated stable conformations of the complexes. These theoretical studies provide guidance and theoretical basis for designing and developing novel BRD4-BD2 inhibitors.

β-keto amyrin isolated from Cryptostegia grandiflora R. br. inhibits inflammation caused by Daboia russellii viper venom: Direct binding of β-keto amyrin to phospholipase A2

The search for mechanism-based anti-inflammatory therapies is of fundamental importance to avoid undesired off-target effects. Phospholipase A2 (PLA2) activity is a potential molecular target for anti-inflammatory drugs because it fuels arachidonic acid needed to synthesize inflammation mediators, such as prostaglandins. Herein, we aim to investigate the molecular mechanism by which β-keto amyrin isolated from a methanolic extract of Cryptostegia grandiflora R. Br. Leaves can inhibit inflammation caused by Daboia russellii viper (DR) venom that mainly contains PLA2. We found that β-keto amyrin neutralizes DR venom-induced paw-edema in a mouse model. Molecular docking of PLA2 with β-keto amyrin complex resulted in a higher binding energy score of -8.86 kcal/mol and an inhibition constant of 611.7 nM. Diclofenac had a binding energy of -7.04 kcal/mol and an IC50 value of 620 nM, which predicts a poorer binding interaction than β-keto amyrin. The higher conformational stability of β-keto amyrin interaction compared to diclofenac is confirmed by molecular dynamics simulation. β-keto amyrin isolated from C. grandiflora inhibits the PLA2 activity contained in Daboia russellii viper venom. The anti-inflammatory property of β-keto amyrin is due to its direct binding into the active site of PLA2, thus inhibiting its enzyme activity.

Tubulin CFEOM mutations both inhibit or activate kinesin motor activity

Kinesin-mediated transport along microtubules is critical for axon development and health. Mutations in the kinesin Kif21a, or the microtubule subunit β-tubulin, inhibit axon growth and/or maintenance resulting in the eye-movement disorder congenital fibrosis of the extraocular muscles (CFEOM). While most examined CFEOM-causing β-tubulin mutations inhibit kinesin-microtubule interactions, Kif21a mutations activate the motor protein. These contrasting observations have led to opposed models of inhibited or hyperactive Kif21a in CFEOM. We show that, contrary to other CFEOM-causing β-tubulin mutations, R380C enhances kinesin activity. Expression of β-tubulin-R380C increases kinesin-mediated peroxisome transport in S2 cells. The binding frequency, percent motile engagements, run length and plus-end dwell time of Kif21a are also elevated on β-tubulin-R380C compared with wildtype microtubules in vitro. This conserved effect persists across tubulins from multiple species and kinesins from different families. The enhanced activity is independent of tail-mediated kinesin autoinhibition and thus utilizes a mechanism distinct from CFEOM-causing Kif21a mutations. Using molecular dynamics, we visualize how β-tubulin-R380C allosterically alters critical structural elements within the kinesin motor domain, suggesting a basis for the enhanced motility. These findings resolve the disparate models and confirm that inhibited or increased kinesin activity can both contribute to CFEOM. They also demonstrate the microtubule's role in regulating kinesins and highlight the importance of balanced transport for cellular and organismal health.

Hexaconazole exposure may lead to Parkinson via disrupting glucocerebrosidase and parkin: molecular interaction, dynamics, MMPBSA and DFT based in-silico predictive toxicology

Hexaconazole is a known fungicide for agricultural purposes. It has bioaccumulation ability which makes it important for its toxicological characterization. There are various neurological impacts of pollutants on human health. Therefore, in this study, we have done predictive analyses of the interaction mechanism of hexaconazole by molecular interaction analysis, molecular dynamics simulation, and Poisson-Boltzmann surface area (MM-PBSA) to assess hexaconazole's potency to disrupt the homeostasis of glucocerebrosidase (-7.9 kcal/mol) and parkin (-5.67 kcal/mol) proteins which have significant roles in the manifestation of Parkinson disease. The findings reveal that hexaconazole has the potency to form stable interactions with glucocerebrosidase and parkin. This research provides a molecular and atomic-level understanding of how hexaconazole exposure may disrupt the homeostasis of glucocerebrosidase and parkin. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration, and hydrogen bonding exhibited the potent molecular interactions of hexaconazole, which may lead to neurological manifestations such as Parkinson disease.

Investigating the role of thymol as a promising inhibitor of pyruvate dehydrogenase kinase 3 for targeted cancer therapy

Protein kinases have emerged as major contributors to various diseases. They are currently exploited as a potential target in drug discovery because they play crucial roles in cell signaling, growth, and regulation. Their dysregulation is associated with inflammatory disorders, cancer, and neurodegenerative diseases. Pyruvate dehydrogenase kinase 3 (PDK3) has become an attractive drug target in cancer therapeutics. In the present study, we investigated the effective role of thymol in PDK3 inhibition due to the high affinity predicted through molecular docking studies. Hence, to better understand this inhibition mechanism, we carried out a 100 ns molecular dynamics (MD) simulation to analyse the dynamics and stability of the PDK3-thymol complex. The PDK3-thymol complex was stable and energetically favourable, with many intramolecular hydrogen bond interactions in the PDK3-thymol complex. Enzyme inhibition assay showed significant inhibition of PDK3 by thymol, revealing potential inhibitory action of thymol towards PDK3 (IC50 = 2.66 μM). In summary, we established thymol as one of the potential inhibitors of PDK3, proposing promising therapeutic implications for severe diseases associated with PDK3 dysregulation. This study further advances our understanding of thymol's therapeutic capabilities and potential role in cancer treatment.

Computational investigation of Moringa oleifera phytochemicals targeting EGFR: molecular docking, molecular dynamics simulation and density functional theory studies

Epidermal growth factor receptor (EGFR) is a prominent target for anticancer therapy due to its role in activating several cell signaling cascades. Clinically approved EGFR inhibitors are reported to show treatment resistance and toxicity, this study, therefore, investigates Moringa oleifera phytochemicals to find potent and safe anti-EGFR compounds. For that, phytochemicals were screened based on drug-likeness and molecular docking analysis followed by molecular dynamics simulation, density functional theory analysis and ADMET analysis to identify the effective inhibitors of EGFR tyrosine kinase (EGFR-TK) domain. Known EGFR-TK inhibitors (1-4 generations) were used as control. Among 146 phytochemicals, 136 compounds showed drug-likeness, of which Delta 7-Avenasterol was the most potential EGFR-TK inhibitor with a binding energy of -9.2 kcal/mol followed by 24-Methylenecholesterol (-9.1 kcal/mol), Campesterol (-9.0 kcal/mol) and Ellagic acid (-9.0 kcal/mol). In comparison, the highest binding affinity from control drugs was displayed by Rociletinib (-9.0 kcal/mol). The molecular dynamics simulation (100 ns) exhibited the structural stability of native EGFR-TK and protein-inhibitor complexes. Further, MM/PBSA computed the binding free energies of protein complex with Delta 7-Avenasterol, 24-Methylenecholesterol, Campesterol and Ellagic acid as -154.559 ± 18.591 kJ/mol, -139.176 ± 19.236 kJ/mol, -136.212 ± 17.598 kJ/mol and -139.513 ± 23.832 kJ/mol, respectively. Non-polar interactions were the major contributors to these energies. The density functional theory analysis also established the stability of these inhibitor compounds. ADMET analysis depicted acceptable outcomes for all top phytochemicals without displaying any toxicity. In conclusion, this report has identified promising EGFR-TK inhibitors to treat several cancers that can be further investigated through laboratory and clinical tests.

Antifungal and Antiproliferative Activity of Pistagremic Acid and Flavonoids Extracted from the Galls of Pistacia chinensis subsp. integerrima

Pistacia chinensis subsp. integerrima (J.L. Stewart) Rech. f. is a plant known for its therapeutic applications in traditional medicine, which are related to its antimicrobial, anticancer, antioxidant, anti-inflammatory, analgesic, antidiarrheal, and muscle relaxant properties. The galls of P. chinensis are rich in triterpenes and flavonoids, and we here report the extraction of pistagremic acid (1), apigenin (2) and sakuranetin (3) from this source. The isolated compounds were tested against Aspergillus flavus, Candida albicans, Candida glabrata, Fusarium solani, Microsporum canis and Trichoderma longibrachiatum. The results highlighted the antimicrobial activity of flavonoids 2 and 3, suggesting that this class of molecules may be responsible for the effect related to the traditional use. On the other hand, when the compounds and the extract were tested for their antiproliferative activity on a panel of 4 human cancer cell lines, the triterpene pistagremic acid (1) showed a higher potential, thus demonstrating a different bioactivity profile. Structure-based docking and molecular dynamics simulations were used to help the interpretation of experimental results. Taken together, the here reported findings pave the way for the rationalization of the use of P. chinensis extracts, highlighting the contributions of the different components of galls to the observed bioactivity.

Computational Approaches to Evaluate the Acetylcholinesterase Binding Interaction with Taxifolin for the Management of Alzheimer's Disease

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) are enzymes that break down and reduce the level of the neurotransmitter acetylcholine (ACh). This can cause a variety of cognitive and neurological problems, including Alzheimer's disease. Taxifolin is a natural phytochemical generally found in yew tree bark and has significant pharmacological properties, such as being anti-cancer, anti-inflammatory, and antioxidant. The binding affinity and inhibitory potency of taxifolin to these enzymes were evaluated through molecular docking and molecular dynamics simulations followed by the MMPBSA approach, and the results were significant. Taxifolin's affinity for binding to the AChE-taxifolin complex was -8.85 kcal/mol, with an inhibition constant of 326.70 nM. It was observed to interact through hydrogen bonds. In contrast, the BChE-taxifolin complex binding energy was observed to be -7.42 kcal/mol, and it was significantly nearly equal to the standard inhibitor donepezil. The molecular dynamics and simulation signified the observed interactions of taxifolin with the studied enzymes. The MMPBSA total free energy of binding for AChE-taxifolin was -24.34 kcal/mol, while BChE-taxifolin was -16.14 kcal/mol. The present research suggests that taxifolin has a strong ability to bind and inhibit AChE and BChE and could be used to manage neuron-associated problems; however, further research is required to explore taxifolin's neurological therapeutic potential using animal models of Alzheimer's disease.

Computational investigation of impact of Pb(II) and Ni(II) ions on hUNG enzyme: insights from molecular dynamics simulations

Human uracil DNA glycosylase (hUNG), a crucial player in the initiation of the base excision repair pathway, is susceptible to alterations in function and conformation induced by the accumulation of toxic metals. Despite the recognized impact of toxic metals on DNA repair enzymes, there exists a notable deficiency in theoretical investigations addressing this phenomenon. This study investigates the impact of toxic heavy metal ions, Pb(II) and Ni(II), on the stability of hUNG through molecular dynamics (MD) simulations. The initial analysis involved the identification of key cavities in the hUNG enzyme. Notably, the active site cavity emerged as a promising site for ligand binding. Subsequently, AutoDockTools software was employed to dock Pb(II) and Ni(II) onto the identified cavities, followed by extensive MD simulations. The MD analysis, encompassing parameters such as root mean square deviation, radius of gyration, solvent accessible surface area, hydrogen bond variations, Ramachandran plot, principal component analysis, and root mean square fluctuations, collectively revealed distinct alterations in the behavior of the enzyme upon complexation with Pb(II) and Ni(II). Interestingly, the enzyme exhibited enhanced structural stability, reduced flexibility, and modified hydrogen bonding patterns in the presence of these toxic metal ions. The observed limitation in structural flexibility implies a more rigid and stable conformation when the enzyme complex with Pb(II) and Ni(II) compared to its free form. This structural alteration may lead to a potential reduction in enzymatic activity, suggesting that toxic metal ions influence the functional dynamics of hUNG. These computational findings offer valuable insights into the molecular interactions between metal ions and enzymes.Communicated by Ramaswamy H. Sarma.

Lead generation of UPPS inhibitors targeting MRSA: Using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking, and molecular dynamic simulations

Undecaprenyl Pyrophosphate Synthase (UPPS) is a vital target enzyme in the early stages of bacterial cell wall biosynthesis. UPPS inhibitors have antibacterial activity against resistant strains such as MRSA and VRE. In this study, we used several consecutive computer-based protocols to identify novel UPPS inhibitors. The 3D QSAR pharmacophore model generation (HypoGen algorithm) protocol was used to generate a valid predictive pharmacophore model using a set of UPPS inhibitors with known reported activity. The developed model consists of four pharmacophoric features: one hydrogen bond acceptor, two hydrophobic, and one aromatic ring. It had a correlation coefficient of 0.86 and a null cost difference of 191.39, reflecting its high predictive power. Hypo1 was proven to be statistically significant using Fischer's randomization at a 95% confidence level. The validated pharmacophore model was used for the virtual screening of several databases. The resulting hits were filtered using SMART and Lipinski filters. The hits were docked into the binding site of the UPPS protein, affording 70 hits with higher docking affinities than the reference compound (6TC, - 21.17 kcal/mol). The top five hits were selected through extensive docking analysis and visual inspection based on docking affinities, fit values, and key residue interactions with the UPPS receptor. Moreover, molecular dynamic simulations of the top hits were performed to confirm the stability of the protein-ligand complexes, yielding five promising novel UPPS inhibitors.

Dissecting the geometric and hydrophobic constraints of stapled peptides

Stapled peptides are a promising class of molecules with potential as highly specific probes of protein-protein interactions and as therapeutics. Hydrocarbon stapling affects the peptide properties through the interplay of two factors: enhancing the overall hydrophobicity and constraining the conformational flexibility. By constructing a series of virtual peptides, we study the role of each factor in modulating the structural properties of a hydrocarbon-stapled peptide PM2, which has been shown to enter cells, engage its target Mouse Double Minute 2 (MDM2), and activate p53. Hamiltonian replica exchange molecular dynamics (HREMD) simulations suggest that hydrocarbon stapling favors helical populations of PM2 through a combination of the geometric constraints and the enhanced hydrophobicity of the peptide. To further understand the conformational landscape of the stapled peptides along the binding pathway, we performed HREMD simulations by restraining the peptide at different distances from MDM2. When the peptide approaches MDM2, the binding pocket undergoes dehydration which appears to be greater in the presence of the stapled peptide compared with the linear peptide. In the binding pocket, the helicity of the stapled peptide is increased due to the favorable interactions between the peptide residues as well as the staple and the microenvironment of the binding pocket, contributing to enhanced affinity. The dissection of the multifaceted mechanism of hydrocarbon stapling into individual factors not only deepens fundamental understanding of peptide stapling, but also provides guidelines for the design of new stapled peptides.

Atomistic molecular simulations of Aβ-Zn conformational ensembles

The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

Identification of common candidate genes and pathways for Spina Bifida and Wilm's Tumor using an integrative bioinformatics analysis

Spina Bifida (SB) and Wilm's Tumor (WT) are conditions, both associated with children. Several studies have shown that WT later develops in SB patients, which led us to elucidate common key genes and linked pathways of both conditions, aimed at their concurrent therapeutic management. For this, integrated bioinformatics analysis was employed. A comprehensive manual curation of genes identified 133 and 139 genes associated with SB and WT, respectively, which were used to construct a single protein-protein interaction (PPI) network. Topological parameters analysis of the network showed its scale-free and hierarchical nature. Centrality-based analysis of the network identified 116 hubs, of which, 6 were called the key genes attributed to being common between SB and WT besides being the hubs. Gene enrichment analysis of the 5 most essential modules, identified important biological processes and pathways possibly linking SB to WT. Additionally, miRNA-key gene-transcription factor (TF) regulatory network elucidated a few important miRNAs and TFs that regulate our key genes. In closing, we put forward TP53, DICER1, NCAM1, PAX3, PTCH1, MTHFR; hsa-mir-107, hsa-mir-137, hsa-mir-122, hsa-let-7d; and YY1, SOX4, MYC, STAT3; key genes, miRNAs and TFs, respectively, as the key regulators. Further, MD simulation studies of wild and Glu429Ala forms of MTHFR proteins showed that there is a slight change in MTHFR protein structure due to Glu429Ala polymorphism. We anticipate that the interplay of these three entities will be an interesting area of research to explore the regulatory mechanism of SB and WT and may serve as candidate target molecules to diagnose, monitor, and treat SB and WT, parallelly.Communicated by Ramaswamy H. Sarma.

Synergistic antimicrobial activity, MD simulation studies and crystal structure of natural alcohol motif containing novel substituted cinnamates

A series of natural alcohols motif containing novel substituted cinnamates were developed and screened against five bacterial strains namely, Enterococcus faecal (E. faecalis), Escherichia coli (E. coli), Bacillus subtilis (B. subtilis), Pseudomonas aeruginosa (P. aeruginosa) and Klebsiella pneumonieae (K. pneumonieae). Among all cinnamates, YS17 was identified with 100% bacterial growth inhibition across the panel, except in E. faecalis with MIC values of 0.25 mg/mL against B. subtilis and P. aeruginosa whereas 0.125, 0.5 and 1 mg/mL against E. coli, K. pneumonieae and E. faecalis, respectively. The growth inhibitory property of YS17 was further validated by disk diffusion, synergistic study and in vitro toxicity assays. Interestingly, YS17 exhibits synergistic effect in combination with the standard drug Ampicillin (AMP). The single crystal structure analysis of YS4 and YS6 was also performed which reconfirmed their proposed structures. Molecular docking visualized significant non-covalent interactions between E. coli MetAP and YS17 and the structural and conformational changes were further analysed using MD simulation studies. Overall, the study provided a suitable core for further synthetic alterations for their optimization as an antibacterial agent. Communicated by Ramaswamy H. Sarma.

Computational Molecular Docking and Simulation-Based Assessment of Anti-Inflammatory Properties of Nyctanthes arbor-tristis Linn Phytochemicals

The leaves, flowers, seeds, and bark of the Nyctanthes arbor-tristis Linn plant have been pharmacologically evaluated to signify the medicinal importance traditionally described for various ailments. We evaluated the anti-inflammatory potentials of 26 natural compounds using AutoDock 4.2 and Molecular Dynamics (MDS) performed with the GROMACS tool. SwissADME evaluated ADME (adsorption, distribution, metabolism, and excretion) parameters. Arb_E and Beta-sito, natural compounds of the plant, showed significant levels of binding affinity against COX-1, COX-2, PDE4, PDE7, IL-17A, IL-17D, TNF-α, IL-1β, prostaglandin E2, and prostaglandin F synthase. The control drug celecoxib exhibited a binding energy of -9.29 kcal/mol, and among the tested compounds, Arb_E was the most significant (docking energy: -10.26 kcal/mol). Beta_sito was also observed with high and considerable docking energy of -8.86 kcal/mol with the COX-2 receptor. COX-2 simulation in the presence of Arb_E and control drug celecoxib, RMSD ranged from 0.15 to 0.25 nm, showing stability until the end of the simulation. Also, MM-PBSA analysis showed that Arb_E bound to COX-2 exhibited the lowest binding energy of -277.602 kJ/mol. Arb_E and Beta_sito showed interesting ADME physico-chemical and drug-like characteristics with significant drug-like effects. Therefore, the studied natural compounds could be potential anti-inflammatory molecules and need further in vitro/in vivo experimentation to develop novel anti-inflammatory drugs.

Structure-based molecular docking and molecular dynamics simulations study for the identification of dipeptidyl peptidase 4 inhibitors in type 2 diabetes

Inhibition of dipeptidyl peptidase-4 (DPP4) activity has emerged as a promising therapeutic approach for the treatment of type 2 diabetes mellitus (T2DM). Bioinformatics-driven approaches have emerged as crucial tools in drug discovery. Molecular docking and molecular dynamics (MD) simulations are effective tools in drug discovery, as they reduce the time and cost associated with experimental screening. In this study, we employed structure-assisted in-silico methods, including molecular docking and MD simulations, to identify SRT2183, a small molecule that may potentially inhibit the activity of DPP4 enzyme. The interaction between the small molecule "SRT2183" and DPP4 exhibited a binding affinity of -9.9 Kcal/Mol, leading to the formation of hydrogen bonds with the amino acid residues MET348, SER376, and THR351 of DPP4. The MD simulations over a period of 100 ns indicated stable protein-ligand interactions, with no significant conformational rearrangements observed within the simulated timeframe. In conclusion, our results suggest that the small molecule SRT2183 may have the potential to inhibit the DPP4 enzyme and pave the way for the therapeutics of T2DM.Communicated by Ramaswamy H. Sarma.

Representing structures of the multiple conformational states of proteins

Biomolecules exhibit dynamic behavior that single-state models of their structures cannot fully capture. We review some recent advances for investigating multiple conformations of biomolecules, including experimental methods, molecular dynamics simulations, and machine learning. We also address the challenges associated with representing single- and multiple-state models in data archives, with a particular focus on NMR structures. Establishing standardized representations and annotations will facilitate effective communication and understanding of these complex models to the broader scientific community.

Macromolecules: Synthesis, antimicrobial, POM analysis and computational approaches of some glucoside derivatives bearing acyl moieties

Macromolecules i.e., carbohydrate derivatives are crucial to biochemical and medical research. Herein, we designed and synthesized eight methyl α-D-glucopyranoside (MGP) derivatives (2-8) in good yields following the regioselective direct acylation method. The structural configurations of the synthesized MGP derivatives were analyzed and verified using multiple physicochemical and spectroscopic techniques. Antimicrobial experiments revealed that almost all derivatives demonstrated noticeable antifungal and antibacterial efficacy. The synthesized derivatives showed minimum inhibitory concentration (MIC) values ranging from 0.75 µg/mL to 1.50 µg/mL and minimum bactericidal concentrations (MBCs) ranging from 8.00 µg/mL to 16.00 µg/mL. Compound 6 inhibited Ehrlich ascites carcinoma (EAC) cell proliferation by 10.36% with an IC50 of 2602.23 μg/mL in the MTT colorimetric assay. The obtained results were further rationalized by docking analysis of the synthesized derivatives against 4URO and 4XE3 receptors to explore the binding affinities and nonbonding interactions of MGP derivatives with target proteins. Compound 6 demonstrated the potential to bind with the target with the highest binding energy. In a stimulating environment, a molecular dynamics study showed that MGP derivatives have a stable conformation and binding pattern. The MGP derivatives were examined using POM (Petra/Osiris/Molinspiration) bioinformatics, and as a result, these derivatives showed good toxicity, bioavailability, and pharmacokinetics. Various antifungal/antiviral pharmacophore (Oδ-, O'δ-) sites were identified by using POM investigations, and compound 6 was further tested against other pathogenic fungi and viruses, such as Micron and Delta mutants of SARS-CoV-2.

Molecular modelling of SdiA protein by selected flavonoid and terpenes compounds to attenuate virulence in Klebsiella pneumoniae

Klebsiella pneumoniae is one of the perturbing multidrug resistant (MDR) and ESKAPE pathogens contributing to the mounting morbidity, mortality and extended rate of hospitalization. Its virulence, often regulated by quorum sensing (QS) reinforces the need to explore alternative and prospective antivirulence agents, relatively from plants secondary metabolites. Computer aided drug discovery using molecular modelling techniques offers advantage to investigate prospective drugs to combat MDR pathogens. Thus, this study employed virtual screening of selected terpenes and flavonoids from medicinal plants to interrupt the QS associated SdiA protein in K. pneumoniae to attenuate its virulence. 4LFU was used as a template to model the structure of SdiA. ProCheck, Verify3D, Ramachandran plot scores, and ProSA-Web all attested to the model's good quality. Since SdiA protein in K. pneumoniae leads to the expression of virulence, 31 prospective bioactive compounds were docked for antagonistic potential. The stability of the protein-ligand complex, atomic motions and inter-atomic interactions were further investigated through molecular dynamics simulations (MDS) at 100 ns production runs. The binding free energy was estimated using the molecular mechanics/poisson-boltzmann surface area (MM/PB-SA). Furthermore, the drug-likeness properties of the studied compounds were validated. Docking studies showed phytol possesses the highest binding affinity (-9.205 kcal/mol) while glycitein had -9.752 kcal/mol highest docking score. The MDS of the protein in complex with the best-docked compounds revealed phytol with the highest binding energy of -44.2625 kcal/mol, a low root-mean-square deviation (RMSD) value of 1.54 Å and root-mean-square fluctuation (RMSF) score of 1.78 Å. Analysis of the drug-likeness properties prediction and bioavailability of these compounds revealed their conformed activity to lipinski's rules with bioavailability scores of 0.55 F. The studied terpenes and flavonoids compounds effectively thwart SdiA protein, therefore regulate inter- or intra cellular communication and associated in virulence Enterobacteriaceae, serving as prospective antivirulence drugs.Communicated by Ramaswamy H. Sarma.

Virtual Screening-Accelerated Discovery of a Phosphodiesterase 9 Inhibitor with Neuroprotective Effects in the Kainate Toxicity In Vitro Model

In the central nervous system, some specific phosphodiesterase (PDE) isoforms modulate pathways involved in neuronal plasticity. Accumulating evidence suggests that PDE9 may be a promising therapeutic target for neurodegenerative diseases. In the current study, computational techniques were used to identify a nature-inspired PDE9 inhibitor bearing the scaffold of an isoflavone, starting from a database of synthetic small molecules using a ligand-based approach. Furthermore, docking studies supported by molecular dynamics investigations allowed us to evaluate the features of the ligand-target complex. In vitro assays confirmed the computational results, showing that the selected compound inhibits the enzyme in the nanomolar range. Additionally, we evaluated the expression of gene and protein levels of PDE9 in organotypic hippocampal slices, observing an increase following exposure to kainate (KA). Importantly, the PDE9 inhibitor reduced CA3 damage induced by KA in a dose-dependent manner in organotypic hippocampal slices. Taken together, these observations strongly support the potential of the identified nature-inspired PDE9 inhibitor and suggest that such a molecule could represent a promising lead compound to develop novel therapeutic tools against neurological diseases..

Engineering chitin deacetylase AsCDA for improving the catalytic efficiency towards crystalline chitin

Chitin deacetylase (CDA) catalyzing the deacetylation of crystal chitin is a crucial step in the biosynthesis of chitosan, and also a scientific problem to be solved, which restricts the high-value utilization of chitin resources. This study aims to improve the catalytic efficiency of AsCDA from Acinetobacter schindleri MCDA01 by a semi-rational design using alanine scanning mutagenesis and saturation mutagenesis. The quadruple mutant M11 displayed a 2.31 and 1.73-fold improvement in kcat/Km and specific activity over AsCDA, which can remove 68 % of the acetyl groups from α-chitin. Furthermore, structural analysis suggested that additional hydrogen bonds, contributing the flexibility of amino acids and increasing the negative charge in M11 increased the catalytic efficiency. The microstructure changes of α-chitin pretreated by the mutant M11 were observed and evaluated using ¹³C CP/MAS NMR spectroscopy, FT-IR spectroscopy, XRD and SEM, and the results showed that M11 more efficiently catalyzed the release of acetyl groups from α-chitin. This study would provide a theoretical basis for the molecular modification of CDAs and accelerate the process of industrial production of chitosan by CDAs.

Hexaconazole exposure disrupt acetylcholinesterase, leading to mental illness

Hexaconazole is widely used in agricultural work, and it has been observed that it has potential to disrupt endocrine function and it has also capacity of bioaccumulation. In this study, we examined how the hexaconazole disrupts the usual balance of acetylcholinesterase. It has been already reported that heavy pesticide exposures may be a reason for several mental illnesses because these pesticides may disrupt normal balance of acetylcholinesterase. In this paper, we have done a complete molecular and dynamics analysis to understand the behavior of hexaconazole with acetylcholinesterase so that its toxicological aspect may be explored. Our findings revealed that hexaconazole has potency to interact with acetylcholinesterase in a stable manner. The binding energy of hexaconazole was found to be -7.95 kcal/mol. However, chlorpyrifos, known inhibitors of acetylcholinesterase, has binding energy of -7.17 kcal/mol. With respect to stability analysis, hexaconazole has similar stability like chlorpyrifos. Root-mean-square deviation, root-mean-square fluctuation, radius of gyration, hydrogen bonding, and solvent accessible surface area were similar to chlorpyrifos. In addition, density functional theory computations analysis reveals that hexaconazole is energetically stable like chlorpyrifos, which is necessary for establishing a stable ligand-protein complex. The result of this complete molecular analysis reveals that hexaconazole may disrupt the acetylcholinesterase balance, which leads to mental illness.

Screening of marine natural products for potential inhibitors targeting biotin biosynthesis pathway in Mycobacterium tuberculosis

Tuberculosis (TB) remains as one of the major public health concerns worldwide. A successful TB control and treatment is very challenging, due to continuing emergence of Mycobacterium tuberculosis strains resistant to known drugs. Therefore, the development of new drugs with different chemical and biological approaches is necessary to obtain more efficient anti-tubercular therapeutics. Biotin is an essential cofactor for lipid biosynthesis and gluconeogenesis in M. tuberculosis. M. tuberculosis relies on de novo biotin biosynthesis to obtain this vital cofactor since it cannot scavenge sufficient biotin from a mammalian host. In this study, comprehensive in silico methods including structure-based virtual screening, molecular docking, and molecular dynamic simulation analysis for ∼8000 marine natural products were performed against two essential enzymes involved in biotin synthesis and ligation of M. tuberculosis namely, pyridoxal 5'-phosphate-dependent transaminase (BioA) and mycobacterial biotin protein ligase (MtBPL). Two compounds; CMNPD10112 and CMNPD10113 are unveiled to bind the enzymes consistently and with high affinities. The binding pattern of compounds is further noticed in very stable binding modes as analyzed by molecular dynamics simulation and the mean RMSD of the complexes is within 4 Å. The intermolecular binding free energies validated complexes are less than -40 kcal/mol, which demonstrates strong and stable complexes formation. The identified hit compounds could be seeds for design of effective anti-mycobacterium therapeutics by inhibition of bacterial growth through blocking the biotin biosynthesis.Communicated by Ramaswamy H. Sarma.

Identification of essential oil phytocompounds as natural inhibitors of Odorant-binding protein to prevent malaria through in silico approach

Malaria predominantly affects millions annually in the African and Asian tropical and subtropical countries. With no effective vaccine, malaria prevention is exclusively dependent on preventing human-vector interaction. Anopheles gambiae, the main vector of the malaria parasite Plasmodium falciparum contains Odorant Binding proteins (OBPs) which are considered an attractive drug target for anti-malarial therapy. To identify a potential anti-malarial compound, we performed a structure-based screening of 876 phytocompounds derived from essential oils against the OBP4 by molecular docking. The compounds having better docking scores were assessed for drug-likeness, toxicity, and molecular interaction analysis. As per the results, strong affinities and high stability were demonstrated by two phytocompounds viz. Alpha-cyperone (-8.1 kcal mol-1) and Humulene oxide (-8.1 kcal mol-1) with OBP4. The hydrophobic interactions involve Phe123, Ala106, Thr57, Ala52, Thr69, and Ile64 within the binding cavities, which may block the OBP4 receptor resulting in disorientation. After that, the potential compounds were subjected to molecular dynamics (MD) simulation to evaluate their structural stability and dynamics at the active site of OBP4. The MM-PBSA result revealed that Alpha-cyperone and Humulene oxide had binding free energy of -92.44 kJ mol-1 and -113.25 kJ mol-1, respectively. Simulation outcomes demonstrate that these phytocompounds displayed considerable significant structural and pharmacological properties. The LD50 value of Alpha-cyperone and Humulene oxide also suggested that both are safe and suitable for use in natural repellent development. We suggest that the use of these compounds can minimize the treatment period and the various side effects associated with the currently available anti-malarial drugs.Communicated by Ramaswamy H. Sarma.

Pharmacoinformatics and Breed-Based De Novo Hybridization Studies to Develop New Neuraminidase Inhibitors as Potential Anti-Influenza Agents

Influenza represents a profoundly transmissible viral ailment primarily afflicting the respiratory system. Neuraminidase inhibitors constitute a class of antiviral therapeutics employed in the management of influenza. These inhibitors impede the liberation of the viral neuraminidase protein, thereby impeding viral dissemination from the infected cell to host cells. As such, neuraminidase has emerged as a pivotal target for mitigating influenza and its associated complications. Here, we apply a de novo hybridization approach based on a breed-centric methodology to elucidate novel neuraminidase inhibitors. The breed technique amalgamates established ligand frameworks with the shared target, neuraminidase, resulting in innovative inhibitor constructs. Molecular docking analysis revealed that the seven synthesized breed molecules (designated Breeds 1-7) formed more robust complexes with the neuraminidase receptor than conventional clinical neuraminidase inhibitors such as zanamivir, oseltamivir, and peramivir. Pharmacokinetic evaluations of the seven breed molecules (Breeds 1-7) demonstrated favorable bioavailability and optimal permeability, all falling within the specified parameters for human application. Molecular dynamics simulations spanning 100 nanoseconds corroborated the stability of these breed molecules within the active site of neuraminidase, shedding light on their structural dynamics. Binding energy assessments, which were conducted through MM-PBSA analysis, substantiated the enduring complexes formed by the seven types of molecules and the neuraminidase receptor. Last, the investigation employed a reaction-based enumeration technique to ascertain the synthetic pathways for the synthesis of the seven breed molecules.

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Discovery of novel potent drugs for influenza by inhibiting the vital function of neuraminidase via fragment-based drug design (FBDD) and molecular dynamics simulation strategies

The current work describes a fragment linking methodology to generate new neuraminidase inhibitors. A total number of 28,977 fragments from Zinc 20 have been obtained and screened for neuraminidase receptor affinity. Using Schrödinger software, the highest-scoring 270 fragment hits (with scores greater than -7.6) were subjected to fragment combining to create 100 new molecules. These 100 novel compounds were studied using XP docking to evaluate the molecular interaction modes and their binding affinity to neuraminidase receptor. The top ten molecules were selected, for ADMET, drug-likeness features. Based on these characteristics, the best four developed molecules and Zanamivir were submitted to a molecular dynamics simulation investigation to estimate their dynamics within the neuraminidase receptor using Gromacs software. All MD simulation findings show that the generated complexes are very stable when compared to the clinical inhibitor (Zanamivir). In addition, the four designed neuraminidase inhibitors formed very stable complexes with neuraminidase receptor (with total binding energies ranging from -83.50 to -107.85 Kj/mol) according to the total binding energy calculated by MM-PBSA. For the objective of developing new influenza medications, these novel molecules have the potential to be further evaluated in vitro and in vivo for influenza drug discovery.Communicated by Ramaswamy H. Sarma.

Thiadiazine thione derivatives as anti-leishmanial agents: synthesis, biological evaluation, structure activity relationship, ADMET, molecular docking and molecular dynamics simulation studies

During last decades, 3,5-disubstituted-tetrahydro-2H-thiadiazine-2-thione scaffold remains the center of interest due to their ease of preparation, diverse range substituents at N-3 and N-5 positions, and profound biological activities. In the current study, a series of 3,5-disubstituted-tetrahydro-2H-thiadiazine-2-thiones were synthesized in good to excellent yield, and the structure of the compounds were confirmed by various spectroscopic techniques such as FTIR, 1H-NMR, 13C-NMR and Mass spectrometry, and finally evaluated against Leishmania major. Whereas, all the evaluated compounds (1-33), demonstrate potential leishmanicidal activities with IC50 values in the range of (1.30- 149.98 uM). Among the evaluated compounds such as 3, 4, 6, and 10 exhibited excellent leishmanicidal activities with IC50 values of (2.17 μM), (2.39 μM), (2.00 μM), and (1.39 μM), respectively even better than the standard amphotericin B (IC50 = 0.50) and pentamidine (IC50 = 7.52). In order to investigate binding interaction of the most active compounds, molecular docking study was conducted with Leishmania major. Further molecular dynamic simulation study was also carried out to assess the stability and correct binding of the most active compound 10, within active site of the Leishamania major. Likewise, the physiochemical properties, drug likeness, and ADMET of the most active compounds were investigated, it was found that none of the compounds violate Lipiniski's rule of five, which show that this class of compounds had enough potential to be used as drug candidate in near future.Communicated by Ramaswamy H. Sarma.

Exploring and targeting potential druggable antimicrobial resistance targets ArgS, SecY, and MurA in Staphylococcus sciuri with TCM inhibitors through a subtractive genomics strategy

Staphylococcus sciuri (also currently Mammaliicoccus sciuri) are anaerobic facultative and non-motile bacteria that cause significant human pathogenesis such as endocarditis, wound infections, peritonitis, UTI, and septic shock. Methicillin-resistant S. sciuri (MRSS) strains also infects animals that include healthy broilers, cattle, dogs, and pigs. The emergence of MRSS strains thereby poses a serious health threat and thrives the scientific community towards novel treatment options. Herein, we investigated the druggable genome of S. sciuri by employing subtractive genomics that resulted in seven genes/proteins where only three of them were predicted as final targets. Further mining the literature showed that the ArgS (WP_058610923), SecY (WP_058611897), and MurA (WP_058612677) are involved in the multi-drug resistance phenomenon. After constructing and verifying the 3D protein homology models, a screening process was carried out using a library of Traditional Chinese Medicine compounds (consisting of 36,043 compounds). The molecular docking and simulation studies revealed the physicochemical stability parameters of the docked TCM inhibitors in the druggable cavities of each protein target by identifying their druggability potential and maximum hydrogen bonding interactions. The simulated receptor-ligand complexes showed the conformational changes and stability index of the secondary structure elements. The root mean square deviation (RMSD) graph showed fluctuations due to structural changes in the helix-coil-helix and beta-turn-beta changes at specific points where the pattern of the RMSD and root mean square fluctuation (RMSF) (< 1.0 Å) support any major domain shifts within the structural framework of the protein-ligand complex and placement of ligand was well complemented within the binding site. The β-factor values demonstrated instability at few points while the radius of gyration for structural compactness as a time function for the 100-ns simulation of protein-ligand complexes showed favorable average values and denoted the stability of all complexes. It is assumed that such findings might facilitate researchers to robustly discover and develop effective therapeutics against S. sciuri alongside other enteric infections.

In vitro antitumor activity, molecular dynamics simulation, DFT study, ADME prediction, and Eg5 binding of enastron analogues

The objective of this study is to evaluate a series of molecules based on cyclosulfamide as potential anticancer agents. Additionally, the study aims to analyze the obtained results through in silico studies; by conducting experiments and utilizing theoretical methods. In this context, we investigated the cytotoxic activity of enastron analogues on three human cell lines PRI (lymphoblastic cell line) derived from B-cell lymphoma. JURKAT (ATCC TIB-152) acute T cell leukaemia and K562 (ATCC CLL-243) is a chronic myelogenous leukaemia. Most of the tested compounds showed good inhibitory activity compared with the reference ligand (chlorambucil). The 5a derivative demonstrated the strongest effect against all cancer cells used. Furthermore, molecular docking simulations of the Eg5-enastron analogue complex revealed that studied molecules have the ability to inhibit the Eg5 enzyme, as evidenced by their calculated docking score. Following the promising results from the molecular docking study, the complex Eg5-4a underwent a 100 ns molecular dynamics simulation using Desmond. During the simulation, the receptor-ligand pairing demonstrated substantial stability after the initial 70 ns. In addition, we used DFT calculations to analyze the electronic and geometric characteristics of the studied compounds. The HOMO and LUMO band gap energies, and the molecular electrostatic potential surface were also deducted for the stable structure of each compound. Also, we studied the prediction of absorption, distribution, metabolism and excretion (ADME) of the compounds.

Identification of Potent hDHODH Inhibitors for Lung Cancer via Virtual Screening of a Rationally Designed Small Combinatorial Library

Cancer is characterized by altered cellular metabolism, and metabolic enzymes are considered as a promising target for anticancer therapy. Pyrimidine metabolism dysregulation is associated with various types of cancer, particularly lung cancer, which is one of the leading causes of cancer-related mortality worldwide. Recent studies have shown that small-cell lung cancer cells are particularly reliant on the pyrimidine biosynthesis pathway and are sensitive to its disruption. DHODH, the rate-limiting enzyme of the de novo pyrimidine production pathway, is essential in the production of RNA and DNA and is overexpressed in malignancies such as AML, skin cancer, breast cancer, and lung cancer, thereby highlighting DHODH as a viable target for developing drugs to combat lung cancer. Herein, rational drug design and computational techniques were used to discover novel DHODH inhibitors. A small combinatorial library was generated, and the top hits were synthesized and tested for anticancer activity against three lung cancer cell lines. Among the tested compounds, compound 5c possessed a stronger cytotoxicity (TC₅₀ of 11 μM) compared to the standard FDA-approved drug (Regorafenib, TC₅₀ of 13 μM) on the A549 cell line. Furthermore, compound 5c demonstrated potent inhibitory activity against hDHODH at a nanomolar level of 421 nM. DFT, molecular docking, molecular dynamic simulations, and free energy calculations were also carried out to understand the inhibitory mechanisms of the synthesized scaffolds. These in silico studies identified key mechanisms and structural features that will be crucial for future studies.

Structural Dynamics of the Precatalytic State of Human Cytochrome c upon T28C, G34C, and A50C Mutations: A Molecular Dynamics Simulation Perspective

The native structure of cytochrome c (cytc) contains hexacoordinate heme iron with His18 and Met80 residues ligated at the axial sites. Mutations of cytc at Ω-loops have been investigated in modulating the peroxidase activity and, hence, related to the initiation of the apoptotic pathway. Our previous experimental data reported on the peroxidase activity of the cysteine-directed mutants at different parts of the Ω-loop of human cytc (hCytc), that is, T28C, G34C, and A50C. In this work, we performed 1 μs molecular dynamics (MD) simulations to elucidate the detailed structural and dynamic changes upon these mutations, particularly at the proximal Ω-loop. The structures of hCytc were modeled in the hexacoordinated form, which was referred to as the "precatalytic state". The results showed that the structural features of the G34C mutant were more distinctive than those of other mutants. G34C mutation caused local destabilization and flexibility at the proximal Ω-loop (residues 12-28) and an extended distance between this Ω-loop region and heme iron. Besides, analysis of the orientation of the Arg38 side chain of the G34C mutant revealed the Arg38 conformer facing away from the heme iron. The obtained MD results also suggested structural diversity of the precatalytic states for the three hCytc mutants, specifically the effect of G34C mutation on the flexibility of the proximal Ω-loops. Therefore, our MD simulations combined with previous experimental data provide detailed insights into the structural basis of hCytc that could contribute to its pro-apoptotic function.

Association of peptidyl prolyl cis/trans isomerase Rrd1 with C terminal domain of RNA polymerase II

The largest subunit of RNAPII extends as the conserved unstructured heptapeptide consensus repeats Y₁S₂P₃T₄S₅P₆S₇ and their posttranslational modification, especially the phosphorylation state at Ser2, Ser5 and Ser7 of CTD recruits different transcription factors involved in transcription. In the current study, fluorescence anisotropy, pull down assay and molecular dynamics simulation studies employed to conclude that peptidyl-prolyl cis/trans-isomerase Rrd1 has strong affinity for unphosphorylated CTD rather than phosphorylated CTD for mRNA transcription. Rrd1 preferentially interacts with unphosphorylated GST-CTD in comparison to hyperphosphorylated GST-CTD in vitro. Fluorescence anisotropy revealed that recombinant Rrd1 prefers to bind unphosphorylated CTD peptide in comparison to phosphorylated CTD peptide. In computational studies, the RMSD of Rrd1-unphosphorylated CTD complex was greater than the RMSD of Rrd1-pCTD complex. During 50 ns MD simulation run Rrd1-pCTD complex get dissociated twice viz. 20 ns to 30 ns and 40 ns to 50 ns, while Rrd1-unpCTD complex remain stable throughout the process. Additionally, the Rrd1-unphosphorylated CTD complexes acquire comparatively higher number of H-bonds, water bridges and hydrophobic interactions occupancy than Rrd1-pCTD complex, concludes that the Rrd1 interacts more strongly with the unphosphorylated CTD than the pCTD.

Identification of Natural Compounds of the Apple as Inhibitors against Cholinesterase for the Treatment of Alzheimer's Disease: An In Silico Molecular Docking Simulation and ADMET Study

Alzheimer's disease (AD), the most common type of dementia in older people, causes neurological problems associated with memory and thinking. The key enzymes involved in Alzheimer's disease pathways are acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). Because of this, there is a lot of interest in finding new AChE inhibitors. Among compounds that are not alkaloids, flavonoids have stood out as good candidates. The apple fruit, Malus domestica (Rosaceae), is second only to cranberries regarding total phenolic compound concentration. Computational tools and biological databases were used to investigate enzymes and natural compounds. Molecular docking techniques were used to analyze the interactions of natural compounds of the apple with enzymes involved in the central nervous system (CNS), acetylcholinesterase, and butyrylcholinesterase, followed by binding affinity calculations using the AutoDock tool. The molecular docking results revealed that CID: 107905 exhibited the best interactions with AChE, with a binding affinity of -12.2 kcal/mol, and CID: 163103561 showed the highest binding affinity with BuChE, i.e., -11.2 kcal/mol. Importantly, it was observed that amino acid residue Trp286 of AChE was involved in hydrogen bond formation, Van Der Walls interactions, and Pi-Sigma/Pi-Pi interactions in the studied complexes. Moreover, the results of the Molecular Dynamics Simulation (MDS) analysis indicated interaction stability. This study shows that CID: 12000657 could be used as an AChE inhibitor and CID: 135398658 as a BuChE inhibitor to treat Alzheimer's disease and other neurological disorders.

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.

Spectroscopic analysis, kinetic mechanism, computational docking, and molecular dynamics of active metabolites from the aerial parts of Astragalus membranaceusBunge as tyrosinase inhibitors

A new isoflavane derivative (2), a new natural isoflavane (6), four new oleanane-type triterpenoid saponins (23, 25, 28, and 29), and twenty three known secondary metabolites (1, 3-5, 7-22, 24, 26, and 27) were isolated from the aerial parts of Astragalus membranaceus Bunge. The chemical structures of these compounds were elucidated through spectroscopic analysis and compared with those identified in previous studies. Tyrosinase inhibition ability of isolated compounds (1-29) was evaluated. Of these, compounds 3, 4, 6, and 14 exhibited inhibitory effects, with IC₅₀ values ranging from 24.6 to 59.2 μM. According to kinetic analysis, compounds 3 and 4 were non-competitive inhibitors of tyrosinase, whereas compounds 6 and 14 inhibited tyrosinase in uncompetitive and competitive modes, respectively. Molecular docking analysis identified that compounds 3, 4, and 6 could bind to allosteric sites and compound 14 could bind to the catalytic site of tyrosinase, which is consistent with the results of kinetic studies. Molecular dynamics behaviors of the active compounds in complex with tyrosinase were investigated via 60 ns simulation which demonstrated their high stability. These findings indicate that the aerial parts of A. membranaceus are a potential source of natural tyrosinase inhibitors.

Inhibition of microtubule affinity regulating kinase 4 by an acetylcholinesterase inhibitor, Huperzine A: Computational and experimental approaches

Microtubule affinity regulating kinase 4 (MARK4), 752 amino acids long, belonging to the AMPK superfamily, plays a vital role in regulating microtubules due to its potential to phosphorylate microtubule-associated proteins (MAP's) and thus, MARK4 plays a key role in Alzheimer's disease (AD) pathology. MARK4 is a druggable target for cancer, neurodegenerative diseases, and metabolic disorders. In this study, we have evaluated the MARK4 inhibitory potential of Huperzine A (HpA), an acetylcholinesterase inhibitor (AChEI), a potential AD drug. Molecular docking revealed the key residues governing the MARK4-HpA complex formation. The structural stability and conformational dynamics of the MARK4-HpA complex was assessed by employing Molecular dynamics (MD) simulation. The results suggested that the binding of HpA with MARK4 leads to minimal structural alterations in the native conformation of MARK4, implying the stability of the MARK4-HpA complex. Isothermal titration calorimetry (ITC) studies deciphered that HpA binds to MARK4 spontaneously. Moreover, the kinase assay depicted significant inhibition of MARK by HpA (IC₅₀ = 4.91 μM), implying it to be a potent MARK4 inhibitor that can be implicated in the treatment of MARK4-directed diseases.

Novel phytochemical inhibitors targeting monkeypox virus thymidine and serine/threonine kinase: integrating computational modeling and molecular dynamics simulation

Due to the rapid spread of the monkeypox virus and rise in the number of cases, there is an urgent need for the development of effective drugs against the infection. Serine/threonine protein kinase (Ser/Thr kinase) and Thymidine Kinase (TK) plays an imperative role in the replication and virulence of monkeypox virus and thus is deliberated as an attractive target in anti-viral drug development. In the present study, the 3D structure of monkeypox virus Ser/Thr kinase and TK was generated via molecular modeling techniques and performed their thorough structural analysis. We have screened potent anti-viral phytochemicals from the literature to inhibit Ser/Thr kinase and TK. As part of the initial screening, the physicochemical properties of the compounds were examined. Following this, a structure-based molecular docking technique was used to select compounds based on their binding affinity towards Ser/Thr kinase and TK. In order to find more potent hits against Ser/Thr kinase and TK, further examinations of ADMET properties, PAINS patterns and blood-brain barrier permeability were conducted. As a result, thalimonine and galanthamine were identified from the screening process bearing appreciable binding affinity towards Ser/Thr kinase and TK respectively, which showed a worthy set of drug-like properties. In the end, molecular dynamics simulations were performed for 100 ns, which showed decent stability of both protein-ligand complex throughout the trajectory. Due to the possibility that both monkeypox virus target proteins may be inhibited by thalimonine and galanthamine, our study highlights the need to investigate in vivo effects of thalimonine and galanthamine.Communicated by Ramaswamy H. Sarma.

Bioremediation potential of laccase for catalysis of glyphosate, isoproturon, lignin, and parathion: Molecular docking, dynamics, and simulation

The present study aimed to investigate the catalytic degradation produced by laccase in the detoxification of glyphosate, isoproturon, lignin polymer, and parathion. We explored laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion using molecular docking (MD) and molecular dynamics simulation (MDS) approaches. The results suggest that laccase interacts well with glyphosate, lignin polymer, isoproturon, and parathion during biodegradation. We calculated the root mean square deviations (RMSD) of laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion as 0.24 ± 0.02, 0.59 ± 0.32, 0.43 ± 0.07, and 0.43 ± 0.06 nm, respectively. In an aqueous solution, the stability of laccase with glyphosate, lignin polymer, isoproturon, and parathion is mediated through the formation of hydrophobic interactions, hydrogen bonds, and van der Waals interactions. The presence of xenobiotic toxic compounds in the active site changed the conformation of laccase. MDS of the laccase-substrate complexes confirmed their stability during catalytic degradation. Laccase assay results confirmed that the degradation of syringol, dihydroconiferyl alcohol, guaiacol, parathion, isoproturon, and glyphosate were 100%, 99.31%, 95.69%, 60.96%, 54.51%, and 48.34% within 2 h, respectively. Taken together, we describe a novel method to understand the molecular-level biodegradation of xenobiotic compounds through laccase and its potential application in contaminant removal.

In silico protein engineering shows that novel mutations affecting NAD+ binding sites may improve phosphite dehydrogenase stability and activity

Pseudomonas stutzeri phosphite dehydrogenase (PTDH) catalyzes the oxidation of phosphite to phosphate in the presence of NAD, resulting in the formation of NADH. The regeneration of NADH by PTDH is greater than any other enzyme due to the substantial change in the free energy of reaction (G°' = - 63.3 kJ/mol). Presently, improving the stability of PTDH is for a great importance to ensure an economically viable reaction process to produce phosphite as a byproduct for agronomic applications. The binding site of NAD⁺ with PTDH includes thirty-four residues; eight of which have been previously mutated and characterized for their roles in catalysis. In the present study, the unexplored twenty-six key residues involved in the binding of NAD⁺ were subjected to in silico mutagenesis based on the physicochemical properties of the amino acids. The effects of these mutations on the structure, stability, activity, and interaction of PTDH with NAD⁺ were investigated using molecular docking, molecular dynamics simulations, free energy calculations, and secondary structure analysis. We identified seven novel mutations, A155I, G157I, L217I, P235A, V262I, I293A, and I293L, that reduce the compactness of the protein while improving PTDH stability and binding to NAD⁺.

In silico investigation of the interactions of certain drugs proposed for the treatment of Covid-19 with the paraoxonase-1

Coronavirus disease 2019 (Covid-19) has caused one of the biggest pandemics of modern times, infected over 240 million people and killed over 4.9 million people, and continues to do so. Although many drugs are widely recommended in the treatment of this disease, the interactions of these drugs with an anti-atherosclerotic enzyme, paraoxonase-1 (PON1), are not well known. In our study, we investigated the interactions of 18 different drugs, which are claimed to be effective against covid-19, with the PON1 enzyme and its genetics variants L55M and Q192R with molecular docking, molecular dynamics simulation and free energy calculation method MM/PBSA. We found that ruxolitinib, dexamethasone, colchicine; dexamethasone, sitagliptin, baricitinib and galidesivir, ruxolitinib, hydroxychloroquine were the most effective compounds in binding PON1-w, PON1L55M and PON1Q192R respectively. Mainly, sitagliptin, galidesivir and hydroxychloroquine have attracted attention by showing very high affinity (<-300 kJ/mol) according to the MM/PBSA method. We concluded that the drug interactions should be considered and more attention should be paid in the use of these drugs.Communicated by Ramaswamy H. Sarma.