Version 9 of the Amber simulation programs includes a new semi-empirical hybrid QM/MM functionality. This includes support for implicit solvent (generalized Born) and for periodic explicit solvent simulations using a newly developed QM/MM implementation of the particle mesh Ewald (PME) method. The code provides sufficiently accurate gradients to run constant energy QM/MM MD simulations for many nanoseconds. The link atom approach used for treating the QM/MM boundary shows improved performance, and the user interface has been rewritten to bring the format into line with classical MD simulations. Support is provided for the PM3, PDDG/PM3, PM3CARB1, AM1, MNDO, and PDDG/MNDO semi-empirical Hamiltonians as well as the self-consistent charge density functional tight binding (SCC-DFTB) method. Performance has been improved to the point where using QM/MM, for a QM system of 71 atoms within an explicitly solvated protein using periodic boundaries and PME requires less than twice the cpu time of the corresponding classical simulation.

Printed with a renewable-source ink on paper containing at least 50% wastepaper, including 10% post consumer waste. Cellulase C. bescii CelA, a highly active and stable enzyme, exhibits a new cellulose digestion paradigm promoting inter-cellulase synergy. C. bescii CelA, a hydrolytic enzyme with multiple functional domains, may have several advantages over other fungal and bacterial cellulases for use in biofuels production: very high specific activity, stability at elevated temperatures , and a novel digestion mechanism.

Glycoside hydrolases (GHs) cleave glycosidic linkages in carbohydrates, typically via inverting or retaining mechanisms, the latter of which proceeds via a two-step mechanism that includes formation of a glycosyl-enzyme intermediate. We present two new structures of the catalytic domain of Hypocrea jecorina GH Family 7 cellobiohydrolase Cel7A, namely a Michaelis complex with a full cellononaose ligand and a glycosyl-enzyme intermediate, that reveal details of the 'static' reaction coordinate. We also employ transition path sampling to determine the 'dynamic' reaction coordinate for the catalytic cycle. The glycosylation reaction coordinate contains components of forming and breaking bonds and a conformational change in the nucleophile. Deglycosylation proceeds via a product-assisted mechanism wherein the glycosylation product, cellobiose, positions a water molecule for nucleophilic attack on the anomeric carbon of the glycosyl-enzyme intermediate. In concert with previous structures, the present results reveal the complete hydrolytic reaction coordinate for this naturally and industrially important enzyme family.

Cellulose is a crystalline polymer of β1,4-D-glucose that is difficult to deconstruct to sugars by enzymes. The recalcitrance of cellulose microfibrils is a function of both the shape of cellulose microfibrils and the intrinsic work required to decrystallize individual chains, the latter of which is calculated here from the surfaces of four crystalline cellulose polymorphs: cellulose Iβ, cellulose Iα, cellulose II, and cellulose III(I). For edge chains, the order of decrystallization work is as follows (from highest to lowest): Iβ, Iα, ΙΙΙ(Ι), and II. For cellulose Iβ, we compare chains from three different locations on the surface and find that an increasing number of intralayer hydrogen bonds (from 0 to 2) increases the intrinsic decrystallization work. From these results, we propose a microkinetic model for the deconstruction of cellulose (and chitin) by processive enzymes, which when taken with a previous study [Horn et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 18089] identifies the thermodynamic and kinetic attributes of enzyme and substrate engineering for enhanced cellulose (or chitin) conversion. Overall, this study provides new insights into the molecular interactions that form the structural basis of cellulose, which is the primary building block of plant cell walls, and highlights the need for experimentally determining microfibril shape at the nanometer length scale when comparing conversion rates of cellulose polymorphs by enzymes.

Background: Lytic polysaccharide monooxygenases (LPMOs) represent a recently discovered enzymatic route to cleave carbohydrates. Results:We report the first basidiomycete LPMO structure and describe enzyme-cellulose interactions with simulation. Conclusion:We characterize the copper-containing active site and identify loops important for substrate recognition and binding. Significance: This structure is the first LPMO from a model basidiomycete fungus that contains many LPMO genes.

