Short-time dynamics of pH-dependent conformation and substrate binding in the active site of beta-glucosidases: A computational study

A mutant GH3 family β-glucosidase from Oenococcus oeni exhibits superior adaptation to wine stresses and potential for improving wine aroma and phenolic profiles

In this study, we conducted a comprehensive investigation into a GH3 family β-glucosidase (BGL) from the wild-type strain of Oenococcus oeni and its mutated counterpart from the acid-tolerant mutant strain. Our analysis revealed the mutant BGL's remarkable capacity to adapt to wine-related stress conditions, including heightened tolerance to low pH, elevated ethanol concentrations, and metal ions. Additionally, the mutant BGL exhibited superior hydrolytic activity towards various substrates. Through de novo modeling, we identified specific amino acid mutations responsible for its resilience to low pH and high ethanol environments. In simulated wine conditions, the mutant BGL outperformed both wild-type and commercial BGLs, efficiently releasing terpene and phenolic aglycones from glycosides in wine grapes. These findings not only expand our understanding of O. oeni BGLs but also highlight their potential in enhancing wine production. The mutant BGL's enhanced adaptation to wine stress conditions opens promising avenue for improving wine quality and flavor.

Comparative modeling and enzymatic affinity of novel haloacid dehalogenase from Bacillus megaterium strain BHS1 isolated from alkaline Blue Lake in Turkey

This study presents the initial structural model of L-haloacid dehalogenase (DehLBHS1) from Bacillus megaterium BHS1, an alkalotolerant bacterium known for its ability to degrade halogenated environmental pollutants. The model provides insights into the structural features of DehLBHS1 and expands our understanding of the enzymatic mechanisms involved in the degradation of these hazardous pollutants. Key amino acid residues (Arg40, Phe59, Asn118, Asn176, and Trp178) in DehLBHS1 were identified to play critical roles in catalysis and molecular recognition of haloalkanoic acid, essential for efficient binding and transformation of haloalkanoic acid molecules. DehLBHS1 was modeled using I-TASSER, yielding a best TM-score of 0.986 and an RMSD of 0.53 Å. Validation of the model using PROCHECK revealed that 89.2% of the residues were located in the most favored region, providing confidence in its structural accuracy. Molecular docking simulations showed that the non-simulated DehLBHS1 preferred 2,2DCP over other substrates, forming one hydrogen bond with Arg40 and exhibiting a minimum energy of -2.5 kJ/mol. The simulated DehLBHS1 exhibited a minimum energy of -4.3 kJ/mol and formed four hydrogen bonds with Arg40, Asn176, Asp9, and Tyr11, further confirming the preference for 2,2DCP. Molecular dynamics simulations supported this preference, based on various metrics, including RMSD, RMSF, gyration, hydrogen bonding, and molecular distance. MM-PBSA calculations showed that the DehLBHS1-2,2-DCP complex had a markedly lower binding energy (-21.363 ± 1.26 kcal/mol) than the DehLBHS1-3CP complex (-14.327 ± 1.738 kcal/mol). This finding has important implications for the substrate specificity and catalytic function of DehLBHS1, particularly in the bioremediation of 2,2-DCP in contaminated alkaline environments. These results provide a detailed view of the molecular interactions between the enzyme and its substrate and may aid in the development of more efficient biocatalytic strategies for the degradation of halogenated compounds.Communicated by Ramaswamy H. Sarma.

Molecular interactions of trichoderma β-1,4-glucosidase (ThBglT12) with mycelial cell wall components of phytopathogenic Macrophomina phaseolina

