Crystal structures of the cellulase Cel48F in complex with inhibitors and substrates give insights into its processive action

Fungicidal activities and molecular docking of the marine alga Ulva lactuca and Punica granatum peel extracts on Alternaria tomato spot disease

In this study, we utilized pomegranate peel and marine algae Ulva lactuca (U. lactuca) as rich and sustained sources of bioactive compounds to combat tomato-black spot disease. n-Hexane extracts from the peel of pomegranate (Punica granatum) (PPE) and the marine algal biomass U. lactuca (ULE) were used alone and in combinations to verify their impact against Alternaria alternata (A. alternata). The applied extracts exhibited severe destructive effects on both fungal growth and structure such as mycelia malformation, underdeveloped conidia, cell wall deformation, and shrinkage. Moreover, increased deformations and protrusions, and notch-like structures, were noticed in A. alternata mycelia treated with mixed extracts (PPE and ULE) compared to all other treatments. The protein and reduced sugar contents in tomato fruits were significantly increased in the infected fruits with A. alternata. The highest enzyme activities of pectinase, cellulase, catalase (CAT), and ascorbate peroxidase (APX) were recorded in infected tomatoes in comparison with the healthy ones. Molecular docking studies showed that each extract is rich with bioactive compounds that have a promising inhibition effect on A. alternata cellulases. Pomegranate and Ulva extract showed promising antifungal activity against A. alternata which revealed their feasibility and applicability as biocontrol agents in postharvest disease management and food preservation against fungal pathogens.

Impact of cellulose properties on enzymatic degradation by bacterial GH48 enzymes: Structural and mechanistic insights from processive Bacillus licheniformis Cel48B cellulase

Processive cellulases are highly efficient molecular engines involved in the cellulose breakdown process. However, the mechanism that processive bacterial enzymes utilize to recruit and retain cellulose strands in the catalytic site remains poorly understood. Here, integrated enzymatic assays, protein crystallography and computational approaches were combined to study the enzymatic properties of the processive BlCel48B cellulase from Bacillus licheniformis. Hydrolytic efficiency, substrate binding affinity, cleavage patterns, and the apparent processivity of bacterial BlCel48B are significantly impacted by the cellulose size and its surface morphology. BlCel48B crystallographic structure was solved with ligands spanning -5 to -2 and +1 to +2 subsites. Statistical coupling analysis and molecular dynamics show that co-evolved residues on active site are critical for stabilizing ligands in the catalytic tunnel. Our results provide mechanistic insights into BlCel48B molecular-level determinants of activity, substrate binding, and processivity on insoluble cellulose, thus shedding light on structure-activity correlations of GH48 family members in general.

Convergent evolution of processivity in bacterial and fungal cellulases

Some cellulases exhibit “processivity”: the ability to degrade crystalline cellulose through successive hydrolytic catalytic reactions without the release of the enzyme from the substrate surface. We previously observed the movement of fungal processive cellulases by high-speed atomic force microscopy, and here, we use the same technique to directly observe the processive movement of bacterial cellobiohydrolases settling a long-standing controversy. Although fungal and bacterial processive cellulases have completely different protein folds, they have evolved to acquire processivity through the same strategy of adding subsites to extend the substrate-binding site and forming a tunnel-like active site by increasing the number of loops covering the active site. This represents an example of protein-level convergent evolution to acquire the same functions from different ancestors.

Cellulose is the most abundant biomass on Earth, and many microorganisms depend on it as a source of energy. It consists mainly of crystalline and amorphous regions, and natural degradation of the crystalline part is highly dependent on the degree of processivity of the degrading enzymes (i.e., the extent of continuous hydrolysis without detachment from the substrate cellulose). Here, we report high-speed atomic force microscopic (HS-AFM) observations of the movement of four types of cellulases derived from the cellulolytic bacteria Cellulomonas fimi on various insoluble cellulose substrates. The HS-AFM images clearly demonstrated that two of them (CfCel6B and CfCel48A) slide on crystalline cellulose. The direction of processive movement of CfCel6B is from the nonreducing to the reducing end of the substrate, which is opposite that of processive cellulase Cel7A of the fungus Trichoderma reesei (TrCel7A), whose movement was first observed by this technique, while CfCel48A moves in the same direction as TrCel7A. When CfCel6B and TrCel7A were mixed on the same substrate, “traffic accidents” were observed, in which the two cellulases blocked each other’s progress. The processivity of CfCel6B was similar to those of fungal family 7 cellulases but considerably higher than those of fungal family 6 cellulases. The results indicate that bacteria utilize family 6 cellulases as high-processivity enzymes for efficient degradation of crystalline cellulose, whereas family 7 enzymes have the same function in fungi. This is consistent with the idea of convergent evolution of processive cellulases in fungi and bacteria to achieve similar functionality using different protein foldings.

Insights into the mechanism(s) of digestion of crystalline cellulose by plant class C GH9 endoglucanases

