- Open Access
Molecular dynamics study of enhanced Man5B enzymatic activity
Biotechnology for Biofuelsvolume 7, Article number: 83 (2014)
Biofuels are a well-known alternative to the largely used fossil-derived fuels, however the competition with food production is an ethical dilemma. Fortunately a solution is offered by second-generation biofuels which can be produced from agricultural waste or, more specifically, from plant cell wall polysaccharides. The conversion process involves typically enzymatic hydrolysis of lignocellulosic biomass and then separation of its constituent sugars that are further fermented to produce ethanol. Over the years several technologies have been developed that allow this conversion process to occur and the objective is now to make this process cost-competitive in today’s markets.
We observe that reduction of enzymatic efficiency in the presence of gluco-oligosaccharides is associated with a loss of the enzyme’s flexibility, the latter being required to bind new substrate, while the presence of manno-oligosaccharides does not pose this problem. Molecular dynamics simulations identify key contacts between substrates and the enzyme catalytic pocket that might be modified through site-directed mutagenesis to prevent loss of enzymatic efficiency.
Based on previous experimental studies and the new molecular dynamics data, we suggest that cellohexaose in the active site pocket slows down or even inhibits Man5B enzymatic activity. The assumption of such a mechanism is reasonable since when the gluco-oligosaccharide substrate is attached to the catalytic pocket it takes much longer to leave the pocket and thus prevents other substrates from reaching the active site. The insight is of crucial importance since the inhibition of enzymes by the enzymatic product or by an unsuitable substrate is a major technological problem in reducing the competitiveness of second-generation biofuel production.
Deconstruction of plant cell walls to fermentable sugar using enzymatic hydrolysis is being pursued for the production of so-called second-generation biofuels. Driven by significant research efforts worldwide, a large number of enzymes that may be used for biofuel production have been identified and biochemically characterized[1–3]. However, further data to elucidate cell wall deconstruction are needed for a fuller understanding of the enzyme reaction and development of enhanced conversion processes[3, 4].
Here we discuss the dynamics of Caldanaerobius polysaccharolyticus Man5B, an enzyme that cleaves both β-1,4 glucosidic and β-1,4 mannosidic linkages. Man5B and Man5A, two glycoside hydrolase family 5 (GH5) enzymes from the same bacterium, were shown to act synergistically and at high temperature on enzymatic conversion of plant cell wall polysaccharides to fermentable sugars[5, 6], a property that is highly desirable in the emerging biofuel industry[5, 7, 8].
A biochemical characterization of the two thermophilic β–mannanases was performed in an earlier report. The results provided insight into the physiological role of these enzymes in mannan degradation. Man5A is anchored to the cell surface of C. polysaccharolyticus through its surface layer homology (SLH) domain and generates oligosaccharides, which are then shuttled into the cytoplasm by the products of a gene cluster within which is also located the gene encoding Man5B. Man5B, a cytoplasmic enzyme, has been shown to cleave the transported oligosaccharides into mono- and disaccharides for subsequent metabolism. In reports on the enzymatic activities of Man5A and Man5B it was demonstrated that Man5B and Man5A show highly specific activities with glucomannan as a substrate. Interestingly, however, in addition to cleaving β-1,4 mannosidic linkages the two enzymes also cleaved β-1,4 glucosidic bonds. It has been reported that GH5 mannanases with known three-dimensional structures act specifically on glucose or mannose, however, due to their absolute specificity for mannose at the -1 sub-site they cleave only mannosidic bonds, as also observed for other mannanases. Therefore, the wider capacity of the C. polysaccharolyticus GH5 enzymes to cleave β-1,4 mannosidic and β-1,4 glucosidic linkages is of great importance. Since the mechanisms underlying the two different enzymatic activities in the two enzymes are unknown, in the present study we subjected Man5B to molecular dynamics simulations to unravel how the substrates dock to the catalytic site of the enzyme and how enzyme dynamics are affected by both mannohexaose and cellohexaose.
Results and discussion
Understanding the dynamics of glycoside hydrolases is key for the development of cost-competitive second-generation biofuels. In our study we docked cellohexaose and mannohexaose as substrates to the Man5B catalytic site and carried out molecular dynamics simulations employing the program NAMD[11, 12] to elucidate the mechanism by which Man5B, a thermophilic enzyme, hydrolyzes cello-oligosaccharide and manno-oligosaccharide substrates (the latter more efficiently). For the docking (see Figure 1A) we used the software VMD; for a template we used GH5 structures with mono- and disaccharides presented as substrates taken from the protein data bank [PDB:1CEN and PDB:3AMG] and reported in[14, 15]. After docking and subsequent equilibration stable positions for the ligands were established, three similar conformations for cellohexaose and three similar conformations for mannohexaose were obtained as shown in Additional file1. As shown in Figure 1B and C, a perfect pocket that accommodates the C6 group of the sugar in position -1 was recognized together with a pocket for the sugar in position 1. As can be seen in Figure 1C, the hexasaccharide substrates with their sugar chains were found to be slightly twisted into a conformation that could favor hydrolysis.
The docking procedure described resulted in placements of the hexasaccharide substrates in which only the disaccharide segments, corresponding to the region of the mono- and disaccharide substrates of structures [PDB:1CEN and PDB:3AMG], are buried in the enzyme’s catalytic pocket. One may wonder in how far the hexasaccharide placement is an artifact of the docking and equilibration procedure adopted and if modeling should rather await the availability of crystal structures for hexasaccharide substrates. Indeed, modeled substrate and/or enzyme complexes can be erroneous. However, three general arguments can be made in favor of employing the docked structures described over employing at a later stage not yet available crystal structures. First, the relaxation time of almost 20 ns (including time used for the equilibration period of the molecular dynamics protocol) is likely long enough for sufficient local relaxation of the initial geometry adopted. Second, multimeric substrates adopt natural, rather disordered, geometries that are well accounted for in molecular dynamics (MD) simulations. Third, crystallographic structures are representative of very closely packed systems characterized by limited solvent accessibility as well as strong crystal contacts and, as a result, are often not representative of enzymes in their functional states.
To analyze how the substrates can interfere with the dynamics of the enzyme we performed a root mean square deviation (RMSD) analysis of all our MD trajectories. As shown in Figure 2A for one of the trajectories for each system, protein RMSD values show that both cellohexaose and mannohexaose substrates reduce the flexibility of Man5B. A similar observation had been reported before for several systems, including both glycoside hydrolases and other enzymes[16, 17]. The ligand RMSD values depicted in Figure 2B show that both cellohexaose and mannohexaose are very flexible in the catalytic site. In the case of the oligosaccharide chain of six sugars, two of the sugars on one side of the chain are always unbound and move, accordingly, freely in the solvent. However, the RMSD analysis shows that on the other end of the chain, close to the catalytic residues (position 1 and -1, indicated in Figure 1C), the substrates are in close contact with the enzyme and, as a result, exhibit only small amplitude thermal motion (data not shown).
A characterization of the protein-substrate interaction was achieved by a map of contacts between protein and substrate throughout the molecular dynamics simulations. For this purpose we tracked the strongest interactions, namely the hydrogen bonds. In Table 1 we present a list of hydrogen bonds and their average prevalence calculated from the three sample simulations carried out for each of the two substrates, mannohexaose and cellohexaose. This analysis shows that both substrates are interacting mostly with the same amino acid residues and that some of the interactions are more prevalent for mannohexaose. As regards the average number of hydrogen bond interactions for each substrate, presented in Figure 3A, molecular dynamics simulations show that the average number of interactions between mannohexaose and Man5B becomes more stable over time when compared to cellohexaose, but is not necessarily smaller than the interactions between cellohexaose and Man5B. The fact that hydrogen bonds are more stable for mannohexaose as a substrate might be an indication of the better enzymatic efficiency for cleavage of this substrate as occurrence of the hydrolysis reaction requires a stable conformation.Analysis of the hydrogen bonds can also point to the amino acids that contribute most to the contact between enzyme and substrate and, therefore, are candidates for mutation seeking altered enzymatic activity. Figure 3B and C shows that the most prevalent interactions are located around the catalytic center. However, mutations on the respective amino acids may not only affect binding strength, but may at the same time interfere unfavorably with optimal enzymatic activity since the overall reaction process depends on the position of the amino acids around the catalytic amino acids. In general, the mechanism of glycoside hydrolases (GH) involves a nucleophilic attack by a deprotonated carboxylic acid nucleophile (GLU258 side chain) on the anomeric carbon, displacing the attached sugar residue and forming a covalent enzyme-sugar adduct. Subsequently, an activated water molecule displaces the enzymic carboxylic acid (GLU137 side chain) resulting in the net retention of stereochemical configuration at the anomeric carbon.
The aforementioned results on the most prevalent hydrogen bonds characterize the interaction of both substrates with the enzyme. However, the small difference observed between Man5B acting on cellohexaose or on mannohexaose does not permit a final conclusion on why the enzyme is more efficient at cleaving the linkages of mannohexaose rather than those of cellohexaose. To shed light on how the two substrates affect the enzyme dynamics differently, we employed principal component analysis (PCA), which permits one to identify large movements of the protein. PCA can describe the conformational flexibility of enzymes commonly associated with enzyme activity. In fact, a great variety of internal motions ranging over a wide range of time scales are involved in the catalytic process; such motions are difficult to capture in straightforward molecular dynamics simulations, but can be captured in a PCA analysis[19, 20]. Figure 4 shows for Man5B how each residue participates in the one PCA normal mode judged by the present authors to be the most essential one for the enzyme dynamics. The mode selected is the one with largest amplitude in the case of both the control simulation (with no substrate) and the sample simulations with mannohexaose as substrate. In the simulations with cellohexaose, the PCA mode representing the same movement was found to be only the mode with the third largest amplitude for two of the sample simulations and the mode with the fourth largest amplitude for the third sample simulation. Figure 4A illustrates that the selected mode corresponds to the opening and closing of a cleft suitable for binding and releasing the substrates as well as in bringing the substrates into contact with the catalytic center. As shown in Figure 4B and in detail for the cleft region in Figure 4C, when substrate is bound a small reduction in the normal mode amplitude is observed for mannohexaose but a very large reduction is seen for cellohexaose.
The behavior revealed by the PCA mode shown in Figure 4 suggests a stronger interaction of Man5B with cellohexaose than with mannohexaose, and one might want to conclude that the enzyme works better on cellohexaose. However, the experimental data tell the opposite. We interpret our PCA results therefore as strongly suggesting that the opening and closing of the clefts around the enzyme’s catalytic site is a much slower process when cellohexaose is bound, and that once the latter substrate is bound opening of the clefts might even not occur. Such an opening and closing movement is likely to be related to the catalysis, allowing the entrance of a new substrate and the release of the products. Based on our MD simulations we suggest then that cellohexaose actually inhibits the enzymatic activity, since upon entering and strongly binding to the catalytic pocket the cleaved substrate takes too long to leave and thereby prevents fresh substrate molecules from reaching the reaction site. Inefficient release of the reaction products, or possibly of a different substrate bound, can lead to enzyme inhibition and can thereby reduce the efficiency of biomass conversion[8, 21, 22].
The search for new and more efficient glycoside hydrolases (GH) has intensified over the last few years due to a need for such enzymes in the biofuel industry. Man5B, a cytoplasmic enzyme of the glycoside hydrolase family 5, has been shown to cleave manno- and gluco-oligosaccharides into mono- and disaccharides for subsequent metabolism. Experimental assays show that Man5B acts more efficiently on manno-oligosaccharides than on gluco-oligosaccharides, however the mechanism for this behavior was not clear. We have performed a molecular dynamics study to identify this mechanism and to elucidate which amino acid residues are controlling enzymatic efficiency. The insights gained from these studies are critical to the development of more efficient enzymes through rational targeting of residues for site-directed mutagenesis.
The molecular dynamics simulations yielded surprising results in that they showed Man5B to bind cellohexaose nearly as tightly as mannohexaose, as shown in Figure 3A. The RMSD shows that the protein is stabilized by both ligands, however, mannohexaose shows to be slightly more flexible in the catalytic pocket as shown in Figure 2B. The same behavior is not observed when one analyzes the RMSD for only the carbohydrates at position -1 and 1 (data not shown). At these positions mannohexaose looks slightly more stable than cellohexaose, as is also evident from Table 1, where hydrogen bond prevalence suggests that mannohexaose hydrogen bonds are more prevalent and, therefore, that the conformation of this substrate close to the catalytic amino acids is more stable. Such behavior, together with a PCA of the trajectories, shown in Figure 4, suggests a stronger binding of Man5B to cellohexaose than to mannohexaose, except for the binding to the amino acids in the catalytic center. The analysis may explain why Man5B is less efficient in cleaving cellohexaose than mannohexaose: the PCA mode corresponding to the opening and closing of a molecular cleft containing the catalytic pockets of Man5B has only a small amplitude in the case of the Man5B-cellohexaose complex, while the amplitude is large enough for both for the Man5B-mannohexaose complex and for Man5B without a substrate. When cellohexaose is present, the much smaller PCA amplitude for a mode representing the same movement indicates that this movement might be much slower, or even completely inhibited. The result together with the experimental findings on enzyme activity implies that cellohexaose prevents the opening of the enzyme to release the reaction product. It is also possible that cellohexaose slightly misaligned in the enzyme cannot be released and rebound, thereby inhibiting the enzyme. The experimental data also indicates that the enzyme is more efficient to degrade mannohexaose than mannotetraose which, according to our RMSD analysis, can be related to the large flexibility of the substrate sugar present at position 5 and 6. This larger flexibility is likely to be involved in the faster mechanism of the opening and closing of the clefts. This suggestion is supported by the observed stronger interaction between the cellohexaose and amino acids outside of the catalytic center.
The inhibition of enzymes by the enzymatic product or by an unsuitable substrate is known to be a key problem in biofuel production[1, 22–25] and the characterization of the mechanism of inhibition is of fundamental importance for the second-generation biofuel industry. Now that an explanation for Man5B inhibition by cellohexaose is suggested, further simulations can assist bioengineers in altering Man5B by mutating amino acids contributing to the loss of flexibility in the presence of gluco-oligosaccharides. In this regard obvious candidates for mutations are the residues in the region of the flaps around the catalytic pocket, namely amino acid residues 90-94 and 208-212, in particular side groups ASN92 and TRP210, that are forming the tunnel where the substrates lies. Smaller side chains in the positions 92 and 210 could inhibit the formation of the catalytic tunnel and enhance the enzymatic activity for gluco-oligosaccharides.However, mutations on the amino acid residues of the cleft region may not only affect binding strength, but may at the same time interfere with optimal enzymatic activity. Some of these residues are important for the stability of the substrate inside the catalytic pocket; in particular side group TRP210, shown in Figure 3B and C, appears in the simulation frequently with its rings in a parallel contact with the substrate’s carbohydrates rings. A mutation, for example destabilizing substrate conformations critical for the catalyzed hydrolysis reaction, could then decrease the enzymatic efficiency of Man5B. The results provided by the present study suggest that the activity of Man5B, and likely that of other glycoside hydrolases, is much more complex than expected for the reaction step alone, involving a complete set of large-scale motions of the enzyme that are much more rate limiting than the reaction itself. Based on the present study and the previous experimental studies, an extensive screening study employing site-directed mutagenesis of the aforementioned side groups needs to be performed to check how the activity of the enzyme is affected by mutations in the cleft region.
The structure of Man5B has been solved by means of X-ray crystallography at 1.60 Å resolution and is available at the Protein Data Bank [PDB:3W0K]. Three different systems were constructed for the molecular dynamics study of the enzyme’s activity: a control system with no substrate, a system with cellohexaose docked to the Man5B catalytic site, and a system with mannohexaose docked. Both substrates were placed in the catalytic site using tools available in the VMD software as well as using information from other structures of the glycoside hydrolase family 5 [PDB:1CEN and PDB:3AMG][14, 15]. The enzymes were crystalized with mono- [chain A of PDB:3AMG] and disaccharides [PDB:1CEN and chain B of PDB:3AMG] in their catalytic pocket. The positions of the mono- and disaccharide substrates in the templates were employed then as a guide to fit the hexasaccharide substrates to Man5B. Three structures were generated, one using [PDB:1CEN] as template and two using the different chains in [PDB:3AMG] as templates. A stable conformation of the hexasaccharide substrates was determined by means of NAMD’s energy minimization protocol, where positions of atoms of the hexasaccharide substrate that were present in the various templates were first restrained to the respective atomic position of the mono- and disaccharide substrates in our templates. The three slightly different conformations were subjected to an equilibration protocol through which a stable conformation of the hexasaccharide substrates was determined. In this protocol the position of the atoms of the substrate that were present in the template were first restrained to the position found in the template as just pointed out. The positions of the atoms of the protein backbone were also restrained and a short (500 ps) MD simulation was performed. Subsequently, equilibration simulations without constraint to any atoms were carried out as detailed below. After the equilibration simulations the substrates for the systems with slightly different initial conformations assumed actually similar conformations as shown in Additional file1. All seven systems, the control system without substrate and three systems each for the enzyme with mannohexaose and cellohexaose as substrate, were then solvated and the net charge of the protein (the substrates have no net charge) was neutralized using three sodium atoms as counter-ions, which were randomly arranged in the solvent.
The simulations in this study were performed employing the NAMD molecular dynamics package[11, 12]. The CHARMM36 force field[27, 28] along with the TIP3 water model was used to describe all systems. The simulations were done assuming periodic boundary conditions in the NpT ensemble with temperature maintained at 65°C (338 K) using Langevin dynamics for pressure, kept at 1 bar, and temperature coupling. A distance cutoff of 11.0 Å was applied to short range non-bonded interactions, whereas long range electrostatic interactions were calculated using the particle-mesh Ewald (PME) method. The equations of motion were integrated using the r-RESPA multiple-time-step scheme to update the van der Waals interactions every two steps and electrostatic interactions every four steps. The time-step of integration was chosen to be 2 fs for all simulations performed in this study. For the control simulation, namely the simulation without substrate, the first 2 ns of the simulations served to equilibrate the system ramping the temperature from 0 K to 338 K. During the first half of the equilibration the position of the atoms of the backbone were restrained. The equilibration protocol for the systems where substrates were present was slightly different: during the first 3 ns the atoms of the backbone of the protein were constrained, followed by 15 ns where atoms of the backbone that were located up to 5 Å of the substrate were constrained. A total of 100 ns of molecular dynamics were performed for each system.
Analyses of MD trajectories were carried out employing VMD and its plugins. We determined the RMSD for both the ligand and the protein. Hydrogen bonds were assigned based on two geometric criteria for every trajectory frame saved: first, distances between acceptor and hydrogen should be less than 3.0 Å; second, the angle between hydrogen-donor-acceptor should be smaller than 30 degrees. Using the ProDy plugin PCA was performed to detect correlation in atomic motion over a molecular dynamics trajectory. PCA is useful in analyzing the slow motion of flexible regions of an enzyme.
Glycoside hydrolase family 5
Principal component analysis
Protein data bank
Root mean square deviation
Surface layer homology.
Himmel ME, Ding S-Y, Johnson DK, Adney WS, Nimlos MR, Brady JW, Foust TD: Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science 2007, 315: 804-807.
Horn SJ, Vaaje-Kolstad G, Westereng B, Eijsink VG: Novel enzymes for the degradation of cellulose. Biotechnol Biofuels 2012, 5: 45.
Beckham GT, Bomble YJ, Bayer EA, Himmel ME, Crowley MF: Applications of computational science for understanding enzymatic deconstruction of cellulose. Curr Opin Biotechnol 2011, 22: 231-238.
Ding S-Y, Liu Y-S, Zeng Y, Himmel ME, Baker JO, Bayer EA: How does plant cell wall nanoscale architecture correlate with enzymatic digestibility? Science 2012, 338: 1055-1060.
Han Y, Dodd D, Hespen CW, Ohene-Adjei S, Schroeder CM, Mackie RI, Cann IKO: Comparative analyses of two thermophilic enzymes exhibiting both beta-1,4 mannosidic and beta-1,4 glucosidic cleavage activities from Caldanaerobius polysaccharolyticus. J Bacteriol 2010, 192: 4111-4121.
Su X, Agarwal V, Dodd D, Bae B, Mackie RI, Nair SK, Cann IKO: Mutational insights into the roles of amino acid residues in ligand binding for two closely related family 16 carbohydrate binding modules. J Biol Chem 2010, 285: 34665-34676.
Batista PR, Costa MGDS, Pascutti PG, Bisch PM, de Souza W: High temperatures enhance cooperative motions between CBM and catalytic domains of a thermostable cellulase: mechanism insights from essential dynamics. Phys Chem Chem Phys 2011, 13: 13709-13720.
Yennamalli RM, Rader AJ, Kenny AJ, Wolt JD, Sen TZ: Endoglucanases: insights into thermostability for biofuel applications. Biotechnol Biofuels 2013, 6: 136.
Cann IKO, Kocherginskaya S, Michael R, White BA, Mackie RI: Molecular cloning, sequencing, and expression of a novel multidomain mannanase gene from thermoanaerobacterium polysaccharolyticum molecular cloning, sequencing, and expression of a novel multidomain mannanase gene from thermoanaerobacterium polysaccharolyticum. J Bacteriol 1999, 181: 1643-1651.
Tailford LE, Ducros VM, Flint JE, Roberts SM, Morland C, Zechel DL, Smith N, Bjørnvad ME, Borchert TV, Wilson KS, Davies GJ, Gilbert HJ: Understanding how diverse beta-mannanases recognize heterogeneous substrates. Biochemistry 2009, 48: 7009-7018.
Kale L, Skeel R, Bhandarkar M, Brunner R, Gursoy A, Krawetz N, Phillips J, Shinozaki A, Varadarajan K, Schulten K: NAMD2: greater scalability for parallel molecular dynamics. J Comput Phys 1998, 151: 283-312.
Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kalé L, Schulten K: Scalable molecular dynamics with NAMD. J Comput Chem 2005, 26: 1781-1802.
Humphrey W: VMD: visual molecular dynamics. J Mol Graph 1996, 14: 33-38.
Domínguez R, Souchon H, Lascombe M, Alzari PM: The crystal structure of a family 5 endoglucanase mutant in complexed and uncomplexed forms reveals an induced fit activation mechanism. J Mol Biol 1996, 257: 1042-1051.
Wu T-H, Huang C-H, Ko T-P, Lai H-L, Ma Y, Chen C-C, Cheng Y-S, Liu J-R, Guo R-T: Diverse substrate recognition mechanism revealed by thermotoga maritima Cel5A structures in complex with cellotetraose, cellobiose and mannotriose. Biochim Biophys Acta 1814, 2011: 1832-1840.
Yang L-Q, Sang P, Tao Y, Fu Y-X, Zhang K-Q, Xie Y-H, Liu S-Q: Protein dynamics and motions in relation to their functions: several case studies and the underlying mechanisms. J Biomol Struct Dyn 2013, 32: 37-41.
Batista PR, Pandey G, Pascutti PG, Bisch PM, Perahia D, Robert CH: Free energy profiles along consensus normal modes provide insight. J Chem Theor Comput 2011, 7: 2348-2352.
Yon JM, Perahia D, Ghélis C: Conformational dynamics and enzyme activity. Biochimie 1998, 80: 33-42.
Balsera MA, Wriggers W, Oono Y, Schulten K: Principal component analysis and long time protein dynamics. J Phys Chem 1996, 100: 2567-2572.
Isin B, Schulten K, Tajkhorshid E, Bahar I: Mechanism of signal propagation upon retinal isomerization: insights from molecular dynamics simulations of rhodopsin restrained by normal modes. Biophys J 2008, 95: 789-803.
Teugjas H, Väljamäe P: Product inhibition of cellulases studied with 14C-labeled cellulose substrates. Biotechnol Biofuels 2013, 6: 104.
Xiao Z, Zhang X, Gregg DJ, Saddler JN: Effects of sugar inhibition on cellulases and beta-glucosidase during enzymatic hydrolysis of softwood substrates. Appl Biochem Biotechnol 2004, 115: 1115-1126.
Jönsson LJ, Alriksson B, Nilvebrant N-O: Bioconversion of lignocellulose: inhibitors and detoxification. Biotechnol Biofuels 2013, 6: 16.
Duan C-J, Feng J-X: Mining metagenomes for novel cellulase genes. Biotechnol Lett 2010, 32: 1765-1775.
Dunlop MJ: Engineering microbes for tolerance to next-generation biofuels. Biotechnol Biofuels 2011, 4: 32.
Oyama T, Schmitz GE, Dodd D, Han Y, Burnett A, Nagasawa N, Mackie RI, Nakamura H, Morikawa K, Cann I: Mutational and structural analyses of caldanaerobius polysaccharolyticus Man5b reveal novel active site residues for family 5 glycoside hydrolases. PLoS One 2013, 8: e80448.
Best RB, Zhu X, Shim J, Lopes PEM, Mittal J, Feig M, Mackerell AD: Optimization of the additive CHARMM all-atom protein force field targeting Improved sampling of the backbone, ψ and side-chain χ 1 and χ 2 dihedral angles. J Chem Theor Comput 2012, 8: 3257-3273.
MacKerell AD, Bashford D, Dunbrack RL, Evanseck JD, Field MJ, Fischer S, Gao J, Guo H, Ha S, Joseph-McCarthy D, Kuchnir L, Kuczera K, Lau FTK, Mattos C, Michnick S, Ngo T, Nguyen DT, Prodhom B, Reiher WE, Roux B, Schlenkrich M, Smith JC, Stote R, Straub J, Watanabe M, Wiórkiewicz-Kuczera J, Yin D, Karplus M: All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B 1998, 102: 3586-3616.
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of simple potential functions for simulating liquid water. J Chem Phys 1983, 79: 926-935.
Darden T, York D, Pedersen L: Particle mesh Ewald: an N log (N) method for Ewald sums in large systems. J Chem Phys 1993, 98: 10089-10092.
Bakan A, Meireles LM, Bahar I: ProDy: protein dynamics inferred from theory and experiments. Bioinformatics 2011, 27: 1575-1577.
This work was supported by grants from the National Institutes of Health (NIH, 9P41GM104601 to KS) and the National Science Foundation (NSF, MCB-1157615 to KS). Simulations made use of the Texas Advanced Computing Center (TACC) as part of the Extreme Science and Engineering Discovery Environment (XSEDE, MCA93S028 to KS and MCB130089 and MCB130115 to IC).
The authors declare that they have no competing interests.
RCB carried out the molecular dynamics studies, participated in its design, performed the analyses, and drafted the manuscript. IC conceived the study, participated in its design and analyses, and helped to draft the manuscript. KS participated in the design and coordination of the study, in the provision of simulation techniques, and in the drafting of the manuscript. All authors read and approved the final manuscript.