Search

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

Family 45 glycoside hydrolases (GH45) are endoglucanases that are integral to cellulolytic secretomes, and their ability to break down cellulose has been successfully exploited in textile and detergent industries. In addition to their industrial relevance, understanding the molecular mechanism of GH45-catalyzed hydrolysis is of fundamental importance because of their structural similarity to cell wall–modifying enzymes such as bacterial lytic transglycosylases (LTs) and expansins present in bacteria, plants, and fungi. Our understanding of the catalytic itinerary of GH45s has been incomplete because a crystal structure with substrate spanning the −1 to +1 subsites is currently lacking. Here we constructed and validated a putative Michaelis complex in silico and used it to elucidate the hydrolytic mechanism in a GH45, Cel45A from the fungus Humicola insolens, via unbiased simulation approaches. These molecular simulations revealed that the solvent-exposed active-site architecture results in lack of coordination for the hydroxymethyl group of the substrate at the −1 subsite. This lack of coordination imparted mobility to the hydroxymethyl group and enabled a crucial hydrogen bond with the catalytic acid during and after the reaction. This suggests the possibility of a nonhydrolytic reaction mechanism when the catalytic base aspartic acid is missing, as is the case in some LTs (murein transglycosylase A) and expansins. We calculated reaction free energies and demonstrate the thermodynamic feasibility of the hydrolytic and nonhydrolytic reaction mechanisms. Our results provide molecular insights into the hydrolysis mechanism in HiCel45A, with possible implications for elucidating the elusive catalytic mechanism in LTs and expansins.

Significant effort has been expended toward the discovery and/or engineering of improved cellulases. An alternative to this approach is utilizing multifunctional enzymes; however, essentially most if not all relevant bacterial enzymes of this type do not express well in fungi. Therefore, developing a systematic understanding of how to construct multifunctional enzymes that are expressible in commercial fungal hosts is crucial to developing next generation enzymes for biomass deconstruction. Multifunctional cellulolytic enzymes, such as CelA from Caldicellulosiruptor bescii, show extremely high cellulolytic activity; however, a systematic understanding of its mechanism of action does not exist and it is not readily expressed in current industrial hosts. CelA is comprised of GH 9 and GH 48 catalytic domains connected by three type III cellulose-binding modules (CBMs). We have engineered several multifunctional enzymes designed to mimic CelA, and successfully expressed them in T. reesei. We then assessed their biophysical and kinetic performance parameters. The CBM3b-containing construct demonstrated increased initial binding rate to cellulose, enhanced digestion of biomass, and was able to increase the activity of a commercial cellulase formulation acting on pretreated biomass. The same construct containing CBM1 also demonstrated enhancement of a commercial cellulase formulation in the digestion of pretreated biomass; however, it had lower initial binding rates and did not demonstrate improved activity on its own. Interestingly, a CBM-less construct decreased the activity of the commercial cellulase formulation slightly. Examination of the mechanism of the CBM3b-containing construct revealed a novel biomass deconstruction behavior similar to, but yet distinct from that of native CelA.

Cellulase enzymes deconstruct recalcitrant cellulose into soluble sugars, making them a biocatalyst of biotechnological interest for use in the nascent lignocellulosic bioeconomy. Cellobiohydrolases (CBHs) are cellulases capable of liberating many sugar molecules in a processive manner without dissociating from the substrate. Within the complete processive cycle of CBHs, dissociation from the cellulose substrate is rate limiting, but the molecular mechanism of this step is unknown. Here, we present a direct comparison of potential molecular mechanisms for dissociation via Hamiltonian replica exchange molecular dynamics of the model fungal CBH, Trichoderma reesei Cel7A. Computational rate estimates indicate that stepwise cellulose dethreading from the binding tunnel is 4 orders of magnitude faster than a clamshell mechanism, in which the substrate-enclosing loops open and release the substrate without reversing. We also present the crystal structure of a disulfide variant that covalently links substrate-enclosing loops on either side of the substrate-binding tunnel, which constitutes a CBH that can only dissociate via stepwise dethreading. Biochemical measurements indicate that this variant has a dissociation rate constant essentially equivalent to the wild type, implying that dethreading is likely the predominant mechanism for dissociation.

