Displaying publications 1 - 20 of 224 in total

Abstract:
Sort:
  1. Zifruddin AN, Mohamad Yusoff MA, Abd Ghani NS, Nor Muhammad NA, Lam KW, Hassan M
    Comput Biol Chem, 2023 Apr;103:107811.
    PMID: 36645937 DOI: 10.1016/j.compbiolchem.2023.107811
    Metisa plana (Lepidoptera: Psychidae) bagworm is a leaf-eater caterpillar ubiquitously found as a damaging pest in oil palm plantations, specifically in Malaysia. Various strategies have been implemented, including the usage of chemical insecticides. However, the main challenges include the development of insecticide resistance and its detrimental effects on the environment and non-target organisms. Therefore, a biorational insecticide is introduced by targeting the juvenile hormone (JH) biosynthetic pathway, which is mainly present in the insect and vital for the insect's growth, diapause, metamorphosis, and adult reproduction. This study aimed to investigate the potential inhibitor for the rate-limiting enzyme involved in the JH pathway known as farnesol dehydrogenase. A 255 amino acids sequence encoded for the putative M. plana farnesol dehydrogenase (MpFolDH) open reading frame had been identified and isolated. The three-dimensional structure of MpFolDH was predicted to have seven β- sheets with α-helices at both sides, showing typical characteristics for classical short-chain dehydrogenase and associated with oxidoreductase activity. Then, the ensemble-based virtual screening was conducted based on the ZINC20 database, in which 43 768 compounds that fulfilled pesticide-likeness criteria were screened by site-specific molecular docking. After a short molecular dynamics simulation (5 ns) was conducted towards 102 compounds, only the top 10 compounds based on their most favourable binding energy were selected for a more extended simulation (100 ns). Based on the protein-ligand stability, protein compactness, residues rigidity, binding interaction, binding energy throughout the 100 ns simulation, and physicochemical analysis, ZINC000408743205 was selected as a potential inhibitor for this enzyme. Amino acids decomposition analysis indicates Ile18, Ala95, Val198 and Val202 were the critical contributor residues for MpFolDH-inhibitors(s) complex.
    Matched MeSH terms: Molecular Dynamics Simulation
  2. 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
  3. Zainal Abidin MH, Abd Halim KB, Huyop F, Tengku Abdul Hamid TH, Abdul Wahab R, Abdul Hamid AA
    J Mol Graph Model, 2019 07;90:219-225.
    PMID: 31103914 DOI: 10.1016/j.jmgm.2019.05.003
    Dehalogenase E (DehE) is a non-stereospecific enzyme produced by the soil bacterium, Rhizobium sp. RC1. Till now, the catalytic mechanism of DehE remains unclear although several literature concerning its structure and function are available. Since DehE is non-stereospecific, the enzyme was hypothesized to follow a 'direct attack mechanism' for the catalytic breakdown of a haloacid. For a molecular insight, the DehE modelled structure was docked in silico with the substrate 2-chloropropionic acid (2CP) in the active site. The ideal position of DehE residues that allowed a direct attack mechanism was then assessed via molecular dynamics (MD) simulation. It was revealed that the essential catalytic water was hydrogen bonded to the 'water-bearer', Asn114, at a relatively constant distance of ∼2.0 Å after 50 ns. The same water molecule was also closely sited to the catalytic Asp189 at an average distance of ∼2.0 Å, signifying the imperative role of the latter to initiate proton abstraction for water activation. This reaction was crucial to promote a direct attack on the α-carbon of 2CP to eject the halide ion. The water molecule was oriented favourably towards the α-carbon of 2CP at an angle of ∼75°, mirrored by the formation of stable enzyme-substrate orientations throughout the simulation. The data therefore substantiated that the degradation of a haloacid by DehE followed a 'direct attack mechanism'. Hence, this study offers valuable information into future advancements in the engineering of haloacid dehalogenases with improved activity and selectivity, as well as functionality in solvents other than water.
    Matched MeSH terms: Molecular Dynamics Simulation
  4. Yusuf M, Mohamed N, Mohamad S, Janezic D, Damodaran KV, Wahab HA
    J Chem Inf Model, 2016 Jan 25;56(1):82-100.
    PMID: 26703840 DOI: 10.1021/acs.jcim.5b00331
    Increased reports of oseltamivir (OTV)-resistant strains of the influenza virus, such as the H274Y mutation on its neuraminidase (NA), have created some cause for concern. Many studies have been conducted in the attempt to uncover the mechanism of OTV resistance in H274Y NA. However, most of the reported studies on H274Y focused only on the drug-bound system, so the direct effects of the mutation on NA itself prior to drug binding still remain unclear. Therefore, molecular dynamics simulations of NA in apo form, followed by principal component analysis and interaction energy calculations, were performed to investigate the structural changes of the NA binding site as a result of the H274Y mutation. It was observed that the disruption of the NA binding site due to the H274Y mutation was initiated by the repulsive effect of Y274 on the 250-loop, which in turn altered the hydrogen-bonding network around residue 274. The rotated W295 side chain caused the upward movement of the 340-loop. Consequently, sliding box docking results suggested that the binding pathway of OTV was compromised because of the disruption of this binding site. This study also highlighted the importance of the functional group at C6 of the sialic acid mimicry. It is hoped that these results will improve the understanding of OTV resistance and shed some light on the design of a novel anti-influenza drug.
    Matched MeSH terms: Molecular Dynamics Simulation*
  5. Yotmanee P, Rungrotmongkol T, Wichapong K, Choi SB, Wahab HA, Kungwan N, et al.
    J Mol Graph Model, 2015 Jul;60:24-33.
    PMID: 26086900 DOI: 10.1016/j.jmgm.2015.05.008
    The pathogenic dengue virus (DV) is a growing global threat, particularly in South East Asia, for which there is no specific treatment available. The virus possesses a two-component (NS2B/NS3) serine protease that cleaves the viral precursor proteins. Here, we performed molecular dynamics simulations of the NS2B/NS3 protease complexes with six peptide substrates (capsid, intNS3, 2A/2B, 4B/5, 3/4A and 2B/3 containing the proteolytic site between P(1) and P(1)' subsites) of DV type 2 to compare the specificity of the protein-substrate binding recognition. Although all substrates were in the active conformation for cleavage reaction by NS2B/NS3 protease, their binding strength was somewhat different. The simulated results of intermolecular hydrogen bonds and decomposition energies suggested that among the ten substrate residues (P(5)-P(5)') the P(1) and P(2) subsites play a major role in the binding with the focused protease. The arginine residue at these two subsites was found to be specific preferential binding at the active site with a stabilization energy of intNS3>2A/2B>4B/5>3/4A>2B/3 in a relative correspondence with previous experimentally derived values.
    Matched MeSH terms: Molecular Dynamics Simulation*
  6. Yeni Y, Supandi S, Dwita LP, Suswandari S, Shaharun MS, Sambudi NS
    J Pharm Bioallied Sci, 2020 Nov;12(Suppl 2):S836-S840.
    PMID: 33828386 DOI: 10.4103/jpbs.JPBS_103_20
    Background: Inflammatory mediators produced by cyclooxygenase (COX) and lipoxygenase (LOX) pathways are responsible for many human diseases, such as cancer, arthritis, and neurological disorders. Flavonoid-containing plants, such as Ipomoea batatas leaves, have shown potential anti-inflammatory activity.

    Objectives: This study aimed to predict the actions of 10 compounds in I. batatas leaves, which are YGM-0a [cyanidin 3-0-sophoroside-5-0-glucosede], YGM-0f [cyanidin 3-O-(2-0-(6-0-(E)-p-coumaroyl-β-D-glucopyranosyl)-β-D-glucopyranoside)-5-0-β-D-glucopyranoside], YGM-1a [cyanidin 3-(6,6'-caffeylp-hydroxybenzoylsophoroside) -5-glucoside], YGM-1b [cyanidin 3-(6,6'-dicaffeylsophor-oside)-5-glucoside], YGM-2 [cyanidin 3-(6-caffeylsophoroside)-5-glucoside], YGM-3 [cyanidin 3-(6,6'-caffeyl-ferulylsophoroside)-5-glucoside], YGM-4b [peonidin 3-(6,6'-dicaffeylsophoroside)-5- glucoside], YGM-5a [peonidin 3-(6,6'-caffeylphydroxybenzo-ylsophoroside)-5-gluco-side], YGM-5b [cyanidin 3-6-caffeylsophoroside)-5-glucosede], and YGM-6 [peonidin 3-(6,6'-caffeylferulylsophoroside)-5-glucoside] as LOX inhibitors, and also predict the stability of ligand-LOX complex.

    Materials and Methods: The compounds were screened through docking studies using PLANTS. Also, the molecular dynamics simulation was conducted using GROMACS at 310K.

    Results: The results showed that the most significant binding affinity toward LOX was shown by YGM-0a and YGM-0a, and the LOX complex in molecular dynamics simulation showed stability for 20 ns.

    Conclusion: Based on Docking Studies and Molecular Dynamics Simulation of I. Batatas Leaves compounds, YGM-0a was shown to be the most probable LOX inhibitor.

    Matched MeSH terms: Molecular Dynamics Simulation
  7. Yau MQ, Emtage AL, Chan NJY, Doughty SW, Loo JSE
    J Comput Aided Mol Des, 2019 05;33(5):487-496.
    PMID: 30989574 DOI: 10.1007/s10822-019-00201-3
    The recent expansion of GPCR crystal structures provides the opportunity to assess the performance of structure-based drug design methods for the GPCR superfamily. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA)-based methods are commonly used for binding affinity prediction, as they provide an intermediate compromise of speed and accuracy between the empirical scoring functions used in docking and more robust free energy perturbation methods. In this study, we systematically assessed the performance of MM/PBSA in predicting experimental binding free energies using twenty Class A GPCR crystal structures and 934 known ligands. Correlations between predicted and experimental binding free energies varied significantly between individual targets, ranging from r = - 0.334 in the inactive-state CB1 cannabinoid receptor to r = 0.781 in the active-state CB1 cannabinoid receptor, while average correlation across all twenty targets was relatively poor (r = 0.183). MM/PBSA provided better predictions of binding free energies compared to docking scores in eight out of the twenty GPCR targets while performing worse for four targets. MM/PBSA binding affinity predictions calculated using a single, energy minimized structure provided comparable predictions to sampling from molecular dynamics simulations and may be more efficient when computational cost becomes restrictive. Additionally, we observed that restricting MM/PBSA calculations to ligands with a high degree of structural similarity to the crystal structure ligands improved performance in several cases. In conclusion, while MM/PBSA remains a valuable tool for GPCR structure-based drug design, its performance in predicting the binding free energies of GPCR ligands remains highly system-specific as demonstrated in a subset of twenty Class A GPCRs, and validation of MM/PBSA-based methods for each individual case is recommended before prospective use.
    Matched MeSH terms: Molecular Dynamics Simulation
  8. Yau MQ, Emtage AL, Loo JSE
    J Comput Aided Mol Des, 2020 Nov;34(11):1133-1145.
    PMID: 32851579 DOI: 10.1007/s10822-020-00339-5
    Recent breakthroughs in G protein-coupled receptor (GPCR) crystallography and the subsequent increase in number of solved GPCR structures has allowed for the unprecedented opportunity to utilize their experimental structures for structure-based drug discovery applications. As virtual screening represents one of the primary computational methods used for the discovery of novel leads, the GPCR-Bench dataset was created to facilitate comparison among various virtual screening protocols. In this study, we have benchmarked the performance of Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) in improving virtual screening enrichment in comparison to docking with Glide, using the entire GPCR-Bench dataset of 24 GPCR targets and 254,646 actives and decoys. Reranking the top 10% of the docked dataset using MM/PBSA resulted in improvements for six targets at EF1% and nine targets at EF5%, with the gains in enrichment being more pronounced at the EF1% level. We additionally assessed the utility of rescoring the top ten poses from docking and the ability of short MD simulations to refine the binding poses prior to MM/PBSA calculations. There was no clear trend of the benefit observed in both cases, suggesting that utilizing a single energy minimized structure for MM/PBSA calculations may be the most computationally efficient approach in virtual screening. Overall, the performance of MM/PBSA rescoring in improving virtual screening enrichment obtained from docking of the GPCR-Bench dataset was found to be relatively modest and target-specific, highlighting the need for validation of MM/PBSA-based protocols prior to prospective use.
    Matched MeSH terms: Molecular Dynamics Simulation
  9. Yap BK, Buckle MJ, Doughty SW
    J Mol Model, 2012 Aug;18(8):3639-55.
    PMID: 22354276 DOI: 10.1007/s00894-012-1368-5
    5-HT(1A) serotonin and D1 dopamine receptor agonists have been postulated to be able to improve negative and cognitive impairment symptoms of schizophrenia, while partial agonists and antagonists of the D2 and 5-HT(2A) receptors have been reported to be effective in reducing positive symptoms. There is therefore a need for well-defined homology models for the design of more selective antipsychotic agents, since no three-dimensional (3D) crystal structures of these receptors are currently available. In this study, homology models were built based on the high-resolution crystal structure of the β(2)-adrenergic receptor (2RH1) and further refined via molecular dynamics simulations in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer system with the GROMOS96 53A6 united atom force field. Docking evaluations with representative agonists and antagonists using AutoDock 4.2 revealed binding modes in agreement with experimentally determined site-directed mutagenesis data and significant correlations between the computed and experimental pK (i) values. The models are also able to distinguish between antipsychotic agents with different selectivities and binding affinities for the four receptors, as well as to differentiate active compounds from decoys. Hence, these human 5-HT(1A), 5-HT(2A), D1 and D2 receptor homology models are capable of predicting the activities of novel ligands, and can be used as 3D templates for antipsychotic drug design and discovery.
    Matched MeSH terms: Molecular Dynamics Simulation*
  10. Yana J, Chiangraeng N, Nimmanpipug P, Lee VS
    J Mol Graph Model, 2021 09;107:107946.
    PMID: 34119952 DOI: 10.1016/j.jmgm.2021.107946
    Conformational search for the most stable geometry connection of 16 sets of polydopamine (PDA) tetramer subunits has been systematically investigated using density functional theory (DFT) calculations. Our results indicated that the more planar subunits are, the more stable they are. This finding is in good agreement with recent experimental observations, which have suggested that PDA are composed of the nearly planar subunits that appear to be stacked together via the π-π interactions to form graphite-like layered aggregates associated with the balance of the intramolecular hydrogen bonds and steric effects from the indole and catechol moieties. Molecular dynamics (MD) simulations of 16 spherical clusters of the tetramer subunits of PDA in the gas and aqueous phase were performed at 298 K and confirmed the stability of supramolecular tetramer aggregates. The complex formation and binding energy of all 16 clusters are very strong although the shapes of the clusters in aqueous solution are not spherical and are very much different from those in the gas phase. The aggregations of all 16 clusters in aqueous solution were also confirmed from the profiles of the Kratky plot and the radius of gyration of all clusters. Our MD results in both gas phase and aqueous solution pointed out that there are high possibilities of aggregations of the 16 kinds of tetramer subunits although the conformations of each tetramer subunit are not flat. In summary, this work brings an insight into the controversial structure of PDA tetramer units and explains some of the important structural features found in the aqueous phase in comparison to the gas phase.
    Matched MeSH terms: Molecular Dynamics Simulation*
  11. Yahaya MAF, Bakar ARA, Stanslas J, Nordin N, Zainol M, Mehat MZ
    BMC Biotechnol, 2021 06 05;21(1):38.
    PMID: 34090414 DOI: 10.1186/s12896-021-00697-4
    BACKGROUND: Neuroinflammation has been identified to be the key player in most neurodegenerative diseases. If neuroinflammation is left to be unresolved, chronic neuroinflammation will be establish. Such situation is due to the overly-activated microglia which have the tendency to secrete an abundance amount of pro-inflammatory cytokines into the neuron microenvironment. The abundance of pro-inflammatory cytokines will later cause toxic and death to neurons. Toll-like receptor 4 (TLR4)/MD-2 complex found on the cell surface of microglia is responsible for the attachment of LPS and activation of nuclear factor-κB (NF-κB) downstream signalling pathway. Albeit vitexin has been shown to possess anti-inflammatory property, however, little is known on its ability to bind at the binding site of TLR4/MD-2 complex of microglia as well as to be an antagonist for LPS.

    RESULTS: The present study reveals that both vitexin and donepezil are able to bind at the close proximity of LPS binding site located at the TLR4/MD-2 complex with the binding energy of - 4.35 and - 9.14 kcal/mol, respectively. During molecular dynamic simulations, both vitexin and donepezil formed stable complex with TLR4/MD-2 throughout the 100 ns time length with the root mean square deviation (RMSD) values of 2.5 Å and 4.0 Å, respectively. The root mean square fluctuation (RMSF) reveals that both compounds are stable. Interestingly, the radius of gyration (rGyr) for donepezil shows notable fluctuations when compare with vitexin. The MM-GBSA results showed that vitexin has higher binding energy in comparison with donepezil.

    CONCLUSIONS: Taken together, the findings suggest that vitexin is able to bind at the binding site of TLR4/MD-2 complex with more stability than donepezil throughout the course of 100 ns simulation. Hence, vitexin has the potential to be an antagonist candidate for LPS.

    Matched MeSH terms: Molecular Dynamics Simulation
  12. Yaacob N, Mohamad Ali MS, Salleh AB, Rahman RNZRA, Leow ATC
    J Mol Graph Model, 2016 07;68:224-235.
    PMID: 27474867 DOI: 10.1016/j.jmgm.2016.07.003
    The utilization of cold active lipases in organic solvents proves an excellent approach for chiral synthesis and modification of fats and oil due to the inherent flexibility of lipases under low water conditions. In order to verify whether this lipase can function as a valuable synthetic catalyst, the mechanism concerning activation of the lid and interacting solvent residues in the presence of organic solvent must be well understood. A new alkaline cold-adapted lipase, AMS8, from Pseudomonas fluorescens was studied for its structural adaptation and flexibility prior to its exposure to non-polar, polar aprotic and protic solvents. Solvents such as ethanol, toluene, DMSO and 2-propanol showed to have good interactions with active sites. Asparagine (Asn) and tyrosine (Tyr) were key residues attracted to solvents because they could form hydrogen bonds. Unlike in other solvents, Phe-18, Tyr-236 and Tyr-318 were predicted to have aromatic-aromatic side-chain interactions with toluene. Non-polar solvent also was found to possess highest energy binding compared to polar solvents. Due to this circumstance, the interaction of toluene and AMS8 lipase was primarily based on hydrophobicity and molecular recognition. The molecular dynamic simulation showed that lid 2 (residues 148-167) was very flexible in toluene and Ca(2+). As a result, lid 2 moves away from the catalytic areas, leaving an opening for better substrate accessibility which promotes protein activation. Only a single lid (lid 2) showed the movement following interactions with toluene, although AMS8 lipase displayed double lids. The secondary conformation of AMS8 lipase that was affected by toluene observed a reduction of helical strands and increased coil structure. Overall, this work shows that cold active lipase, AMS8 exhibits distinguish interfacial activation and stability in the presence of polar and non-polar solvents.
    Matched MeSH terms: Molecular Dynamics Simulation
  13. Yaacob N, Ahmad Kamarudin NH, Leow ATC, Salleh AB, Raja Abd Rahman RNZ, Mohamad Ali MS
    Molecules, 2017 Aug 12;22(8).
    PMID: 28805665 DOI: 10.3390/molecules22081312
    The alkaline cold-active lipase from Pseudomonas fluorescens AMS8 undergoes major structural changes when reacted with hydrophobic organic solvents. In toluene, the AMS8 lipase catalytic region is exposed by the moving hydrophobic lid 2 (Glu-148 to Gly-167). Solvent-accessible surface area analysis revealed that Leu-208, which is located next to the nucleophilic Ser-207 has a focal function in influencing substrate accessibility and flexibility of the catalytic pocket. Based on molecular dynamic simulations, it was found that Leu-208 strongly facilitates the lid 2 opening via its side-chain. The KM and Kcat/KM of L208A mutant were substrate dependent as it preferred a smaller-chain ester (pNP-caprylate) as compared to medium (pNP-laurate) or long-chain (pNP-palmitate) esters. In esterification of ethyl hexanoate, L208A promotes a higher ester conversion rate at 20 °C but not at 30 °C, as a 27% decline was observed. Interestingly, the wild-type (WT) lipase's conversion rate was found to increase with a higher temperature. WT lipase AMS8 esterification was higher in toluene as compared to L208A. Hence, the results showed that Leu-208 of AMS8 lipase plays an important role in steering a broad range of substrates into its active site region by regulating the flexibility of this region. Leu-208 is therefore predicted to be crucial for its role in interfacial activation and catalysis in toluene.
    Matched MeSH terms: Molecular Dynamics Simulation*
  14. Yaacob MA, Hasan WA, Ali MS, Rahman RN, Salleh AB, Basri M, et al.
    Acta Biochim. Pol., 2014;61(4):745-52.
    PMID: 25337608
    Genome mining revealed a 1011 nucleotide-long fragment encoding a type I L-asparaginase (J15 asparaginase) from the halo-tolerant Photobacterium sp. strain J15. The gene was overexpressed in pET-32b (+) vector in E. coli strain Rosetta-gami B (DE3) pLysS and purified using two-step chromatographic methods: Ni(2+)-Sepharose affinity chromatography and Q-Sepharose anion exchange chromatography. The final specific activity and yield of the enzyme achieved from these steps were 20 U/mg and 49.2%, respectively. The functional dimeric form of J15-asparaginase was characterised with a molecular weight of ~70 kDa. The optimum temperature and pH were 25°C and pH 7.0, respectively. This protein was stable in the presence of 1 mM Ni(2+) and Mg(2+), but it was inhibited by Mn(2+), Fe(3+) and Zn(2+) at the same concentration. J15 asparaginase actively hydrolysed its native substrate, l-asparagine, but had low activity towards l-glutamine. The melting temperature of J15 asparaginase was ~51°C, which was determined using denatured protein analysis of CD spectra. The Km, Kcat, Kcat/Km of J15 asparaginase were 0.76 mM, 3.2 s(-1), and 4.21 s(-1) mM(-1), respectively. Conformational changes of the J15 asparaginase 3D structure at different temperatures (25°C, 45°C, and 65°C) were analysed using Molecular Dynamic simulations. From the analysis, residues Tyr₂₄ , His₂₂, Gly₂₃, Val₂₅ and Pro₂₆ may be directly involved in the 'open' and 'closed' lid-loop conformation, facilitating the conversion of substrates during enzymatic reactions. The properties of J15 asparaginase, which can work at physiological pH and has low glutaminase activity, suggest that this could be a good candidate for reducing toxic effects during cancer treatment.
    Matched MeSH terms: Molecular Dynamics Simulation
  15. Xing M, Akowuah GA, Gautam V, Gaurav A
    J Biomol Struct Dyn, 2017 Oct;35(13):2910-2924.
    PMID: 27608741 DOI: 10.1080/07391102.2016.1234417
    Phosphodiesterase 4 (PDE4) has been established as a drug target for inflammatory diseases of respiratory tract like asthma and chronic obstructive pulmonary disease. The selective inhibitors of PDE4B, a subtype of PDE4, are devoid of adverse effects like nausea and vomiting commonly associated with non-selective PDE4B inhibitors. This makes the development of PDE4B subtype selective inhibitors a desirable research goal. Thus, in the present study, molecular docking, molecular dynamic simulations and binding free energy were performed to explore potential selective PDE4B inhibitors based on ginger phenolic compounds. The results of docking studies indicate that some of the ginger phenolic compounds demonstrate higher selective PDE4B inhibition than existing selective PDE4B inhibitors. Additionally, 6-gingerol showed the highest PDE4B inhibitory activity as well as selectivity. The comparison of binding mode of PDE4B/6-gingerol and PDE4D/6-gingerol complexes revealed that 6-gingerol formed additional hydrogen bond and hydrophobic interactions with active site and control region 3 (CR3) residues in PDE4B, which were primarily responsible for its PDE4B selectivity. The results of binding free energy demonstrated that electrostatic energy is the primary factor in elucidating the mechanism of PDE4B inhibition by 6-gingerol. Dynamic cross-correlation studies also supported the results of docking and molecular dynamics simulation. Finally, a small library of molecules were designed based on the identified structural features, majority of designed molecules showed higher PDE4B selectivity than 6-gingerol. These results provide important structural features for designing new selective PDE4B inhibitors as anti-inflammatory drugs and promising candidates for synthesis and pre-clinical pharmacological investigations.
    Matched MeSH terms: Molecular Dynamics Simulation
  16. Woon KL, Mustapa SAS, Mohd Jamel NS, Lee VS, Zakaria MZ, Ariffin A
    Chemphyschem, 2020 Sep 17.
    PMID: 32940952 DOI: 10.1002/cphc.202000612
    Material designs that use donor and acceptor units are often found in organic optoelectronic devices. Molecular level insight into the interactions between donors and acceptors are crucial for understanding how such interactions can modify the optical properties of the organic optoelectronic materials. In this paper, tris(4-(tert-butyl)phenyl)amine (pTPA) was synthesized as a donor in order to compare with unmodified triphenylamine (TPA) in a donor-acceptor system by having 2,4,6-triphenyl-1,3,5-triazine (TRZ) as an acceptor. Dimerization of donors and acceptors occurred in solvent when the concentration of solute is high. At 0 K, using a polarizable continuum model, the nitrogen atom of TPA is found to stack on top of the center of triazine of TRZ, whereas such alignment is offset in pTPA and TRZ. We attributed such alignment in TPA-TRZ as the result of attractive interactions between partial localization of 2pz electrons at the nitrogen atom of TPA and the π deficiency of triazine in TPA-TRZ. By taking into account random motions of the solvent effect at 300 K in quantum molecular dynamics and classical molecular dynamics simulations to interpret the marked difference in emission spectra between TPA-TRZ and pTPA-TRZ, it was revealed that the attractive interaction between pTPA and TRZ in toluene is weaker than TPA and TRZ. Because of the weaker attractive interaction between pTPA and TRZ in toluene, the dimers adopted numerous ground state conformations resulting in broad emission bands superimposed with multiple small Gaussian peaks. This is in contrast to TPA-TRZ which has only one dominant dimer conformation. This study demonstrates that the strength of intermolecular interactions between donors and acceptors should be taken into consideration in designing supramolecular structures.
    Matched MeSH terms: Molecular Dynamics Simulation
  17. Wongrattanakamon P, Lee VS, Nimmanpipug P, Sirithunyalug B, Chansakaow S, Jiranusornkul S
    Toxicol. Mech. Methods, 2017 May;27(4):253-271.
    PMID: 27996361 DOI: 10.1080/15376516.2016.1273428
    In this work, molecular docking, pharmacophore modeling and molecular dynamics (MD) simulation were rendered for the mouse P-glycoprotein (P-gp) (code: 4Q9H) and bioflavonoids; amorphigenin, chrysin, epigallocatechin, formononetin and rotenone including a positive control; verapamil to identify protein-ligand interaction features including binding affinities, interaction characteristics, hot-spot amino acid residues and complex stabilities. These flavonoids occupied the same binding site with high binding affinities and shared the same key residues for their binding interactions and the binding region of the flavonoids was revealed that overlapped the ATP binding region with hydrophobic and hydrophilic interactions suggesting a competitive inhibition mechanism of the compounds. Root mean square deviations (RMSDs) analysis of MD trajectories of the protein-ligand complexes and NBD2 residues, and ligands pointed out these residues were stable throughout the duration of MD simulations. Thus, the applied preliminary structure-based molecular modeling approach of interactions between NBD2 and flavonoids may be gainful to realize the intimate inhibition mechanism of P-gp at NBD2 level and on the basis of the obtained data, it can be concluded that these bioflavonoids have the potential to cause herb-drug interactions or be used as lead molecules for the inhibition of P-gp (as anti-multidrug resistance agents) via the NBD2 blocking mechanism in future.
    Matched MeSH terms: Molecular Dynamics Simulation
  18. Wang Y, Wei DQ, Wang JF
    J Chem Inf Model, 2010 May 24;50(5):875-8.
    PMID: 20443585 DOI: 10.1021/ci900458u
    T1 lipase is isolated from the palm Geobacillus zalihae strain T1 in Malaysia, functioning as a secreted protein responsible for the catalyzing hydrolysis of long-chain triglycerides into fatty acids and glycerol at high temperatures. In the current study, using 30 ns molecular dynamics simulations at different temperatures, an aqueous activation was detected for T1 lipase. This aqueous activation in T1 lipase was mainly caused by a double-flap movement mechanism. The double flaps were constituted by the hydrophobic helices 6 and 9. Helix 6 employed two major components with the hydrophilic part at the surface and the hydrophobic part inside. In the aqueous solution, the hydrophobic part could provide enough power for helix 6 to move away, driving the protein into an open configuration and exposing the catalytic triad. Our findings could provide structural evidence to support the double-flap movement, revealing the catalytic mechanism for T1 lipase.
    Matched MeSH terms: Molecular Dynamics Simulation*
  19. Wan Jusoh WZA, Abdul Rahman S, Ahmad AL, Mohd Mokhtar N
    Data Brief, 2019 Jun;24:103910.
    PMID: 31193576 DOI: 10.1016/j.dib.2019.103910
    This paper focus to examine the best molecular interaction between Polyamide Thin Film Composite (PA TFC) layers with different properties of the support membrane. The support membrane of Nylon 66 (N66) and Polyvinylidene fluoride (PVDF) was chosen to represent the hydrophilic and hydrophobic model respectively in the Molecular Dynamic (MD) simulation. The Condensed-Phase Optimized Molecular Potential for Atomistic Simulation Studies (COMPASS) force field was used with the total simulation runs were set 1000 picoseconds run production ensembles. The temperature and pressure set for both ensembles were 298 K and 1 atm respectively. The validity of our model densities data was check and calculated where the deviation must be less than 6%. The comparison between hydrophobic and hydrophilic of the support membrane data was examined by the distance and magnitude of intensity of the Radial Distribution Function (RDF's) trends.
    Matched MeSH terms: Molecular Dynamics Simulation
  20. 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
Filters
Contact Us

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

External Links