Gaucher disease (GD) is a rare recessively inherited disorder caused by deficiency of a lysosomal enzyme, glucocerebrosidase. Accumulation of glucosylceramide or glucosylsphingosine in macrophages leads to increased production of ferritin and chitotriosidase and to decreases in hemoglobin concentration and platelet count, which are used as blood biomarkers. GD is treated by enzyme replacement therapy (ERT) or, sometimes by substrate reduction therapy. However, no physiological model for analysis of biomarkers change during ERT has been proposed. We aimed to develop a pathophysiological model to analyze biomarker's response to ERT and several covariates impact.