Biofuels such as γ-valerolactone, bioethanol, and biodiesel are derived from potentially fermentable cellulose and vegetable oils. Plant class C GH9 endoglucanases are CBM49-encompassing hydrolases that cleave the β (1 → 4) glycosidic linkage of contiguous D-glucopyranose residues of crystalline cellulose. Here, I analyse 3D-homology models of characterised and putative class C enzymes to glean insights into the contribution of the GH9, linker, and CBM49 to the mechanism(s) of crystalline cellulose digestion. Crystalline cellulose may be accommodated in a surface groove which is imperfectly bounded by the GH9_CBM49, GH9_linker, and linker_CBM49 surfaces and thence digested in a solvent accessible subsurface cavity. The physical dimensions and distortions thereof, of the groove, are mediated in part by the bulky side chains of aromatic amino acids that comprise it and may also result in a strained geometry of the bound cellulose polymer. These data along with an almost complete absence of measurable cavities, along with poorly conserved, hydrophobic, and heterogeneous amino acid composition, increased atomic motion of the CBM49_linker junction, and docking experiements with ligands of lower degrees of polymerization suggests a modulatory rather than direct role for CBM49 in catalysis. Crystalline cellulose is the de facto substrate for CBM-containing plant and non-plant GH9 enzymes, a finding supported by exceptional sequence- and structural-homology. However, despite the implied similarity in general acid-base catalysis of crystalline cellulose, this study also highlights qualitative differences in substrate binding and glycosidic bond cleavage amongst class C members. Results presented may aid the development of novel plant-based GH9 endoglucanases that could extract and utilise potential fermentable carbohydrates from biomass.

Graphical Abstract

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.

Crystal Structures of the C-Terminally Truncated Endoglucanase Cel9Q from Clostridium thermocellum Complexed with Cellodextrins and Tris

Endoglucanase CtCel9Q is one of the enzyme components of the cellulosome, which is an active cellulase system in the thermophile Clostridium thermocellum. The precursor form of CtCel9Q comprises a signal peptide, a glycoside hydrolase family 9 catalytic domain, a type 3c carbohydrate-binding module (CBM), and a type I dockerin domain. Here, we report the crystal structures of C-terminally truncated CtCel9Q (CtCel9QΔc) complexed with Tris, Tris+cellobiose, cellobiose+cellotriose, cellotriose, and cellotetraose at resolutions of 1.50, 1.70, 2.05, 2.05 and 1.75 Å, respectively. CtCel9QΔc forms a V-shaped homodimer through residues Lys529-Glu542 on the type 3c CBM, which pairs two β-strands (β4 and β5 of the CBM). In addition, a disulfide bond was formed between the two Cys535 residues of the protein monomers in the asymmetric unit. The structures allow the identification of four minus (-) subsites and two plus (+) subsites; this is important for further understanding the structural basis of cellulose binding and hydrolysis. In the oligosaccharide-free and cellobiose-bound CtCel9QΔc structures, a Tris molecule was found to be bound to three catalytic residues of CtCel9Q and occupied subsite -1 of the CtCel9Q active-site cleft. Moreover, the enzyme activity assay in the presence of 100 mm Tris showed that the Tris almost completely suppressed CtCel9Q hydrolase activity.

Effects of Tyr555 and Trp678 on the processivity of cellobiohydrolase A from Ruminiclostridium thermocellum: A simulation study

Cellobiohydrolase A from Ruminiclostridium thermocellum (Cbh9A) is a processive exoglucanase from family 9 and is an important cellobiohydrolase that hydrolyzes cello-oligosaccharide into cellobiose. Residues Tyr555 and Trp678 considerably affect catalytic activity, but their mechanisms are still unknown. To investigate how the Tyr555 and Trp678 affect the processivity of Cbh9A, conventional molecular dynamics, steered molecular dynamics, and free energy calculation were performed to simulate the processive process of wild type (WT)-Cbh9A, Y555S mutant, and W678G mutant. Analysis of simulation results suggests that the binding free energies between the substrate and WT-Cbh9A are lower than those of Y555S and W678G mutants. The pull forces and energy barrier in Y555S and W678G mutants also reduced significantly during the steered molecular dynamics (SMD) simulation compared with that of the WT-Cbh9A. And the potential mean force calculations showed that the pulling energy barrier of Y555S and W678G mutants is much lower than that of WT-Cbh9A.

Origin, evolution, and divergence of plant class C GH9 endoglucanases

BACKGROUND

Glycoside hydrolases of the GH9 family encode cellulases that predominantly function as endoglucanases and have wide applications in the food, paper, pharmaceutical, and biofuel industries. The partitioning of plant GH9 endoglucanases, into classes A, B, and C, is based on the differential presence of transmembrane, signal peptide, and the carbohydrate binding module (CBM49). There is considerable debate on the distribution and the functions of these enzymes which may vary in different organisms. In light of these findings we examined the origin, emergence, and subsequent divergence of plant GH9 endoglucanases, with an emphasis on elucidating the role of CBM49 in the digestion of crystalline cellulose by class C members.