In the terrestrial biosphere, biomass deconstruction is conducted by microbes employing a variety of complementary strategies, many of which remain to be discovered. Moreover, the biofuels industry seeks more efficient (and less costly) cellulase formulations upon which to launch the nascent sustainable bioenergy economy. The glycan decoration of fungal cellulases has been shown to protect these enzymes from protease action and to enhance binding to cellulose. We show here that thermal tolerant bacterial cellulases are glycosylated as well, although the types and extents of decoration differ from their Eukaryotic counterparts. Our major findings are that glycosylation of CelA is uniform across its three linker peptides and composed of mainly galactose disaccharides (which is unique) and that this glycosylation dramatically impacts the hydrolysis of insoluble substrates, proteolytic and thermal stability, and substrate binding and changes the dynamics of the enzyme. This study suggests that the glycosylation of CelA is crucial for its exceptionally high cellulolytic activity on biomass and provides the robustness needed for this enzyme to function in harsh environments including industrial settings.

Biological funneling of lignin-derived aromatic compounds is a promising approach for valorizing its catalytic depolymerization products. Industrial processes for aromatic bioconversion will require efficient enzymes for key reactions, including demethylation of O-methoxy-aryl groups, an essential and often rate-limiting step. The recently characterized GcoAB cytochrome P450 system comprises a coupled monoxygenase (GcoA) and reductase (GcoB) that catalyzes oxidative demethylation of the O-methoxy-aryl group in guaiacol. Here, we evaluate a series of engineered GcoA variants for their ability to demethylate o-and p-vanillin, which are abundant lignin depolymerization products. Two rationally designed, single amino acid substitutions, F169S and T296S, are required to convert GcoA into an efficient catalyst toward the o- and p-isomers of vanillin, respectively. Gain-of-function in each case is explained in light of an extensive series of enzyme-ligand structures, kinetic data, and molecular dynamics simulations. Using strains of Pseudomonas putida KT2440 already optimized for p-vanillin production from ferulate, we demonstrate demethylation by the T296S variant in vivo. This work expands the known aromatic O-demethylation capacity of cytochrome P450 enzymes toward important lignin-derived aromatic monomers.

Significance Deconstruction of recalcitrant polymers, such as cellulose or chitin, is accomplished in nature by synergistic enzyme cocktails that evolved over millions of years. In these systems, soluble dimeric or oligomeric intermediates are typically released via interfacial biocatalysis, and additional enzymes often process the soluble intermediates into monomers for microbial uptake. The recent discovery of a two-enzyme system for polyethylene terephthalate (PET) deconstruction, which employs one enzyme to convert the polymer into soluble intermediates and another enzyme to produce the constituent PET monomers (MHETase), suggests that nature may be evolving similar deconstruction strategies for synthetic plastics. This study on the characterization of the MHETase enzyme and synergy of the two-enzyme PET depolymerization system may inform enzyme cocktail-based strategies for plastics upcycling., Plastics pollution represents a global environmental crisis. In response, microbes are evolving the capacity to utilize synthetic polymers as carbon and energy sources. Recently, Ideonella sakaiensis was reported to secrete a two-enzyme system to deconstruct polyethylene terephthalate (PET) to its constituent monomers. Specifically, the I. sakaiensis PETase depolymerizes PET, liberating soluble products, including mono(2-hydroxyethyl) terephthalate (MHET), which is cleaved to terephthalic acid and ethylene glycol by MHETase. Here, we report a 1.6 Å resolution MHETase structure, illustrating that the MHETase core domain is similar to PETase, capped by a lid domain. Simulations of the catalytic itinerary predict that MHETase follows the canonical two-step serine hydrolase mechanism. Bioinformatics analysis suggests that MHETase evolved from ferulic acid esterases, and two homologous enzymes are shown to exhibit MHET turnover. Analysis of the two homologous enzymes and the MHETase S131G mutant demonstrates the importance of this residue for accommodation of MHET in the active site. We also demonstrate that the MHETase lid is crucial for hydrolysis of MHET and, furthermore, that MHETase does not turnover mono(2-hydroxyethyl)-furanoate or mono(2-hydroxyethyl)-isophthalate. A highly synergistic relationship between PETase and MHETase was observed for the conversion of amorphous PET film to monomers across all nonzero MHETase concentrations tested. Finally, we compare the performance of MHETase:PETase chimeric proteins of varying linker lengths, which all exhibit improved PET and MHET turnover relative to the free enzymes. Together, these results offer insights into the two-enzyme PET depolymerization system and will inform future efforts in the biological deconstruction and upcycling of mixed plastics.

