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.
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.
DNA polymerase β is a 39kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31kDa domain responsible for the polymerase activity and an 8kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S44/55 and S55 phosphorylation were less dramatic than S44 and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is likely the major contributor to structural fluctuations that lead to loss of enzymatic activity.
DNA polymerase β is a 39 kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31 kDa domain responsible for the polymerase activity and an 8 kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics simulations. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S55/44 and S55 phosphorylation were less dramatic and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is the major contributor to structural fluctuations that lead to loss of enzymatic activity.
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.