RESULTS

Since, the digestion of crystalline cellulose mandates the presence of a well-defined set of aromatic and polar amino acids and/or an attributable domain that can mediate this conversion, we hypothesize a vertical mode of transfer of genes that could favour the emergence of class C like GH9 endoglucanase activity in land plants from potentially ancestral non plant taxa. We demonstrated the concomitant occurrence of a GH9 domain with CBM49 and other homologous carbohydrate binding modules, in putative endoglucanase sequences from several non-plant taxa. In the absence of comparable full length CBMs, we have characterized several low strength patterns that could approximate the CBM49, thereby, extending support for digestion of crystalline cellulose to other segments of the protein. We also provide data suggestive of the ancestral role of putative class C GH9 endoglucanases in land plants, which includes detailed phylogenetics and the presence and subsequent loss of CBM49, transmembrane, and signal peptide regions in certain populations of early land plants. These findings suggest that classes A and B of modern vascular land plants may have emerged by diverging directly from CBM49 encompassing putative class C enzymes.

CONCLUSION

Our detailed phylogenetic and bioinformatics analysis of putative GH9 endoglucanase sequences across major taxa suggests that plant class C enzymes, despite their recent discovery, could function as the last common ancestor of classes A and B. Additionally, research into their ability to digest or inter-convert crystalline and amorphous forms of cellulose could make them lucrative candidates for engineering biofuel feedstock.

Conformational diversity analysis reveals three functional mechanisms in proteins

Protein motions are a key feature to understand biological function. Recently, a large-scale analysis of protein conformational diversity showed a positively skewed distribution with a peak at 0.5 Å C-alpha root-mean-square-deviation (RMSD). To understand this distribution in terms of structure-function relationships, we studied a well curated and large dataset of ~5,000 proteins with experimentally determined conformational diversity. We searched for global behaviour patterns studying how structure-based features change among the available conformer population for each protein. This procedure allowed us to describe the RMSD distribution in terms of three main protein classes sharing given properties. The largest of these protein subsets (~60%), which we call "rigid" (average RMSD = 0.83 Å), has no disordered regions, shows low conformational diversity, the largest tunnels and smaller and buried cavities. The two additional subsets contain disordered regions, but with differential sequence composition and behaviour. Partially disordered proteins have on average 67% of their conformers with disordered regions, average RMSD = 1.1 Å, the highest number of hinges and the longest disordered regions. In contrast, malleable proteins have on average only 25% of disordered conformers and average RMSD = 1.3 Å, flexible cavities affected in size by the presence of disordered regions and show the highest diversity of cognate ligands. Proteins in each set are mostly non-homologous to each other, share no given fold class, nor functional similarity but do share features derived from their conformer population. These shared features could represent conformational mechanisms related with biological functions.

Investigation of an "alternate water supply system" in enzymatic hydrolysis in the processive endocellulase Cel7A from Rasamsonia emersonii by molecular dynamics simulation

Cel7A from Rasamsonia emersonii is one of the processive endocellulases classified under family 7 glycoside hydrolase. Molecular dynamics simulations were carried out to obtain the optimized sliding and hydrolyzing conformations, in which the reducing ends of sugar chains are located on different sites. Hydrogen bonds are investigated to clarify the interactions between protein and substrate in either conformation. Nine hydrogen bonding interactions are identified in the sliding conformation, and six similar interactions are also found correspondingly in the hydrolyzing conformation. In addition, four strong hydrophobic interactions are also determined. The domain cross-correlation map analysis shows movement correlation of protein including autocorrelation between residues. The root mean square fluctuations analysis represents the various flexibilities of different fragment in the two conformations. Comparing the two conformations reveals the water-supply mechanism of selective hydrolysis of cellulose in Cel7A. The mechanism can be described as follow. When the reducing end of substrate slides from the unhydrolyzing site (sliding conformation) to the hydrolyzing site (hydrolyzing conformation), His225 is pushed down and rotated, the rotation leads to the movement of Glu209 with the interstrand hydrogen bonding in β-sheet. It further makes Asp211 close to the hydrolysis center and provides a water molecule bounding on its carboxyl in the previous unhydrolyzing site. After the hydrolysis takes place and the product is excluded from the enzyme, the Asp211 comes back to its initial position. In summary, Asp211 acts as an elevator to transport outer water molecules into the hydrolysis site for every other glycosidic bond.

Using Carbohydrate Interaction Assays to Reveal Novel Binding Sites in Carbohydrate Active Enzymes

