Fungal pretreatment, using lignin-degrading microorganisms to improve lignocellulosic feedstocks with minimal energy input, is a potential alternative to physiochemical pretreatment methods. Identifying the kinetics for fungal pretreatment during solid substrate cultivation is needed to help establish the processing conditions for effective scale up of this technology. In this study, a set of mathematical models were proposed for describing the interactions between holocellulose consumption, lignin degradation, cellulase, ligninolytic enzyme, and the growth of Phanerochaete chrysosporium during a 14 day fungal pretreatment process. Model parameters were estimated and validated by the System Biology Toolbox in MatLab. Developed models provided sufficiently accurate predictions for fungal growth (R (2) = 0.97), holocellulose consumption (R (2) = 0.97), lignin degradation (R (2) = 0.93) and ligninolytic enzyme production (R (2) = 0.92), and fair prediction for cellulase production (R (2) = 0.61). The models provide valuable information for understanding the interactive mechanisms in biological systems as well as for fungal pretreatment process scale up and improvement.