Efficacy of a β-1,4-glucosidase from Trichoderma harzianum T12 (ThBglT12) in disrupting the cell wall of the phytopathogenic fungus M. phaseolina (Macrophomina phaseolina) was studied, as the underlying molecular mechanisms of cell wall recognition remains elusive. In this study, the binding location identified by a consensus of residues predicted by COACH tool, blind docking, and multiple sequence alignment revealed that molecular recognition by ThBglT12 occurred through interactions between the α-1,3-glucan, β-1,3-glucan, β-1,3/1,4-glucan, and chitin components of M. phaseolina, with corresponding binding energies of -7.4, -7.6, -7.5 and -7.8 kcal/mol. The residue consensus verified the participation of Glu172, Tyr304, Trp345, Glu373, Glu430, and Trp431 in the active site pocket of ThBglT12 to bind the ligands, of which Trp345 was the common interacting residue. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), total energy, and minimum distance calculation from molecular dynamics (MD) simulation further confirmed the stability and the closeness of the binding ligands into the ThBglT12 active site pocket. The h-bond occupancy by Glu373 and Trp431 instated the role of the nucleophile for substrate recognition and specificity, crucial for cleaving the β-1,4 linkage. Further investigation showed that the proximity of Glu373 to the anomeric carbon of β-1,3/1,4-glucan (3.5 Å) and chitin (5.5 Å) indicates the nucleophiles' readiness to form enzyme-substrate intermediates. Plus, the neighboring water molecule appeared to be correctly positioned and oriented towards the anomeric carbon to hydrolyze the β-1,3/1,4-glucan and chitin, in less than 4.0 Å. In a nutshell, the study verified that the ThBglT12 is a good alternative fungicide to inhibit the growth of M. phaseolina.Communicated by Ramaswamy H. Sarma.

Hierarchical Reactivity of Enzyme-Mediated Phosphorus Recycling from Organic Mixtures by Aspergillus niger Phytase

Biological recycling of inorganic phosphorus (P_i) from organic phosphorus (P_o) compounds by phosphatase-type enzymes, including phytases, is an important contributor to the pool of bioavailable P to plants and microorganisms. However, studies of mixed-substrate reactions with these enzymes are lacking. Here, we explore the reactivity of a phytase extract from the fungus Aspergillus niger toward a heterogeneous mixture containing, in addition to phytate, different structures of environmentally relevant P_o compounds such as ribonucleotides and sugar phosphates. Using a high-resolution liquid chromatography-mass spectrometry method to monitor simultaneously the parent P_o compounds and their by-products, we captured sequential substrate-specific evolution of P_i from the mixture, with faster hydrolysis of multiphosphorylated compounds (phytate, diphosphorylated sugars, and di- and tri-phosphorylated ribonucleotides) than hydrolysis of monophosphorylated compounds (monophosphorylated sugars and monophosphorylated ribonucleotides). The interaction mechanisms and energies revealed by molecular docking simulations of each P_o compound within the enzyme's active site explained the substrate hierarchy observed experimentally. Specifically, the favorable orientation for binding of the negatively charged phosphate moieties with respect to the positive potential surface of the active site was important. Collectively, our findings provide mechanistic insights about the broad but hierarchical role of phytase-type enzymes in P_i recycling from the heterogeneous assembly of P_o compounds in agricultural soils or wastes.

Probing the Effect of Glucose on the Activity and Stability of β-Glucosidase: An All-Atom Molecular Dynamics Simulation Investigation

β-Glucosidase (EC 3.2.1.21) plays an essential role in the removal of glycosyl residues from disaccharide cellobiose to produce glucose during the hydrolysis of lignocellulosic biomass. Although there exist a few β-glucosidase that are tolerant to large concentrations of glucose, these enzymes are typically prone to glucose inhibition. Understanding the basis of this inhibition is important for the production of cheaper biofuels from lignocellulose. In this study, all-atom molecular dynamics simulation at different temperatures and glucose concentrations was used to understand the molecular basis of glucose inhibition of GH1 β-glucosidase (B8CYA8) from Halothermothrix orenii. Our results show that glucose induces a broadening of the active site tunnel through residues lining the tunnel and facilitates the accumulation of glucose. In particular, we observed that glucose accumulates at the tunnel entrance and near the catalytic sites to block substrate accessibility and inhibit enzyme activity. The reduction of enzyme activity was also confirmed experimentally through specific activity measurements in the presence of 0-2.5 M glucose. We also show that the increase in glucose concentrations leads to a decrease in the number of water molecules inside the tunnel to affect substrate hydrolysis. Overall, the results help in understanding the role of residues along the active site tunnel for the engineering of glucose-tolerant β-glucosidase.