Carbohydrate active enzymes often contain auxiliary binding sites located either on independent domains termed carbohydrate binding modules (CBMs) or as so-called surface binding sites (SBSs) on the catalytic module at a certain distance from the active site. The SBSs are usually critical for the activity of their cognate enzyme, though they are not readily detected in the sequence of a protein, but normally require a crystal structure of a complex for their identification. A variety of methods, including affinity electrophoresis (AE), insoluble polysaccharide pulldown (IPP) and surface plasmon resonance (SPR) have been used to study auxiliary binding sites. These techniques are complementary as AE allows monitoring of binding to soluble polysaccharides, IPP to insoluble polysaccharides and SPR to oligosaccharides. Here we show that these methods are useful not only for analyzing known binding sites, but also for identifying new ones, even without structural data available. We further verify the chosen assays discriminate between known SBS/CBM containing enzymes and negative controls. Altogether 35 enzymes are screened for the presence of SBSs or CBMs and several novel binding sites are identified, including the first SBS ever reported in a cellulase. This work demonstrates that combinations of these methods can be used as a part of routine enzyme characterization to identify new binding sites and advance the study of SBSs and CBMs, allowing them to be detected in the absence of structural data.

Strategies to reduce end-product inhibition in family 48 glycoside hydrolases

Family 48 cellobiohydrolases are some of the most abundant glycoside hydrolases in nature. They are able to degrade cellulosic biomass and therefore serve as good enzyme candidates for biofuel production. Family 48 cellulases hydrolyze cellulose chains via a processive mechanism, and produce end products composed primarily of cellobiose as well as other cellooligomers (dp ≤ 4). The challenge of utilizing cellulases in biofuel production lies in their extremely slow turnover rate. A factor contributing to the low enzyme activity is suggested to be product binding to enzyme and the resulting performance inhibition. In this study, we quantitatively evaluated the product inhibitory effect of four family 48 glycoside hydrolases using molecular dynamics simulations and product expulsion free-energy calculations. We also suggested a series of single mutants of the four family 48 glycoside hydrolases with theoretically reduced level of product inhibition. The theoretical calculations provide a guide for future experimental studies designed to produce mutant cellulases with enhanced activity.

Structures of exoglucanase from Clostridium cellulovorans: cellotetraose binding and cleavage

Exoglucanase/cellobiohydrolase (EC 3.2.1.176) hydrolyzes a β-1,4-glycosidic bond from the reducing end of cellulose and releases cellobiose as the major product. Three complex crystal structures of the glycosyl hydrolase 48 (GH48) cellobiohydrolase S (ExgS) from Clostridium cellulovorans with cellobiose, cellotetraose and triethylene glycol molecules were solved. The product cellobiose occupies subsites +1 and +2 in the open active-site cleft of the enzyme-cellotetraose complex structure, indicating an enzymatic hydrolysis function. Moreover, three triethylene glycol molecules and one pentaethylene glycol molecule are located at active-site subsites -2 to -6 in the structure of the ExgS-triethylene glycol complex shown here. Modelling of glucose into subsite -1 in the active site of the ExgS-cellobiose structure revealed that Glu50 acts as a proton donor and Asp222 plays a nucleophilic role.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Kinetics of cellobiohydrolase (Cel7A) variants with lowered substrate affinity

Cellobiohydrolases are exo-active glycosyl hydrolases that processively convert cellulose to soluble sugars, typically cellobiose. They effectively break down crystalline cellulose and make up a major component in industrial enzyme mixtures used for deconstruction of lignocellulosic biomass. Identification of the rate-limiting step for cellobiohydrolases remains controversial, and recent reports have alternately suggested either association (on-rate) or dissociation (off-rate) as the overall bottleneck. Obviously, this uncertainty hampers both fundamental mechanistic understanding and rational design of enzymes with improved industrial applicability. To elucidate the role of on- and off-rates, respectively, on the overall kinetics, we have expressed a variant in which a tryptophan residue (Trp-38) in the middle of the active tunnel has been replaced with an alanine. This mutation weakens complex formation, and the population of substrate-bound W38A was only about half of the wild type. Nevertheless, the maximal, steady-state rate was twice as high for the variant enzyme. It is argued that these opposite effects on binding and activity can be reconciled if the rate-limiting step is after the catalysis (i.e. in the dissociation process).

Molecular Dynamics and Metadynamics Simulations of the Cellulase Cel48F

Molecular dynamics (MD) and metadynamics techniques were used to study the cellulase Cel48F-sugar. Cellulase is enzyme that breaks cellulose fibers into small sugar units and is potentially useful in second generation alcohol production. In MD simulations, the overall structure of equilibrated Cel48F did not significantly change along the trajectory, retaining root mean square deviation below 0.15 nm. A set of 15 residues interacting with the sugar chains via hydrogen bonding throughout the simulation was observed. The free energy of dissociation (ΔG_diss.) of the chains in the catalytic tunnel of Cel48F was determined by metadynamics. The ΔG_diss. values of the chains entering and leaving the wild-type Cel48F cavity were 13.9 and 62.1 kcal/mol, respectively. We also mutated the E542 and Q543 to alanine residue and obtained ΔG_diss. of 41.8 and 45.9 kcal/mol, respectively. These mutations were found to facilitate smooth dissociation of the sugar chain across the Cel48F tunnel. At the entry of the Cel48F tunnel, three residues were mutated to alanine: T110, T213, and L274. Contrary to the T110A-Cel48F, the mutants T213-Cel48F and L274-Cel48F prevented the sugar chain from passing across the leaving site. The present results can be a guideline in mutagenesis studies to improve processing by Cel48F.

