Displaying publications 1 - 20 of 60 in total

Abstract:
Sort:
  1. Al-Nema M, Gaurav A, Lee VS
    Heliyon, 2020 Sep;6(9):e04856.
    PMID: 32984588 DOI: 10.1016/j.heliyon.2020.e04856
    Inhibition of phosphodiesterase 4 (PDE4) is a promising therapeutic approach for the treatment of inflammatory pulmonary disorders, i.e. asthma and chronic obstructive pulmonary disease. However, the treatment with non-selective PDE4 inhibitors is associated with side effects such as nausea and vomiting. Among the subtypes of PDE4 inhibited by these inhibitors, PDE4B is expressed in immune, inflammatory and airway smooth muscle cells, whereas, PDE4D is expressed in the area postrema and nucleus of the solitary tract. Thus, PDE4D inhibition is responsible for the emetic response. In this regard, a selective PDE4B inhibitor is expected to be a potential drug candidate for the treatment of inflammatory pulmonary disorders. Therefore, a shared feature pharmacophore model was developed and used as a query for the virtual screening of Maybridge and SPECS databases. A number of filters were applied to ensure only compounds with drug-like properties were selected. Accordingly, nine compounds have been identified as final hits, where HTS04529 showed the highest affinity and selectivity for PDE4B over PDE4D in molecular docking. The docked complexes of HTS04529 with PDE4B and PDE4D were subjected to molecular dynamics simulations for 100ns to assess their binding stability. The results showed that HTS04529 was bound tightly to PDE4B and formed a more stable complex with it than with PDE4D.
  2. Mahalapbutr P, Lee VS, Rungrotmongkol T
    J Agric Food Chem, 2020 Jul 29;68(30):7974-7983.
    PMID: 32551626 DOI: 10.1021/acs.jafc.0c02580
    Human sweet taste receptor (hSTR) recognizes a wide array of sweeteners, resulting in sweet taste perception. Maltitol and lactitol have been extensively used in place of sucrose due to their capability to prevent dental caries. Herein, several molecular modeling approaches were applied to investigate the structural and energetic properties of these two polyols/hSTR complexes. Triplicate 500 ns molecular dynamics (MD) simulations and molecular mechanics/generalized Born surface area (MM/GBSA)-based free energy calculations revealed that the TAS1R2 monomer is the preferential binding site for maltitol and lactitol rather than the TAS1R3 region. Several polar residues (D142, S144, Y215, D278, E302, R383, and especially N143) were involved in polyols binding through electrostatic attractions and H-bond formations. The molecular complexation process not only induced the stable form of ligands but also stimulated the conformational adaptation of the TAS1R2 monomer to become a close-packed structure through an induced-fit mechanism. Notably, the binding affinity of the maltitol/TAS1R2 complex (ΔGbind of -17.93 ± 1.49 kcal/mol) was significantly higher than that of the lactitol/TAS1R2 system (-8.53 ± 1.78 kcal/mol), in line with the experimental relative sweetness. These findings provide an in-depth understanding of the differences in the sweetness response between maltitol and lactitol, which could be helpful to design novel polyol derivatives with higher sweet taste perception.
  3. Waris KH, Lee VS, Mohamad S
    Environ Sci Pollut Res Int, 2021 Sep;28(35):47785-47799.
    PMID: 34296410 DOI: 10.1007/s11356-021-15434-9
    The aim of this review is to highlight and provide an update on the current development of pesticide remediation methods, focusing on the utilization of different cyclodextrin (CD) molecules. Because of less environmental impact and non-toxic nature, CDs are beneficial for pesticide remediation, reducing environmental risk and health hazards. They are advantageous for the removal of pesticides from contaminated areas, as well as for better pesticide formulation and, posing significant effects on the hydrolysis or degradation of pesticides. The review focuses on the current trend and innovations regarding the methods and strategies employed for using CDs in designing pesticide remediation. Nowadays, in addition to the conventional experimental techniques, molecular simulation approaches are significantly contributing to the study of such phenomena and hence are recognized as a widely used tool.
  4. Al-Nema M, Gaurav A, Lee VS
    Comput Biol Med, 2023 Apr 03;159:106869.
    PMID: 37071939 DOI: 10.1016/j.compbiomed.2023.106869
    In recent years, the PDE1B enzyme has become a desirable drug target for the treatment of psychological and neurological disorders, particularly schizophrenia disorder, due to the expression of PDE1B in brain regions involved in volitional behaviour, learning and memory. Although several inhibitors of PDE1 have been identified using different methods, none of these inhibitors has reached the market yet. Thus, searching for novel PDE1B inhibitors is considered a major scientific challenge. In this study, pharmacophore-based screening, ensemble docking and molecular dynamics simulations have been performed to identify a lead inhibitor of PDE1B with a new chemical scaffold. Five PDE1B crystal structures have been utilised in the docking study to improve the possibility of identifying an active compound compared to the use of a single crystal structure. Finally, the structure-activity- relationship was studied, and the structure of the lead molecule was modified to design novel inhibitors with a high affinity for PDE1B. As a result, two novel compounds have been designed that exhibited a higher affinity to PDE1B compared to the lead compound and the other designed compounds.
  5. Pirojsirikul T, Lee VS, Nimmanpipug P
    Mol Biotechnol, 2024 Feb 19.
    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.
  6. Aznan AM, Abdullah Z, Lee VS, Tiekink ER
    Acta Crystallogr Sect E Struct Rep Online, 2014 Aug 1;70(Pt 8):58-61.
    PMID: 25249854 DOI: 10.1107/S1600536814012227
    The title compound, C12H11N3O2, is a second monoclinic polymorph (P21, with Z' = 4) of the previously reported monoclinic (P21/c, with Z' = 2) form [Akhmad Aznan et al. (2010 ▶). Acta Cryst. E66, o2400]. Four independent mol-ecules comprise the asymmetric unit, which have the common features of a syn disposition of the pyridine N atom and the toluene ring, and an intra-molecular amine-nitro N-H⋯O hydrogen bond. The differences between mol-ecules relate to the dihedral angles between the rings which range from 2.92 (19) to 26.24 (19)°. The geometry-optimized structure [B3LYP level of theory and 6-311 g+(d,p) basis set] has the same features except that the entire mol-ecule is planar. In the crystal, the three-dimensional architecture is consolidated by a combination of C-H⋯O, C-H⋯π, nitro-N-O⋯π and π-π inter-actions [inter-centroid distances = 3.649 (2)-3.916 (2) Å].
  7. Nokinsee D, Shank L, Lee VS, Nimmanpipug P
    Enzyme Res, 2015;2015:262364.
    PMID: 26788364 DOI: 10.1155/2015/262364
    Tyrosinase is a key enzyme in melanogenesis. Generally, mushroom tyrosinase from A. bisporus had been used as a model in skin-whitening agent tests employed in the cosmetic industry. The recently obtained crystal structure of bacterial tyrosinase from B. megaterium has high similarity (33.5%) to the human enzyme and thus it was used as a template for constructing of the human model. Binding of tyrosinase to a series of its inhibitors was simulated by automated docking calculations. Docking and MD simulation results suggested that N81, N260, H263, and M280 are involved in the binding of inhibitors to mushroom tyrosinase. E195 and H208 are important residues in bacterial tyrosinase, while E230, S245, N249, H252, V262, and S265 bind to inhibitors and are important in forming pi interaction in human tyrosinase.
  8. 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.
  9. Wongrattanakamon P, Lee VS, Nimmanpipug P, Jiranusornkul S
    Data Brief, 2016 Dec;9:35-42.
    PMID: 27626051 DOI: 10.1016/j.dib.2016.08.004
    The data is obtained from exploring the modulatory activities of bioflavonoids on P-glycoprotein function by ligand-based approaches. Multivariate Linear-QSAR models for predicting the induced/inhibitory activities of the flavonoids were created. Molecular descriptors were initially used as independent variables and a dependent variable was expressed as pFAR. The variables were then used in MLR analysis by stepwise regression calculation to build the linear QSAR data. The entire dataset consisted of 23 bioflavonoids was used as a training set. Regarding the obtained MLR QSAR model, R of 0.963, R (2)=0.927, [Formula: see text], SEE=0.197, F=33.849 and q (2)=0.927 were achieved. The true predictabilities of QSAR model were justified by evaluation with the external dataset (Table 4). The pFARs of representative flavonoids were predicted by MLR QSAR modelling. The data showed that internal and external validations may generate the same conclusion.
  10. Wisitponchai T, Shoombuatong W, Lee VS, Kitidee K, Tayapiwatana C
    BMC Bioinformatics, 2017 Apr 19;18(1):220.
    PMID: 28424069 DOI: 10.1186/s12859-017-1628-6
    BACKGROUND: Computational analysis of protein-protein interaction provided the crucial information to increase the binding affinity without a change in basic conformation. Several docking programs were used to predict the near-native poses of the protein-protein complex in 10 top-rankings. The universal criteria for discriminating the near-native pose are not available since there are several classes of recognition protein. Currently, the explicit criteria for identifying the near-native pose of ankyrin-protein complexes (APKs) have not been reported yet.

    RESULTS: In this study, we established an ensemble computational model for discriminating the near-native docking pose of APKs named "AnkPlex". A dataset of APKs was generated from seven X-ray APKs, which consisted of 3 internal domains, using the reliable docking tool ZDOCK. The dataset was composed of 669 and 44,334 near-native and non-near-native poses, respectively, and it was used to generate eleven informative features. Subsequently, a re-scoring rank was generated by AnkPlex using a combination of a decision tree algorithm and logistic regression. AnkPlex achieved superior efficiency with ≥1 near-native complexes in the 10 top-rankings for nine X-ray complexes compared to ZDOCK, which only obtained six X-ray complexes. In addition, feature analysis demonstrated that the van der Waals feature was the dominant near-native pose out of the potential ankyrin-protein docking poses.

    CONCLUSION: The AnkPlex model achieved a success at predicting near-native docking poses and led to the discovery of informative characteristics that could further improve our understanding of the ankyrin-protein complex. Our computational study could be useful for predicting the near-native poses of binding proteins and desired targets, especially for ankyrin-protein complexes. The AnkPlex web server is freely accessible at http://ankplex.ams.cmu.ac.th .

  11. Yeap BK, Othman AS, Lee VS, Lee CY
    J Econ Entomol, 2007 Apr;100(2):467-74.
    PMID: 17461072
    The phylogenetic relationship of Coptotermes gestroi (Wasmann) and Coptotermes vastator Light (Isoptera: Rhinotermitidae) was determined using DNA sequence comparisons of mitochondrial genes. Partial sequences of the ribosomal RNA small subunit 12S, ribosomal RNA large subunit 16S, and mitochondrial COII were obtained from nine populations of C. gestroi from South East Asia (Malaysia, Singapore, Thailand, and Indonesia) and four populations of C. vastator from the Philippines and Hawaii. In addition, four populations of Coptotermes formosanus Shiraki and Globitermes sulphureus (Haviland) were used as the outgroups. Consensus sequences were obtained and aligned. C. vastator and C. gestroi are synonymous, based on high sequence homology across the 12S, 16S, and COII genes. The interspecific pairwise sequence divergence, based on Kimura 2-parameter model between C. gestroi and C. vastator, varied only up to 0.80%. Morphometric measurements of 16 characteristics revealed numerous overlaps between the examined individuals of both species. Based on the molecular phylogenetics and morphometric data, it is proposed that C. vastator is a junior synonym of C. gestroi.
  12. Chong WL, Saparpakorn P, Sangma C, Lee VS, Hannongbua S
    Heliyon, 2023 Jan;9(1):e12667.
    PMID: 36618128 DOI: 10.1016/j.heliyon.2022.e12667
    SARS-CoV-2 virus continues to evolve and mutate causing most of the mutated variants resist to many of the therapeutic monoclonal antibodies (mAbs). Despite several mAbs retained neutralizing capability for Omicron BA.1 and BA.2, reduction in neutralization potency was reported. Hence, effort of searching for mAb that is broader in neutralization breadth without losing the neutralizing ability is continued. MW06 was reported with capability in neutralizing most of the variants of concern (VOC) and it binds to the conserved region (left flank) near epitope mAb sotrovimab (S309). In this study, binding affinity of mAb MW06 and its cocktail formulation with MW05 for receptor binding domain (RBD) SARS-CoV-2 virus was investigated under molecular dynamics simulations (MDs). Binding free energies computed by Molecular Mechanics Generalised Born Surface Area (MM-GBSA) algorithm predicted the binding affinity of MW06 for RBD BA.1 (-53 kcal/mol) as strong as RBD wildtype (-58 kcal/mol) while deterioration was observed for RBD BA.2 (-43 kcal/mol). Alike S309 and MW06, simulated cocktail mAb (MW05 and MW06)-RBD interactions suggested the neutralizing capability of the cocktail formulation for RBD BA.1 and BA.2 reduced. Meanwhile, residue pairs that favour the communication between the mAb and RBD have been identified by decomposing the free energy per pairwise residue basis. Apart from understanding the effects of mutation occurred in the RBD region on human angiotensin-converting enzyme 2 (hACE2) binding, impact of heavily mutated RBD on mAb-RBD interactions was investigated in this study as well. In addition to energetic profile obtained from MDs, plotting the dynamics cross-correlation map of the mAb-RBD complex under elastic network model (ENM) was aimed to understand the cross-correlations between residue fluctuations. It allows simple and rapid analysis on the motions or dynamics of the protein residues of mAbs and RBD in complex. Protein residues having correlated motions are normally part of the structural domains of the protein and their respective motions and protein function are related. Motion of mutated RBD residues and mAb residues was less correlated while their respective interactions energy computed to be higher. The combined techniques of MDs and ENM offered simplicity in understanding dynamics and energy contribution that explain binding affinity of mAb-RBD complexes.
  13. Charoenkwan P, Chotpatiwetchkul W, Lee VS, Nantasenamat C, Shoombuatong W
    Sci Rep, 2021 Dec 10;11(1):23782.
    PMID: 34893688 DOI: 10.1038/s41598-021-03293-w
    Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906-0.910) and 2-17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs.
  14. Rosli AN, Ahmad MR, Alias Y, Zain SM, Lee VS, Woi PM
    J Mol Model, 2014 Dec;20(12):2533.
    PMID: 25433601 DOI: 10.1007/s00894-014-2533-9
    Design of neutral receptor molecules (ionophores) for beryllium(II) using unsaturated carbonitrile models has been carried out via density functional theory, G3, and G4 calculations. The first part of this work focuses on gas phase binding energies between beryllium(II) and 2-cyano butadiene (2-CN BD), 3-cyano propene (3-CN P), and simpler models with two separate fragments; acrylonitrile and ethylene. Interactions between beryllium(II) and cyano nitrogen and terminal olefin in the models have been examined in terms of geometrical changes, distribution of charge over the entire π-system, and rehybridization of vinyl carbon orbitals. NMR shieldings and vibrational frequencies probed charge centers and strength of interactions. The six-membered cyclic complexes have planar structures with the rehybridized carbon slightly out of plane (16° in 2-CN BD). G3 results show that in 2-CN BD complex participation of vinyl carbon further stabilizes the cyclic adduct by 16.3 kcal mol(-1), whereas, in simpler models, interaction between beryllium(II) and acetonitrile is favorable by 46.4 kcal mol(-1) compared with that of ethylene. The terminal vinyl carbon in 2-CN BD rehybridizes to sp (3) with an increase of 7 % of s character to allow interaction with beryllium(II). G4 calculations show that the Be(II) and 2-CN BD complex is more strongly bound than those with Mg(II) and Ca(II) by 98.5 and 139.2 kcal mol(-1) (-1), respectively. QST2 method shows that the cyclic and acyclic forms of Be(II)-2-CN BD complexes are separated by 12.3 kcal mol(-1) barrier height. Overlap population analysis reveals that Ca(II) can be discriminated based on its tendency to form ionic interaction with the receptor models.
  15. Fatima A, Abdul AB, Abdullah R, Karjiban RA, Lee VS
    Int J Mol Sci, 2015 Jan 26;16(2):2747-66.
    PMID: 25629232 DOI: 10.3390/ijms16022747
    Breast cancer is the second most common cancer among women worldwide. Several signaling pathways have been implicated as causative and progression agents. The tumor necrosis factor (TNF) α protein plays a dual role in promoting and inhibiting cancer depending largely on the pathway initiated by the binding of the protein to its receptor. Zerumbone, an active constituent of Zingiber zerumbet, Smith, is known to act on the tumor necrosis factor pathway upregulating tumour necrosis factor related apoptosis inducing ligand (TRAIL) death receptors and inducing apoptosis in cancer cells. Zerumbone is a sesquiterpene that is able to penetrate into the hydrophobic pockets of proteins to exert its inhibiting activity with several proteins. We found a good binding with the tumor necrosis factor, kinase κB (IKKβ) and the Nuclear factor κB (NF-κB) component proteins along the TNF pathway. Our results suggest that zerumbone can exert its apoptotic activities by inhibiting the cytoplasmic proteins. It inhibits the IKKβ kinase that activates the NF-κB and also binds to the NF-κB complex in the TNF pathway. Blocking both proteins can lead to inhibition of cell proliferating proteins to be downregulated and possibly ultimate induction of apoptosis.
  16. Woi PM, Bakar MA, Rosli AN, Lee VS, Ahmad MR, Zain S, et al.
    J Mol Model, 2014 May;20(5):2219.
    PMID: 24770548 DOI: 10.1007/s00894-014-2219-3
    DFT and G4 results reveal that cations display the following trends in imparting its positive charge to acrylonitrile; H⁺ > Li⁺ > Na⁺ > K⁺ for group I and Be²⁺ > Mg²⁺ > Ca²⁺ for group II. Solvation by water molecules and interaction with cation make the cyano bond more polarized and exhibits ketene-imine character. Bond order in nitrile-cation complexes has been predicted based on the s character of the covalent bond orbitals. Mulliken, CHELPG, and NPA charges are in good agreement in predicting positive charge buildup and GIAO nuclear deshileding on C1. G4 enthalpies show that Mg²⁺ is more strongly bound to acrylonitrile than to acetonitrile by 3 kcal mol⁻¹, and the proton affinity of the former is higher by 0.8 kcal mol⁻¹. G4 enthalpies of reductions support prior experimental observation that metalated conjugated nitriles show enhanced reactivity toward weak nucleophiles to afford Michael addition products.
  17. Abdallah HH, Mavri J, Repič M, Lee VS, Wahab HA
    Int J Mol Sci, 2012;13(2):1269-83.
    PMID: 22408390 DOI: 10.3390/ijms13021269
    Genistein, daidzein, glycitein and quercetin are flavonoids present in soybean and other vegetables in high amounts. These flavonoids can be metabolically converted to more active forms, which may react with guanine in the DNA to form complexes and can lead to DNA depurination. We assumed two ultimate carcinogen forms of each of these flavonoids, diol epoxide form and diketone form. Density functional theory (DFT) and Hartree-Fock (HF) methods were used to study the reaction thermodynamics between active forms of flavonoids and DNA guanine. Solvent reaction field method of Tomasi and co-workers and the Langevin dipoles method of Florian and Warshel were used to calculate the hydration free energies. Activation free energy for each reaction was estimated using the linear free energy relation. Our calculations show that diol epoxide forms of flavonoids are more reactive than the corresponding diketone forms and are hence more likely flavonoid ultimate carcinogens. Genistein, daidzein and glycitein show comparable reactivity while quercetin is less reactive toward DNA.
  18. Abdullah NH, Thomas NF, Sivasothy Y, Lee VS, Liew SY, Noorbatcha IA, et al.
    Int J Mol Sci, 2016 Feb 14;17(2):143.
    PMID: 26907251 DOI: 10.3390/ijms17020143
    The mammalian hyaluronidase degrades hyaluronic acid by the cleavage of the β-1,4-glycosidic bond furnishing a tetrasaccharide molecule as the main product which is a highly angiogenic and potent inducer of inflammatory cytokines. Ursolic acid 1, isolated from Prismatomeris tetrandra, was identified as having the potential to develop inhibitors of hyaluronidase. A series of ursolic acid analogues were either synthesized via structure modification of ursolic acid 1 or commercially obtained. The evaluation of the inhibitory activity of these compounds on the hyaluronidase enzyme was conducted. Several structural, topological and quantum chemical descriptors for these compounds were calculated using semi empirical quantum chemical methods. A quantitative structure activity relationship study (QSAR) was performed to correlate these descriptors with the hyaluronidase inhibitory activity. The statistical characteristics provided by the best multi linear model (BML) (R² = 0.9717, R²cv = 0.9506) indicated satisfactory stability and predictive ability of the developed model. The in silico molecular docking study which was used to determine the binding interactions revealed that the ursolic acid analog 22 had a strong affinity towards human hyaluronidase.
Related Terms
Filters
Contact Us

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

External Links