Displaying publications 1 - 20 of 224 in total

Abstract:
Sort:
  1. Effects of anthocyanidins on the conformational transition of Aβ(1-42) peptide: Insights from molecular docking and molecular dynamics simulations
    Zakaria N, Wan Harun WMRS, Mohammad Latif MA, Azaman SNA, Abdul Rahman MB, Faujan NH
    J Mol Graph Model, 2024 Jun;129:108732.
    PMID: 38412813 DOI: 10.1016/j.jmgm.2024.108732
    Recent evidence from in vitro and in vivo studies has shown that anthocyanins and anthocyanidins can reduce and inhibit the amyloid beta (Aβ) species, one of the hallmarks of Alzheimer's disease (AD). However, their inhibition mechanisms on Aβ species at molecular details remain elusive. Therefore, in the present study, molecular modelling methods were employed to investigate their inhibitory mechanisms on Aβ(1-42) peptide. The results highlighted that anthocyanidins effectively inhibited the conformational transitions of helices into beta-sheet (β-sheet) conformation within Aβ(1-42) peptide by two different mechanisms: 1) the obstruction of two terminals from coming into contact due to the binding of anthocyanidins with residues of N- and second hydrophobic core (SHC)-C-terminals, and 2) the prevention of the folding process due to the binding of anthocyanidin with the central polar (Asp23 and Lys28) and native helix (Asp23, Lys28, and Leu34) residues. These new findings on the inhibition of β-sheet formation by targeting both N- and SHC-C-terminals, and the long-established target, D23-K28 salt bridge residues, not with the conventional central hydrophobic core (CHC) as reported in the literature, might aid in designing more potent inhibitors for AD treatment.
    Matched MeSH terms: Molecular Dynamics Simulation
  2. Probing the anti-Aβ42 aggregation and protective effects of prenylated xanthone against Aβ42-induced toxicity in transgenic Caenorhabditis elegans model
    Thew HY, Boon Keat K, Tan YC, Ong YS, Parat MO, Murugaiyah V, et al.
    Chem Biol Interact, 2024 May 01;394:110978.
    PMID: 38552766 DOI: 10.1016/j.cbi.2024.110978
    Alzheimer's disease (AD) is a progressive neurodegenerative disorder characterized by the accumulation of amyloid-β (Aβ) protein aggregates, leading to synaptic dysfunction and neuronal cell death. In this study, we used a comprehensive approach encompassing in vitro assays, computational analyses, and an in vivo Caenorhabditis elegans model to evaluate the inhibitory effects of various xanthones, focusing on Garcinone D (GD), on Aβ42 oligomer formation. Dot blot analysis revealed concentration-dependent responses among xanthones, with GD consistently inhibiting Aβ42 oligomer formation at low concentrations (0.1 and 0.5 μM, inhibitions of 84.66 ± 2.25% and 85.06 ± 6.57%, respectively). Molecular docking and dynamics simulations provided insights into the molecular interactions between xanthones and Aβ42, highlighting the disruption of key residues involved in Aβ42 aggregation. The neuroprotective potential of GD was established using transgenic C. elegans GMC101, with substantial delays in paralysis reported at higher concentrations. Our findings show that GD is a potent suppressor of Aβ42 oligomer formation, suggesting its potential as a therapeutic candidate for AD. The concentration-dependent effects observed in both in vitro and in vivo models underscore the need for nuanced dose-response assessments. These findings contribute novel insights into the therapeutic landscape of xanthones against AD, emphasizing the multifaceted potential of GD for further translational endeavors in neurodegenerative disorder research.
    Matched MeSH terms: Molecular Dynamics Simulation
  3. Preparation and Preliminary Structure-Activity Relationship Studies of Schwarzinicine A Analogs as Vasorelaxant Agents
    Lee FK, Chan NJ, Krishnan P, Datu Abdul Salam DS, Chee XW, Muhamad A, et al.
    J Nat Prod, 2024 Apr 26;87(4):675-691.
    PMID: 38442031 DOI: 10.1021/acs.jnatprod.3c00707
    Schwarzinicines A-D, a series of alkaloids recently discovered from Ficus schwarzii, exhibit pronounced vasorelaxant activity in rat isolated aorta. Building on this finding, a concise synthesis of schwarzinicines A and B has been reported, allowing further investigations into their biological properties. Herein, a preliminary exploration of the chemical space surrounding the structure of schwarzinicine A (1) was carried out aiming to identify structural features that are essential for vasorelaxant activity. A total of 57 analogs were synthesized and tested for vasorelaxant activity in rat isolated aorta. Both efficacy (Emax) and potency (EC50) of these analogs were compared. In addition to identifying structural features that are required for activity or associated with potency enhancement effect, four analogs showed significant potency improvements of up to 40.2-fold when compared to 1. Molecular dynamics simulation of a tetrameric 44-bound transient receptor potential canonical-6 (TRPC6) protein indicated that 44 could potentially form important interactions with the residues Glu509, Asp530, Lys748, Arg758, and Tyr521. These results may serve as a foundation for guiding further structural optimization of the schwarzinicine A scaffold, aiming to discover even more potent analogs.
    Matched MeSH terms: Molecular Dynamics Simulation
  4. Proposed molecular mechanism of non-competitive inhibition using molecular dynamics simulations between α-glucosidase enzyme and mangostin compound as antidiabetic
    Maulana AF, Maksum IP, Sriwidodo S, Rukayadi Y
    J Mol Model, 2024 Apr 18;30(5):136.
    PMID: 38634946 DOI: 10.1007/s00894-024-05934-z
    CONTEXT: Further understanding of the molecular mechanisms is necessary since it is important for designing new drugs. This study aimed to understand the molecular mechanisms involved in the design of drugs that are inhibitors of the α-glucosidase enzyme. This research aims to gain further understanding of the molecular mechanisms underlying antidiabetic drug design. The molecular docking process yielded 4 compounds with the best affinity energy, including γ-Mangostin, 1,6-dimethyl-ester-3-isomangostin, 1,3,6-trimethyl-ester-α-mangostin, and 3,6,7-trimethyl-ester-γ-mangostin. Free energy calculation with molecular mechanics with generalized born and surface area solvation indicated that the 3,6,7-trimethyl-γ-mangostin had a better free energy value compared to acarbose and simulated maltose together with 3,6,7-trimethyl-γ-mangostin compound. Based on the analysis of electrostatic, van der Waals, and intermolecular hydrogen interactions, 3,6,7-trimethyl-γ-mangostin adopts a noncompetitive inhibition mechanism, whereas acarbose adopts a competitive inhibition mechanism. Consequently, 3,6,7-trimethyl-ester-γ-mangostin, which is a derivative of γ-mangostin, can provide better activity in silico with molecular docking approaches and molecular dynamics simulations.

    METHOD: This research commenced with retrieving protein structures from the RCSB database, generating the formation of ligands using the ChemDraw Professional software, conducting molecular docking with the Autodock Vina software, and performing molecular dynamics simulations using the Amber software, along with the evaluation of RMSD values and intermolecular hydrogen bonds. Free energy, electrostatic interactions, and Van der Waals interaction were calculated using MM/GBSA. Acarbose, used as a positive control, and maltose are simulated together with test compound that has the best free energy. The forcefields used for molecular dynamics simulations are ff19SB, gaff2, and tip3p.

    Matched MeSH terms: Molecular Dynamics Simulation
  5. Functional expression, purification, biochemical and biophysical characterizations, and molecular dynamics simulation of a histidine acid phosphatase from Saccharomyces cerevisiae
    Nezhad NG, Jamaludin SZB, Rahman RNZRA, Yahaya NM, Oslan SN, Shariff FM, et al.
    World J Microbiol Biotechnol, 2024 Apr 17;40(6):171.
    PMID: 38630327 DOI: 10.1007/s11274-024-03970-8
    A histidine acid phosphatase (HAP) (PhySc) with 99.50% protein sequence similarity with PHO5 from Saccharomyces cerevisiae was expressed functionally with the molecular mass of ∼110 kDa through co-expression along with the set of molecular chaperones dnaK, dnaJ, GroESL. The purified HAP illustrated the optimum activity of 28.75 ± 0.39 U/mg at pH 5.5 and 40 ˚C. The Km and Kcat values towards calcium phytate were 0.608 ± 0.09 mM and 650.89 ± 3.6 s- 1. The half-lives (T1/2) at 55 and 60 ˚C were 2.75 min and 55 s, respectively. The circular dichroism (CD) demonstrated that PhySc includes 30.5, 28.1, 21.3, and 20.1% of random coils, α-Helix, β-Turns, and β-Sheet, respectively. The Tm recorded by CD for PhySc was 56.5 ± 0.34˚C. The molecular docking illustrated that His59 and Asp322 act as catalytic residues in the PhySc. MD simulation showed that PhySc at 40 ˚C has higher structural stability over those of the temperatures 60 and 80 ˚C that support the thermodynamic in vitro investigations. Secondary structure content results obtained from MD simulation indicated that PhySc consists of 34.03, 33.09, 17.5, 12.31, and 3.05% of coil, helix, turn, sheet, and helix310, respectively, which is almost consistent with the experimental results.
    Matched MeSH terms: Molecular Dynamics Simulation*
  6. Exploring the binding mechanisms of thermally and ultrasonically induced molten globule-like β-lactoglobulin with heptanal as revealed by multi-spectroscopic techniques and molecular simulation
    Han C, Zheng Y, Huang S, Xu L, Zhou C, Sun Y, et al.
    Int J Biol Macromol, 2024 Apr;263(Pt 1):130300.
    PMID: 38395276 DOI: 10.1016/j.ijbiomac.2024.130300
    This work employed the model protein β-lactoglobulin (BLG) to investigate the contribution of microstructural changes to regulating the interaction patterns between protein and flavor compounds through employing computer simulation and multi-spectroscopic techniques. The formation of molten globule (MG) state-like protein during the conformational evolution of BLG, in response to ultrasonic (UC) and heat (HT) treatments, was revealed through multi-spectroscopic characterization. Differential MG structures were distinguished by variations in surface hydrophobicity and the microenvironment of tryptophan residues. Fluorescence quenching measurements indicated that the formation of MG enhanced the binding affinity of heptanal to protein. LC-MS/MS and NMR revealed the covalent bonding between heptanal and BLG formed by Michael addition and Schiff-base reactions, and MG-like BLG exhibited fewer chemical shift residues. Molecular docking and molecular dynamics simulation confirmed the synergistic involvement of hydrophobic interactions and hydrogen bonds in shaping BLG-heptanal complexes thus promoting the stability of BLG structures. These findings indicated that the production of BLG-heptanal complexes was driven synergistically by non-covalent and covalent bonds, and their interaction processes were influenced by processes-induced formation of MG potentially tuning the release and retention behaviors of flavor compounds.
    Matched MeSH terms: Molecular Dynamics Simulation
  7. Structural vaccinology, molecular simulation and immune simulation approaches to design multi-epitopes vaccine against John Cunningham virus
    Suleman M, Khan TA, Ejaz H, Maroof S, Alshammari A, Albekairi NA, et al.
    Microb Pathog, 2024 Apr;189:106572.
    PMID: 38354987 DOI: 10.1016/j.micpath.2024.106572
    The JCV (John Cunningham Virus) is known to cause progressive multifocal leukoencephalopathy, a condition that results in the formation of tumors. Symptoms of this condition such as sensory defects, cognitive dysfunction, muscle weakness, homonosapobia, difficulties with coordination, and aphasia. To date, there is no specific and effective treatment to completely cure or prevent John Cunningham polyomavirus infections. Since the best way to control the disease is vaccination. In this study, the immunoinformatic tools were used to predict the high immunogenic and non-allergenic B cells, helper T cells (HTL), and cytotoxic T cells (CTL) epitopes from capsid, major capsid, and T antigen proteins of JC virus to design the highly efficient subunit vaccines. The specific immunogenic linkers were used to link together the predicted epitopes and subjected to 3D modeling by using the Robetta server. MD simulation was used to confirm that the newly constructed vaccines are stable and properly fold. Additionally, the molecular docking approach revealed that the vaccines have a strong binding affinity with human TLR-7. The codon adaptation index (CAI) and GC content values verified that the constructed vaccines would be highly expressed in E. coli pET28a (+) plasmid. The immune simulation analysis indicated that the human immune system would have a strong response to the vaccines, with a high titer of IgM and IgG antibodies being produced. In conclusion, this study will provide a pre-clinical concept to construct an effective, highly antigenic, non-allergenic, and thermostable vaccine to combat the infection of the John Cunningham virus.
    Matched MeSH terms: Molecular Dynamics Simulation
  8. Fulltext Identification of Pyrazole Derivatives of Usnic Acid as Novel Inhibitor of SARS-CoV-2 Main Protease Through Virtual Screening Approaches
    Roney M, Singh G, Huq AKMM, Forid MS, Ishak WMBW, Rullah K, et al.
    Mol Biotechnol, 2024 Apr;66(4):696-706.
    PMID: 36752937 DOI: 10.1007/s12033-023-00667-5
    The infection produced by the SARS-CoV-2 virus remains a significant health crisis worldwide. The lack of specific medications for COVID-19 necessitates a concerted effort to find the much-desired therapies for this condition. The main protease (Mpro) of SARS-CoV-2 is a promising target, vital for virus replication and transcription. In this study, fifty pyrazole derivatives were tested for their pharmacokinetics and drugability, resulting in eight hit compounds. Subsequent molecular docking simulations on SARS-CoV-2 main protease afforded two lead compounds with strong affinity at the active site. Additionally, the molecular dynamics (MD) simulations of lead compounds (17 and 39), along with binding free energy calculations, were accomplished to validate the stability of the docked complexes and the binding poses achieved in docking experiments. Based on these findings, compound 17 and 39, with their favorable projected pharmacokinetics and pharmacological characteristics, are the proposed potential antiviral candidates which require further investigation to be used as anti-SARS-CoV-2 medication.
    Matched MeSH terms: Molecular Dynamics Simulation*
  9. Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes
    Pirojsirikul T, Lee VS, Nimmanpipug P
    Mol Biotechnol, 2024 Apr;66(4):582-591.
    PMID: 38374320 DOI: 10.1007/s12033-024-01082-0
    We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
    Matched MeSH terms: Molecular Dynamics Simulation*
  10. Special issue: Biomolecular Modelling and Simulation
    Choong YS, Mancera R, Lee VS
    Mol Biotechnol, 2024 Apr;66(4):567.
    PMID: 38337130 DOI: 10.1007/s12033-024-01073-1
    Matched MeSH terms: Molecular Dynamics Simulation*
  11. Fulltext Molecular Dynamics Simulations Reveal Novel Interacting Regions of Human Prion Protein to Brucella abortus Hsp60 Protein
    Le-Dao HA, Dinh TT, Tran TL, Lee VS, Tran-Van H
    Mol Biotechnol, 2024 Apr;66(4):687-695.
    PMID: 36633832 DOI: 10.1007/s12033-023-00655-9
    The distinctive morphology characteristics of microfold cells (M cells) allow the vaccine antigen not only to interact with immune cells directly, but also to effectively stimulate mucosal immune responses via receptors on its apical surface. Human prion protein, a transmembrane receptor for Brucella abortus Hsp60, is highly expressed on the M cell surface. Nonetheless, this protein tends to express in inclusion body in prokaryotic hosts. In this study, the shorter interacting regions of human prion protein were identified via computational methods such as docking and molecular dynamics simulations to minimize its aggregation tendency. The computational calculations revealed three novel human prion protein-interacting regions, namely PrP125, PrP174, and PrP180. In accordance with in silico prediction, the biologically synthesized peptides fusing with GST tag demonstrated their specific binding to Hsp60 protein via pull-down assay. Hence, this finding laid the groundwork for M-cell targeting candidate validation through these newly identified interacting regions.
    Matched MeSH terms: Molecular Dynamics Simulation*
  12. Immunoinformatics-driven In silico vaccine design for Nipah virus (NPV): Integrating machine learning and computational epitope prediction
    Shahab M, Iqbal MW, Ahmad A, Alshabrmi FM, Wei DQ, Khan A, et al.
    Comput Biol Med, 2024 Mar;170:108056.
    PMID: 38301512 DOI: 10.1016/j.compbiomed.2024.108056
    The Nipah virus (NPV) is a highly lethal virus, known for its significant fatality rate. The virus initially originated in Malaysia in 1998 and later led to outbreaks in nearby countries such as Bangladesh, Singapore, and India. Currently, there are no specific vaccines available for this virus. The current work employed the reverse vaccinology method to conduct a comprehensive analysis of the entire proteome of the NPV virus. The aim was to identify and choose the most promising antigenic proteins that could serve as potential candidates for vaccine development. We have also designed B and T cell epitopes-based vaccine candidate using immunoinformatics approach. We have identified a total of 5 novel Cytotoxic T Lymphocytes (CTL), 5 Helper T Lymphocytes (HTL), and 6 linear B-cell potential antigenic epitopes which are novel and can be used for further vaccine development against Nipah virus. Then we performed the physicochemical properties, antigenic, immunogenic and allergenicity prediction of the designed vaccine candidate against NPV. Further, Computational analysis indicated that these epitopes possessed highly antigenic properties and were capable of interacting with immune receptors. The designed vaccine were then docked with the human immune receptors, namely TLR-2 and TLR-4 showed robust interaction with the immune receptor. Molecular dynamics simulations demonstrated robust binding and good dynamics. After numerous dosages at varied intervals, computational immune response modeling showed that the immunogenic construct might elicit a significant immune response. In conclusion, the immunogenic construct shows promise in providing protection against NPV, However, further experimental validation is required before moving to clinical trials.
    Matched MeSH terms: Molecular Dynamics Simulation
  13. QSAR modeling approaches to identify a novel ACE2 inhibitor that selectively bind with the C and N terminals of the ectodomain
    Jawarkar RD, Zaki MEA, Al-Hussain SA, Al-Mutairi AA, Samad A, Mukerjee N, et al.
    J Biomol Struct Dyn, 2024 Mar;42(5):2550-2569.
    PMID: 37144753 DOI: 10.1080/07391102.2023.2205948
    Due to the high rates of drug development failure and the massive expenses associated with drug discovery, repurposing existing drugs has become more popular. As a result, we have used QSAR modelling on a large and varied dataset of 657 compounds in an effort to discover both explicit and subtle structural features requisite for ACE2 inhibitory activity, with the goal of identifying novel hit molecules. The QSAR modelling yielded a statistically robust QSAR model with high predictivity (R2tr=0.84, R2ex=0.79), previously undisclosed features, and novel mechanistic interpretations. The developed QSAR model predicted the ACE2 inhibitory activity (PIC50) of 1615 ZINC FDA compounds. This led to the detection of a PIC50 of 8.604 M for the hit molecule (ZINC000027990463). The hit molecule's docking score is -9.67 kcal/mol (RMSD 1.4). The hit molecule revealed 25 interactions with the residue ASP40, which defines the N and C termini of the ectodomain of ACE2. The HIT molecule conducted more than thirty contacts with water molecules and exhibited polar interaction with the ARG522 residue coupled with the second chloride ion, which is 10.4 nm away from the zinc ion. Both molecular docking and QSAR produced comparable findings. Moreover, MD simulation and MMGBSA studies verified docking analysis. The MD simulation showed that the hit molecule-ACE2 receptor complex is stable for 400 ns, suggesting that repurposed hit molecule 3 is a viable ACE2 inhibitor.
    Matched MeSH terms: Molecular Dynamics Simulation
  14. Interaction of Asn297-Linked Glycan Ligands with the Fc Fragment of the Immunoglobulin Class G1: A Molecular Dynamics Simulation Study
    Tong M, Liu P, Li C, Zhang Z, Sun W, Dong P, et al.
    J Chem Inf Model, 2024 Feb 12;64(3):785-798.
    PMID: 38262973 DOI: 10.1021/acs.jcim.3c01584
    The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.
    Matched MeSH terms: Molecular Dynamics Simulation*
  15. Carbazole derivatives as promising competitive and allosteric inhibitors of human serotonin transporter: computational pharmacology
    Ayipo YO, Ahmad I, Chong CF, Zainurin NA, Najib SY, Patel H, et al.
    J Biomol Struct Dyn, 2024;42(2):993-1014.
    PMID: 37021485 DOI: 10.1080/07391102.2023.2198016
    The human serotonin transporters (hSERTs) are neurotransmitter sodium symporters of the aminergic G protein-coupled receptors, regulating the synaptic serotonin and neuropharmacological processes related to neuropsychiatric disorders, notably, depression. Selective serotonin reuptake inhibitors (SSRIs) such as fluoxetine and (S)-citalopram are competitive inhibitors of hSERTs and are commonly the first-line medications for major depressive disorder (MDD). However, treatment-resistance and unpleasant aftereffects constitute their clinical drawbacks. Interestingly, vilazodone emerged with polypharmacological (competitive and allosteric) inhibitions on hSERTs, amenable to improved efficacy. However, its application usually warrants adjuvant/combination therapy, another subject of critical adverse events. Thus, the discovery of alternatives with polypharmacological potentials (one-drug-multiple-target) and improved safety remains essential. In this study, carbazole analogues from chemical libraries were explored using docking and molecular dynamics (MD) simulation. Selectively, two IBScreen ligands, STOCK3S-30866 and STOCK1N-37454 predictively bound to the active pockets and expanded boundaries (extracellular vestibules) of the hSERTs more potently than vilazodone and (S)-citalopram. For instance, the two ligands showed docking scores of -9.52 and -9.59 kcal/mol and MM-GBSA scores of -92.96 and -65.66 kcal/mol respectively compared to vilazodone's respective scores of -7.828 and -59.27 against the central active site of the hSERT (PDB 7LWD). Similarly, the two ligands also docked to the allosteric pocket (PDB 5I73) with scores of -8.15 and -8.40 kcal/mol and MM-GBSA of -96.14 and -68.46 kcal/mol whereas (S)-citalopram has -6.90 and -69.39 kcal/mol respectively. The ligands also conferred conformational stability on the receptors during 100 ns MD simulations and displayed interesting ADMET profiles, representing promising hSERT modulators for MDD upon experimental validation.Communicated by Ramaswamy H. Sarma.
    Matched MeSH terms: Molecular Dynamics Simulation
  16. Comparative modeling and enzymatic affinity of novel haloacid dehalogenase from Bacillus megaterium strain BHS1 isolated from alkaline Blue Lake in Turkey
    Wahhab BH, Oyewusi HA, Wahab RA, Mohammad Hood MH, Abdul Hamid AA, Al-Nimer MS, et al.
    J Biomol Struct Dyn, 2024;42(3):1429-1442.
    PMID: 37038649 DOI: 10.1080/07391102.2023.2199870
    This study presents the initial structural model of L-haloacid dehalogenase (DehLBHS1) from Bacillus megaterium BHS1, an alkalotolerant bacterium known for its ability to degrade halogenated environmental pollutants. The model provides insights into the structural features of DehLBHS1 and expands our understanding of the enzymatic mechanisms involved in the degradation of these hazardous pollutants. Key amino acid residues (Arg40, Phe59, Asn118, Asn176, and Trp178) in DehLBHS1 were identified to play critical roles in catalysis and molecular recognition of haloalkanoic acid, essential for efficient binding and transformation of haloalkanoic acid molecules. DehLBHS1 was modeled using I-TASSER, yielding a best TM-score of 0.986 and an RMSD of 0.53 Å. Validation of the model using PROCHECK revealed that 89.2% of the residues were located in the most favored region, providing confidence in its structural accuracy. Molecular docking simulations showed that the non-simulated DehLBHS1 preferred 2,2DCP over other substrates, forming one hydrogen bond with Arg40 and exhibiting a minimum energy of -2.5 kJ/mol. The simulated DehLBHS1 exhibited a minimum energy of -4.3 kJ/mol and formed four hydrogen bonds with Arg40, Asn176, Asp9, and Tyr11, further confirming the preference for 2,2DCP. Molecular dynamics simulations supported this preference, based on various metrics, including RMSD, RMSF, gyration, hydrogen bonding, and molecular distance. MM-PBSA calculations showed that the DehLBHS1-2,2-DCP complex had a markedly lower binding energy (-21.363 ± 1.26 kcal/mol) than the DehLBHS1-3CP complex (-14.327 ± 1.738 kcal/mol). This finding has important implications for the substrate specificity and catalytic function of DehLBHS1, particularly in the bioremediation of 2,2-DCP in contaminated alkaline environments. These results provide a detailed view of the molecular interactions between the enzyme and its substrate and may aid in the development of more efficient biocatalytic strategies for the degradation of halogenated compounds.Communicated by Ramaswamy H. Sarma.
    Matched MeSH terms: Molecular Dynamics Simulation
  17. Fulltext Phenolic compounds of Theobroma cacao L. show potential against dengue RdRp protease enzyme inhibition by In-silico docking, DFT study, MD simulation and MMGBSA calculation
    Huq AKMM, Roney M, Dubey A, Nasir MH, Tufail A, Aluwi MFFM, et al.
    PLoS One, 2024;19(3):e0299238.
    PMID: 38483871 DOI: 10.1371/journal.pone.0299238
    BACKGROUND: Currently, there is no antiviral medication for dengue, a potentially fatal tropical infectious illness spread by two mosquito species, Aedes aegypti and Aedes albopictus. The RdRp protease of dengue virus is a potential therapeutic target. This study focused on the in silico drug discovery of RdRp protease inhibitors.

    METHODS: To assess the potential inhibitory activity of 29 phenolic acids from Theobroma cacao L. against DENV3-NS5 RdRp, a range of computational methods were employed. These included docking, drug-likeness analysis, ADMET prediction, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. The aim of these studies was to confirm the stability of the ligand-protein complex and the binding pose identified during the docking experiment.

    RESULTS: Twenty-one compounds were found to have possible inhibitory activities against DENV according to the docking data, and they had a binding affinity of ≥-37.417 kcal/mol for DENV3- enzyme as compared to the reference compound panduratin A. Additionally, the drug-likeness investigation produced four hit compounds that were subjected to ADMET screening to obtain the lead compound, catechin. Based on ELUMO, EHOMO, and band energy gap, the DFT calculations showed strong electronegetivity, favouravle global softness and chemical reactivity with considerable intra-molecular charge transfer between electron-donor to electron-acceptor groups for catechin. The MD simulation result also demonstrated favourable RMSD, RMSF, SASA and H-bonds in at the binding pocket of DENV3-NS5 RdRp for catechin as compared to panduratin A.

    CONCLUSION: According to the present findings, catechin showed high binding affinity and sufficient drug-like properties with the appropriate ADMET profiles. Moreover, DFT and MD studies further supported the drug-like action of catechin as a potential therapeutic candidate. Therefore, further in vitro and in vivo research on cocoa and its phytochemical catechin should be taken into consideration to develop as a potential DENV inhibitor.

    Matched MeSH terms: Molecular Dynamics Simulation
  18. Fulltext Unravelling the link between SARS-CoV-2 mutation frequencies, patient comorbidities, and structural dynamics
    Azzeri A, Mohamed NA, Wan Rosli SH, Abdul Samat MN, Rashid ZZ, Mohamad Jamali MA, et al.
    PLoS One, 2024;19(3):e0291892.
    PMID: 38483913 DOI: 10.1371/journal.pone.0291892
    Genomic surveillance is crucial for tracking emergence and spread of novel variants of pathogens, such as SARS-CoV-2, to inform public health interventions and to enforce control measures. However, in some settings especially in low- and middle- income counties, where sequencing platforms are limited, only certain patients get to be selected for sequencing surveillance. Here, we show that patients with multiple comorbidities potentially harbour SARS-CoV-2 with higher mutation rates and thus deserve more attention for genomic surveillance. The relationship between the patient comorbidities, and type of amino acid mutations was assessed. Correlation analysis showed that there was a significant tendency for mutations to occur within the ORF1a region for patients with higher number of comorbidities. Frequency analysis of the amino acid substitution within ORF1a showed that nsp3 P822L of the PLpro protease was one of the highest occurring mutations. Using molecular dynamics, we simulated that the P822L mutation in PLpro represents a system with lower Root Mean Square Deviation (RMSD) fluctuations, and consistent Radius of gyration (Rg), Solvent Accessible Surface Area (SASA) values-indicate a much stabler protein than the wildtype. The outcome of this study will help determine the relationship between the clinical status of a patient and the mutations of the infecting SARS-CoV-2 virus.
    Matched MeSH terms: Molecular Dynamics Simulation
  19. Immunoinformatics-based potential multi-peptide vaccine designing against Jamestown Canyon Virus (JCV) capable of eliciting cellular and humoral immune responses
    Shahab M, Aiman S, Alshammari A, Alasmari AF, Alharbi M, Khan A, et al.
    Int J Biol Macromol, 2023 Dec 31;253(Pt 2):126678.
    PMID: 37666399 DOI: 10.1016/j.ijbiomac.2023.126678
    Jamestown Canyon virus (JCV) is a deadly viral infection transmitted by various mosquito species. This mosquito-borne virus belongs to Bunyaviridae family, posing a high public health threat in the in tropical regions of the United States causing encephalitis in humans. Common symptoms of JCV include fever, headache, stiff neck, photophobia, nausea, vomiting, and seizures. Despite the availability of resources, there is currently no vaccine or drug available to combat JCV. The purpose of this study was to develop an epitope-based vaccine using immunoinformatics approaches. The vaccine aimed to be secure, efficient, bio-compatible, and capable of stimulating both innate and adaptive immune responses. In this study, the protein sequence of JCV was obtained from the NCBI database. Various bioinformatics methods, including toxicity evaluation, antigenicity testing, conservancy analysis, and allergenicity assessment were utilized to identify the most promising epitopes. Suitable linkers and adjuvant sequences were used in the design of vaccine construct. 50s ribosomal protein sequence was used as an adjuvant at the N-terminus of the construct. A total of 5 CTL, 5 HTL, and 5 linear B cell epitopes were selected based on non-allergenicity, immunological potential, and antigenicity scores to design a highly immunogenic multi-peptide vaccine construct. Strong interactions between the proposed vaccine and human immune receptors, i.e., TLR-2 and TLR-4, were revealed in a docking study using ClusPro software, suggesting their possible relevance in the immunological response to the vaccine. Immunological and physicochemical properties assessment ensured that the proposed vaccine demonstrated high immunogenicity, solubility and thermostability. Molecular dynamics simulations confirmed the strong binding affinities, as well as dynamic and structural stability of the proposed vaccine. Immune simulation suggest that the vaccine has the potential to effectively stimulate cellular and humoral immune responses to combat JCV infection. Experimental and clinical assays are required to validate the results of this study.
    Matched MeSH terms: Molecular Dynamics Simulation
  20. Fulltext Pharmacotherapeutic Potential of Natural Products to Target the SARS-CoV-2 PLpro Using Molecular Screening and Simulation Approaches
    Sayaf AM, Ahmad H, Aslam MA, Ghani SA, Bano S, Yousafi Q, et al.
    Appl Biochem Biotechnol, 2023 Nov;195(11):6959-6978.
    PMID: 36961512 DOI: 10.1007/s12010-023-04466-1
    Because of the essential role of PLpro in the regulation of replication and dysregulation of the host immune sensing, it is considered a therapeutic target for novel drug development. To reduce the risk of immune evasion and vaccine effectiveness, small molecular therapeutics are the best complementary approach. Hence, we used a structure-based drug-designing approach to identify potential small molecular inhibitors for PLpro of SARS-CoV-2. Initial scoring and re-scoring of the best hits revealed that three compounds NPC320891 (2,2-Dihydroxyindene-1,3-Dione), NPC474594 (Isonarciclasine), and NPC474595 (7-Deoxyisonarciclasine) exhibit higher docking scores than the control GRL0617. Investigation of the binding modes revealed that alongside the essential contacts, i.e., Asp164, Glu167, Tyr264, and Gln269, these molecules also target Lys157 and Tyr268 residues in the active site. Moreover, molecular simulation demonstrated that the reported top hits also possess stable dynamics and structural packing. Furthermore, the residues' flexibility revealed that all the complexes demonstrated higher flexibility in the regions 120-140, 160-180, and 205-215. The 120-140 and 160-180 lie in the finger region of PLpro, which may open/close during the simulation to cover the active site and push the ligand inside. In addition, the total binding free energy was reported to be - 32.65 ± 0.17 kcal/mol for the GRL0617-PLpro, for the NPC320891-PLpro complex, the TBE was - 35.58 ± 0.14 kcal/mol, for the NPC474594-PLpro, the TBE was - 43.72 ± 0.22 kcal/mol, while for NPC474595-PLpro complex, the TBE was calculated to be - 41.61 ± 0.20 kcal/mol, respectively. Clustering of the protein's motion and FEL further revealed that in NPC474594 and NPC474595 complexes, the drug was seen to have moved inside the binding cavity along with the loop in the palm region harboring the catalytic triad, thus justifying the higher binding of these two molecules particularly. In conclusion, the overall results reflect favorable binding of the identified hits strongly than the control drug, thus demanding in vitro and in vivo validation for clinical purposes.
    Matched MeSH terms: Molecular Dynamics Simulation
Related Terms
Filters
Contact Us

Please provide feedback to Administrator (afdal@afpm.org.my)

External Links