Molecular and biochemical analyses of CbCel9A/Cel48A, a highly secreted multi-modular cellulase by Caldicellulosiruptor bescii during growth on crystalline cellulose

During growth on crystalline cellulose, the thermophilic bacterium Caldicellulosiruptor bescii secretes several cellulose-degrading enzymes. Among these enzymes is CelA (CbCel9A/Cel48A), which is reported as the most highly secreted cellulolytic enzyme in this bacterium. CbCel9A/Cel48A is a large multi-modular polypeptide, composed of an N-terminal catalytic glycoside hydrolase family 9 (GH9) module and a C-terminal GH48 catalytic module that are separated by a family 3c carbohydrate-binding module (CBM3c) and two identical CBM3bs. The wild-type CbCel9A/Cel48A and its truncational mutants were expressed in Bacillus megaterium and Escherichia coli, respectively. The wild-type polypeptide released twice the amount of glucose equivalents from Avicel than its truncational mutant that lacks the GH48 catalytic module. The truncational mutant harboring the GH9 module and the CBM3c was more thermostable than the wild-type protein, likely due to its compact structure. The main hydrolytic activity was present in the GH9 catalytic module, while the truncational mutant containing the GH48 module and the three CBMs was ineffective in degradation of either crystalline or amorphous cellulose. Interestingly, the GH9 and/or GH48 catalytic modules containing the CBM3bs form low-density particles during hydrolysis of crystalline cellulose. Moreover, TM3 (GH9/CBM3c) and TM2 (GH48 with three CBM3 modules) synergistically hydrolyze crystalline cellulose. Deletion of the CBM3bs or mutations that compromised their binding activity suggested that these CBMs are important during hydrolysis of crystalline cellulose. In agreement with this observation, seven of nine genes in a C. bescii gene cluster predicted to encode cellulose-degrading enzymes harbor CBM3bs. Based on our results, we hypothesize that C. bescii uses the GH48 module and the CBM3bs in CbCel9A/Cel48A to destabilize certain regions of crystalline cellulose for attack by the highly active GH9 module and other endoglucanases produced by this hyperthermophilic bacterium.

Sequence, structure, and evolution of cellulases in glycoside hydrolase family 48

Currently, the cost of cellulase enzymes remains a key economic impediment to commercialization of biofuels. Enzymes from glycoside hydrolase family 48 (GH48) are a critical component of numerous natural lignocellulose-degrading systems. Although computational mining of large genomic data sets is a promising new approach for identifying novel cellulolytic activities, current computational methods are unable to distinguish between cellulases and enzymes with different substrate specificities that belong to the same protein family. We show that by using a robust computational approach supported by experimental studies, cellulases and non-cellulases can be effectively identified within a given protein family. Phylogenetic analysis of GH48 showed non-monophyletic distribution, an indication of horizontal gene transfer. Enzymatic function of GH48 proteins coded by horizontally transferred genes was verified experimentally, which confirmed that these proteins are cellulases. Computational and structural studies of GH48 enzymes identified structural elements that define cellulases and can be used to computationally distinguish them from non-cellulases. We propose that the structural element that can be used for in silico discrimination between cellulases and non-cellulases belonging to GH48 is an ω-loop located on the surface of the molecule and characterized by highly conserved rare amino acids. These markers were used to screen metagenomics data for "true" cellulases.

The crystallization and structural analysis of cellulases (and other glycoside hydrolases): strategies and tactics

