Displaying publications 1 - 20 of 245 in total

Abstract:
Sort:
  1. Unraveling the Structural Dynamics of an Enzyme Encapsulated within a Metal-Organic Framework
    Tuan Kob TNA, Ismail MF, Abdul Rahman MB, Cordova KE, Mohammad Latif MA
    J Phys Chem B, 2020 05 07;124(18):3678-3685.
    PMID: 32275422 DOI: 10.1021/acs.jpcb.0c02145
    Herein, we detail an atomic-level investigation of the cutinase enzyme encapsulated within a model metal-organic framework (MOF) platform using quantum mechanics calculations and molecular dynamics simulations. Cutinase, when encapsulated in an isoreticularly expanded MOF-74 (cutinase@IRMOF-74-VI), was proven to maintain its structural stability at temperatures that would otherwise denature the enzyme in its unprotected native state. Hydrogen bonding and salt bridge interactions, most notably involving arginine residues at the surface of the enzyme, were critical for stabilizing cutinase within the pore channels of IRMOF-74-VI. The findings reported support the viability of enzyme encapsulation in a porous material by demonstrating that a model enzyme not only retains its structural integrity but also remains accessible and active under extreme and foreign conditions.
    Matched MeSH terms: Molecular Dynamics Simulation
  2. 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*
  3. Molecular docking and molecular dynamics simulations of a mutant Acinetobacter haemolyticus alkaline-stable lipase against tributyrin
    Anuar NFSK, Wahab RA, Huyop F, Amran SI, Hamid AAA, Halim KBA, et al.
    J Biomol Struct Dyn, 2021 Apr;39(6):2079-2091.
    PMID: 32174260 DOI: 10.1080/07391102.2020.1743364
    We previously reported on a mutant lipase KV1 (Mut-LipKV1) from Acinetobacter haemolyticus which optimal pH was raised from 8.0 to 11.0 after triple substitutions of surface aspartic acid (Asp) with lysine (Lys). Herein, this study further examined the Mut-LipKV1 by molecular docking, molecular dynamics (MD) simulations and molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) calculations to explore the structural requirements that participated in the effective binding of tributyrin and its catalytic triad (Ser165, Asp259 and His289) and identify detailed changes that occurred post mutation. Mut-LipKV1 bound favorably with tributyrin (-4.1 kcal/mol) and formed a single hydrogen bond with His289, at pH 9.0. Despite the incongruent docking analysis data, results of MD simulations showed configurations of both the tributyrin-Mut-LipKV1 (RMSD 0.3 nm; RMSF 0.05 - 0.3 nm) and the tributyrin-wildtype lipase KV1 (tributyrin-LipKV1) complexes (RMSD 0.35 nm; RMSF 0.05 - 0.4 nm) being comparably stable at pH 8.0. MM-PBSA analysis indicated that van der Waals interactions made the most contribution during the molecular binding process, with the Mut-LipKV1-tributyrin complex (-44.04 kcal/mol) showing relatively lower binding energy than LipKV1-tributyrin (-43.83 kcal/mol), at pH 12.0. All tributyrin-Mut-LipKV1 complexes displayed improved binding free energies over a broader pH range from 8.0 - 12.0, as compared to LipKV1-tributyrin. Future empirical works are thus, important to validate the improved alkaline-stability of Mut-LipKV1. In a nutshell, our research offered a considerable insight for further improving the alkaline tolerance of lipases.Communicated by Ramaswamy H. Sarma.
    Matched MeSH terms: Molecular Dynamics Simulation*
  4. Fulltext Molecular Dynamics Simulations Reveal Novel Interacting Regions of Human Prion Protein to Brucella abortus Hsp60 Protein
    Le-Dao HA, Dinh TT, Tran TL, Lee VS, Tran-Van H
    Mol Biotechnol, 2024 Apr;66(4):687-695.
    PMID: 36633832 DOI: 10.1007/s12033-023-00655-9
    The distinctive morphology characteristics of microfold cells (M cells) allow the vaccine antigen not only to interact with immune cells directly, but also to effectively stimulate mucosal immune responses via receptors on its apical surface. Human prion protein, a transmembrane receptor for Brucella abortus Hsp60, is highly expressed on the M cell surface. Nonetheless, this protein tends to express in inclusion body in prokaryotic hosts. In this study, the shorter interacting regions of human prion protein were identified via computational methods such as docking and molecular dynamics simulations to minimize its aggregation tendency. The computational calculations revealed three novel human prion protein-interacting regions, namely PrP125, PrP174, and PrP180. In accordance with in silico prediction, the biologically synthesized peptides fusing with GST tag demonstrated their specific binding to Hsp60 protein via pull-down assay. Hence, this finding laid the groundwork for M-cell targeting candidate validation through these newly identified interacting regions.
    Matched MeSH terms: Molecular Dynamics Simulation*
  5. Interaction of Asn297-Linked Glycan Ligands with the Fc Fragment of the Immunoglobulin Class G1: A Molecular Dynamics Simulation Study
    Tong M, Liu P, Li C, Zhang Z, Sun W, Dong P, et al.
    J Chem Inf Model, 2024 Feb 12;64(3):785-798.
    PMID: 38262973 DOI: 10.1021/acs.jcim.3c01584
    The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.
    Matched MeSH terms: Molecular Dynamics Simulation*
  6. Structural investigation, computational analysis, and theoretical cryoprotectant approach of antifreeze protein type IV mutants
    Eskandari A, Leow TC, Rahman MBA, Oslan SN
    Eur Biophys J, 2024 Nov;53(7-8):385-403.
    PMID: 39327310 DOI: 10.1007/s00249-024-01719-7
    Antifreeze proteins (AFPs) have unique features to sustain life in sub-zero environments due to ice recrystallization inhibition (IRI) and thermal hysteresis (TH). AFPs are in demand as agents in cryopreservation, but some antifreeze proteins have low levels of activity. This research aims to improve the cryopreservation activity of an AFPIV. In this in silico study, the helical peptide afp1m from an Antarctic yeast AFP was modeled into a sculpin AFPIV, to replace each of its four α-helices in turn, using various computational tools. Additionally, a new linker between the first two helices of AFPIV was designed, based on a flounder AFPI, to boost the ice interaction activity of the mutants. Bioinformatics tools such as ExPASy Prot-Param, Pep-Wheel, SOPMA, GOR IV, Swiss-Model, Phyre2, MODFOLD, MolPropity, and ProQ were used to validate and analyze the structural and functional properties of the model proteins. Furthermore, to evaluate the AFP/ice interaction, molecular dynamics (MD) simulations were executed for 20, 100, and 500 ns at various temperatures using GROMACS software. The primary, secondary, and 3D modeling analysis showed the best model for a redesigned antifreeze protein (AFP1mb, with afp1m in place of the fourth AFPIV helix) with a QMEAN (Swiss-Model) Z score value of 0.36, a confidence of 99.5%, a coverage score of 22%, and a p value of 0.01. The results of the MD simulations illustrated that AFP1mb had more rigidity and better ice interactions as a potential cryoprotectant than the other models; it also displayed enhanced activity in limiting ice growth at different temperatures.
    Matched MeSH terms: Molecular Dynamics Simulation*
  7. Nonequilibrium Langevin dynamics: A demonstration study of shear flow fluctuations in a simple fluid
    Belousov R, Cohen EGD, Rondoni L
    Phys Rev E, 2017 Aug;96(2-1):022125.
    PMID: 28950465 DOI: 10.1103/PhysRevE.96.022125
    The present paper is based on a recent success of the second-order stochastic fluctuation theory in describing time autocorrelations of equilibrium and nonequilibrium physical systems. In particular, it was shown to yield values of the related deterministic parameters of the Langevin equation for a Couette flow in a microscopic molecular dynamics model of a simple fluid. In this paper we find all the remaining constants of the stochastic dynamics, which then is simulated numerically and compared directly with the original physical system. By using these data, we study in detail the accuracy and precision of a second-order Langevin model for nonequilibrium physical systems theoretically and computationally. We find an intriguing relation between an applied external force and cumulants of the resulting flow fluctuations. This is characterized by a linear dependence of an athermal cumulant ratio, an apposite quantity introduced here. In addition, we discuss how the order of a given Langevin dynamics can be raised systematically by introducing colored noise.
    Matched MeSH terms: Molecular Dynamics Simulation
  8. Molecular dynamics revealed the effect of epoxy group on triglyceride digestion
    Li X, Liu YJ, Nian BB, Cao XY, Tan CP, Liu YF, et al.
    Food Chem, 2022 Mar 30;373(Pt B):131285.
    PMID: 34740049 DOI: 10.1016/j.foodchem.2021.131285
    The digestion behavior of epoxy triglyceride, the main cytotoxic product of deep-frying oil, remains unknown, which may affect its biosafety. In this study, epoxy triglyceride (EGT) and triglyceride (GT) were used to reveal the effect of epoxy group on digestion. Digestibility rate analysis showed that the free fatty acids release rate of EGT was slower. To clarify this phenomenon, binding ability with salt ions in digestive juice and particle size were also been studied. Cluster size analysis indicated that epoxy group increased triglyceride particle size, resulting in smaller contact area between EGT and lipase. Interface behaviors displayed EGT decreased binding ability with salt ions in digestive juice. Spectroscopic analysis showed EGT caused the red shift of lipase peak, indicating that epoxy group changed lipase structure. Molecular dynamics simulation suggested EGT leads to loosen lipase structure. In conclusion, this study highlights that epoxy group could weaken the triglyceride digestion.
    Matched MeSH terms: Molecular Dynamics Simulation*
  9. Special issue: Biomolecular Modelling and Simulation
    Choong YS, Mancera R, Lee VS
    Mol Biotechnol, 2024 Apr;66(4):567.
    PMID: 38337130 DOI: 10.1007/s12033-024-01073-1
    Matched MeSH terms: Molecular Dynamics Simulation*
  10. Molecular dynamics studies on T1 lipase: insight into a double-flap mechanism
    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*
  11. 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*
  12. Fulltext birgHPC: creating instant computing clusters for bioinformatics and molecular dynamics
    Chew TH, Joyce-Tan KH, Akma F, Shamsir MS
    Bioinformatics, 2011 May 1;27(9):1320-1.
    PMID: 21398666 DOI: 10.1093/bioinformatics/btr109
    birgHPC, a bootable Linux Live CD has been developed to create high-performance clusters for bioinformatics and molecular dynamics studies using any Local Area Network (LAN)-networked computers. birgHPC features automated hardware and slots detection as well as provides a simple job submission interface. The latest versions of GROMACS, NAMD, mpiBLAST and ClustalW-MPI can be run in parallel by simply booting the birgHPC CD or flash drive from the head node, which immediately positions the rest of the PCs on the network as computing nodes. Thus, a temporary, affordable, scalable and high-performance computing environment can be built by non-computing-based researchers using low-cost commodity hardware.
    Matched MeSH terms: Molecular Dynamics Simulation*
  13. In silico analysis of a putative dehalogenase from the genome of halophilic bacterium Halomonas smyrnensis AAD6T
    Oyewusi HA, Akinyede KA, Abdul Wahab R, Huyop F
    J Biomol Struct Dyn, 2023 Jan;41(1):319-335.
    PMID: 34854349 DOI: 10.1080/07391102.2021.2006085
    Microbial-assisted removal of natural or synthetic pollutants is the prevailing green, low-cost technology to treat polluted environments. However, the challenge with enzyme-assisted bioremediation is the laborious nature of dehalogenase-producing microorganisms' bioprospecting. This bottleneck could be circumvented by in-silico analysis of certain microorganisms' whole-genome sequences to predict their protein functions and enzyme versatility for improved biotechnological applications. Herein, this study performed structural analysis on a dehalogenase (DehHsAAD6) from the genome of Halomonas smyrnensis AAD6 by molecular docking and molecular dynamic (MD) simulations. Other bioinformatics tools were also employed to identify substrate preference (haloacids and haloacetates) of the DehHsAAD6. The DehHsAAD6 preferentially degraded haloacids and haloacetates (-3.2-4.8 kcal/mol) and which formed three hydrogen bonds with Tyr12, Lys46, and Asp182. MD simulations data revealed the higher stability of DehHsAAD6-haloacid- (RMSD 0.22-0.3 nm) and DehHsAAD6-haloacetates (RMSF 0.05-0.14 nm) complexes, with the DehHsAAD6-L-2CP complex being the most stable. The detail of molecular docking calculations ranked complexes with the lowest binding free energies as: DehHsAAD6-L-2CP complex (-4.8 kcal/mol) = DehHsAAD6-MCA (-4.8 kcal/mol) < DehHsAAD6-TCA (-4.5 kcal/mol) < DehHsAAD6-2,3-DCP (-4.1 kcal/mol) < DehHsAAD6-D-2CP (-3.9 kcal/mol) < DehHsAAD6-2,2-DCP (-3.5 kcal/mol) < DehHsAAD6-3CP (-3.2 kcal/mol). In a nutshell, the study findings offer valuable perceptions into the elucidation of possible reaction mechanisms of dehalogenases for extended substrate specificity and higher catalytic activity.Communicated by Ramaswamy H. Sarma.
    Matched MeSH terms: Molecular Dynamics Simulation*
  14. Molecular dynamics simulation and structural analysis of aquaporin Z from an Antarctic Pseudomonas sp. strain AMS3
    Balakrishnan S, Rahman RNZRA, Noor NDM, Latip W, Ali MSM
    J Biomol Struct Dyn, 2023;41(21):11498-11509.
    PMID: 36598349 DOI: 10.1080/07391102.2022.2164519
    Aquaporin is a water channel protein that facilitates the movement of water across the cell membrane. Aquaporin from the Antarctic region has been noted for its psychrophilic properties and its ability to perform at a lower temperature but there remains limited understanding of the water mechanism of Antarctic Pseudomonas sp. strain AMS3 However, studies regarding aquaporin isolated from psychrophilic Pseudomonas sp. are still scattered. Recently, the genome sequence of an Antarctic Pseudomonas sp. strain AMS3 revealed a gene sequence encoding for a putative aquaporin designated as AqpZ1 AMS3. In this study, structure analysis and a molecular dynamics (MD) simulation of a predicted model of a fully hydrated aquaporin tetramer embedded in a lipid bilayer was performed at different temperatures for structural flexibility and stability analysis. The MD simulation results revealed that the structures were able to remain stable at low to medium temperatures. The protein was observed to have high flexibility in the loop region as compared to the helices region throughout the simulated temperatures. The selectivity filter and NPA motifs play a major role in solute selectivity and the pore radius of the protein. The structural and functional characterization of this psychrophilic aquaporin provides new insights for the future applications of this protein.Communicated by Ramaswamy H. Sarma.
    Matched MeSH terms: Molecular Dynamics Simulation*
  15. Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes
    Pirojsirikul T, Lee VS, Nimmanpipug P
    Mol Biotechnol, 2024 Apr;66(4):582-591.
    PMID: 38374320 DOI: 10.1007/s12033-024-01082-0
    We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
    Matched MeSH terms: Molecular Dynamics Simulation*
  16. Study of the antidiabetic mechanism of berberine compound on FOXO1 transcription factor through molecular docking and molecular dynamics simulations
    Maksum IP, Rustaman R, Deawati Y, Rukayadi Y, Utami AR, Nafisa ZK
    J Mol Model, 2024 Jul 09;30(8):260.
    PMID: 38981921 DOI: 10.1007/s00894-024-06060-6
    CONTEXT: Diabetes mellitus (DM) is a metabolic disorder disease that causes hyperglycemia conditions and associated with various chronic complications leading to mortality. Due to high toxicity of conventional diabetic drugs, the exploration of natural compounds as alternative diabetes treatments has been widely carried out. Previous in silico studies have highlighted berberine, a natural compound, as a promising alternative in antidiabetic therapy, potentially acting through various pathways, including the inhibition of the FOXO1 transcription factor in the gluconeogenesis pathway. However, the specific mechanism by which berberine interacts with FOXO1 remains unclear, and research in this area is relatively limited. Therefore, this study aims to determine the stability of berberine structure with FOXO1 based on RMSD, RMSF, binding energy, and trajectory analysis to determine the potential of berberine to inhibit the gluconeogenesis pathway. This research was conducted by in silico method with molecular docking using AutoDock4.2 and molecular dynamics study using Amber20, then visualized by VMD.

    METHODS: Docking between ligand and FOXO1 receptor was carried out with Autodock4.2. For molecular dynamics simulations, the force fields of DNA.OL15, protein.ff14SB, gaff2, and tip3p were used.

    Matched MeSH terms: Molecular Dynamics Simulation*
  17. N-substituted tetrahydro-beta-carboline as mu-opioid receptors ligands: in silico study; molecular docking, ADMET and molecular dynamics approach
    Alananzeh WA, Al-Qattan MN, Ayipo YO, Mordi MN
    Mol Divers, 2024 Jun;28(3):1273-1289.
    PMID: 37133710 DOI: 10.1007/s11030-023-10655-1
    Manipulating intracellular signals by interaction with transmembranal G-protein-coupled receptors (GPCRs) is the way of action of more than 30% of available medicines. Designing molecules against GPCRs is most challenging due to their flexible binding orthosteric and allosteric pockets, a property that lead to different mode and extent of activation of intracellular mediators. Here, in the current study we aimed to design N-substituted tetrahydro-beta-carbolines (THβC's) targeting Mu Opioid Receptors (MORs). We performed ligand docking study for reference and designed compounds against active and inactive states of MOR, as well as the active state bound to intracellular mediator of Gi. The reference compounds include 40 known agonists and antagonists, while the designed compounds include 25,227 N-substituted THβC analogues. Out of the designed compounds, 15 compounds were comparatively having better extra precision (XP) Gscore and were analyzed for absorption, distribution, metabolism, and excretion-toxicity (ADMET) properties, drug-likness, and molecular dynamic (MD) simulation. The results showed that N-substituted tetrahydro-beta-carbolines with and without C6-methoxy group substitutions (THBC/6MTHBC) analogues of A1/B1 and A9/B9 have relatively acceptable affinity and within pocket-stability toward MOR compared to the reference compounds of morphine (agonist) and naloxone (antagonist). Moreover, the designed analogues interact with key residue within the binding pocket of Asp 147 that is reported to be involved in receptor activation. In conclusion, the designed THBC analogues represent a good starting point for designing opioid receptor ligands other than morphinan scaffold, that have good synthetic accessibility which promotes feasible structural manipulation to tailor pharmacological effects with minimal side effects.
    Matched MeSH terms: Molecular Dynamics Simulation*
  18. Investigation of laccase activity in cholinium-based ionic liquids using experimental and molecular dynamics techniques
    Chan KK, Pereira AF, Valente AI, Tavares APM, Coutinho JAP, Ooi CW
    Int J Biol Macromol, 2024 Oct;277(Pt 4):134443.
    PMID: 39217672 DOI: 10.1016/j.ijbiomac.2024.134443
    Laccases hold great potential for biotechnological applications, particularly in environmental pollutant remediation. Laccase activity is governed by the solvent environment, and ionic liquids (ILs) emerge as a versatile solvent for activation or stabilization of enzymes. Herein, effects of cholinium-based ILs formulated with carboxylic acids, inorganic acid, and amino acids as anionic species, on the catalytic activity of laccase from Trametes versicolor were investigated by experimental and computational approaches. Experimental results showed that laccase activity was enhanced by 21.39 % in 0.5 M cholinium dihydrogen citrate ([Cho][DHC]), in relation to the laccase activity in phosphate buffer medium. However, cholinium aminoate ILs negatively affected laccase activity, as evidenced by the partial deactivation of laccase in both cholinium glycinate and cholinium phenylalaninate, at concentrations of 0.1 M and 0.5 M, respectively. Molecular dynamics studies revealed that the enhancement of laccase activity in [Cho][DHC] might be attributed to the highly stabilized and compact structure of laccase, facilitating a better internal electron transfer during the laccase-substrate interactions. Enhanced catalytic performance of laccase in [Cho][DHC] was postulated to be driven by the high accumulation level of dihydrogen citrate anions around laccase's surface. [Cho][DHC] holds great promise as a cosolvent in laccase-catalyzed biochemical reactions.
    Matched MeSH terms: Molecular Dynamics Simulation*
  19. 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
  20. Fulltext New reaction kinetic model derived from molecular dynamics simulation data on non-linearity of rate coefficients and their relation to activity coefficients
    Jesudason, C.G.
    ASM Science Journal, 2007;1(1):7-18.
    MyJurnal
    Molecular dynamics reaction simulation showed that the rate constant is not constant over the concentration profile of reactants and products over a fixed temperature regime, and this variation is expressed in terms of the defined reactivity coefficients. The ratio of these coefficients for the forward and backward reactions were found to equal that of the activity coefficient ratio for the product and reactant species. A theory was developed to explain kinetics in general based on these observations. Several other theorems had first to be developed, most striking of all was the inference that the excess Helmholtz free energy was the thermodynamical function which had a direct relation to these activity factors than the Gibbs free energy. The theory is applied to a class of ionic reactions which could not be rationalized using the standard Bjørn-Bjerrum theory of ionic reactions.
    Matched MeSH terms: Molecular Dynamics Simulation
Related Terms
Filters
Contact Us

Please provide feedback to Administrator ([email protected])

External Links