Stochastic model of lignocellulosic material saccharification

The processing of agricultural wastes towards extraction of renewable resources is recently being considered as a promising alternative to conventional biofuel production. The degradation of agricultural residues is a complex chemical process that is currently time intensive and costly. Various pre-treatment methods are being investigated to determine the subsequent modification of the material and the main obstacles in increasing the enzymatic saccharification. In this study, we present a computational model that complements the experimental approaches. We decipher how the three-dimensional structure of the substrate impacts the saccharification dynamics. We model a cell wall microfibril composed of cellulose and surrounded by hemicellulose and lignin, with various relative abundances and arrangements. This substrate is subjected to digestion by different cocktails of well characterized enzymes. The saccharification dynamics is simulated in silico using a stochastic procedure based on a Gillespie algorithm. As we additionally implement a fitting procedure that optimizes the parameters of the simulation runs, we are able to reproduce experimental saccharification time courses for corn stover. Our model highlights the synergistic action of enzymes, and confirms the linear decrease of sugar conversion when either lignin content or crystallinity of the substrate increases. Importantly, we show that considering the crystallinity of cellulose in addition to the substrate composition is essential to interpret experimental saccharification data. Finally, our findings support the hypothesis of xylan being partially crystalline.