Displaying publications 1 - 20 of 216 in total

Abstract:
Sort:
  1. 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. 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. 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. 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*
  5. 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
  6. 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*
  7. 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*
  8. 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*
  9. 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*
  10. 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*
  11. 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*
  12. 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
  13. 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
  14. Ma L, Liu X, Soh AK, He L, Wu C, Ni Y
    Soft Matter, 2019 May 15.
    PMID: 31090782 DOI: 10.1039/c9sm00507b
    Topological defect nucleation and boundary branching in crystal growth on a curved surface are two typical elastic instabilities driven by curvature induced stress, and have usually been discussed separately in the past. In this work they are simultaneously considered during crystal growth on a sphere. Phase diagrams with respect to sphere radius, size, edge energy and stiffness of the crystal for the equilibrium crystal morphologies are achieved by theoretical analysis and validated by Brownian dynamics simulations. The simulation results further demonstrate the detail of morphological evolution governed by these two different stress relaxation modes. Topological defect nucleation and boundary branching not only compete with each other but also coexist in a range of combinations of factors. Clarification of the interaction mechanism provides a better understanding of various curved crystal morphologies for their potential applications.
    Matched MeSH terms: Molecular Dynamics Simulation
  15. Khan A, Hussain S, Ahmad S, Suleman M, Bukhari I, Khan T, et al.
    Comput Biol Med, 2022 02;141:105163.
    PMID: 34979405 DOI: 10.1016/j.compbiomed.2021.105163
    The spike protein of SARS-CoV-2 and the host ACE2 receptor plays a vital role in the entry to the cell. Among which the hotspot residue 501 is continuously subjected to positive selection pressure and induces unusual virulence. Keeping in view the importance of the hot spot residue 501, we predicted the potentially emerging structural variants of 501 residue. We analyzed the binding pattern of wild type and mutants (Spike RBD) to the ACE2 receptor by deciphering variations in the amino acids' interaction networks by graph kernels along with evolutionary, network metrics, and energetic information. Our analysis revealed that N501I, N501T, and N501V increase the binding affinity and alter the intra and inter-residue bonding networks. The N501T has shown strong positive selection and fitness in other animals. Docking results and repeated simulations (three times) confirmed the structural stability and tighter binding of these three variants, correlated with the previous results following the global stability trend. Consequently, we reported three variants N501I, N501T, and N501V could worsen the situation further if they emerged. The relations between the viral fitness and binding affinity is a complicated game thus the emergence of high affinity mutations in the SARS-CoV-2 RBD brings up the question of whether or not positive selection favours these mutations or not?
    Matched MeSH terms: Molecular Dynamics Simulation
  16. Sayaf AM, Ahmad H, Aslam MA, Ghani SA, Bano S, Yousafi Q, et al.
    Appl Biochem Biotechnol, 2023 Nov;195(11):6959-6978.
    PMID: 36961512 DOI: 10.1007/s12010-023-04466-1
    Because of the essential role of PLpro in the regulation of replication and dysregulation of the host immune sensing, it is considered a therapeutic target for novel drug development. To reduce the risk of immune evasion and vaccine effectiveness, small molecular therapeutics are the best complementary approach. Hence, we used a structure-based drug-designing approach to identify potential small molecular inhibitors for PLpro of SARS-CoV-2. Initial scoring and re-scoring of the best hits revealed that three compounds NPC320891 (2,2-Dihydroxyindene-1,3-Dione), NPC474594 (Isonarciclasine), and NPC474595 (7-Deoxyisonarciclasine) exhibit higher docking scores than the control GRL0617. Investigation of the binding modes revealed that alongside the essential contacts, i.e., Asp164, Glu167, Tyr264, and Gln269, these molecules also target Lys157 and Tyr268 residues in the active site. Moreover, molecular simulation demonstrated that the reported top hits also possess stable dynamics and structural packing. Furthermore, the residues' flexibility revealed that all the complexes demonstrated higher flexibility in the regions 120-140, 160-180, and 205-215. The 120-140 and 160-180 lie in the finger region of PLpro, which may open/close during the simulation to cover the active site and push the ligand inside. In addition, the total binding free energy was reported to be - 32.65 ± 0.17 kcal/mol for the GRL0617-PLpro, for the NPC320891-PLpro complex, the TBE was - 35.58 ± 0.14 kcal/mol, for the NPC474594-PLpro, the TBE was - 43.72 ± 0.22 kcal/mol, while for NPC474595-PLpro complex, the TBE was calculated to be - 41.61 ± 0.20 kcal/mol, respectively. Clustering of the protein's motion and FEL further revealed that in NPC474594 and NPC474595 complexes, the drug was seen to have moved inside the binding cavity along with the loop in the palm region harboring the catalytic triad, thus justifying the higher binding of these two molecules particularly. In conclusion, the overall results reflect favorable binding of the identified hits strongly than the control drug, thus demanding in vitro and in vivo validation for clinical purposes.
    Matched MeSH terms: Molecular Dynamics Simulation
  17. Roney M, Huq AKMM, Issahaku AR, Soliman MES, Hossain MS, Mustafa AH, et al.
    J Biomol Struct Dyn, 2023;41(21):12186-12203.
    PMID: 36645141 DOI: 10.1080/07391102.2023.2166123
    Dengue fever is a significant public health concern throughout the world, causing an estimated 500,000 hospitalizations and 20,000 deaths each year, despite the lack of effective therapies. The DENV-2 RdRp has been identified as a potential target for the development of new and effective dengue therapies. This research's primary objective was to discover an anti-DENV inhibitor using in silico ligand- and structure-based approaches. To begin, a ligand-based pharmacophore model was developed, and 130 distinct natural products (NPs) were screened. Docking of the pharmacophore-matched compounds were performed to the active site of DENV-2 RdRp protease . Eleven compounds were identified as potential DENV-2 RdRp inhibitors based on docking energy and binding interactions. ADMET and drug-likeness were done to predict their pharmacologic, pharmacokinetic, and drug-likeproperties . Compounds ranked highest in terms of pharmacokinetics and drug-like appearances were then subjected to additional toxicity testing to determine the leading compound. Additionally, MD simulation of the lead compound was performed to confirm the docked complex's stability and the binding site determined by docking. As a result, the lead compound (compound-108) demonstrated an excellent match to the pharmacophore, a strong binding contact and affinity for the RdRp enzyme, favourable pharmacokinetics, and drug-like characteristics. In summary, the lead compound identified in this study could be a possible DENV-2 RdRp inhibitor that may be further studied on in vitro and in vivo models to develop as a drug candidate.Communicated by Ramaswamy H. Sarma.
    Matched MeSH terms: Molecular Dynamics Simulation
  18. Nguan H, Ahmadi S, Hashim R
    Phys Chem Chem Phys, 2014 Jan 7;16(1):324-34.
    PMID: 24257208 DOI: 10.1039/c3cp52385c
    Through atomistic molecular dynamic simulations using a GROMOS53a6 force field for the carbohydrate, we studied the lyotropic reverse hexagonal phase HII from a glycolipid, namely the Guerbet branched-chain β-d-glucoside, at 14% and 22% water concentrations. Our simulations showed that at low water concentration (14%) the sugar head group overlapped extensively and protruded into the water channel. In contrast, in the 22% concentration system a water column free from the sugar headgroup ('free' water) was formed as expected for the system close to the limit of maximum hydration. In both concentrations, we found anomalous water diffusion in the xy-plane, i.e. the two-dimensional space confined by the surface of the cylinder. On the other hand, along the z-axis, the water diffusion obeyed the Einstein relation for the 22% system, while for the 14% system it was slightly anomalous. For the 22% system, the diffusion along the z-axis of the 'free' water obeyed the Einstein relation, while that of the 'bound' water is slightly anomalous. The xy-plane displacement of the 'bound' water was higher than that for the 'free' water at times longer than 200 ps, as a consequence of the exchange of water molecules between the two regions. Based on our findings, we proposed an alternative explanation to the observed spatial heterogeneity in the HII phase from probe diffusion by Penaloza et al. (Phys. Chem. Chem. Phys., 2012, 14(15), 5247-5250). We found the extent of contact with water was different at different oxygen atoms within the sugar ring. Generally, a higher probability of hydrogen bonding but a shorter lifetime was found in 22% water compared to the case of 14% water. Finally, we examined the extension and compression of the alkyl chain of a columnar.
    Matched MeSH terms: Molecular Dynamics Simulation*
  19. Mohamad Ali MS, Mohd Fuzi SF, Ganasen M, Abdul Rahman RN, Basri M, Salleh AB
    Biomed Res Int, 2013;2013:925373.
    PMID: 23738333 DOI: 10.1155/2013/925373
    The psychrophilic enzyme is an interesting subject to study due to its special ability to adapt to extreme temperatures, unlike typical enzymes. Utilizing computer-aided software, the predicted structure and function of the enzyme lipase AMS8 (LipAMS8) (isolated from the psychrophilic Pseudomonas sp., obtained from the Antarctic soil) are studied. The enzyme shows significant sequence similarities with lipases from Pseudomonas sp. MIS38 and Serratia marcescens. These similarities aid in the prediction of the 3D molecular structure of the enzyme. In this study, 12 ns MD simulation is performed at different temperatures for structural flexibility and stability analysis. The results show that the enzyme is most stable at 0°C and 5°C. In terms of stability and flexibility, the catalytic domain (N-terminus) maintained its stability more than the noncatalytic domain (C-terminus), but the non-catalytic domain showed higher flexibility than the catalytic domain. The analysis of the structure and function of LipAMS8 provides new insights into the structural adaptation of this protein at low temperatures. The information obtained could be a useful tool for low temperature industrial applications and molecular engineering purposes, in the near future.
    Matched MeSH terms: Molecular Dynamics Simulation*
  20. Achari VM, Nguan HS, Heidelberg T, Bryce RA, Hashim R
    J Phys Chem B, 2012 Sep 27;116(38):11626-34.
    PMID: 22967067
    Glycolipids form materials of considerable potential for a wide range of surfactant and thin film applications. Understanding the effect of glycolipid covalent structure on the properties of their thermotropic and lyotropic assemblies is a key step toward rational design of new glycolipid-based materials. Here, we perform molecular dynamics simulations of anhydrous bilayers of dodecyl β-maltoside, dodecyl β-cellobioside, dodecyl β-isomaltoside, and a C(12)C(10) branched β-maltoside. Specifically, we examine the consequences of chain branching and headgroup identity on the structure and dynamics of the lamellar assemblies. Chain branching of the glycolipid leads to measurable differences in the dimensions and interactions of the lamellar assembly, as well as a more fluid-like hydrophobic chain region. Substitution of the maltosyl headgroup of βMal-C(12) by an isomaltosyl moiety leads to a significant decrease in bilayer spacing as well as a markedly altered pattern of inter-headgroup hydrogen bonding. The distinctive simulated structures of the two regioisomers provide insight into the difference of ~90 °C in their observed clearing temperatures. For all four simulated glycolipid systems, with the exception of the sn-2 chain of the branched maltoside, the alkyl chains are ordered and exhibit a distinct tilt, consistent with recent crystallographic analysis of a branched chain Guerbet glycoside. These insights into structure-property relationships from simulation provide an important molecular basis for future design of synthetic glycolipid materials.
    Matched MeSH terms: Molecular Dynamics Simulation*
Filters
Contact Us

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

External Links