Modeling changes in biomarkers in Gaucher disease patients receiving enzyme replacement therapy using a pathophysiological model
Orphanet Journal of Rare Diseases volume 9, Article number: 95 (2014)
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.
Changes in blood ferritin, chitotriosidase, hemoglobin and platelets were analyzed in French GD Registry patients receiving imiglucerase/alglucerase as ERT. We used simplified exponential pathophysiological model, with initial concentration, biomarkers amplitude of variation and rate constant of normalization during ERT. Changes in four biomarkers were analyzed separately and then all four together from initiation to discontinuation of ERT, or until the end of follow-up. Several covariates were tested, including age at ERT initiation, splenectomy, sex, genotype (N370S/N370S), and ERT dose.
An exponential model gave a good data fit. The four biomarkers analysis showed that the rate of nomalization was the same for all biomarkers, with a half-life of 0.5 years. Predicted values of biomarkers at ERT’s steady state were 40% and 10% of initial concentrations, for ferritin and chitotriosidase, respectively, and 120% and 200% for hemoglobin and platelets, respectively. We found that 3 covariates had an effect on initial concentration or on amplitude of variation in ferritin, hemoglobin and platelets: women and patients under 15 years of age had lower ferritin and hemoglobin concentrations, and patients under 15 years of age had higher platelet count. Splenectomized patients had higher ferritin concentrations and platelet count and lower amplitude of variation of hemoglobin.
We report the first dynamic model of biomarker changes in GD. It enabled us to estimate that 95% of biomarker response to ERT was achieved in 2 years, but with high inter-patient variability. We also found that with the current treatment, normalization of chitotriosidase and ferritin will occur in about 65% of patients.
Gaucher disease (GD)  is a recessively inherited lysosomal storage disorder caused by deficiency of a lysosomal enzyme, glucocerebrosidase (EC 188.8.131.52), which leads to insufficient clearance of the enzyme’s substrate, cellular glucosylceramide. Pathologic accumulations of glucosylceramide (or other substrates, such as glucosylsphingosine) in the lysosomes of tissue macrophages (Gaucher cells) results in splenomegaly, hepatomegaly and multiple forms of skeletal disease . Three clinical phenotypes have been described: type 1, the prevalent form usually defined by the absence of central nervous system impairment; and types 2 and 3, both rare and severe, have central neurological involvement . GD diagnosis is confirmed by the detection of low glucocerebrosidase activity, usually less than 30% of the normal value in peripheral leukocytes. Genotyping can sometimes provide prognostic information . More than 250 mutations of the GBA1 gene encoding lysosomal glucocerebrosidase have been reported as being associated with GD, but the predominant mutation in type 1 GD is called N370S (or c.1226A > G) . The N370S mutation is usually protective against neuronopathic disease. GD can be treated by enzyme replacement therapy (ERT). The first enzyme preparation used to treat GD consisted of placenta-derived glucocerebrosidase (alglucerase available in 1991, Genzyme Corporation) with modified mannose-terminated glycans, allowing more selective uptake by tissue macrophages, the prominent storage cells in GD [6–8]. This preparation was replaced in 1996 by recombinant enzyme (imiglucerase), which was therapeutically equivalent in terms both of safety and efficacy . Two new biosimilar agents are now available: velaglucerase-alfa (Shire)  and taliglucerase-alfa (Pfizer) . ERT reduces macrophagic substrate accumulation, but no routine substrate assay is currently available. A substrate reduction therapy (miglustat, Actelion) that can be prescribed in a few indications has been available since 2002.
The levels of several biomarkers (e.g., in our article, ferritin and chitotriosidase, but also tartrate-resistant acid phosphatase or angiotensin-converting enzyme, not analyzed in this study) change during the clinical course of GD [12–14] due to macrophagic activation: their concentrations rise with disease progression and generally decrease during ERT [12, 15, 16]. These variations can predict bone complications . High blood ferritin during the course of type 1 GD may reflect macrophage activation triggered by substrate accumulation, demonstrated by increases in CCL18 and macrophage inflammatory protein-1α or 1β [18, 19]. ERT is associated with a dramatic decrease of blood ferritin [16, 17], which is more pronounced in patients with an intact spleen . Chitotriosidase is massively produced by storage cells and there is a linear relationship between chitotriosidase and glucosylceramide levels, as shown in spleen sections from patients with GD . Chitotriosidase values drop sharply during ERT, when substrate accumulation decreases, coinciding with clinical improvements .
Hematological abnormalities (anemia and thrombocytopenia) are common in GD, because Gaucher cell infiltration leads to hypersplenism (increased destruction or sequestration of red blood cells or platelets) and bone-marrow insufficiency (decreased production) . Splenectomy, performed essentially before 1991, increases baseline platelet count and decreases the slope of platelet clearance during ERT . Anemia and thrombocytopenia can be used as biomarkers to manage GD patients.
Grabowski et al.  developed an Emax model to describe changes in hemoglobin and platelets and in splenic and hepatic volume during ERT of patients in the International Collaborative Gaucher Group Registry. Biomarkers of French GD patients from a single center was modeled before and during ERT by Stirnemann et al. , but no physiological model was proposed to analyze changes in biomarkers levels during ERT. Nonlinear mixed effects models  are widely used to analyze biological processes described by repeated longitudinal data. They allow estimation of the mean value of the parameters and their inter-individual variability. These models allow a sparse sampling design with few data points per individual in a large set of individuals.
The aim of this study was to develop a pathophysiological model explaining the response of biomarkers (ferritin, chitotriosidase, hemoglobin, platelets) to ERT and to analyze the influence of several covariates.
Patients and methods
Patients and data
We analyzed changes in four biomarkers, ferritin, chitotriosidase, hemoglobin and platelets, in patients receiving ERT with alglucerase from 1991 to 1996 and thereafter imiglucerase, from the French GD Registry (FGDR) . The French Data Protection Commission (CNIL) approval of the FGDR required oral or written informed consent from patients or their parents. The local Institutional Review Board of Northern Paris Hospitals, Paris–Diderot University, AP–HP (Ethics Committee) reviewed and approved the initial research project (220–08). These patients were included in the study if the last follow-up visit was between 2009 and 2010. This investigation was undertaken to describe changes in biomarkers from initiation of ERT until 31 December 2010 corresponding to the limit of collected retrospective data in the FGDR. These data correspond also to those used in the first description of patients of the FGDR . All measurements of biomarkers from 2 years before the initiation of ERT to discontinuation of ERT (interruption of more than 6 months) or the end of follow-up (December 2010) were used. When the chitotriosidase activity was undetectable (patients with probable homozygous chitotriosidase-gene deficiency), this biomarker was not retested  and the data were not included in the analysis. When ERT was interrupted for less than 6 months, patients were considered to be still under treatment. Several covariates were tested including age at initiation of ERT, splenectomy, genotype (N370S/N370S or others), sex, dose at initiation of treatment (divided into 3 classes: lower than 90 IU/kg/month, between 90 and 120 IU/kg/month or higher than 120 IU/kg/month) and average dose in the third year of treatment (divided into the same 3 classes) for patients followed up for at least 3 years.
Model for biomarker changes during ERT
We defined a pathophysiological model  of GD and treatment by ERT (Figure 1). GD is caused by a glucocerebrosidase deficiency that leads to an accumulation of glucosylceramide, which is no longer degraded. As a consequence, there is an increase in ferritin and chitotriosidase production. With ERT, patients are supplied with glucocerebrosidase which leads to the degradation of the glucosylceramide and then improvements of biomarkers levels. The differential equation of this model for chitotriosidase is presented in Supplementary material (Additional file 1). As only biomarker measurements were available, we made a further simplification. Since the rate constant of biomarker elimination is very rapid compared with the rate constant of normalization of glucosylceramide under treatment (k), it was neglected so that we have an exponential increase for chitotriosidase. For chitotriosidase activity, the model is:
where C0 is the initial concentration, r is the amplitude of variation and k is the rate constant of normalization during ERT. For serum ferritin, a similar reasoning leads to the same equation. From the model, we can derive the corrected values of biomarkers at steady state of ERT as C0r (Figure 1) and the normalization half-life as .
Even though two different mechanisms explain hematological abnormalities (destruction by hypersplenism and decrease hematopoiesis), we proposed to simplify the model with only one predominant mechanism (destruction). By similar reasoning, for hemoglobin and platelets, the model is:
which increases from C0 to C0r (Figure 1).
We used a nonlinear mixed effects model to analyze measurements of the four biomarkers. First, we performed a separate analysis of each biomarker. We assumed an exponential random effect on all parameters and that these random effects have a normal distribution with a mean of 0 and a standard deviation of ω. The residual error model was supposed to be proportional with a standard deviation, σ. We also tested the correlation between random effects. Model selection was performed using the Bayesian information criterion (BIC) which adds to the log-likelihood penalty increases with the number of model parameters and the sample size used (BIC = −2 × logL + k log (N)), where L is the maximum likelihood, k is the number of model parameters and N is the number of patients. BIC is a parsimonious criterion where the lowest value corresponds to the best model. We compared our exponential model with the Emax model of Grabowski et al.  using the BIC.
We tested the impact of the covariates mentioned above to explain part of the variability . As only discrete covariates were considered, the effect of the covariate is to change the parameter by exp(β) where β is the estimated regression coefficient. The final model was built using a two-step approach. In the first step, individual empirical Bayesian estimates of parameters were generated from the basic model without covariates. We searched for univariate associations with covariates using Wilcoxon or Spearman tests. In the second step, all covariates with a p-value of < 0.2 were entered into a multivariate model. We then performed a multivariate analysis using a forward selection. P-values of the Wald test were assessed at the 0.05 level. Splenectomy has an impact on changes in platelet count . So we considered (and tested) that splenectomy has an impact on the initial platelet count and on its amplitude of variation.
We then analyzed the 4 biomarkers jointly and evaluated whether one or more parameters were correlated or if several parameters had similar values, to simplify the model. We evaluated the final model with various goodness-of-fit plots : individual fits, graphs of individual weighted residuals versus time and visual predictive check. The significant covariates obtained previously in models of each biomarker were added. Then, we performed a backward selection and computed the P-values of the Wald test. Following estimation of population parameters, we could estimate individual parameters. Individual empirical Bayesian estimates were obtained as maximum a posteriori. From the individual parameters, we estimated individual half-life, time to 95% of response (as 4.3 half-lives), values of biomarkers at steady state of treatment and whether this value is normal or not.
Estimations were performed using the Stochastic Approximation Expectation Maximization (SAEM) algorithm in MONOLIX 4.2.0 (Lixoft, Orsay, France, available at http://www.lixoft.com) . BIC and log-likelihood were estimated by importance sampling.
In the FGDR, 233 patients receiving imiglucerase had a medical visit in 2009 or 2010. The diagnosis of GD was confirmed biochemically in all patients, demonstrating a deficiency of glucocerebrosidase activity either in leukocytes or in a fibroblast cell line. Their median (range) follow-up with ERT was 9 (0–19) years. Median age at the initiation of ERT was 34.3 (1.0-76.0) years, with 18% < 15 years of age, 115 (49.4%) were male, and 61 (26.2%) had a splenectomy before the initiation of the ERT. Twenty-six (11.1%) were N370S/N370S homozygotes. Other characteristics of the patients are presented in Table 1. The median (range) dose at initiation of ERT was 120 IU/kg/month (28 IU/kg/month - 240 IU/kg/month) with 169 (79.0%) patients with 120 IU/kg/month. The average dose during the third year was lower than 90 IU/kg/month for 40 (17.4%) patients, higher than 120 IU/kg/month for 13 (5.7%) and between 90 and 120 IU/kg/month for 176 (76.9%). In some patients there was no monitoring of biochemical parameters during follow-up, or monitoring of only one of the four biomarkers. Therefore, data were available for 129, 142, 191 and 191 patients with a total number of observations of 586, 596, 1287, and 1336 for ferritin, chitotriosidase, hemoglobin and platelets, respectively. Median length of follow-up during ERT was 5 years for ferritin and 6 years for the other biomarkers. Over 65% of the patients were followed up for more than 3 years.
Changes in biomarkers
Figure 2 shows the individual changes in each biomarker during ERT in a spaghetti plot. The exponential model fitted the data. Comparison of our exponential model with the Emax model gave a gain BIC of 3 for ferritin, 12 for chitotriosidase, and 22 for platelets for our model and a loss of 7 for hemoglobin. We tested the impact of covariates in the model of separate biomarkers and found that 3 covariates (initiation of ERT < 15 years of age, splenectomy and sex) had an effect on the initial concentration and/or on the amplitude of variation in ferritin, hemoglobin and platelets. Genotype and dose had no significant impact on the parameters and none of the covariates had a significant impact on chitotriosidase.
The best model of analysis of the four biomarkers together was the one with the same parameter k for all biomarkers, so we considered a rate constant of normalization during ERT common to the four biomarkers. The model with the same parameter k had a smaller BIC of 22 than the model with 4 different values for parameters k. Values of parameters without covariates (only splenectomy for platelets) are given in Table 2 and Table 3: Predicted values of biomarkers at steady state of ERT were 38% and 10% of initial concentrations, for ferritin and chitotriosidase, respectively, and 117% and 200% for hemoglobin and platelets, respectively. Estimated corrected values at steady state of ERT were 179.4 ng/L, 792 nmol/h.mL and 13.9 g/dL for ferritin, chitotriosidase and hemoglobin, respectively. The estimated corrected platelet counts at steady state of ERT were 156,030/mm3 and 195,037/mm3 for non-splenectomized and splenectomized patients, respectively. Individual predicted values show that, under current ERT, normalization will occur for 58% of the patients for ferritin, 66% for chitotriosidase, 88% for hemoglobin and 64% for platelets. Figure 3 shows individual fits for one splenectomized patient and one non-splenectomized. Relative standard errors were lower than 20% for both fixed effects and variabilities. Diagnostic plots (Additional file 2) show that the model describes the data adequately.
Values of parameters in the joint model with covariates are given in Table 4. The covariate effects were the same as in the models considering the biomarkers separately. Figure 4 shows the response of the biomarkers to ERT, according to significant covariates. Non-splenectomized (p-value = 10−6) and age <15 years (p-value < 10−10) patients had a lower initial concentration of ferritin. Women (p-value = 0.002) and non-splenectomized (p-value = 10−4) patients had a lower amplitude of variation in ferritin. Women (p-value = 0.01) and age <15 years (p-value = 10−8) patients had a lower initial concentrations of hemoglobin. Splenectomized (p-value = 10−5) patients had a lower amplitude of variation in hemoglobin. Splenectomized (p-value < 10−10) and age <15 years (p-value = 10−6) patients had a higher initial platelet count. Splenectomized (p-value = 10−8) patients had a lower amplitude of variation in platelets.
In the final model, we estimated that during ERT the half-life of normalization was 0.5 years (95% CI = [0.4-0.6]), with a variability of 97%. Half-life was 0.5 (0.1-3.1) years and 95% of response was obtained after 2 (0.4-13.3) years under ERT. Variabilities were reasonable for ferritin, hemoglobin and platelet count. Greater variability was found for chitotriosidase, doubtless because of variability in measurement error (a sample dilution is usually needed when values are high).
Our pathophysiological model predicts changes in biomarkers on ERT and estimates the rate constant of normalization. For the final model, we estimated a normalization half-life of 0.5 years during ERT for all four biomarkers. Only 2 studies have modeled the changes in biomarkers and they used different models: Emax and linear mixed models. In the study of Grabowski et al. , the ERT dose effect had a significant impact on biomarkers with a large sample size. Emax is an empirical model with a hyperbolic function; it is not the result of a physiological model. These models are used in pharmacometrics  to study the link between dose and concentration. However, the drug often acts on a biological quantity, modifying its production or elimination. We modeled this quantity by a model with a production rate and a rate constant of elimination. Using physiological knowledge, we obtained differential equations to explain changes in biomarkers over time. For instance, 90% response is obtained after 9 T50 for an Emax model, whereas it is reached after only 3.3 half-lives (similar to T50) for an exponential model.
Our results show improvement in all biomarkers under ERT (decrease in ferritin and chitotriosidase and increase in hemoglobin and platelets). Stein et al.  also highlighted an increase in ferritin and a renormalization under ERT. Hollak et al.  reported a decrease of chitotriosidase of 32% in 1 year, and we found a 95% response in 2 years and a 36% response in 1 year. De Fost et al.  reported that 53% of patients had anemia at baseline and 58% had thrombocytopenia, with renormalization under ERT, and showed a similar pattern of response after 1 year under ERT.
Patients <15 years of age have lower initial concentrations of ferritin and hemoglobin but higher platelet counts; at initiation of ERT, women have a lower concentration of hemoglobin and splenectomized patients have a higher platelet count and ferritin. For GD children from the International Collaborative Gaucher Group Registry, Kaplan P et al.  noted that 50% had platelet counts less than 120 × 103/mm3 and 40% had anemia at the time of diagnosis. Hemoglobin and ferritin tend to be lower in GD women, as for other women, probably because of menstruation, with no link to GD. In our model, no covariates had a significant impact on the chitotriosidase changes. In contrast, Stirnemann et al.  found that splenectomy and GD genotype affected chitotriosidase activity, but their study was limited by a small number of samples, and they used a different model.
Our model allows estimation of the rate constant of biomarker improvement using a pathophysiological model. Bone events are the most debilitating and disabling complication of GD. With substrate overload, Gaucher cells activate and induce proinflammatory cytokine synthesis which can modify the activity of the osteoblast-osteoclast system and promote lytic phenomena and intraosseous vascular complications [34, 35]. Further analysis of the interaction between biomarkers and bone events is needed.
Recently, plasma glucosylsphingosine has been proposed as a biomarker for GD [36, 37] and could be used as a reflection of intracellular glucosylceramide. Our model may also predict changes in intracellular glucosylceramide and/or glucosylsphingosine.
Our model could be used to study the effect of biomarker changes on complications such as the occurrence of bone events. Using a few measurements of biomarkers post-treatment, we can estimate the rate constant of normalization and individual amplitude of variation in biomarkers in order to predict further response to treatment. Our model could help to define an individual risk for complications and to refine the best ERT regimens.
Parameters of this pathological model could be used to predict value at steady state with the dosages in the first months of ERT. These predicted values could be used as an individual objective of treatment to modify dosage for each patient. Clinical studies are needed to confirm this hypothesis.
We developed a model which could in the future be used to manage patients, but further studies are needed to confirm clinical applications. This model with 2 parameters (in addition to the baseline value) predicts the level of improvement of biomarkers and the time to 95% response, using only a few measurements per patient.
This is to our knowledge the first study of changes in biomarkers in Gaucher disease using a pathophysiological model. With the model we estimate that 95% of biomarker response to ERT is achieved in 2 years, but with high inter-patient variability. We also found that with the current treatment, normalization of chitotriosidase and ferritin will occur in about 65% of patients.
Bayesian information criterion
French Gaucher disease registry
Grabowski GA: Advances in Gaucher Disease: Basic and Clinical Perspectives. London: Future Medicine Ltd; 2013.
Brady RO, Kanfer JN, Shapiro D: Metabolism of glucocerebrosides. II. Evidence of an enzymatic deficiency in Gaucher’s disease. Biochem Biophys Res Commun. 1965, 18: 221-225. 10.1016/0006-291X(65)90743-6.
Beutler E: Gaucher disease. Arch Intern Med. 1999, 159: 881-882. 10.1001/archinte.159.8.881.
Grabowski GA: Recent clinical progress in Gaucher disease. Curr Opin Pediatr. 2005, 17: 519-524. 10.1097/01.mop.0000172702.33128.19.
Jmoudiak M, Futerman AH: Gaucher disease: pathological mechanisms and modern management. Br J Haematol. 2005, 129: 178-188. 10.1111/j.1365-2141.2004.05351.x.
Barton NW, Brady RO, Dambrosia JM, Di Bisceglie AM, Doppelt SH, Hill SC, Mankin HJ, Murray GJ, Parker RI, Argoff CE: Replacement therapy for inherited enzyme deficiency–macrophage-targeted glucocerebrosidase for Gaucher’s disease. N Engl J Med. 1991, 324: 1464-1470. 10.1056/NEJM199105233242104.
Hollak CE, Aerts JM, Goudsmit R, Phoa SS, Ek M, Van Weely S, Von dem Borne AE, Van Oers MH: Individualised low-dose alglucerase therapy for type 1 Gaucher’s disease. Lancet. 1995, 345: 1474-1478. 10.1016/S0140-6736(95)91037-9.
Bijsterbosch MK, Donker W, Van de Bilt H, Van Weely S, Van Berkel TJ, Aerts JM: Quantitative analysis of the targeting of mannose-terminal glucocerebrosidase: predominant uptake by liver endothelial cells. Eur J Biochem. 1996, 237: 344-349. 10.1111/j.1432-1033.1996.00344.x.
Grabowski GA, Barton NW, Pastores G, Dambrosia JM, Banerjee TK, McKee MA, Parker C, Schiffmann R, Hill SC, Brady RO: Enzyme therapy in type 1 Gaucher disease: comparative efficacy of mannose-terminated glucocerebrosidase from natural and recombinant sources. Ann Intern Med. 1995, 122: 33-39. 10.7326/0003-4819-122-1-199501010-00005.
Zimran A, Altarescu G, Philips M, Attias D, Jmoudiak M, Deeb M, Wang N, Bhirangi K, Cohn GM, Elstein D: Phase 1/2 and extension study of velaglucerase alfa replacement therapy in adults with type 1 Gaucher disease: 48-month experience. Blood. 2010, 115: 4651-4656. 10.1182/blood-2010-02-268649.
Zimran A, Brill-Almon E, Chertkoff R, Petakov M, Blanco-Favela F, Muñoz ET, Solorio-Meza SE, Amato D, Duran G, Giona F, Heitner R, Rosenbaum H, Giraldo P, Mehta A, Park G, Phillips M, Elstein D, Altarescu G, Szleifer M, Hashmueli S, Aviezer D: Pivotal trial with plant cell-expressed recombinant glucocerebrosidase, taliglucerase alfa, a novel enzyme replacement therapy for Gaucher disease. Blood. 2011, 118: 5767-5773. 10.1182/blood-2011-07-366955.
Hollak CE, Van Weely S, Van Oers MH, Aerts JM: Marked elevation of plasma chitotriosidase activity: a novel hallmark of Gaucher disease. J Clin Invest. 1994, 93: 1288-1292. 10.1172/JCI117084.
Morgan MA, Hoffbrand AV, Laulicht M, Luck W, Knowles S: Serum ferritin concentration in Gaucher’s disease. Br Med J (Clin Res Ed). 1983, 286: 1864. 10.1136/bmj.286.6381.1864.
Cohen LA, Gutierrez L, Weiss A, Leichtmann-Bardoogo Y, Zhang D, Crooks DR, Sougrat R, Morgenstern A, Galy B, Hentze MW, Lazaro FJ, Rouault TA, Meyron-Holtz EG: Serum ferritin is derived primarily from macrophages through a nonclassical secretory pathway. Blood. 2010, 116: 1574-1584. 10.1182/blood-2009-11-253815.
Cabrera-Salazar MA, O’Rourke E, Henderson N, Wessel H, Barranger JA: Correlation of surrogate markers of Gaucher disease: implications for long-term follow up of enzyme replacement therapy. Clin Chim Acta. 2004, 344: 101-107. 10.1016/j.cccn.2004.02.018.
Mekinian A, Stirnemann J, Belmatoug N, Heraoui D, Fantin B, Fain O, Charpentier A, Rose C: Ferritinemia during type 1 Gaucher disease: mechanisms and progression under treatment. Blood Cells Mol Dis. 2012, 49: 53-57. 10.1016/j.bcmd.2012.04.002.
Stirnemann J, Belmatoug N, Vincent C, Fain O, Fantin B, Mentré F: Bone events and evolution of biologic markers in Gaucher disease before and during treatment. Arthritis Res Ther. 2010, 12: R156. 10.1186/ar3111.
Boot RG, Verhoek M, De Fost M, Hollak CEM, Maas M, Bleijlevens B, Van Breemen MJ, Van Meurs M, Boven LA, Laman JD, Moran MT, Cox TM, Aerts JMFG: Marked elevation of the chemokine CCL18/PARC in Gaucher disease: a novel surrogate marker for assessing therapeutic intervention. Blood. 2004, 103: 33-39. 10.1182/blood-2003-05-1612.
Van Breemen MJ, De Fost M, Voerman JSA, Laman JD, Boot RG, Maas M, Hollak CEM, Aerts JM, Rezaee F: Increased plasma macrophage inflammatory protein (MIP)-1alpha and MIP-1beta levels in type 1 Gaucher disease. Biochim Biophys Acta. 2007, 1772: 788-796. 10.1016/j.bbadis.2007.04.002.
Stein P, Yu H, Jain D, Mistry PK: Hyperferritinemia and iron overload in type 1 Gaucher disease. Am J Hematol. 2010, 85: 472-476. 10.1002/ajh.21721.
Bussink AP, Van Eijk M, Renkema GH, Aerts JM, Boot RG: The biology of the Gaucher cell: the cradle of human chitinases. Int Rev Cytol. 2006, 252: 71-128.
Guggenbuhl P, Grosbois B, Chalès G: Gaucher disease. Joint Bone Spine. 2008, 75: 116-124. 10.1016/j.jbspin.2007.06.006.
Grabowski GA, Kacena K, Cole JA, Hollak CEM, Zhang L, Yee J, Mistry PK, Zimran A, Charrow J, Vom Dahl S: Dose–response relationships for enzyme replacement therapy with imiglucerase/alglucerase in patients with Gaucher disease type 1. Genet Med. 2009, 11: 92-100. 10.1097/GIM.0b013e31818e2c19.
Lindstrom MJ, Bates DM: Nonlinear mixed effects models for repeated measures data. Biometrics. 1990, 46: 673. 10.2307/2532087.
Stirnemann J, Vigan M, Hamroun D, Heraoui D, Rossi-Semerano L, Berger MG, Rose C, Camou F, De Roux-Serratrice C, Grosbois B, Kaminsky P, Robert A, Caillaud C, Froissart R, Levade T, Masseau A, Mignot C, Sedel F, Dobbelaere D, Vanier MT, Valayanopoulos V, Fain O, Fantin B, Billette T, Mentré F, Belmatoug N: The French Gaucher’s disease registry: clinical characteristics, complications and treatment of 562 patients. Orphanet J Rare Dis. 2012, 7: 77. 10.1186/1750-1172-7-77.
Aerts JM, Hollak CE: Plasma and metabolic abnormalities in Gaucher’s disease. Baillieres Clin Haematol. 1997, 10: 691-709. 10.1016/S0950-3536(97)80034-0.
Upton RN, Mould DR: Basic concepts in population modeling, simulation, and model-based drug development: part 3—introduction to pharmacodynamic modeling methods. CPT: Pharmacomet Syst Pharmacol. 2014, 3: e88. 10.1038/psp.2013.71.
Mould DR, Upton RN: Basic concepts in population modeling, simulation, and model-based drug development-part 2: introduction to pharmacokinetic modeling methods. CPT: Pharmacomet Syst Pharmacol. 2013, 2: e38-10.1038/psp.2013.14.
Karlsson MO, Savic RM: Diagnosing model diagnostics. Clin Pharmacol Ther. 2007, 82: 17-20. 10.1038/sj.clpt.6100241.
Kuhn E, Lavielle M: Maximum likelihood estimation in nonlinear mixed effects models. Comput Stat Data Anal. 2005, 49: 1020-1038. 10.1016/j.csda.2004.07.002.
Hollak CE, Maas M, Aerts JM: Clinically relevant therapeutic endpoints in type I Gaucher disease. J Inherit Metab Dis. 2001, 24 (Suppl 2): 97-105. discussion 87–88.
De Fost M, Hollak CEM, Groener JEM, Aerts JMFG, Maas M, Poll LW, Wiersma MG, Häussinger D, Brett S, Brill N, Vom Dahl S: Superior effects of high-dose enzyme replacement therapy in type 1 Gaucher disease on bone marrow involvement and chitotriosidase levels: a 2-center retrospective analysis. Blood. 2006, 108: 830-835. 10.1182/blood-2005-12-5072.
Kaplan P, Andersson HC, Kacena KA, Yee JD: The clinical and demographic characteristics of nonneuronopathic Gaucher disease in 887 children at diagnosis. Arch Pediatr Adolesc Med. 2006, 160: 603-608. 10.1001/archpedi.160.6.603.
Pavlova EV, Deegan PB, Tindall J, McFarlane I, Mehta A, Hughes D, Wraith JE, Cox TM: Potential biomarkers of osteonecrosis in Gaucher disease. Blood Cells Mol Dis. 2011, 46: 27-33. 10.1016/j.bcmd.2010.10.010.
Mucci JM, Suqueli García F, De Francesco PN, Ceci R, Di Genaro S, Fossati CA, Delpino MV, Rozenfeld PA: Uncoupling of osteoblast-osteoclast regulation in a chemical murine model of Gaucher disease. Gene. 2013, 532: 186-191. 10.1016/j.gene.2013.09.072.
Dekker N, Van Dussen L, Hollak CEM, Overkleeft H, Scheij S, Ghauharali K, Van Breemen MJ, Ferraz MJ, Groener JEM, Maas M, Wijburg FA, Speijer D, Tylki-Szymanska A, Mistry PK, Boot RG, Aerts JM: Elevated plasma glucosylsphingosine in Gaucher disease: relation to phenotype, storage cell markers, and therapeutic response. Blood. 2011, 118: e118-127. 10.1182/blood-2011-05-352971.
Rolfs A, Giese A-K, Grittner U, Mascher D, Elstein D, Zimran A, Böttcher T, Lukas J, Hübner R, Gölnitz U, Röhle A, Dudesek A, Meyer W, Wittstock M, Mascher H: Glucosylsphingosine is a highly sensitive and specific biomarker for primary diagnostic and follow-up monitoring in Gaucher disease in a non-jewish, caucasian cohort of Gaucher disease patients. PLoS One. 2013, 8: e79732. 10.1371/journal.pone.0079732.
The authors wish to thank the RCLD assistant of the RCLD, Monia Bengherbia, the RCLD secretary, Samira Zebiche, the GD patients and their physicians, who contributed to the FDGR. The development of the original software for the French Gaucher Disease Registry (FGDR) was funded by a grant from the nonprofit organization VML (Vaincre les Maladies Lysosomales). The FGDR was funded, in part, by INSERM and InVS (Institut national de Veille Sanitaire).
FM had consulting fees and INSERM UMR1137 received a grant from Sanofi. JS had travel fees from Sanofi-Genzyme. NB had travel fees, consulting fees, fees for speaking and received grants (Genzyme, Shire, Actelion, Pfizer) donated to the Department of Clinical Research of the Assistance-Publique Hôpitaux de Paris. The other authors declared no financial disclosures.
MV, JS and FM designed the research and analyzed and interpreted the data. MV, JS and FM wrote the first draft of the paper, which was then corrected and approved by all authors.
About this article
Cite this article
Vigan, M., Stirnemann, J., Caillaud, C. et al. Modeling changes in biomarkers in Gaucher disease patients receiving enzyme replacement therapy using a pathophysiological model. Orphanet J Rare Dis 9, 95 (2014). https://doi.org/10.1186/1750-1172-9-95
- Gaucher disease
- French registry
- Enzyme replacement therapy