Fungi of the Trichoderma species are valued industrial enzymes in support of the 'zero-waste' technology to convert agro-industrial biomass into valuable products, i.e. nanocellulose (NC). In this study, an in silico approach using substrate docking and molecular dynamic (MD) simulation was used to predict the order of which the multilayers of cellulosic polymers, i.e. lignin, hemicellulose and cellulose in oil palm leaves (OPL) are degraded by fungal enzymes, endocellulase and exocellulase. The study aimed to establish the catalytic tendencies of the enzymes to optimally degrade the cellulosic components of OPL for high yield production of NC. Energy minimized endocellulase and exocellulase models revealed satisfactory scores of PROCHECK (90.0% and 91.2%), Verify3D (97.23% and 98.85%) and ERRAT (95.24% and 91.00%) assessments. Active site prediction by blind docking, COACH meta-server and multiple sequence alignment indicated the catalytic triads for endocellulase and exocellulase were Ser116-His205-Glu249 and Ser382-Arg124-Asp385, respectively. Binding energy of endocellulase docked with hemicellulose (-6.0 kcal mol-1) was the most favourable followed by lignin (-5.6 kcal mol-1) and cellulose (-4.4 kcal mol-1). Exocellulase, contrarily, bonded favorably with lignin (-8.7 kcal mol-1), closely followed by cellulose (-8.5 kcal mol-1) and hemicellulose (-8.4 kcal mol-1). MDs simulations showed that interactions of complexes, endocellulase-hemicellulose and the exocellulase-cellulose being the most stable. Thus, the findings of the study successfully identified the specific actions of sugar-acting enzymes for NC production. Communicated by Ramaswamy H. Sarma.
Alkaline-stable lipases are highly valuable biocatalysts that catalyze reactions under highly basic conditions. Herein, computational predictions of lipase from Acinetobacter haemolyticus and its mutant, Mut-LipKV1 was performed to identify functionally relevant mutations that enhance pH performance under increasing basicity. Mut-LipKV1 was constructed by in silico site directed mutagenesis of several outer loop acidic residues, aspartic acid (Asp) into basic ones, lysine (Lys) at positions 51, 122 and 247, followed by simulation under extreme pH conditions (pH 8.0-pH 12.0). The energy minimized Mut-LipKV1 model exhibited good quality as shown by PROCHECK, ERRAT and Verify3D data that corresponded to 79.2, 88.82 and 89.42% in comparison to 75.2, 86.15, and 95.19% in the wild-type. Electrostatic surface potentials and charge distributions of the Mut-LipKV1 model was more stable and better adapted to conditions of elevated pHs (pH 8.0 - 10.0). Mut-LipKV1 exhibited a mixture of neutral and positive surface charge distribution compared to the predominantly negative charge in the wild-type lipase at pH 8.0. Data of molecular dynamics simulations also supported the increased alkaline-stability of Mut-LipKV1, wherein the lipase was more stable at a higher pH 9.0 (RMSD = ∼0.3 nm, RMSF = ∼0.05-0.2 nm), over the optimal pH 8.0 of the wild-type lipase (RMSD = 0.3 nm, RMSF = 0.05-0.20 nm). Thus, the adaptive strategy of replacing surface aspartic acid to lysine in lipase was successful in yielding a more alkaline-stable Mut-LipKV1 under elevated basic conditions.Communicated by Ramaswamy H. Sarma.
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.
Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzyme deficiency disorder affecting over 400 million individuals worldwide. G6PD protects red blood cells (RBC) from the harmful effects of oxidative substances. There are more than 400 G6PD mutations, of which 186 variants have shown to be linked to G6PD deficiency by decreasing the activity or stability of the enzyme. Different variants manifest different clinical phenotypes which complicate comprehending the mechanism of the disease. In order to carry out computational approaches to elucidate the structural changes of different G6PD variants that are common to the Asian population, a complete G6PD monomer-ligand complex was constructed using AutoDock 4.2, and the molecular dynamics simulation package GROMACS 4.6.7 was used to study the protein dynamics. The G410D and V291M variants were chosen to represent classes I and II respectively and were created by in silico site-directed mutagenesis. Results from the Root mean square deviation (RMSD), Root mean square fluctuation (RMSF) and Radius of gyration (Rg) analyses provided insights on the structure - function relationship for the variants. G410D indicated impaired dimerization and structural NADP binding while the impaired catalytic activity for V291M was indicated by a conformational change at its mutation site.
Less sedimentation and convection in a microgravity environment has become a well-suited condition for growing high quality protein crystals. Thermostable T1 lipase derived from bacterium Geobacilluszalihae has been crystallized using the counter diffusion method under space and earth conditions. Preliminary study using YASARA molecular modeling structure program for both structures showed differences in number of hydrogen bond, ionic interaction, and conformation. The space-grown crystal structure contains more hydrogen bonds as compared with the earth-grown crystal structure. A molecular dynamics simulation study was used to provide insight on the fluctuations and conformational changes of both T1 lipase structures. The analysis of root mean square deviation (RMSD), radius of gyration, and root mean square fluctuation (RMSF) showed that space-grown structure is more stable than the earth-grown structure. Space-structure also showed more hydrogen bonds and ion interactions compared to the earth-grown structure. Further analysis also revealed that the space-grown structure has long-lived interactions, hence it is considered as the more stable structure. This study provides the conformational dynamics of T1 lipase crystal structure grown in space and earth condition.