In nature, many microbes secrete mixtures of glycoside hydrolases, oxidoreductases, and accessory enzymes to deconstruct polysaccharides and lignin in plants. These enzymes are often decorated with N- and O-glycosylation, the roles of which have been broadly attributed to protection from proteolysis, as the extracellular milieu is an aggressive environment. Glycosylation has been shown to sometimes affect activity, but these effects are not fully understood. Here, we examine N- and O-glycosylation on a model, multimodular glycoside hydrolase family 7 cellobiohydrolase (Cel7A), which exhibits an O-glycosylated carbohydrate-binding module (CBM) and an O-glycosylated linker connected to an N- and O-glycosylated catalytic domain (CD)-a domain architecture common to many biomass-degrading enzymes. We report consensus maps for Cel7A glycosylation that include glycan sites and motifs. Additionally, we examine the roles of glycans on activity, substrate binding, and thermal and proteolytic stability. N-glycan knockouts on the CD demonstrate that N-glycosylation has little impact on cellulose conversion or binding, but does have major stability impacts. O-glycans on the CBM have little impact on binding, proteolysis, or activity in the whole-enzyme context. However, linker O-glycans greatly impact cellulose conversion via their contribution to proteolysis resistance. Molecular simulations predict an additional role for linker O-glycans, namely that they are responsible for maintaining separation between ordered domains when Cel7A is engaged on cellulose, as models predict α-helix formation and decreased cellulose interaction for the nonglycosylated linker. Overall, this study reveals key roles for N- and O-glycosylation that are likely broadly applicable to other plant cell-wall-degrading enzymes.

In the upgrading of biomass pyrolysis vapors to hydrocarbons, dehydration accomplishes a primary objective of removing oxygen, and acidic zeolites represent promising catalysts for the dehydration reaction. Here, we utilized density functional theory calculations to estimate adsorption energetics and intrinsic kinetics of alcohol dehydration over H-ZSM-5, H-BEA, and H-AEL zeolites. The ONIOM (our Own N-layered Integrated molecular Orbital and molecular Mechanics) calculations of adsorption energies were observed to be inconsistent when benchmarked against QM (Quantum Mechanical)/Hartree–Fock and periodic boundary condition calculations. However, reaction coordinate calculations of adsorbed species and transition states were consistent across all levels considered. Comparison of ethanol, isopropanol (IPA), and tert-amyl alcohol (TAA) over these three zeolites allowed for a detailed examination of how confinement impacts on reaction mechanisms and kinetics. The TAA, seen to proceed via a carbocationic mechanism, was found to have the lowest activation barrier, followed by IPA and then ethanol, both of which dehydrate via a concerted mechanism. Barriers in H-BEA were consistently found to be lower than in H-ZSM-5 and H-AEL, attributed to late transition states and either elevated strain or inaccurately estimating long-range electrostatic interactions in H-AEL, respectively. Molecular dynamics simulations revealed that the diffusivity of these three alcohols in H-ZSM-5 were significantly overestimated by Knudsen diffusion, which will complicate experimental efforts to develop a kinetic model for catalytic fast pyrolysis.

Significance Lignin is an abundant but underutilized heterogeneous polymer found in terrestrial plants. In current lignocellulosic biorefinery paradigms, lignin is primarily slated for incineration, but for a nonfood plant-based bioeconomy to be successful, lignin valorization is critical. An emerging concept to valorize lignin uses aromatic–catabolic pathways and microbes to funnel heterogeneous lignin-derived aromatic compounds to single high-value products. For this approach to be viable, the discovery and engineering of enzymes to conduct key reactions is critical. In this work, we have engineered a two-component cytochrome P450 enzyme system to conduct one of the most important reactions in biological lignin conversion: aromatic O-demethylation of syringol, the base aromatic unit of S-lignin, which is highly abundant in hardwoods and grasses., Microbial conversion of aromatic compounds is an emerging and promising strategy for valorization of the plant biopolymer lignin. A critical and often rate-limiting reaction in aromatic catabolism is O-aryl-demethylation of the abundant aromatic methoxy groups in lignin to form diols, which enables subsequent oxidative aromatic ring-opening. Recently, a cytochrome P450 system, GcoAB, was discovered to demethylate guaiacol (2-methoxyphenol), which can be produced from coniferyl alcohol-derived lignin, to form catechol. However, native GcoAB has minimal ability to demethylate syringol (2,6-dimethoxyphenol), the analogous compound that can be produced from sinapyl alcohol-derived lignin. Despite the abundance of sinapyl alcohol-based lignin in plants, no pathway for syringol catabolism has been reported to date. Here we used structure-guided protein engineering to enable microbial syringol utilization with GcoAB. Specifically, a phenylalanine residue (GcoA-F169) interferes with the binding of syringol in the active site, and on mutation to smaller amino acids, efficient syringol O-demethylation is achieved. Crystallography indicates that syringol adopts a productive binding pose in the variant, which molecular dynamics simulations trace to the elimination of steric clash between the highly flexible side chain of GcoA-F169 and the additional methoxy group of syringol. Finally, we demonstrate in vivo syringol turnover in Pseudomonas putida KT2440 with the GcoA-F169A variant. Taken together, our findings highlight the significant potential and plasticity of cytochrome P450 aromatic O-demethylases in the biological conversion of lignin-derived aromatic compounds.

