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.
The stability of biocatalysts is an important criterion for a sustainable industrial operation economically. T1 lipase is a thermoalkalophilic enzyme derived from Geobacillus zalihae strain T1 (T1 lipase) that was isolated from palm oil mill effluent (POME) in Malaysia. We report here the results of high temperatures molecular dynamics (MD) simulations of T1 lipase in explicit solvent. We found that the N-terminal moiety of this enzyme was accompanied by a large flexibility and dynamics during temperature-induced unfolding simulations which preceded and followed by clear structural changes in two specific regions; the small domain (consisting of helices alpha3 and alpha5, strands beta1 and beta2, and connecting loops) and the main catalytic domain or core domain (consisting of helices alpha6- alpha9 and connecting loops which located above the active site) of the enzyme. The results suggest that the small domain of model enzyme is a critical region to the thermostability of this organism.
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.
An all-atom level MD simulation in explicit solvent at high temperature is a powerful technique to increase our knowledge about the structurally important regions modulating thermal stability in thermenzymes. In this respect, two large-sized thermoalkalophilic enzymes from Bacillus stearothermophilus L1 (L1 lipase) and Geobacillus zalihae strain T1 (T1 lipase) are well-established representatives. In this paper, comparative results from temperature-induced MD simulations of both model systems at 300 K, 400 K and 500 K are presented and discussed with respect to identification of highly flexible regions critical to thermostability. From our MD simulation results, specific regions along the L1 lipase and T1 lipase polypeptide chain including the small domain and the main catalytic domain or core domain of both enzymes show a marked increase in fluctuations and dynamics followed by clear structural changes. Overall, the N-terminal moiety of both enzymes and their small domains exhibit hyper-sensitivity to thermal stress. The results appear to propose that these regions are critical in determining of the overall thermal stability of both organisms.
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.
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.
Previously, the hypothetical protein, KPN00728 from Klebsiella pneumoniae MGH78578 was the Succinate dehydrogenase (SDH) chain C subunit via structural prediction and molecular docking simulation studies. However, due to limitation in docking simulation, an in-depth understanding of how SDH interaction occurs across the transmembrane of mitochondria could not be provided.
Klebsiella pneumoniae causes neonatal sepsis and nosocomial infections. One of the strains, K. pneumoniae MGH 78578, shows high level of resistance to multiple microbial agents. In this study, domain family, amino acid sequence and topology analyses were performed on one of its hypothetical protein, YggG (KPN_03358). Structural bioinformatics approaches were used to predict the structure and functionality of YggG protein. The open reading frame (ORF) of yggG, which was a putative metalloprotease gene, was also cloned, expressed and characterized. The ORF was PCR amplified from K. pneumoniae MGH 78578 genomic DNA and cloned into a pET14-b vector for heterologous expression in Escherichia coli. The purified YggG protein was subsequently assayed for casein hydrolysis under different conditions. This protein was classified as peptidase M48 family and subclan gluzincin. It was predicted to contain one transmembrane domain by TMpred. Optimal protein expression was achieved by induction with 0.6 mM isopropyl thiogalactoside (IPTG) at 25 °C for six hours. YggG was purified as soluble protein and confirmed to be proteolytically active under the presence of 1.25 mM zinc acetate and showed optimum activity at 37 °C and pH 7.4. We confirmed for the first time that the yggG gene product is a zinc-dependent metalloprotease.
The activation of lipases has been postulated to proceed by interfacial activation, temperature switch activation, or aqueous activation. Recently, based on molecular dynamics (MD) simulation experiments, the T1 lipase activation mechanism was proposed to involve aqueous activation in addition to a double-flap mechanism. Because the open conformation structure is still unavailable, it is difficult to validate the proposed theory unambiguously to understand the behavior of the enzyme. In this study, we try to validate the previous reports and uncover the mystery behind the activation process using structural analysis and MD simulations. To investigate the effects of temperature and environmental conditions on the activation process, MD simulations in different solvent environments (water and water-octane interface) and temperatures (20, 50, 70, 80, and 100°C) were performed. Based on the structural analysis of the lipases in the same family of T1 lipase (I.5 lipase family), we proposed that the lid domain comprises α6 and α7 helices connected by a loop, thus forming a helix-loop-helix motif involved in interfacial activation. Throughout the MD simulations experiments, lid displacements were only observed in the water-octane interface, not in the aqueous environment with respect to the temperature effect, suggesting that the activation process is governed by interfacial activation coupled with temperature switch activation. Examining the activation process in detail revealed that the large structural rearrangement of the lid domain was caused by the interaction between the hydrophobic residues of the lid with octane, a nonpolar solvent, and this conformation was found to be thermodynamically favorable.
5-HT(1A) serotonin and D1 dopamine receptor agonists have been postulated to be able to improve negative and cognitive impairment symptoms of schizophrenia, while partial agonists and antagonists of the D2 and 5-HT(2A) receptors have been reported to be effective in reducing positive symptoms. There is therefore a need for well-defined homology models for the design of more selective antipsychotic agents, since no three-dimensional (3D) crystal structures of these receptors are currently available. In this study, homology models were built based on the high-resolution crystal structure of the β(2)-adrenergic receptor (2RH1) and further refined via molecular dynamics simulations in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer system with the GROMOS96 53A6 united atom force field. Docking evaluations with representative agonists and antagonists using AutoDock 4.2 revealed binding modes in agreement with experimentally determined site-directed mutagenesis data and significant correlations between the computed and experimental pK (i) values. The models are also able to distinguish between antipsychotic agents with different selectivities and binding affinities for the four receptors, as well as to differentiate active compounds from decoys. Hence, these human 5-HT(1A), 5-HT(2A), D1 and D2 receptor homology models are capable of predicting the activities of novel ligands, and can be used as 3D templates for antipsychotic drug design and discovery.
The structure of psychrophilic chitinase (CHI II) from Glaciozyma antarctica PI12 has yet to be studied in detail. Due to its low sequence identity (<30 %), the structural prediction of CHI II is a challenge. A 3D model of CHI II was built by first using a threading approach to search for a suitable template and to generate an optimum target-template alignment, followed by model building using MODELLER9v7. Analysis of the catalytic insertion domain structure in CHI II revealed an increase in the number of aromatic residues and longer loops compared to mesophilic and thermophilic chitinases. A molecular dynamics simulation was used to examine the stability of the CHI II structure at 273, 288 and 300 K. Structural analysis of the substrate-binding cleft revealed a few exposed aromatic residues. Substitutions of certain amino acids in the surface and loop regions of CHI II conferred an increased flexibility to the enzyme, allowing for an adaptation to cold temperatures. A substrate binding comparison of CHI II with the mesophilic chitinase from Coccidioides immitis, 1D2K, suggested that the psychrophilic adaptation and catalytic activity at low temperatures were achieved through a reduction in the number of salt bridges, fewer hydrogen bonds and an increase in the exposure of the hydrophobic side chains to the solvent.
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.
Exotic functions of antifreeze proteins (AFP) and antifreeze glycopeptides (AFGP) have recently been attracted with much interest to develop them as commercial products. AFPs and AFGPs inhibit ice crystal growth by lowering the water freezing point without changing the water melting point. Our group isolated the Antarctic yeast Glaciozyma antarctica that expresses antifreeze protein to assist it in its survival mechanism at sub-zero temperatures. The protein is unique and novel, indicated by its low sequence homology compared to those of other AFPs. We explore the structure-function relationship of G. antarctica AFP using various approaches ranging from protein structure prediction, peptide design and antifreeze activity assays, nuclear magnetic resonance (NMR) studies and molecular dynamics simulation. The predicted secondary structure of G. antarctica AFP shows several α-helices, assumed to be responsible for its antifreeze activity. We designed several peptide fragments derived from the amino acid sequences of α-helical regions of the parent AFP and they also showed substantial antifreeze activities, below that of the original AFP. The relationship between peptide structure and activity was explored by NMR spectroscopy and molecular dynamics simulation. NMR results show that the antifreeze activity of the peptides correlates with their helicity and geometrical straightforwardness. Furthermore, molecular dynamics simulation also suggests that the activity of the designed peptides can be explained in terms of the structural rigidity/flexibility, i.e., the most active peptide demonstrates higher structural stability, lower flexibility than that of the other peptides with lower activities, and of lower rigidity. This report represents the first detailed report of downsizing a yeast AFP into its peptide fragments with measurable antifreeze activities.
Currently, there is no three-dimensional structure of D-specific dehalogenase (DehD) in the protein database. We modeled DehD using ab initio technique, performed molecular dynamics (MD) simulation and docking of D-2-chloropropionate (D-2CP), D-2-bromopropionate (D-2BP), monochloroacetate (MCA), monobromoacetate (MBA), 2,2-dichloropropionate (2,2-DCP), d,l-2,3-dichloropropionate (d,l-2,3-DCP), and 3-chloropropionate (3-CP) into the DehD active site. The sequences of DehD and D-2-haloacid dehalogenase (HadD) from Pseudomonas putida AJ1 have 15% sequence similarity. The model had 80% of the amino acid residues in the most favored region when compared to the crystal structure of DehI from Pseudomonas putida PP3. Docking analysis revealed that Arg107, Arg134 and Tyr135 interacted with D-2CP, and Glu20 activated the water molecule for hydrolytic dehalogenation. Single residue substitutions at 25-30 °C showed that polar residues of DehD were stable when substituted with nonpolar residues and showed a decrease in activity within the same temperature range. The molecular dynamics simulation of DehD and its variants showed that in R134A variant, Arg107 interacted with D-2CP, while in Y135A, Gln221 and Arg231 interacted with D-2CP. It is our emphatic belief that the new model will be useful for the rational design of DehDs with enhanced potentials.
Palm oil-based esters (POEs) are unsaturated and non-ionic esters with a great potential to act as chemical penetration enhancers and drug carriers for transdermal drug nano-delivery. A ratio of palmitate ester and nonionic Tween80 with and without diclofenac acid was chosen from an experimentally determined phase diagram. Molecular dynamics simulations were performed for selected compositions over a period of 15 ns. Both micelles showed a prolate-like shape, while adding the drug produced a more compact micellar structure. Our results proposed that the drug could behave as a co-surfactant in our simulated model.
Tuberculosis (TB), an epidemic disease, affects the world with death rate of two million people every year. The bacterium Mycobacterium tuberculosis was found to be a more potent and disease-prolonged bacterium among the world due to multi-drug resistance. Emergence of new drug targets is needed to overcome the bacterial resistance that leads to control epidemic tuberculosis. The pathway thiamine biosynthesis was targeting M. tuberculosis due to its role in intracellular growth of the bacterium. The screening of enzymes involved in thiamin biosynthesis showed novel target thiazole synthase (ThiG) involved in catalysis of rearrangement of 1-deoxy-D-xylulose 5-phosphate (DXP) to produce the thiazole phosphate moiety of thiamine. We carried out homology modeling for ThiG to understand the structure-function relationship, and the model was refined with MD simulations. The results showed that the model predicted with (α + β)8-fold of synthase family proteins. Molecular docking of ThiG model with substrate DXP showed binding mode and key residues ARG46, ASN69, THR41, and LYS96 involved in the catalysis. First-line anti-tuberculosis drugs were docked with ThiG to identify the inhibition. The report showed the anti-tuberculosis drugs interact well with ThiG which may lead to block thiamin biosynthesis pathway.
Molecular dynamics simulation was used to study the dynamic differences between native Aspergillus niger PhyA phytase and a mutant with 20 % greater thermostability. Atomic root mean square deviation, radius of gyration, and number of hydrogen bonds and salt bridges are examined to determine thermostability factors. The results suggest that, among secondary structure elements, loops have the most impact on the thermal stability of A. niger phytase. In addition, the location rather than the number of hydrogen bonds is found to have an important contribution to thermostability. The results also show that salt bridges may have stabilizing or destabilizing effect on the enzyme and influence its thermostability accordingly.
Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.
We describe a server that allows the interrogation of the Protein Data Bank for hypothetical 3D side chain patterns that are not limited to known patterns from existing 3D structures. A minimal side chain description allows a variety of side chain orientations to exist within the pattern, and generic side chain types such as acid, base and hydroxyl-containing can be additionally deployed in the search query. Moreover, only a subset of distances between the side chains need be specified. We illustrate these capabilities in case studies involving arginine stacks, serine-acid group arrangements and multiple catalytic triad-like configurations. The IMAAAGINE server can be accessed at http://mfrlab.org/grafss/imaaagine/.
Immobilized lipases were applied to the enzymatic conversion of oils from spent coffee ground into biodiesel. Two lipases were selected for the study because of their conformational behavior analysed by Molecular Dynamics (MD) simulations taking into account that immobilization conditions affect conformational behavior of the lipases and ultimately, their efficiency upon immobilization. The enzymatic synthesis of biodiesel was initially carried out on a model substrate (triolein) in order to select the most promising immobilized biocatalysts. The results indicate that oils can be converted quantitatively within hours. The role of the nature of the immobilization support emerged as a key factor affecting reaction rate, most probably because of partition and mass transfer barriers occurring with hydrophilic solid supports. Finally, oil from spent coffee ground was transformed into biodiesel with yields ranging from 55% to 72%. The synthesis is of particular interest in the perspective of developing sustainable processes for the production of bio-fuels from food wastes and renewable materials. The enzymatic synthesis of biodiesel is carried out under mild conditions, with stoichiometric amounts of substrates (oil and methanol) and the removal of free fatty acids is not required.