Substrate binding versus escape dynamics in a pH-affected fungal beta-glucosidase revealed by molecular dynamics simulations

The cellulolytic ability of fungal species is important to both natural and engineered biocycling of plant matter. One essential step is the conversion of cellobiose into glucose catalyzed by beta-glucosidases. Mutagenesis studies have implicated altering the substrate binding pocket to influence the pH-activity profile of this enzyme. However, structural understanding of the pH-affected substrate binding environment is lacking. Here we conducted molecular dynamics simulations of fully hydrated TrBgl2, a beta-glucosidase of Trichoderma reesei, equilibrated at its optimal pH (pH 6) and two unfavorable pHs (pH 5 and pH 7.5). We identified structural arrangement of specific residues that facilitated substrate escape from the catalytic site at pH 5 but locked the bound substrate in an unfavorable orientation at pH 7.5. For comparative analysis, we also performed simulations of a mutated TrBgl2 with previously demonstrated improved catalysis as a function of pH. We captured the responsible conformational changes in the engineered substrate binding pocket.

Insights into the structural biology of Gaucher disease

Gaucher disease, the most common lysosomal storage disorder, is caused by mutations in the gene encoding the acid-β-glucosidase lysosomal hydrolase enzyme that cleaves glucocerebroside into glucose and ceramide. Reduced enzyme activity and impaired structural stability arise due to >300 known disease-causing mutations. Several of these mutations have also been associated with an increased risk of Parkinson disease (PD). Since the discovery of the acid-β-glucosidase X-ray structure, there have been major advances in our understanding of the structural properties of the protein. Analysis of specific residues has provided insight into their functional and structural importance and provided insight into the pathogenesis of Gaucher disease and the contribution to PD. Disease-causing mutations are positioned throughout the acid-β-glucosidase structure, with many located far from the active site and thus retaining some enzymatic activity however, thus far no clear relationship between mutation location and disease severity has been established. Here, we review the crystal structure of acid-β-glucosidase, while highlighting important structural aspects of the protein in detail. This review discusses the structural stability of acid-β-glucosidase, which can be altered by pH and glycosylation, and explores the relationship between known Gaucher disease and PD mutations, structural stability and disease severity.

Molecular Structural Basis for the Cold Adaptedness of the Psychrophilic β-Glucosidase BglU in Micrococcus antarcticus

Psychrophilic enzymes play crucial roles in cold adaptation of microbes and provide useful models for studies of protein evolution, folding, and dynamic properties. We examined the crystal structure (2.2-Å resolution) of the psychrophilic β-glucosidase BglU, a member of the glycosyl hydrolase 1 (GH1) enzyme family found in the cold-adapted bacterium Micrococcus antarcticus. Structural comparison and sequence alignment between BglU and its mesophilic and thermophilic counterpart enzymes (BglB and GlyTn, respectively) revealed two notable features distinct to BglU: (i) a unique long-loop L3 (35 versus 7 amino acids in others) involved in substrate binding and (ii) a unique amino acid, His299 (Tyr in others), involved in the stabilization of an ordered water molecule chain. Shortening of loop L3 to 25 amino acids reduced low-temperature catalytic activity, substrate-binding ability, the optimal temperature, and the melting temperature (Tm). Mutation of His299 to Tyr increased the optimal temperature, the Tm, and the catalytic activity. Conversely, mutation of Tyr301 to His in BglB caused a reduction in catalytic activity, thermostability, and the optimal temperature (45 to 35°C). Loop L3 shortening and H299Y substitution jointly restored enzyme activity to the level of BglU, but at moderate temperatures. Our findings indicate that loop L3 controls the level of catalytic activity at low temperatures, residue His299 is responsible for thermolability (particularly heat lability of the active center), and long-loop L3 and His299 are jointly responsible for the psychrophilic properties. The described structural basis for the cold adaptedness of BglU will be helpful for structure-based engineering of new cold-adapted enzymes and for the production of mutants useful in a variety of industrial processes at different temperatures.