Microbial aromatic catabolism offers a promising approach to convert lignin, a vast source of renewable carbon, into useful products. Aryl-O-demethylation is an essential biochemical reaction to ultimately catabolize coniferyl and sinapyl lignin-derived aromatic compounds, and is often a key bottleneck for both native and engineered bioconversion pathways. Here, we report the comprehensive characterization of a promiscuous P450 aryl-O-demethylase, consisting of a cytochrome P450 protein from the family CYP255A (GcoA) and a three-domain reductase (GcoB) that together represent a new two-component P450 class. Though originally described as converting guaiacol to catechol, we show that this system efficiently demethylates both guaiacol and an unexpectedly wide variety of lignin-relevant monomers. Structural, biochemical, and computational studies of this novel two-component system elucidate the mechanism of its broad substrate specificity, presenting it as a new tool for a critical step in biological lignin conversion.

Plant cell-wall polysaccharides represent a vast source of food in nature. To depolymerize polysaccharides to soluble sugars, many organisms use multifunctional enzyme mixtures consisting of glycoside hydrolases, lytic polysaccharide mono-oxygenases, polysaccharide lyases, and carbohydrate esterases, as well as accessory, redox-active enzymes for lignin depolymerization. Many of these enzymes that degrade lignocellulose are multimodular with carbohydrate-binding modules (CBMs) and catalytic domains connected by flexible, glycosylated linkers. These linkers have long been thought to simply serve as a tether between structured domains or to act in an inchworm-like fashion during catalytic action. To examine linker function, we performed molecular dynamics (MD) simulations of the Trichoderma reesei Family 6 and Family 7 cellobiohydrolases (TrCel6A and TrCel7A, respectively) bound to cellulose. During these simulations, the glycosylated linkers bind directly to cellulose, suggesting a previously unknown role in enzyme action. The prediction from the MD simulations was examined experimentally by measuring the binding affinity of the Cel7A CBM and the natively glycosylated Cel7A CBM-linker. On crystalline cellulose, the glycosylated linker enhances the binding affinity over the CBM alone by an order of magnitude. The MD simulations before and after binding of the linker also suggest that the bound linker may affect enzyme action due to significant damping in the enzyme fluctuations. Together, these results suggest that glycosylated linkers in carbohydrate-active enzymes, which are intrinsically disordered proteins in solution, aid in dynamic binding during the enzymatic deconstruction of plant cell walls. biofuels | cellulase | post-translational modification | carbohydrate recognition

Translocation of carbohydrate polymers through protein tunnels and clefts is a ubiquitous biochemical phenomenon in proteins such as polysaccharide synthases, glycoside hydrolases, and carbohydrate-binding modules. Although static snapshots of carbohydrate polymer binding in proteins have long been studied via crystallography and spectroscopy, the molecular details of polysaccharide chain processivity have not been elucidated. Here, we employ simulation to examine how a cellulose chain translocates by a disaccharide unit during the processive cycle of a glycoside hydrolase family 7 cellobiohydrolase. Our results demonstrate that these biologically and industrially important enzymes employ a two-step mechanism for chain threading to form a Michaelis complex and that the free energy barrier to chain threading is significantly lower than the hydrolysis barrier. Taken with previous studies, our findings suggest that the rate-limiting step in enzymatic cellulose degradation is the glycosylation reaction, not chain processivity. Based on the simulations, we find that strong electrostatic interactions with polar residues that are conserved in GH7 cellobiohydrolases, but not in GH7 endoglucanases, at the leading glucosyl ring provide the thermodynamic driving force for polysaccharide chain translocation. Also, we consider the role of aromatic-carbohydrate interactions, which are widespread in carbohydrate-active enzymes and have long been associated with processivity. Our analysis suggests that the primary role for these aromatic residues is to provide tunnel shape and guide the carbohydrate chain to the active site. More broadly, this work elucidates the role of common protein motifs found in carbohydrate-active enzymes that synthesize or depolymerize polysaccharides by chain translocation mechanisms coupled to catalysis.