A monte carlo simulation of the depolymerization of linear homopolymers by endo-enzymes exhibiting random-attack probability and single-attack mechanism: application to the (1-->3), (1-->4)-beta-D-glucan/endo-(1-->3),(1-->4)-beta-D-glucanase system

Potential of time-stepping stochastic models as tools for guiding the design and operation of processes for the enzymatic hydrolysis of polysaccharides - A review

Processes for the enzymatic hydrolysis of polysaccharides in biorefineries are becoming increasingly important. The complex network of reactions involved in polysaccharide hydrolysis can be described by stochastic models that advance in steps of time. Such models have the potential to be important tools for guiding process design and operation, and several have been developed over the last two decades. We evaluate these models. Many of the current stochastic models for the hydrolysis of colloidal polysaccharides use empirical parameters that have no recognized biological meaning. Only one model uses classical parameters of enzyme kinetics, namely specificity constants and saturation constants. Recent stochastic models for the hydrolysis of insoluble cellulose give valuable insights into the molecular-level phenomenon that limit hydrolysis rates. We conclude that, if stochastic models of enzymatic polysaccharide hydrolysis are to become widely used tools for guiding process development, then further improvements are required.

Development and application of a model for chitosan hydrolysis by a family 18 chitinase

Hydrolysis of partially deacetylated chitosans by ChitinaseB (ChiBeta) from Serratia marcescens results in mixtures of oligosaccharides typically between 2 and 20 sugar residues. The amounts of different oligomer fractions depend on the degree of acetylation of the starting chitosans. We have used experimentally determined distributions of hydrolysis products to develop a model for chitosan hydrolysis by ChiB. Important elements of the model include interaction parameters for acetylated/deacetylated units in each of the six subsites in the active cleft and degree of processivity (multiple attack). The hydrolysis reaction is described as a chemical reaction with an activation barrier that depends on the substrate sequence presented to the enzyme subsites. Using a Monte Carlo approach, the interaction parameters were refined by minimizing the difference between observed and predicted amounts of hydrolysis products obtained upon degradation of chitosan with a degree of acetylation of 65%. The final model can accurately predict complex patterns of oligosaccharides produced in the hydrolysis of chitosans with various degrees of acetylation, as well as patterns observed during reactions with chito-oligosaccharides. The behavior of a ChiB mutant with a mutation in subsite +2 (Gly188Asp), which reduces the affinity for an acetylated sugar, could be predicted correctly by introducing one single change in the model parameters (the interaction energy for an acetylated unit in the +2 subsite). The proposed model may be used to explore degradation products for different enzyme-substrates combinations and to optimize conditions for preparation of specific oligosaccharides. In addition, the model provides insight into subsite interaction parameters and the degree of processivity, which complements previous experimental studies on the mode of action of ChiB.

Mode of action of recombinant Azotobacter vinelandii mannuronan C-5 epimerases AlgE2 and AlgE4

The enzymes mannuronan C-5 epimerases catalyze conversion of beta-D-mannuronic acid to alpha-L-guluronic acid in alginates at the polymer level and thereby introduce sequences that have functional properties relevant to gelation. The enzymatic conversion by recombinant mannuronan C-5 epimerases AlgE4 and AlgE2 on alginate type substrates with different degree of polymerization and initial low fraction of alpha-L-guluronic acid was investigated. Essentially no enzymatic activity was found for fractionated mannuronan oligomer substrates with an average degree of polymerization, DP(n), less than or equal 6, whereas increasing the DP(n) yielded increased epimerization activity. This indicates that these enzymes have an active site consisting of binding domains for consecutive residues that requires interaction with 7 or more consecutive residues to show enzymatic activity. The experimentally determined kinetics of the reaction, and the residue sequence arrangement introduced by the epimerization, were modeled using Monte Carlo simulation accounting for the various competing intrachain substrates and assuming either a processive mode of action or preferred attack. The comparison between experimental data and simulation results suggests that epimerization by AlgE4 is best described by a processive mode of action, whereas the mode of action of AlgE2 appears to be more difficult to determine.

Iteration model of starch hydrolysis by amylolytic enzymes

An elaborate computer program to simulate the process of starch hydrolysis by amylolytic enzymes was been developed. It is based on the Monte Carlo method and iteration kinetic model, which predict productive and non-productive amylase complexes with substrates. It describes both multienzymatic and multisubstrate reactions simulating the "real" concentrations of all components versus the time of the depolymerization reaction the number of substrates, intermediate products, and final products are limited only by computer memory. In this work, it is assumed that the "proper" substrate for amylases is the glucoside linkages in starch molecules. Dynamic changes of substrate during the simulation adequately influence the increase or decrease of reaction velocity, as well as the kinetics of depolymerization. The presented kinetic model, can be adapted to describe most enzymatic degradations of a polymer. This computer program has been tested on experimental data obtained for alpha- and beta-amylases.