The three-dimensional (3-D) structures of cellulases, and other glycoside hydrolases, are a central feature of research in carbohydrate chemistry and biochemistry. 3-D structure is used to inform protein engineering campaigns, both academic and industrial, which are typically used to improve the stability or activity of an enzyme. Examples of classical protein engineering goals include higher thermal stability, reduced metal-ion dependency, detergent and protease resistance, decreased product inhibition, and altered specificity. 3-D structure may also be used to interpret the behavior of enzyme variants that are derived from screening or random mutagenesis approaches, with a view to establishing an iterative design process. In other areas, 3-D structure is used as one of the many tools to probe enzymatic catalysis, typically dovetailing with physical organic chemistry approaches to provide complete reaction mechanisms for enzymes by visualizing catalytic site interactions at different stages of the reaction. Such mechanistic insight is not only fundamentally important, impacting on inhibitor and drug design approaches with ramifications way beyond cellulose hydrolysis, but also provides the framework for the design of enzyme variants to use as biocatalysts for the synthesis of bespoke oligosaccharides. Here we review some of the strategies and tactics that may be applied to the X-ray structure solution of cellulases (and other carbohydrate-active enzymes). The general approach is first to decide why you are doing the work, then to establish correct domain boundaries for truncated constructs (typically the catalytic domain only), and finally to pursue crystallization of pure, homogeneous, and monodisperse protein with appropriate ligand and additive combinations. Cellulase-specific strategies are important for the delineation of domain boundaries, while glycoside hydrolases generally also present challenges and opportunities for the selection and optimization of ligands to both aid crystallization, and also provide structural and mechanistic insight. As the many roles for plant cell wall degrading enzymes increase, so does the need for rapid high-quality structure determination to provide a sound structural foundation for understanding mechanism and specificity, and for future protein engineering strategies.

Measuring processivity

Natural cellulolytic enzyme systems as well as leading commercial cellulase cocktails are dominated by enzymes that degrade cellulose chains in a processive manner. Despite the abundance of processivity among natural cellulases, the molecular basis as well as the biotechnological implications of this mechanism are only partly understood. One of the major limitations lies in the fact that it is not straightforward to measure and quantify processivity in what essentially are biphasic experimental systems. Here, we describe and discuss both well-established methods and newer methods for measuring cellulase processivity. In addition, we discuss recent insights from studies on chitinases that may help direct further studies on processivity in cellulases.

Determination of the catalytic base in family 48 glycosyl hydrolases

The catalytic base in family 48 glycosyl hydrolases has not been previously established experimentally. Based on structural and modeling data published to date, we used site-directed mutagenesis and azide rescue activity assays to show definitively that the catalytic base in Thermobifida fusca Cel48A is aspartic acid 225. Of the tested mutants, only Cel48A with the D225E mutation retained partial activity on soluble and insoluble substrates. In azide rescue experiments, only the D225G mutation, in the smallest residue tested, showed an increase in activity with added azide.

Catalytic mechanism of cellulose degradation by a cellobiohydrolase, CelS

The hydrolysis of cellulose is the bottleneck in cellulosic ethanol production. The cellobiohydrolase CelS from Clostridium thermocellum catalyzes the hydrolysis of cello-oligosaccharides via inversion of the anomeric carbon. Here, to examine key features of the CelS-catalyzed reaction, QM/MM (SCCDFTB/MM) simulations are performed. The calculated free energy profile for the reaction possesses a 19 kcal/mol barrier. The results confirm the role of active site residue Glu87 as the general acid catalyst in the cleavage reaction and show that Asp255 may act as the general base. A feasible position in the reactant state of the water molecule responsible for nucleophilic attack is identified. Sugar ring distortion as the reaction progresses is quantified. The results provide a computational approach that may complement the experimental design of more efficient enzymes for biofuel production.

A highly accurate statistical approach for the prediction of transmembrane beta-barrels

MOTIVATION

Transmembrane beta-barrels (TMBBs) belong to a special structural class of proteins predominately found in the outer membranes of Gram-negative bacteria, mitochondria and chloroplasts. TMBBs are surface-exposed proteins that perform a variety of functions ranging from nutrient acquisition to osmotic regulation. These properties suggest that TMBBs have great potential for use in vaccine or drug therapy development. However, membrane proteins, such as TMBBs, are notoriously difficult to identify and characterize using traditional experimental approaches and current prediction methods are still unreliable.

RESULTS

A prediction method based on the physicochemical properties of experimentally characterized TMBB structures was developed to predict TMBB-encoding genes from genomic databases. The Freeman-Wimley prediction algorithm developed in this study has an accuracy of 99% and MCC of 0.748 when using the most efficient prediction criteria, which is better than any previously published algorithm.

AVAILABILITY

The MS Windows-compatible application is available for download at http://www.tulane.edu/~biochem/WW/apps.html.

Crystal structure of glycoside hydrolase family 55 {beta}-1,3-glucanase from the basidiomycete Phanerochaete chrysosporium

Glycoside hydrolase family 55 consists of beta-1,3-glucanases mainly from filamentous fungi. A beta-1,3-glucanase (Lam55A) from the Basidiomycete Phanerochaete chrysosporium hydrolyzes beta-1,3-glucans in the exo-mode with inversion of anomeric configuration and produces gentiobiose in addition to glucose from beta-1,3/1,6-glucans. Here we report the crystal structure of Lam55A, establishing the three-dimensional structure of a member of glycoside hydrolase 55 for the first time. Lam55A has two beta-helical domains in a single polypeptide chain. These two domains are separated by a long linker region but are positioned side by side, and the overall structure resembles a rib cage. In the complex, a gluconolactone molecule is bound at the bottom of a pocket between the two beta-helical domains. Based on the position of the gluconolactone molecule, Glu-633 appears to be the catalytic acid, whereas the catalytic base residue could not be identified. The substrate binding pocket appears to be able to accept a gentiobiose unit near the cleavage site, and a long cleft runs from the pocket, in accordance with the activity of this enzyme toward various beta-1,3-glucan oligosaccharides. In conclusion, we provide important features of the substrate-binding site at the interface of the two beta-helical domains, demonstrating an unexpected variety of carbohydrate binding modes.