Biomass pyrolysis offers a promising means to rapidly depolymerize lignocellulosic biomass for subsequent catalytic upgrading to renewable fuels. Substantial efforts are currently ongoing to optimize pyrolysis processes including various fast pyrolysis and catalytic fast pyrolysis schemes. In all cases, complex aqueous streams are generated containing solubilized organic compounds that are not converted to target fuels or chemicals and are often slated for wastewater treatment, in turn creating an economic burden on the biorefinery. Valorization of the species in these aqueous streams, however, offers significant potential for substantially improving the economics and sustainability of thermochemical biorefineries. To that end, here we provide a thorough characterization of the aqueous streams from four pilot-scale pyrolysis processes: namely, from fast pyrolysis, fast pyrolysis with downstream fractionation, in situ catalytic fast pyrolysis, and ex situ catalytic fast pyrolysis. These configurations and ...

Glycoside hydrolase family 7 (GH7) cellobiohydrolases (CBHs) are enzymes commonly employed in plant cell wall degradation across eukaryotic kingdoms of life, as they provide significant hydrolytic potential in cellulose turnover. To date, many fungal GH7 CBHs have been examined, yet many questions regarding structure-activity relationships in these important natural and commercial enzymes remain. Here, we present the crystal structures and a biochemical analysis of two GH7 CBHs from social amoeba: Dictyostelium discoideum Cel7A ( Ddi Cel7A) and Dictyostelium purpureum Cel7A ( Dpu Cel7A). Ddi Cel7A and Dpu Cel7A natively consist of a catalytic domain and do not exhibit a carbohydrate-binding module (CBM). The structures of Ddi Cel7A and Dpu Cel7A, resolved to 2.1 Å and 2.7 Å, respectively, are homologous to those of other GH7 CBHs with an enclosed active-site tunnel. Two primary differences between the Dictyostelium CBHs and the archetypal model GH7 CBH, Trichoderma reesei Cel7A ( Tre Cel7A), occur near the hydrolytic active site and the product-binding sites. To compare the activities of these enzymes with the activity of Tre Cel7A, the family 1 Tre Cel7A CBM and linker were added to the C terminus of each of the Dictyostelium enzymes, creating Ddi Cel7A CBM and Dpu Cel7A CBM , which were recombinantly expressed in T. reesei . Ddi Cel7A CBM and Dpu Cel7A CBM hydrolyzed Avicel, pretreated corn stover, and phosphoric acid-swollen cellulose as efficiently as Tre Cel7A when hydrolysis was compared at their temperature optima. The K i of cellobiose was significantly higher for Ddi Cel7A CBM and Dpu Cel7A CBM than for Tre Cel7A: 205, 130, and 29 μM, respectively. Taken together, the present study highlights the remarkable degree of conservation of the activity of these key natural and industrial enzymes across quite distant phylogenetic trees of life. IMPORTANCE GH7 CBHs are among the most important cellulolytic enzymes both in nature and for emerging industrial applications for cellulose breakdown. Understanding the diversity of these key industrial enzymes is critical to engineering them for higher levels of activity and greater stability. The present work demonstrates that two GH7 CBHs from social amoeba are surprisingly quite similar in structure and activity to the canonical GH7 CBH from the model biomass-degrading fungus T. reesei when tested under equivalent conditions (with added CBM-linker domains) on an industrially relevant substrate.

