Displaying publications 1 - 20 of 224 in total

Abstract:
Sort:
  1. A molecular dynamics approach to explore the structural characterization of cataract causing mutation R58H on human γD crystallin
    Karunakaran R, Srikumar PS
    Mol Cell Biochem, 2018 Dec;449(1-2):55-62.
    PMID: 29532225 DOI: 10.1007/s11010-018-3342-8
    The crystallins are a family of monomeric proteins present in the mammalian lens and mutations in these proteins cause various forms of cataracts. The aim of our current study is to emphasize the structural characterization of aggregation propensity of mutation R58H on γD crystallin using molecular dynamics (MD) approach. MD result revealed that difference in the sequence level display a wide variation in the backbone atomic position, and thus exhibits rigid conformational dynamics. Changes in the flexibility of residues favoured to increase the number of intra-molecular hydrogen bonds in mutant R58H. Moreover, notable changes in the hydrogen bonding interaction resulted to cause the misfolding of mutant R58H by introducing α-helix. Principal component analysis (PCA) result suggested that mutant R58H showed unusual conformational dynamics along the two principal components when compared to the wild-type (WT)-γD crystallin. In a nutshell, the increased surface hydrophobicity could be the cause of self-aggregation of mutant R58H leading to aculeiform cataract.
    Matched MeSH terms: Molecular Dynamics Simulation*
  2. A molecular dynamics study of Beta-Glucosidase B upon small substrate binding
    Mazlan NS, Ahmad Khairudin NB
    J Biomol Struct Dyn, 2016 Jul;34(7):1486-94.
    PMID: 26261863 DOI: 10.1080/07391102.2015.1081570
    Paenibacillus polymyxa β-glucosidase B (BglB), belongs to a GH family 1, is a monomeric enzyme that acts as an exo-β-glucosidase hydrolysing cellobiose and cellodextrins of higher degree of polymerization using retaining mechanism. A molecular dynamics (MD) simulation was performed at 300 K under periodic boundary condition for 5 ns using the complexes structure obtained from previous docking study, namely BglB-Beta-d-glucose and BglB-Cellobiose. From the root-mean-square deviation analysis, both enzyme complexes were reported to deviate from the initial structure in the early part of the simulation but it was stable afterwards. The root-mean-square fluctuation analysis revealed that the most flexible regions comprised of the residues from 26 to 29, 43 to 53, 272 to 276, 306 to 325 and 364 to 367. The radius of gyration analysis had shown the structure of BglB without substrate became more compact towards the end of the simulation compare to other two complexes. The residues His122 and Trp410 were observed to form stable hydrogen bond with occupancy higher than 10%. In conclusion, the behaviour of BglB enzyme towards the substrate binding was successfully explored via MD simulation approaches.
    Matched MeSH terms: Molecular Dynamics Simulation*
  3. A theoretical study of supramolecular aggregation of polydopamine tetramer subunits in aqueous solution
    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*
  4. Active Site Flexibility of Mycobacterium tuberculosis Isocitrate Lyase in Dimer Form
    Lee YV, Choi SB, Wahab HA, Choong YS
    J Chem Inf Model, 2017 09 25;57(9):2351-2357.
    PMID: 28820943 DOI: 10.1021/acs.jcim.7b00265
    Tuberculosis (TB) still remains a global threat due to the emergence of a drug-resistant strain. Instead of focusing on the drug target of active stage TB, we are highlighting the isocitrate lyase (ICL) at the dormant stage TB. ICL is one of the persistent factors for Mycobacterium tuberculosis (MTB) to survive during the dormant phase. In addition, the absence of ICL in human has made ICL a potential drug target for TB therapy. However, the dynamic details of ICL which could give insights to the ICL-ligand interaction have yet to be solved. Therefore, a series of ICL dimer dynamics studies through molecular dynamics simulation were performed in this work. The ICL active site entrance gate closure is contributed to by hydrogen bonding and electrostatic interactions with the C-terminal. Analysis suggested that the open-closed behavior of the ICL active site entrance depends on the type of ligand present in the active site. We also observed four residues (Ser91, Asp108, Asp153, and Cys191) which could possibly be the nucleophiles for nucleophilic attack on the cleavage of isocitrate at the C2-C3bond. We hope that the elucidation of ICL dynamics can benefit future works such as lead identification or antibody design against ICL for TB therapeutics.
    Matched MeSH terms: Molecular Dynamics Simulation*
  5. Allosteric pocket of the dengue virus (serotype 2) NS2B/NS3 protease: In silico ligand screening and molecular dynamics studies of inhibition
    Mukhametov A, Newhouse EI, Aziz NA, Saito JA, Alam M
    J Mol Graph Model, 2014 Jul;52:103-13.
    PMID: 25023665 DOI: 10.1016/j.jmgm.2014.06.008
    The allosteric pocket of the Dengue virus (DENV2) NS2B/NS3 protease, which is proximal to its catalytic triad, represents a promising drug target (Othman et al., 2008). We have explored this binding site through large-scale virtual screening and molecular dynamics simulations followed by calculations of binding free energy. We propose two mechanisms for enzyme inhibition. A ligand may either destabilize electronic density or create steric effects relating to the catalytic triad residues NS3-HIS51, NS3-ASP75, and NS3-SER135. A ligand may also disrupt movement of the C-terminal of NS2B required for inter-conversion between the "open" and "closed" conformations. We found that chalcone and adenosine derivatives had the top potential for drug discovery hits, acting through both inhibitory mechanisms. Studying the molecular mechanisms of these compounds might be helpful in further investigations of the allosteric pocket and its potential for drug discovery.
    Matched MeSH terms: Molecular Dynamics Simulation*
  6. Fulltext An S188V mutation alters substrate specificity of non-stereospecific α-haloalkanoic acid dehalogenase E (DehE)
    Hamid AA, Hamid TH, Wahab RA, Omar MS, Huyop F
    PLoS One, 2015;10(3):e0121687.
    PMID: 25816329 DOI: 10.1371/journal.pone.0121687
    The non-stereospecific α-haloalkanoic acid dehalogenase E (DehE) degrades many halogenated compounds but is ineffective against β-halogenated compounds such as 3-chloropropionic acid (3CP). Using molecular dynamics (MD) simulations and site-directed mutagenesis we show here that introducing the mutation S188V into DehE improves substrate specificity towards 3CP. MD simulations showed that residues W34, F37, and S188 of DehE were crucial for substrate binding. DehE showed strong binding ability for D-2-chloropropionic acid (D-2CP) and L-2-chloropropionic acid (L-2CP) but less affinity for 3CP. This reduced affinity was attributed to weak hydrogen bonding between 3CP and residue S188, as the carboxylate of 3CP forms rapidly interconverting hydrogen bonds with the backbone amide and side chain hydroxyl group of S188. By replacing S188 with a valine residue, we reduced the inter-molecular distance and stabilised bonding of the carboxylate of 3CP to hydrogens of the substrate-binding residues. Therefore, the S188V can act on 3CP, although its affinity is less strong than for D-2CP and L-2CP as assessed by Km. This successful alteration of DehE substrate specificity may promote the application of protein engineering strategies to other dehalogenases, thereby generating valuable tools for future bioremediation technologies.
    Matched MeSH terms: Molecular Dynamics Simulation
  7. An in silico approach towards the identification of novel inhibitors of the TLR-4 signaling pathway
    Mahita J, Harini K, Rao Pichika M, Sowdhamini R
    J Biomol Struct Dyn, 2016 Jun;34(6):1345-62.
    PMID: 26264972 DOI: 10.1080/07391102.2015.1079243
    Precise functioning and fine-tuning of Toll-like receptor 4 (TLR4) signaling is a critical requirement for the smooth functioning of the innate immune system, since aberrant TLR4 activation causes excessive production of pro-inflammatory cytokines and interferons. This can result in life threatening conditions such as septic shock and other inflammatory disorders. The TRIF-related adaptor molecule (TRAM) adaptor protein is unique to the TLR4 signaling pathway and abrogation of TRAM-mediated TLR4 signaling is a promising strategy for developing therapeutics aimed at disrupting TRAM interactions with other components of the TLR4 signaling complex. The VIPER motif from the vaccinia virus-producing protein, A46 has been reported to disrupt TRAM-TLR4 interactions. We have exploited this information, in combination with homology modeling and docking approaches, to identify a potential binding site on TRAM lined by the BB loop and αC helix. Virtual screening of commercially available small molecules targeting the binding site enabled to short-list 12 small molecules to abrogate TRAM-mediated TLR4 signaling. Molecular dynamics and molecular mechanics calculations have been performed for the analysis of these receptor-ligand interactions.
    Matched MeSH terms: Molecular Dynamics Simulation
  8. Antitumor effects of berberine against EGFR, ERK1/2, P38 and AKT in MDA-MB231 and MCF-7 breast cancer cells using molecular modelling and in vitro study
    Jabbarzadeh Kaboli P, Leong MP, Ismail P, Ling KH
    Pharmacol Rep, 2019 Feb;71(1):13-23.
    PMID: 30343043 DOI: 10.1016/j.pharep.2018.07.005
    BACKGROUND: Berberine is an alkaloid plant-based DNA intercalator that affects gene regulation, particularly expression of oncogenic and tumor suppressor proteins. The effects of berberine on different signaling proteins remains to be elucidated. The present study aimed to identify the effects of berberine against key oncogenic proteins in breast cancer cells.

    METHODS: Molecular docking and molecular dynamics simulations were used for EGFR, p38, ERK1/2, and AKT. The effects of berberine and lapatinib on MAPK and PI3K pathways in MDA-MB231 and MCF-7 cells were evaluated using immunoflorescence assays, and the amounts of phosphorylated kinases were compared to total kinases after treating with different concentrations of berberine.

    RESULTS: Simulations showed berberine accurately interacted with EGFR, AKT, P38, and ERK1/2 active sites in silico (scores = -7.57 to -7.92 Kcal/mol) and decreased the levels of active forms of corresponding enzymes in both cell lines; however, berberine binding to p38 showed less stability. Cytotoxicity analysis indicated that MDA-MB231 cells were resistant to berberine compared to MCF-7 cells [72 h IC50 = 50 versus 15 μM, respectively). Also, lapatinib strongly activated AKT but suppressed EGFR in MDA-MB231 cells. The activity of EGFR, AKT, P38, and ERK1/2 were affected by berberine; however, berberine dramatically reduced EGFR and AKT phosphorylation.

    CONCLUSION: By way of its multikinase inhibitory effects, berberine might be a useful replacement for lapatinib, an EGFR inhibitor which can cause acquired drug resistance in patients.

    Matched MeSH terms: Molecular Dynamics Simulation*
  9. Fulltext Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE
    Leong SW, Lim TS, Tye GJ, Ismail A, Aziah I, Choong YS
    J Biol Phys, 2014 Sep;40(4):387-400.
    PMID: 25011632 DOI: 10.1007/s10867-014-9357-9
    In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein's β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.
    Matched MeSH terms: Molecular Dynamics Simulation*
  10. Assessing GPCR homology models constructed from templates of various transmembrane sequence identities: Binding mode prediction and docking enrichment
    Loo JSE, Emtage AL, Ng KW, Yong ASJ, Doughty SW
    J Mol Graph Model, 2017 Dec 29;80:38-47.
    PMID: 29306746 DOI: 10.1016/j.jmgm.2017.12.017
    GPCR crystal structures have become more readily accessible in recent years. However, homology models of GPCRs continue to play an important role as many GPCR structures remain unsolved. The new crystal structures now available provide not only additional templates for homology modelling but also the opportunity to assess the performance of homology models against their respective crystal structures and gain insight into the performance of such models. In this study we have constructed homology models from templates of various transmembrane sequence identities for eight GPCR targets to better understand the relationship between transmembrane sequence identity and model quality. Model quality was assessed relative to the crystal structure in terms of structural accuracy as well as performance in two typical structure-based drug design applications: ligand binding pose prediction and docking enrichment in virtual screening. Crystal structures significantly outperformed homology models in both assessments. Accurate ligand binding pose prediction was possible but difficult to achieve using homology models, even with the use of induced fit docking. In virtual screening using homology models still conferred significant enrichment compared to random selection, with a clear benefit also observed in using models optimized through induced fit docking. Our results indicate that while homology models that are reasonably accurate structurally can be constructed, without significant refinement homology models will be outperformed by crystal structures in ligand binding pose prediction and docking enrichment regardless of the template used, primarily due to the extremely high level of structural accuracy needed for such applications.
    Matched MeSH terms: Molecular Dynamics Simulation
  11. Atomistic simulation studies of the α/β-glucoside and galactoside in anhydrous bilayers: effect of the anomeric and epimeric configurations
    Ahmadi S, Manickam Achari V, Nguan H, Hashim R
    J Mol Model, 2014 Mar;20(3):2165.
    PMID: 24623320 DOI: 10.1007/s00894-014-2165-0
    Fully atomistic molecular dynamics simulation studies of thermotropic bilayers were performed using a set of glycosides namely n-octyl-β-D-glucopyranoside (β-C8Glc), n-octyl-α-D-glucopyranoside (α-C8Glc), n-octyl-β-D-galactopyranoside (β-C8Gal), and n-octyl-α-D-galactopyranoside (α-C8Gal) to investigate the stereochemical relationship of the epimeric/anomeric quartet liner glycolipids with the same octyl chain group. The results showed that, the anomeric stereochemistry or the axial/equatorial orientation of C1-O1 (α/β) is an important factor controlling the area and d-spacing of glycolipid bilayer systems in the thermotropic phase. The head group tilt angle and the chain ordering properties are affected by the anomeric effect. In addition, the L(C) phase of β-C8Gal, is tilting less compared to those in the fluid L(α). The stereochemistry of the C4-epimeric (axial/equatorial) and anomeric (α/β) centers simultaneously influence the inter-molecular hydrogen bond. Thus, the trend in the values of the hydrogen bond for these glycosides is β-C8Gal > α-C8Glc > β-C8Glc > α-C8Gal. The four bilayer systems showed anomalous diffusion behavior with an observed trend for the diffusion coefficients; and this trend is β-C8Gal > β-C8Glc > α-C8Gal > α-C8Glc. The "bent" configuration of the α-anomer results in an increase of the hydrophobic area, chain vibration and chain disorganization. Since thermal energy is dispensed more entropically for the chain region, the overall molecular diffusion decreases.
    Matched MeSH terms: Molecular Dynamics Simulation*
  12. Benchmarking the performance of MM/PBSA in virtual screening enrichment using the GPCR-Bench dataset
    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
  13. Binding specificity of polypeptide substrates in NS2B/NS3pro serine protease of dengue virus type 2: A molecular dynamics Study
    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*
  14. Bioassay-guided identification of an anti-inflammatory prenylated acylphloroglucinol from Melicope ptelefolia and molecular insights into its interaction with 5-lipoxygenase
    Shaari K, Suppaiah V, Wai LK, Stanslas J, Tejo BA, Israf DA, et al.
    Bioorg Med Chem, 2011 Nov 1;19(21):6340-7.
    PMID: 21958738 DOI: 10.1016/j.bmc.2011.09.001
    A bioassay-guided investigation of Melicope ptelefolia Champ ex Benth (Rutaceae) resulted in the identification of an acyphloroglucinol, 2,4,6-trihydroxy-3-geranylacetophenone or tHGA, as the active principle inhibiting soybean 15-LOX. The anti-inflammatory action was also demonstrated on human leukocytes, where the compound showed prominent inhibitory activity against human PBML 5-LOX, with an IC(50) value of 0.42 μM, very close to the effect produced by the commonly used standard, NDGA. The compound concentration-dependently inhibited 5-LOX product synthesis, specifically inhibiting cysteinyl leukotriene LTC(4) with an IC(50) value of 1.80 μM, and showed no cell toxicity effects. The anti-inflammatory action does not seem to proceed via redox or metal chelating mechanism since the compound tested negative for these bioactivities. Further tests on cyclooxygenases indicated that the compound acts via a dual LOX/COX inhibitory mechanism, with greater selectivity for 5-LOX and COX-2 (IC(50) value of 0.40 μM). The molecular features that govern the 5-LOX inhibitory activity was thus explored using in silico docking experiments. The residues Ile 553 and Hie 252 were the most important residues in the interaction, each contributing significant energy values of -13.45 (electrostatic) and -5.40 kcal/mol (electrostatic and Van der Waals), respectively. The hydroxyl group of the phloroglucinol core of the compound forms a 2.56Å hydrogen bond with the side chain of the carboxylate group of Ile 553. Both Ile 553 and Hie 252 are crucial amino acid residues which chelate with the metal ion in the active site. Distorting the geometry of these ligands could be the reason for the inhibition activity shown by tHGA. The molecular simulation studies supported the bioassay results and served as a good model for understanding the way tHGA binds in the active site of human 5-LOX enzyme.
    Matched MeSH terms: Molecular Dynamics Simulation
  15. Biophysical and computational characterization of vandetanib-lysozyme interaction
    Kabir MZ, Hamzah NAB, Ghani H, Mohamad SB, Alias Z, Tayyab S
    Spectrochim Acta A Mol Biomol Spectrosc, 2018 Jan 15;189:485-494.
    PMID: 28843881 DOI: 10.1016/j.saa.2017.08.051
    Interaction of an anticancer drug, vandetanib (VDB) with a ligand transporter, lysozyme (LYZ) was explored using multispectroscopic techniques, such as fluorescence, absorption and circular dichroism along with computational analysis. Fluorescence data and absorption results confirmed VDB-LYZ complexation. VDB-induced quenching was characterized as static quenching based on inverse correlation of KSV with temperature as well as kq values. The complex was characterized by the weak binding constant (Ka=4.96-3.14×103M-1). Thermodynamic data (ΔS=+12.82Jmol-1K-1; ΔH=-16.73kJmol-1) of VDB-LYZ interaction revealed participation of hydrophobic and van der Waals forces along with hydrogen bonds in VDB-LYZ complexation. Microenvironmental perturbations around tryptophan and tyrosine residues as well as secondary and tertiary structural alterations in LYZ upon addition of VDB were evident from the 3-D fluorescence, far- and near-UV CD spectral analyses, respectively. Interestingly, addition of VDB to LYZ significantly increased protein's thermostability. Molecular docking results suggested the location of VDB binding site near the LYZ active site while molecular dynamics simulation results suggested stability of VDB-LYZ complex. Presence of Mg2+, Ba2+ and Zn2+ was found to interfere with VDB-LYZ interaction.
    Matched MeSH terms: Molecular Dynamics Simulation*
  16. Biophysical changes of ATP binding pocket may explain loss of kinase activity in mutant DAPK3 in cancer: A molecular dynamic simulation analysis
    Agarwal T, Annamalai N, Maiti TK, Arsad H
    Gene, 2016 Apr 10;580(1):17-25.
    PMID: 26748242 DOI: 10.1016/j.gene.2015.12.066
    DAPK3 belongs to family of DAPK (death-associated protein kinases) and is involved in the regulation of progression of the cell cycle, cell proliferation, apoptosis and autophagy. It is considered as a tumor suppressor kinase, suggesting the loss of its function in case of certain specific mutations. The T112M, D161N and P216S mutations in DAPK3 have been observed in cancer patients. These DAPK3 mutants have been associated with very low kinase activity, which results in the cellular progression towards cancer. However, a clear understanding of the structural and biophysical variations that occur in DAPK3 with these mutations, resulting in the decreased kinase activity has yet not been deciphered. We performed a molecular dynamic simulation study to investigate such structural variations. Our results revealed that mutations caused a significant structural variation in DAPK3, majorly concentrated in the flexible loops that form part of the ATP binding pocket. Interestingly, D161N and P216S mutations collapsed the ATP binding pocket through flexible loops invasion, hindering ATP binding which resulted in very low kinase activity. On the contrary, T112M mutant DAPK3 reduces ATP binding potential through outward distortion of flexible loops. In addition, the mutant lacked characteristic features of the active protein kinase including proper interaction between HR/FD and DFG motifs, well structured hydrophobic spine and Lys42-Glu64 salt bridge interaction. These observations could possibly explain the underlying mechanism associated with the loss of kinase activity with T112M, D161N and P216S mutation in DAPK3.
    Matched MeSH terms: Molecular Dynamics Simulation
  17. Born-Oppenheimer Molecular Dynamics Study on Proton Dynamics of Strong Hydrogen Bonds in Aspirin Crystals, with Emphasis on Differences between Two Crystal Forms
    Brela MZ, Wójcik MJ, Witek ŁJ, Boczar M, Wrona E, Hashim R, et al.
    J Phys Chem B, 2016 04 28;120(16):3854-62.
    PMID: 27045959 DOI: 10.1021/acs.jpcb.6b01601
    In this study, the proton dynamics of hydrogen bonds for two forms of crystalline aspirin was investigated by the Born-Oppenheimer molecular dynamics (BOMD) method. Analysis of the geometrical parameters of hydrogen bonds using BOMD reveals significant differences in hydrogen bonding between the two crystalline forms of aspirin, Form I and Form II. Analysis of the trajectory for Form I shows spontaneous proton transfer in cyclic dimers, which is absent in Form II. Quantization of the O-H stretching modes allows a detailed discussion on the strength of hydrogen-bonding interactions. The focal point of our study is examination of the hydrogen bond characteristics in the crystal structure and clarification of the influence of hydrogen bonding on the presence of the two crystalline forms of aspirin. In the BOMD method, thermal motions were taken into account. Solving the Schrödinger equation for the snapshots of 2D proton potentials, extracted from MD, gives the best agreement with IR spectra. The character of medium-strong hydrogen bonds in Form I of aspirin was compared with that of weaker hydrogen bonds in aspirin Form II. Two proton minima are present in the potential function for the hydrogen bonds in Form I. The band contours, calculated by using one- and two-dimensional O-H quantization, reflect the differences in the hydrogen bond strengths between the two crystalline forms of aspirin, as well as the strong hydrogen bonding in the cyclic dimers of Form I and the medium-strong hydrogen bonding in Form II.
    Matched MeSH terms: Molecular Dynamics Simulation*
  18. 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
  19. Cell-Penetrating Peptides: Correlation between Peptide-Lipid Interaction and Penetration Efficiency
    Her Choong F, Keat Yap B
    Chemphyschem, 2021 03 03;22(5):493-498.
    PMID: 33377300 DOI: 10.1002/cphc.202000873
    Cell-penetrating peptides are used in the delivery of peptides and biologics, with some cell-penetrating peptides found to be more efficient than others. The exact mechanism of how they interact with the cell membrane and penetrate it, however, remains unclear. This study attempts to investigate the difference in free energy profiles of three cell-penetrating peptides (TAT, CPP1 and CPP9) with a model lipid bilayer (DOPC) using molecular dynamics pulling simulations with umbrella sampling. Potential mean force (PMF) and free energy barrier between the peptides and DOPC are determined using WHAM analysis and MM-PBSA analysis, respectively. CPP9 is found to have the smallest PMF value, followed by CPP1 and TAT, consistent with the experimental data. YDEGE peptide, however, does not give the highest PMF value, although it is a non-cell-permeable peptide. YDEGE is also found to form water pores, alongside with TAT and CPP9, suggesting that it is difficult to distinguish true water pore formation from artefacts arising from pulling simulations. On the contrary, free energy analysis of the peptide-DOPC complex at the lipid-water interface with MM-PBSA provides results consistent with experimental data with CPP9 having the least interaction with DOPC and lowest free energy barrier, followed by CPP1, TAT and YDEGE. These findings suggest that peptide-lipid interaction at the lipid-water interface has a direct correlation with the penetration efficiency of peptides across the lipid bilayer.
    Matched MeSH terms: Molecular Dynamics Simulation*
  20. Characterisation and molecular dynamic simulations of J15 asparaginase from Photobacterium sp. strain J15
    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
Related Terms
Filters
Contact Us

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

External Links