X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution

Crystallography supplies unparalleled detail on structural information critical for mechanistic analyses; however, it is restricted to describing low energy conformations of macromolecules within crystal lattices. Small angle X-ray scattering (SAXS) offers complementary information about macromolecular folding, unfolding, aggregation, extended conformations, flexibly linked domains, shape, conformation, and assembly state in solution, albeit at the lower resolution range of about 50 A to 10 A resolution, but without the size limitations inherent in NMR and electron microscopy studies. Together these techniques can allow multi-scale modeling to create complete and accurate images of macromolecules for modeling allosteric mechanisms, supramolecular complexes, and dynamic molecular machines acting in diverse processes ranging from eukaryotic DNA replication, recombination and repair to microbial membrane secretion and assembly systems. This review addresses both theoretical and practical concepts, concerns and considerations for using these techniques in conjunction with computational methods to productively combine solution scattering data with high-resolution structures. Detailed aspects of SAXS experimental results are considered with a focus on data interpretation tools suitable to model protein and nucleic acid macromolecular structures, including membrane protein, RNA, DNA, and protein-nucleic acid complexes. The methods discussed provide the basis to examine molecular interactions in solution and to study macromolecular flexibility and conformational changes that have become increasingly relevant for accurate understanding, simulation, and prediction of mechanisms in structural cell biology and nanotechnology.

Structures of mutants of cellulase Cel48F of Clostridium cellulolyticum in complex with long hemithiocellooligosaccharides give rise to a new view of the substrate pathway during processive action

An efficient breakdown of lignocellulosic biomass is a prerequisite for the production of second-generation biofuels. Cellulases are key enzymes in this process. We crystallized complexes between hemithio-cello-deca and dodecaoses and the inactive mutants E44Q and E55Q of the endo-processive cellulase Cel48F, one of the most abundant cellulases in cellulosomes from Clostridium cellulolyticum, to elucidate its processive mechanism. In both complexes, the cellooligosaccharides occupy similar positions in the tunnel part of the active site but are more or less buried into the cleft, which hosts the active site. In the E44Q complex, it proceeds along the upper part of the cavity, while it occupies in the E55Q complex the same productive binding subsites in the lower part of the cavity that have previously been reported in Cel48F/cellooligosaccharide complexes. In both cases, the sugar moieties are stabilized by stacking interactions with aromatic side chains and H bonds. The upper pathway is gated by Tyr403, which blocks its access in the E55Q complex and offers a new stacking interaction in the E44Q complex. The new structural data give rise to the hypothesis of a two-step mechanism in which processive action and chain disruption occupy different subsites at the end of their trajectory. In the first part of the mechanism, the chain may smoothly slide up to the leaving group site along the upper pathway, while in the second part, the chain is cleaved in the already described productive binding position located in the lower pathway. The solved native structure of Cel48F without any bound sugar in the active site confirms the two side-chain orientations of the proton donor Glu55 as observed in the complex structures.

Structural Basis of the Catalytic Reaction Mechanism of Novel 1,2-α-L-Fucosidase from Bifidobacterium bifidum

1,2-alpha-L-fucosidase (AfcA), which hydrolyzes the glycosidic linkage of Fucalpha1-2Gal via an inverting mechanism, was recently isolated from Bifidobacterium bifidum and classified as the first member of the novel glycoside hydrolase family 95. To better understand the molecular mechanism of this enzyme, we determined the x-ray crystal structures of the AfcA catalytic (Fuc) domain in unliganded and complexed forms with deoxyfuconojirimycin (inhibitor), 2'-fucosyllactose (substrate), and L-fucose and lactose (products) at 1.12-2.10 A resolution. The AfcA Fuc domain is composed of four regions, an N-terminal beta region, a helical linker, an (alpha/alpha)6 helical barrel domain, and a C-terminal beta region, and this arrangement is similar to bacterial phosphorylases. In the complex structures, the ligands were buried in the central cavity of the helical barrel domain. Structural analyses in combination with mutational experiments revealed that the highly conserved Glu566 probably acts as a general acid catalyst. However, no carboxylic acid residue is found at the appropriate position for a general base catalyst. Instead, a water molecule stabilized by Asn423 in the substrate-bound complex is suitably located to perform a nucleophilic attack on the C1 atom of L-fucose moiety in 2'-fucosyllactose, and its location is nearly identical near the O1 atom of beta-L-fucose in the products-bound complex. Based on these data, we propose and discuss a novel catalytic reaction mechanism of AfcA.