The processive cycle of the bacterial cellulose synthase (Bcs) includes the addition of a single glucose moiety to the end of a growing cellulose chain followed by the translocation of the nascent chain across the plasma membrane. The mechanism of this translocation and its precise location within the processive cycle are not well understood. In particular, the molecular details of how a polymer (cellulose) whose basic structural unit is a dimer (cellobiose) can be constructed by adding one monomer (glucose) at a time are yet to be elucidated. Here, we have utilized molecular dynamics simulations and free energy calculations to the shed light on these questions. We find that translocation forward by one glucose unit is quite favorable energetically, giving a free energy stabilization of greater than 10 kcal/mol. In addition, there is only a small barrier to translocation, implying that translocation is not rate limiting within the Bcs processive cycle (given experimental rates for cellulose synthesis in vitro). Perhaps most significantly, our results also indicate that steric constraints at the transmembrane tunnel entrance regulate the dimeric structure of cellulose. Namely, when a glucose molecule is added to the cellulose chain in the same orientation as the acceptor glucose, the terminal glucose freely rotates upon forward motion, thus suggesting a regulatory mechanism for the dimeric structure of cellulose. We characterize both the conserved and non-conserved enzyme-polysaccharide interactions that drive translocation, and find that 20 of the 25 residues that strongly interact with the translocating cellulose chain in the simulations are well conserved, mostly with polar or aromatic side chains. Our results also allow for a dynamical analysis of the role of the so-called `finger helix' in cellulose translocation that has been observed structurally. Taken together, these findings aid in the elucidation of the translocation steps of the Bcs processive cycle and may be widely relevant to polysaccharide synthesizing or degrading enzymes that couple catalysis with chain translocation.

Glycoside Hydrolase Family 7 cellobiohydrolases (GH7 CBHs) catalyze cellulose depolymerization in cellulolytic eukaryotes, making them key discovery and engineering targets. However, there remains a lack of robust structure–activity relationships for these industrially important cellulases. Here, we compare CBHs from Trichoderma reesei (TrCel7A) and Penicillium funiculosum (PfCel7A), which exhibit a multi-modular architecture consisting of catalytic domain (CD), carbohydrate-binding module, and linker. We show that PfCel7A exhibits 60% greater performance on biomass than TrCel7A. To understand the contribution of each domain to this improvement, we measure enzymatic activity for a library of CBH chimeras with swapped subdomains, demonstrating that the enhancement is mainly caused by PfCel7A CD. We solve the crystal structure of PfCel7A CD and use this information to create a second library of TrCel7A CD mutants, identifying a TrCel7A double mutant with near-equivalent activity to wild-type PfCel7A. Overall, these results reveal CBH regions that enable targeted activity improvements., Cellobiohydrolases (CBHs) are critical for natural and industrial biomass degradation but their structure–activity relationships are not fully understood. Here, the authors present the biochemical and structural characterization of two CBHs, identifying protein regions that confer enhanced CBH activity.

Nucleation from solution is a ubiquitous phenomenon with relevance to myriad scientific disciplines, including pharmaceuticals, biomineralization, and disease. One prominent example is the nucleation of clathrate hydrates, multicomponent crystalline inclusion compounds relevant to the energy industry where they block pipelines and also constitute a potential vast energy resource. Despite their importance, the molecular mechanism of incipient hydrate formation remains unknown. Herein, we employ advanced molecular simulation tools (pB histogram, equilibrium path sampling) to provide a statistical-mechanical basis for extracting physical insight into the molecular steps by which clathrates form. Through testing the Mutually Coordinated Guest (MCG) order parameter, we demonstrate that both guest (methane) and host (water) structuring are crucial to accurately describe the nucleation of hydrates and determine a critical nucleus size of MCG-1 = 16 at 255 K and 500 bar. Equipped with a validated (and novel) reaction coordinate, subsequent equilibrium path sampling simulations yield the free energy barrier and nucleation rate. The resulting quantitative nucleation process is described by the MCG clustering mechanism. This constitutes a significant advance in the field of hydrates research, as the fitness of a molecular descriptor has never been statistically verified. More broadly, this work has significance to a wide range of multicomponent nucleation contexts wherein the formation mechanism depends on contributions from both solute and solvent.

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.

Clathrate hydrates are solid crystalline structures most commonly formed from solutions that have nucleated to form a mixed solid composed of water and gas. Understanding the mechanism of clathrate hydrate nucleation is essential to grasp the fundamental chemistry of these complex structures and their applications. Molecular dynamics (MD) simulation is an ideal method to study nucleation at the molecular level because the size of the critical nucleus and formation rate occur on the nano scale. Various analysis methods for nucleation have been developed through MD to analyze nucleation. In particular, the mean first-passage time (MFPT) and survival probability (SP) methods have proven to be effective in procuring the nucleation rate and critical nucleus size for monatomic systems. This study assesses the MFPT and SP methods, previously used for monatomic systems, when applied to analyzing clathrate hydrate nucleation. Because clathrate hydrate nucleation is relatively difficult to observe in MD simulations (due to its high free energy barrier), these methods have yet to be applied to clathrate hydrate systems. In this study, we have analyzed the nucleation rate and critical nucleus size of methane hydrate using MFPT and SP methods from data generated by MD simulations at 255 K and 50 MPa. MFPT was modified for clathrate hydrate from the original version by adding the maximum likelihood estimate and growth effect term. The nucleation rates calculated by MFPT and SP methods are within 5%, and the critical nucleus size estimated by the MFPT method was 50% higher, than values obtained through other more rigorous but computationally expensive estimates. These methods can also be extended to the analysis of other clathrate hydrates.

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.

Polysaccharide depolymerization in nature is primarily accomplished by processive glycoside hydrolases (GHs), which abstract single carbohydrate chains from polymer crystals and cleave glycosidic linkages without dissociating after each catalytic event. Understanding the molecular-level features and structural aspects of processivity is of importance due to the prevalence of processive GHs in biomass-degrading enzyme cocktails. Here, we describe recent advances towards the development of a molecular-level theory of processivity for cellulolytic and chitinolytic enzymes, including the development of novel methods for measuring rates of key steps in processive action and insights gained from structural and computational studies. Overall, we present a framework for developing structure-function relationships in processive GHs and outline additional progress towards developing a fundamental understanding of these industrially important enzymes.

Methane hydrates are ice-like inclusion compounds with importance to the oil and natural gas industry, global climate change, and gas transportation and storage. The molecular mechanism by which these compounds form under conditions relevant to industry and nature remains mysterious. To understand the mechanism of methane hydrate nucleation from supersaturated aqueous solutions, we performed simulations at controlled and realistic supersaturation. We found that critical nuclei are extremely large and that homogeneous nucleation rates are extremely low. Our findings suggest that nucleation of methane hydrates under these realistic conditions cannot occur by a homogeneous mechanism.

Recent experiments have demonstrated that intense, nanosecond laser pulses can induce crystal nucleation from supersaturated solutions that are transparent at the incident wavelengths, a phenomenon termed nonphotochemical laser-induced nucleation (NPLIN). Previous work has proposed that this effect is due to the alignment of solute molecules in solution due to the electric field of the applied laser light, promoting crystalline order. We have used simulations of NPLIN to examine how an orientational bias in solution affects nucleation with Monte Carlo simulations of a Potts lattice gas model. We examine this effect within both a classical, one-step nucleation framework as well as in the context of two-step nucleation. Our results indicate that an orientational bias can reduce the free energy barrier to nucleation within the one-step picture as well as promote the crystallization of amorphous precritical nuclei (the rate-determining step in the two-step picture). However, these effects are only present with field strengths that are much greater than those used in experiments.

Short intense laser pulses of visible and infrared light can dramatically accelerate crystal nucleation from transparent solutions; previous studies invoke mechanisms that are only applicable for nucleation of ordered phases or high dielectric phases. However, we show that similar laser pulses induce CO(2) bubble nucleation in carbonated water. Additionally, in water that is cosupersaturated with argon and glycine, argon bubbles escaping from the water can induce crystal nucleation without a laser. Our findings suggest a possible link between laser-induced nucleation of bubbles and crystals.

In the high friction limit of Kramers' theory, the diffusion coefficient for motion along the reaction coordinate is a crucial parameter in determining reaction rates from mean first passage times. The Einstein relation between mean squared displacement, time, and diffusivity is inaccurate at short times because of ballistic motion and inaccurate at long times because trajectories drift away from maxima in the potential of mean force. Starting from the Smoluchowski equation for a downward parabolic barrier, we show how drift induced by the potential of mean force can be included in estimating the diffusivity. A modified relation between mean squared displacement, time, and diffusivity now also includes a dependence on the barrier curvature. The new relation provides the diffusivity at the top of the barrier from a linear regression that is analogous to the procedure commonly used with Einstein's relation. The new approach has particular advantages over previous approaches when evaluations of the reaction coordinate are costly or when the reaction coordinate cannot be differentiated to compute restraining forces or velocities. We use the new method to study the dynamics of barrier crossing in a Potts lattice gas model of nucleation from solution. Our analysis shows that some current hypotheses about laser-induced nucleation mechanisms lead to a nonzero threshold laser pulse duration below which a laser pulse will not affect nucleation. We therefore propose experiments that might be used to test these hypotheses.