Tunnel engineering to accelerate product release for better biomass-degrading abilities in lignocellulolytic enzymes

Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase

[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O2 and CO. One strategy to engineer [NiFe] hydrogenases and achieve O2- and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site. In this work, we characterized the unbinding pathways of CO in the complex with the wild-type and 10 different mutants of [NiFe] hydrogenase from Desulfovibrio fructosovorans using τ-random accelerated molecular dynamics (τRAMD) to enhance the sampling of unbinding events. The ranking provided by the relative residence times computed with τRAMD is in agreement with experiments. Extensive data analysis of the simulations revealed that from the two bottlenecks proposed in previous studies for the transit of gas molecules (residues 74 and 122 and residues 74 and 476), only one of them (residues 74 and 122) effectively modulates diffusion and residence times for CO. We also computed pathway probabilities for the unbinding of CO, O2, and H2 from the wild-type [NiFe] hydrogenase, and we observed that while the most probable pathways are the same, the secondary pathways are different. We propose that introducing mutations to block the most probable paths, in combination with mutations to open the main secondary path used by H2, can be a feasible strategy to achieve CO and O2 resistance in the [NiFe] hydrogenase from Desulfovibrio fructosovorans.

Computational investigation of cis-1,4-polyisoprene binding to the latex-clearing protein LcpK30

Latex clearing proteins (Lcps) catalyze the oxidative cleavage of the C = C bonds in cis-1,4-polyisoprene (natural rubber), producing oligomeric compounds that can be repurposed to other materials. The active catalytic site of Lcps is buried inside the protein structure, thus raising the question of how the large hydrophobic rubber chains can access the catalytic center. To improve our understanding of hydrophobic polymeric substrate binding to Lcps and subsequent catalysis, we investigated the interaction of a substrate model containing ten carbon-carbon double bonds with the structurally characterized LcpK30, using multiple computational tools. Prediction of the putative tunnels and cavities in the LcpK30 structure, using CAVER-Pymol plugin 3.0.3, fpocket and Molecular Dynamic (MD) simulations provided valuable insights on how substrate enters from the surface to the buried active site. Two dominant tunnels were discovered that provided feasible routes for substrate binding, and the presence of two hydrophobic pockets was predicted near the heme cofactor. The larger of these pockets is likely to accommodate the substrate and to determine the size distribution of the oligomers. Protein-ligand docking was carried out using GOLD software to predict the conformations and interactions of the substrate within the protein active site. Deeper insight into the protein-substrate interactions, including close-contacts, binding energies and potential cleavage sites in the cis-1,4-polyisoprene, were obtained from MD simulations. Our findings provide further justification that the protein-substrate complexation in LcpK30 is mainly driven by the hydrophobic interactions accompanied by mutual conformational changes of both molecules. Two potential binding modes were identified, with the substrate in either extended or folded conformations. Whilst binding in the extended conformation was most favorable, the folded conformation suggested a preference for cleavage of a central double bond, leading to a preference for oligomers with 5 to 6 C = C bonds. The results provide insight into further enzyme engineering studies to improve catalytic activity and diversify the substrate and product scope of Lcps.

Engineering polysaccharide hydrolases in the product-releasing cleft to alter their product profiles

Polysaccharide hydrolases are enzymes capable of hydrolyzing polysaccharides to generate oligosaccharides that have diverse applications in the food, feed and pharmaceutical industries. However, the detailed mechanisms governing the compositions of their hydrolysates remain poorly understood. Previously, we identified a novel neopullulase Amy117, which exclusively converts pullulan to panose by specifically cleaving α-1,4-glycosidic bonds. Yet, several enzymes with high homology to Amy117 produce a mixture of glucose, maltose and panose during pullulan hydrolysis. To explore this particular phenomenon, we compared the sequences and structures between Amy117 and the maltose amylase ThMA, and identified a specific residue Thr299 in Amy117 (equivalent to His294 in ThMA) within the product-releasing cleft of Amy117, which might be responsible for this characteristic feature. Using structure-based rational design, we have successfully converted the product profiles of pullulan hydrolysates between Amy117 and ThMA by simply altering this key residue. Molecular docking analysis indicated that the key residue at the product-releasing outlet altered the product profile by affecting the panose release rate. Moreover, we modeled the long-chain pullulan substrate G8 to examine its potential conformations and found that G8 might undergo a conformational change in the narrow cleft that allows the Amy117 variant to specifically recognize α-1,6-glycosidic bonds.

Biochemical Characterization of a GH46 Chitosanase Provides Insights into the Novel Digestion Specificity

Endo-chitosanases (EC 3.2.1.132) are generally considered to selectively release functional chito-oligosaccharides (COSs) with degrees of polymerization (DPs) ≥ 2. Although numerous endo-chitosanases have been characterized, the digestion specificity of endo-chitosanases needs to be further explored. In this study, a GH46 endo-chitosanase OUC-CsnPa was cloned, expressed, and characterized from Paenibacillus sp. 1-18. The digestion pattern analysis indicated that OUC-CsnPa could produce monosaccharides from chitotetraose [(GlcN)₄], the smallest recognized substrate, in a random endo-acting manner. Especially, the enzyme specificities during chitosan digestion including the regulation of product abundance through a transglycosylation reaction were also evaluated. It was hypothesized that an insertion region in OUC-CsnPa may form a strong force to be involved in stabilizing (GlcN)₄ at its negative subsite for efficient hydrolysis. This is the first comprehensive report to reveal the digestion specificity and subsite specificity of monosaccharide production by endo-chitosanases. Overall, OUC-CsnPa described here highlights the previously unknown digestion properties of the endo-acting chitosanases and provides a unique example of possible structure-function relationships.

Evaluation of lipase access tunnels and analysis of substance transport in comparison with experimental data

Lipases (E.C. 3.1.1.3) have buried active sites and used access tunnels in the transport of substrates and products for biotransformation processes. Computational methods are used to predict the trajectory and energy profile of ligands through these tunnels, and they complement the experimental methodologies because they filter data, optimizing laboratory time and experimental costs. Access tunnels of Burkholderia cepacia lipase (BCL), Candida rugosa lipase (CRL), and porcine pancreas lipase (PPL) and the transport of fatty acids, alcohols and esters through the tunnels were evaluated using the online server CaverWeb V1.0, and server calculation results were compared with experimental data (productivity). BCL showed higher productivity with palmitic acid-C16:0 (4029.95 µmol/h mg); CRL obtained productivity for oleic acid-C18:1 (380.80 µmol/h mg), and PPL achieved productivity for lauric acid-C12:0 (71.27 µmol/h mg). The highest probability of transport for BCL is through the tunnels 1 and 2, for CRL through the tunnel 1, and for PPL through the tunnels 1, 2, 3 and 4. Thus, the best in silico result was the transport of the substrates palmitic acid and ethanol and product ethyl palmitate in tunnel 1 of BCL. This result corroborates with the best result for the productivity data (higher productivity for BCL with palmitic acid-4029.95 µmol/h mg). The combination of in silico evaluation and experimental data gave similar results, demonstrating that in silico approaches are a promising alternative for reducing screening tests and minimizing laboratory time in the bio-catalysis area by identifying the lipases with the greatest reaction potential, as in the case of this proposal.

Advances in Rational Protein Engineering toward Functional Architectures and Their Applications in Food Science

Protein biomolecules including enzymes, cagelike proteins, and specific peptides have been continuously exploited as functional biomaterials applied in catalysis, nutrient delivery, and food preservation in food-related areas. However, natural proteins usually function well in physiological conditions, not industrial conditions, or may possess undesirable physical and chemical properties. Currently, rational protein design as a valuable technology has attracted extensive attention for the rational engineering or fabrication of ideal protein biomaterials with novel properties and functionality. This article starts with the underlying knowledge of protein folding and assembly and is followed by the introduction of the principles and strategies for rational protein design. Basic strategies for rational protein engineering involving experienced protein tailoring, computational prediction, computation redesign, and de novo protein design are summarized. Then, we focus on the recent progress of rational protein engineering or design in the application of food science, and a comprehensive summary ranging from enzyme manufacturing to cagelike protein nanocarriers engineering and antimicrobial peptides preparation is given. Overall, this review highlights the importance of rational protein engineering in food biomaterial preparation which could be beneficial for food science.

Improvement of catalytic performance of endoglucanase CgEndo from Colletotrichum graminicola by site-directed mutagenesis

In order to improve the catalytic efficiency of cellulase for more effective utilization of lignocellulose, a novel endoglucanase (CgEndo) from Colletotrichum graminicola was expressed by Pichia pastoris X33 and modified by site-directed mutagenesis. Two mutants, Y63S and N20D/S113T, with 62.31% and 57.14% increased enzyme activities were obtained, respectively. On this basis, their biochemical properties, kinetic parameters, structural information as well as the application in biomass degradation were investigated and compared with the wild-type CgEngo. The results indicated that the mutation Y63S and N20D/S113T resulted in an improvement of proximity between enzyme and substrate through conformational changes of the catalytic region, which might contribute to the higher enzyme activities and catalysis efficiency (K_cat/K_m) of Y63S and N20D/S113T. These findings laid important foundation for the further engineering of this endoglucanase and practical application in efficient degradation of cellulosic biomass in nature.

Substrate Access Mechanism in a Novel Membrane-Bound Phospholipase A of Pseudomonas aeruginosa Concordant with Specificity and Regioselectivity

PlaF is a cytoplasmic membrane-bound phospholipase A₁ from Pseudomonas aeruginosa that alters the membrane glycerophospholipid (GPL) composition and fosters the virulence of this human pathogen. PlaF activity is regulated by a dimer-to-monomer transition followed by tilting of the monomer in the membrane. However, how substrates reach the active site and how the characteristics of the active site tunnels determine the activity, specificity, and regioselectivity of PlaF for natural GPL substrates have remained elusive. Here, we combined unbiased and biased all-atom molecular dynamics (MD) simulations and configurational free-energy computations to identify access pathways of GPL substrates to the catalytic center of PlaF. Our results map out a distinct tunnel through which substrates access the catalytic center. PlaF variants with bulky tryptophan residues in this tunnel revealed decreased catalysis rates due to tunnel blockage. The MD simulations suggest that GPLs preferably enter the active site with the sn-1 acyl chain first, which agrees with the experimentally demonstrated PLA₁ activity of PlaF. We propose that the acyl chain-length specificity of PlaF is determined by the structural features of the access tunnel, which results in favorable free energy of binding of medium-chain GPLs. The suggested egress route conveys fatty acid (FA) products to the dimerization interface and, thus, contributes to understanding the product feedback regulation of PlaF by FA-triggered dimerization. These findings open up opportunities for developing potential PlaF inhibitors, which may act as antibiotics against P. aeruginosa.

Probing the function of a ligand-modulated dynamic tunnel in bifunctional proline utilization A (PutA)

In many bacteria, the reactions of proline catabolism are catalyzed by the bifunctional enzyme known as proline utilization A (PutA). PutA catalyzes the two-step oxidation of l-proline to l-glutamate using distinct proline dehydrogenase (PRODH) and l-glutamate-γ-semialdehyde dehydrogenase (GSALDH) active sites, which are separated by over 40 Å and connected by a complex tunnel system. The tunnel system consists of a main tunnel that connects the two active sites and functions in substrate channeling, plus six ancillary tunnels whose functions are unknown. Here we used tunnel-blocking mutagenesis to probe the role of a dynamic ancillary tunnel (tunnel 2a) whose shape is modulated by ligand binding to the PRODH active site. The 1.90 Å resolution crystal structure of Geobacter sulfurreducens PutA variant A206W verified that the side chain of Trp206 cleanly blocks tunnel 2a without perturbing the surrounding structure. Steady-state kinetic measurements indicate the mutation impaired PRODH activity without affecting the GSALDH activity. Single-turnover experiments corroborated a severe impairment of PRODH activity with flavin reduction decreased by nearly 600-fold in A206W relative to wild-type. Substrate channeling is also significantly impacted as A206W exhibited a 3000-fold lower catalytic efficiency in coupled PRODH-GSALDH activity assays, which measure NADH formation as a function of proline. The structure suggests that Trp206 inhibits binding of the substrate l-proline by preventing the formation of a conserved glutamate-arginine ion pair and closure of the PRODH active site. Our data are consistent with tunnel 2a serving as an open space through which the glutamate of the ion pair travels during the opening and closing of the active site in response to binding l-proline. These results confirm the essentiality of the conserved ion pair in binding l-proline and support the hypothesis that the ion pair functions as a gate that controls access to the PRODH active site.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

The penicillin binding protein 1A of Helicobacter pylori, its amoxicillin binding site and access routes

BACKGROUND

Amoxicillin-resistant H. pylori strains are increasing worldwide. To explore the potential resistance mechanisms involved, the 3D structure modeling and access tunnel prediction for penicillin-binding proteins (PBP1A) was performed, based on the Streptococcus pneumoniae, PBP 3D structure. Molecular covalent docking was used to determine the interactions between amoxicillin (AMX) and PBP1A.

RESULTS

The AMX-Ser368 covalent complex interacts with the binding site residues (Gly367, Ala369, ILE370, Lys371, Tyr416, Ser433, Thr541, Thr556, Gly557, Thr558, and Asn560) of PBP1A, non-covalently. Six tunnel-like structures, accessing the PBP1A binding site, were characterized, using the CAVER algorithm. Tunnel-1 was the ultimate access route, leading to the drug catalytic binding residue (Ser368). This tunnel comprises of eighteen amino acid residues, 8 of which are shared with the drug binding site. Subsequently, to screen the presence of PBP1A mutations, in the binding site and tunnel residues, in our clinical strains, in vitro assays were performed. H. pylori strains, isolated under gastroscopy, underwent AMX susceptibility testing by E-test. Of the 100 clinical strains tested, 4 were AMX-resistant. The transpeptidase domain of the pbp1a gene of these resistant, plus 10 randomly selected AMX-susceptible strains, were amplified and sequenced. Of the amino acids lining the tunnel-1 and binding site residues, three (Ser414Arg, Val469Met and Thr556Ser) substitutions, were detected in 2 of the 4 resistant and none of the sequenced susceptible strains, respectively.

CONCLUSIONS

We hypothesize that mutations in amino acid residues lining the binding site and/or tunnel-1, resulting in conformational/spatial changes, may block drug binding to PBP1A and cause AMX resistance.

Tunnel engineering for modulating the substrate preference in cytochrome P450BsβHI

An active site is normally located inside enzymes, hence substrates should go through a tunnel to access the active site. Tunnel engineering is a powerful strategy for refining the catalytic properties of enzymes. Here, P450BsβHI (Q85H/V170I) derived from hydroxylase P450Bsβ from Bacillus subtilis was chosen as the study model, which is reported as a potential decarboxylase. However, this enzyme showed low decarboxylase activity towards long-chain fatty acids. Here, a tunnel engineering campaign was performed for modulating the substrate preference and improving the decarboxylation activity of P450BsβHI. The finally obtained BsβHI-F79A variant had a 15.2-fold improved conversion for palmitic acid; BsβHI-F173V variant had a 3.9-fold improved conversion for pentadecanoic acid. The study demonstrates how the substrate preference can be modulated by tunnel engineering strategy.

Structural insights into xylanase mutant 254RL1 for improved activity and lower pH optimum

Xylanases degrade xylan to valuable end products. In our previous study, the alkaline xylanase S7-xyl from Bacillus halodurans S7 was engineered by rational design and the best mutant xylanase 254RL1 exhibited 3.4-fold improvements in specific activity at pH 9.0. Further research found that the enzyme activity at pH 6.0 was almost 2-fold than that at pH 9.0. To elucidate the reason of enhanced performance of 254RL1 at decreased pH optimum, we determined the X-ray crystal structure of 254RL1 at 2.21 Å resolution. The structural analysis revealed that the mutations enlarged the opening of the access tunnel and shortened the tunnel. Moreover, the mutations changed the hydrogen bond network around the catalytic residue and decreased the pK_a value of acid-base catalyst E159 which reduced the pH optimum of the xylanase. The result provided the basis for the acid-alkaline engineering of the glycoside hydrolases.

Progress in Ameliorating Beneficial Characteristics of Microbial Cellulases by Genetic Engineering Approaches for Cellulose Saccharification

Lignocellulosic biomass is a renewable and sustainable energy source. Cellulases are the enzymes that cleave β-1, 4-glycosidic linkages in cellulose to liberate sugars that can be fermented to ethanol, butanol, and other products. Low enzyme activity and yield, and thermostability are, however, some of the limitations posing hurdles in saccharification of lignocellulosic residues. Recent advancements in synthetic and systems biology have generated immense interest in metabolic and genetic engineering that has led to the development of sustainable technology for saccharification of lignocellulosics in the last couple of decades. There have been several attempts in applying genetic engineering in the production of a repertoire of cellulases at a low cost with a high biomass saccharification. A diverse range of cellulases are produced by different microbes, some of which are being engineered to evolve robust cellulases. This review summarizes various successful genetic engineering strategies employed for improving cellulase kinetics and cellulolytic efficiency.