Unravelling carbon metabolism in anaerobic cellulolytic bacteria

Carbon metabolism in anaerobic cellulolytic bacteria has been investigated essentially in Clostridium thermocellum, Clostridium cellulolyticum, Fibrobacter succinogenes, Ruminococcus flavefaciens, and Ruminococcus albus. While cellulose depolymerization into soluble sugars by various cellulases is undoubtedly the first step in bacterial metabolisation of cellulose, it is not the only one to consider. Among anaerobic cellulolytic bacteria, C. cellulolyticum has been investigated metabolically the most in the past few years. Summarizing metabolic flux analyses in continuous culture using either cellobiose (a soluble cellodextrin resulting from cellulose hydrolysis) or cellulose (an insoluble biopolymer), this review aims to stress the importance of the insoluble nature of a carbon source on bacterial metabolism. Furthermore, some general and specific traits of anaerobic cellulolytic bacteria trends, namely, the importance and benefits of (i) cellodextrins with degree of polymerization higher than 2, (ii) intracellular phosphorolytic cleavage, (iii) glycogen cycling on cell bioenergetics, and (iv) carbon overflows in regulation of carbon metabolism, as well as detrimental effects of (i) soluble sugars and (ii) acidic environment on bacterial growth. Future directions for improving bacterial cellulose degradation are discussed.

The Reaction of Phosphohexomutase from Pseudomonas aeruginosa

The enzyme phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa catalyzes the reversible conversion of 1-phospho to 6-phospho-sugars. The reaction entails two phosphoryl transfers, with an intervening 180 degrees reorientation of the reaction intermediate (e.g. glucose 1,6-bisphosphate) during catalysis. Reorientation of the intermediate occurs without dissociation from the active site of the enzyme and is, thus, a simple example of processivity, as defined by multiple rounds of catalysis without release of substrate. Structural characterization of two PMM/PGM-intermediate complexes with glucose 1,6-bisphosphate provides new insights into the reaction catalyzed by the enzyme, including the reorientation of the intermediate. Kinetic analyses of site-directed mutants prompted by the structural studies reveal active site residues critical for maintaining association with glucose 1,6-bisphosphate during its unique dynamic reorientation in the active site of PMM/PGM.

Quantum mechanics studies of cellobiose conformations

Three regions of the conformation space that describes the relative orientations of the two glucose residues of cellobiose were analyzed with quantum mechanics. A central region, in which most crystal structures are found, was covered by a 9 × 9 grid of 20° increments of the linkage torsion angles ϕ and ψ. Besides these 81 constrained minimizations, we studied two central subregions and two regions at the edges of our maps of complete ϕ,ψ space with unconstrained minimization, for a total of 85 target geometries. HF/6-31G(d) and single-point HF/6-311+G(d) calculations were used to find the lowest energies for each geometry. B3LYP/6-31G+G(d) and single point B3LYP/6-11+G(d) calculations were also used for all unconstrained minimizations. For each target, 181 starting geometries were tried (155 for the unconstrained targets). Numerous different starting geometries resulted in the lowest energies for the various target structures. The starting geometries came from five different sets that were based on molecular mechanics energies. Although all five sets contributed to the adiabatic map, use of any single set resulted in discrepancies of 3–7 kcal/mol (1 cal = 4.184 J) with the final map. For most of the targets, the starting geometry that gave the lowest energy depended on the basis set and whether the HF or B3LYP method was used. However, each of the above four calculations gave the same overall lowest energy structure that was found previously by Strati et al. This global minimum, stabilized by highly cooperative hydrogen bonds, is in a region that is essentially not populated by crystal structures. HF/6-31G(d) energy contours of the mapped central region were compatible with the observed crystal structures. Observed structures that lacked O3···O5′ hydrogen bonds were about 1 kcal/mol above the map's minimum, and observed structures that have a pseudo twofold screw axis ranged from about 0.4 to 1.0 kcal/mol. The HF/6-311+G(d) map accommodated the observed structures nearly as well.Key words: cellulose, carbohydrate, conformation, energy, flexibility, folding, helix, shape.

Cellulase, clostridia, and ethanol

Biomass conversion to ethanol as a liquid fuel by the thermophilic and anaerobic clostridia offers a potential partial solution to the problem of the world's dependence on petroleum for energy. Coculture of a cellulolytic strain and a saccharolytic strain of Clostridium on agricultural resources, as well as on urban and industrial cellulosic wastes, is a promising approach to an alternate energy source from an economic viewpoint. This review discusses the need for such a process, the cellulases of clostridia, their presence in extracellular complexes or organelles (the cellulosomes), the binding of the cellulosomes to cellulose and to the cell surface, cellulase genetics, regulation of their synthesis, cocultures, ethanol tolerance, and metabolic pathway engineering for maximizing ethanol yield.