Next Article in Journal
Primary Care Physician Experiences with Integrated Population-Scale Genetic Testing: A Mixed-Methods Assessment
Next Article in Special Issue
Modelling Neuromuscular Diseases in the Age of Precision Medicine
Previous Article in Journal
Convergent and Concurrent Validity between Clinical Recovery and Personal-Civic Recovery in Mental Health
Previous Article in Special Issue
Clinical and Laboratory Associations with Methotrexate Metabolism Gene Polymorphisms in Rheumatoid Arthritis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Serum Pharmacodynamic Biomarkers in Prednisone-Versus Deflazacort-Treated Duchenne Muscular Dystrophy Boys

by
Shefa Tawalbeh
1,2,
Alison Samsel
2,
Heather Gordish-Dressman
3,
Yetrib Hathout
2,*,
CINRG-DNHS Investigators
and
Utkarsh J. Dang
4,*
1
Department of Biomedical Engineering, Binghamton University, Binghamton, NY 13902, USA
2
Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, Binghamton University, Binghamton, NY 13902, USA
3
Division of Biostatistics, Children’s National Hospital, Washington, DC 20010, USA
4
Department of Health Outcomes and Administrative Sciences, School of Pharmacy and Pharmaceutical Sciences, Binghamton University, Binghamton, NY 13902, USA
*
Authors to whom correspondence should be addressed.
Names and affiliations of the CINRG-DNHS Investigators are provided in the Acknowledgments.
J. Pers. Med. 2020, 10(4), 164; https://doi.org/10.3390/jpm10040164
Submission received: 7 August 2020 / Revised: 6 October 2020 / Accepted: 10 October 2020 / Published: 12 October 2020

Abstract

:
Prednisone (Pred) and Deflazacort (Dfz) are commonly used glucocorticoids (GCs) for Duchenne muscular dystrophy (DMD) treatment and management. While GCs are known to delay the loss of ambulation and motor abilities, chronic use can result in onerous side effects, e.g., weight gain, growth stunting, loss of bone density, etc. Here, we use the CINRG Duchenne natural history study to gain insight into comparative safety of Pred versus Dfz treatment through GC-responsive pharmacodynamic (PD) biomarkers. Longitudinal trajectories of SOMAscan® protein data obtained on serum of DMD boys aged 4 to 10 (Pred: n = 7; Dfz: n = 8) were analyzed after accounting for age and time on treatment. Out of the pre-specified biomarkers, seventeen candidate proteins were differentially altered between the two drugs (p < 0.05). These include IGFBP-2 and AGER associated with diabetes complications, and MMP-3 associated with extracellular remodeling. As a follow-up, IGFBP-2, MMP-3, and IGF-I were quantified with an ELISA using a larger sample size of DMD biosamples (Dfz: n = 17, Pred: n = 12; up to 76 sera samples) over a longer treatment duration. MMP-3 and IGFBP-2 validated the SOMAscan® signal, however, IGF-I did not. This study identified GC-responsive biomarkers, some associated with safety, that highlight differential PD response between Dfz and Pred.

Graphical Abstract

1. Introduction

Duchenne muscular dystrophy (DMD) is an X-linked recessive disorder affecting the expression of dystrophin protein, an essential protein that maintains muscle fiber integrity and function [1]. The lack of dystrophin expression leads to muscle inflammation at an early stage of the disease, followed by progressive muscle degeneration and wasting [2]. While there are promising gene therapy and exon-skipping treatments, some of which have received conditional approval from the EMA [3] and FDA [4,5,6,7], these are mutation-specific and only restore a partial amount of truncated dystrophin protein [8]. Hence, DMD patients will continue to need combination therapy using the current standard of care: glucocorticoids (GCs). GC use helps reduce muscle inflammation and delay the loss of motor abilities [9], delay loss of independent ambulation, improve pulmonary function, and delay the onset of cardiomyopathy [10,11,12]. However, a GC regimen can have many side effects, such as weight gain, growth stunting, loss of bone density, hirsutism, Cushingoid features, osteoporosis, hypertension, diabetes, behavioral disturbances, and difficulty in sleeping [12,13,14]. Commonly used steroidal drugs for DMD are prednisone (Pred) and deflazacort (Dfz). Prednisone has been used to treat DMD patients in the USA since the seventies [15]. In February 2017, the FDA approved deflazacort to treat DMD patients aged five years and older [16] although this same drug has been widely used in Europe to treat DMD patients for years. Figure 1 shows structures of prednisone and deflazacort along with their respective active metabolites: prednisolone and 21-desacetyl deflazacort, respectively. Both drugs are given to patients in their prodrug forms, which are then metabolized to their active forms in the liver. The active metabolite prednisolone has a hydroxyl group at carbon 17, while the active 21-desacertyl deflazacort has an oxazoline structure at that same position (red arrow in Figure 1).
Many efforts have been undertaken to investigate the use of biomarkers in DMD [17,18,19] and the efficacy and safety of Pred versus Dfz [10,11,12,20,21]. Both drugs are known to be efficacious in prolonging ambulation [10,11,14,22]. For our purposes here, we provide a small review on safety comparisons from the literature. These studies were compared to clinical outcome measures like changes in height, weight, Cushingoid features, and erythema, among others to evaluate differences in safety. A comparison of high dose Dfz and Pred in 70 systemic lupus patients concluded that Pred was associated with a significant increase in weight gain, Cushingoid severity index, and hirsutism compared to Dfz [23]. In a double-blind, randomized study on 196 DMD subjects, Pred was found to be associated with more weight gain than Dfz [10]. In another randomized study (18 DMD patients), patients treated with Pred showed higher weight gain compared to Dfz-treated patients [13]. As compared to other studies [10,14], a comparison of prednisone/prednisolone versus Dfz treated patients using 340 participants from the Cooperative International Neuromuscular Research Group Duchenne Natural History Study (CINRG-DNHS) showed that participants treated with Dfz had increased frequency of Cushingoid appearance, cataracts, and growth delay. Another study [13] reported that the time of initiating, dosing, and duration of treatments were associated with side effects; longer duration and increased Dfz dosage predicted growth stunting and Dfz was reported to be associated with lighter weight and shorter heights compared to Pred [13]. To summarize, these studies suggest that treatment with Pred is associated with relatively more weight gain, whereas Dfz treatment is associated with relatively more growth stunting. While these studies provide important comparisons between Dfz and Pred, the underlying molecular mechanism leading to these differences, remain unknown. Furthermore, predictive outcome measures, especially for adverse effects are highly desirable.
There are ongoing debates among families, clinicians, and regulatory agencies over which GC drug (Pred or Dfz) and regimen is better for DMD boys. GC responsive biomarkers (for both safety and efficacy) have been defined and confirmed in different diseases and multiple cohorts [17,24,25], using data with pooled corticosteroid drugs not differentiating between Pred and Dfz and different regimens [22]. However, much remains unknown about the comparative effects of Pred vs. Dfz on blood accessible biomarkers and how they can inform clinical decision-making and drug development. Adding pharmacodynamic (PD) biomarkers to the evidence from the aforementioned clinical outcome studies could bring insights into differences between Pred and Dfz at the molecular level. To the best of our knowledge, we are the first to compare serum accessible PD biomarkers in DMD patients (4–16 years old) treated with Dfz or Pred to gain insight into differences at the serum protein level. For this, we use data generated by a high-throughput SOMAScan® technique, followed by confirmation of a specific set of PD biomarkers by ELISAs. The objective here is to define PD biomarkers that differ in their response to Pred and Dfz, and investigate if these biomarkers are known to be associated with reported differences in term of side effects between these two drugs.

2. Materials and Methods

The study protocol was approved by Institutional Review Boards at all participating institutions that provided serum samples. These included, the Office of Research IRB administrations at the University of California Davis, Davis, CA, the University of Pittsburgh, Pittsburgh, PA, Children’s National Health System, Washington DC, the Conjoint Health Research Ethics Board, University of Calgary and the Human Subject Research Review Committee at Binghamton University, NY. Informed written consent was obtained from the parents of the participants or their legal guardians for biomarker studies at each site.

2.1. SOMAscan® Dataset

The CINRG-DNHS is a natural history study of 440 DMD boys. The CINRG-DNHS investigators enrolled these subjects at 20 centers in nine countries [11]. In this prospective cohort study, subjects genetically confirmed to have DMD were followed for up to 10 years. Details about the study, including informed consent and others, have previously been published [26]. At entry, some DMD subjects were steroid-naïve and were then treated during follow up visits, some remained steroid-naïve throughout, while others were already on GCs [26]. Characteristics (height, weight, age, time on GC treatment, and type of GC) were recorded on follow-up visits as well as clinical outcome data on time to run/walk, time to stand, and time to climb until loss of ability [26]. Note that six-minute walk distance and the NorthStar Ambulatory Assessment score were measured for a small subset of patients enrolled later during the study [26]. Blood samples for biomarker investigations were collected by the CINRG-DNHS investigators during some visits on a subset of participants. GC-treated patients were on Pred or Dfz with variation in dose and dose schedules (intermittent or daily). In a recent study [17], we generated biomarker data on a subset of these DMD samples (n = 31 boys, 4 to 10 years old) and used steroid-naïve and steroid-treated paired samples to define 107 (false discovery rate-adjusted p-values < 0.05) PD biomarkers from 1310 serum proteins (SOMAscan® array).
Based on the above-mentioned study [17], we pre-specified the 107 PD biomarkers of interest and focused on an analysis of 35 longitudinal samples from 8 Pred-treated and 7 Dfz-treated ambulatory patients only (screening dataset A) to identify potential candidate serum proteins differentially altered between Pred and Dfz. These subjects were age-matched (see Table 1 for a summary of their characteristics). For these longitudinal data, only subjects with at least 2 visits were used. The name of the study, type of treatment, duration of treatment, age, and some clinical measurements (height, weight, BMI) are tabulated in Supplementary Table S1.

2.2. Validation of Key Pharmacodynamic Biomarkers Using ELISA

A larger data set (confirmation dataset B) contained up to 76 longitudinal samples on 29 ambulatory subjects (17 deflazacort-treated and 12 prednisone-treated) from CINRG DNHS. Note that all samples/subjects from screening dataset A were used in confirmation dataset B as well. For both datasets, we focused on ambulatory patients to minimize any possible confounding issues with the stage of disease (after loss of ambulation). Again, these subjects were age-matched. Summary characteristics of the subjects are provided in Table 2.
A subset of PD serum protein biomarkers, including MMP3, IGFBP-2, and IGF-I, were selected for confirmation analysis using ELISA assays. An ELISA kit from Meso Scale Diagnostics (MSD) was used to measure levels of MMP3. Similarly, ELISA kits from Thermo Fisher Scientific and R&D Systems, Inc. were used to measure levels of IGFBP2 and IGF-I, respectively. These candidates were chosen based on the availability of a validated ELISA assay kit. All three assays were sandwich ELISA, which were highly specific for antigen detection. Dilution factors for serum samples were 1:10, 1:400, and 1:100 for MMP-3, IGFBP-2, and IGF-I, respectively. The three assays were performed following the manufacturer’s instructions. Details about serum samples used in ELISA validation are shown in Table 2.

2.3. Data Analysis and Statistical Methods

As depicted in the schematic (Figure 2), we first pre-specified 107 steroidal PD biomarkers previously identified as being significantly altered over time in their levels between GC-naive and paired GC-treated DMD samples from the same subjects [17]. To compare the serum levels and trajectories of these pre-specified 107 PD biomarkers between DMD patients treated with Dfz and Pred, linear mixed effect models were used to analyze longitudinal measurements for 7 boys on Pred and 8 on Dfz (daily regimen). All statistical analyses were performed using R [27]. Linear mixed effect models were run using the lme4 and lmerTest packages [28,29]. Serum protein levels, which had previously been hybridization control and median signal normalized, were log-transformed (see [17] for details). Random intercept linear mixed effect models were used to investigate the association of (mean centered) time on drug, type of GC (Pred or Dfz), and interaction between time on drug and type of GC on log transformed protein RFUs. The samples were age-matched at baseline (average age of Pred-treated subject at baseline = 5.7 (min age = 4.3, max age = 7.3) years; average age of Dfz-treated subject at baseline = 6 (min age = 4.7, max age = 7.4) years) and over time (see Table 1 for age range at sample collection). Some models did not converge (7 out of 107) due to numerical optimization issues (likely due to small sample size). Here, the interaction coefficient represents the difference in estimated slopes of biomarker trajectory over time between Pred-treated and Dfz-treated DMD patients. If the interaction coefficient is significant, this indicates that treatment type (Pred or Dfz) is associated with different trajectories of biomarker over time. The coefficient for the type of GC represents the difference in mean protein levels between Pred-treated and Dfz-treated subjects for a patient with an average treatment duration with similar biomarker trajectories over time (this is considered after the interaction effect). The same model was also fit to log-transformed ELISA data. To study the effect of time on GC on weight, height, and BMI, data analyses were performed using linear mixed effect models on confirmation dataset B. For this, we used a random intercept linear mixed model adjusting for age at baseline, time on GC, and interaction between type of GC and time on GC. Concordance of SOMAscan® and ELISA measurements was investigated on the subset of samples common to screening dataset A and confirmation dataset B.

3. Results

3.1. Longitudinal Trajectory of Serum PD Biomarkers in Prednisone vs. Deflazacort Treated DMD boys

We examined differences in the longitudinal trajectories of pre-specified 107 PD biomarkers in Pred-treated patients (n = 7) versus Dfz-treated DMD patients (n = 8) accounting for duration of treatment (time on GC) and interaction between treatment duration and type of GC (Pred or Dfz). In general, Pred and Dfz seem to engage the PD biomarkers similarly and altered their trajectory in the same direction in screening dataset A. This was not the case for 17 PD biomarkers that were found differentially altered in their average levels and/or longitudinal trajectories (unadjusted p value < 0.05) in Pred-treated vs. Dfz-treated DMD patients. These 17 differentially altered PD biomarkers are listed in Table 3 with fold change between non-GC-treated DMD subjects and healthy controls (from [17]), fold change between Pred and Dfz-treated subjects, p values for the difference in the mean levels and p values for the difference in longitudinal trajectories between Pred and Dfz treated DMD patients, along with potential significance and references to published literature.
Two major groups of differentially affected steroidal PD biomarkers were thus identified. The first major group consisted of PD biomarkers that were repressed by GC in general, but exhibited a significantly lower mean level in Pred treated group relative to Dfz treated group. These included leukocyte immunoglobulin-like receptor subfamily B member 1 (LILRB1), tumor necrosis factor receptor superfamily member 21 (TNFRSF21), chordin-like protein 1 (CHRDL1), soluble advanced glycosylation end product-specific receptor (sRAGE), annexin A2 (ANXA2), CD166 antigen (CD166), scavenger receptor cysteine-rich type 1 protein M130 (sCD163), induced myeloid leukemia cell differentiation protein (MCL-1), transmembrane glycoprotein NMB (GPNMB), mitogen-activated protein kinase 14 (MAPK14), and neural cell adhesion molecule L1 (NCAM-L1). The second (minor) group consisted of PD biomarkers that increased in their levels following GC treatment and consisted of two subgroups. The first subgroup consisted of those that were significantly elevated in serum samples of Dfz relative to the Pred-treated group such as hemojuvelin (HJV), stromelysin-1 (MMP3) while the second group consists of those biomarkers that were significantly elevated in serum samples of Pred relative to Dfz-treated group such as insulin-like growth factor I (IGF-I), ficolin-1 (FCN1), and cGMP-inhibited 3′,5′-cyclic phosphodiesterase A (PDE3A). One interesting PD biomarker is insulin-like growth factor-binding protein 2 (IGFBP-2); this showed diverging longitudinal trajectories (p = 0.040) in the Dfz-treated group, as compared to the Pred-treated group. Similarly, longitudinal trajectories of Ficolin-1 (FCN1; p = 0.03) and Annexin A2 (ANXA2; 0.054) also seemingly diverged in addition to their differential levels between the two drugs. Figure 3 shows selected examples of these differentially affected steroidal PD biomarkers. For example, IGFBP-2 sharply decreased over time in Dfz-treated DMD patients, while it remained unchanged or slightly increased over time in Pred treated group (p = 0.040). While TNFRSF21 (also known as DR6) was not found to be differentially affected by the two different drugs in terms of longitudinal trajectory slopes but its mean sera levels was significantly lower in the Pred-treated group compared to the Dfz-treated group (p = 0.005). Similarly, mean levels of FCN-1 were lower in the Dfz-treated group relative to the Pred-treated group (p = 0.038), while MMP3 mean levels were significantly elevated in the Dfz-treated group relative to the Pred-treated group.

3.2. Effect of Dfz and Pred on Height and Weight of DMD Boys

We also investigated differences in growth stunting in Dfz treated patients relative to Pred treated patients in confirmation dataset B. DMD patients treated with Dfz had lower height growth rates (p = 0.006 for difference in trajectory slopes over time) compared to those treated with Pred (Figure 4). We did not find any significant difference in weight (p = 0.112) and BMI (p = 0.08) over time between patients treated with Dfz and Pred (Figure 4). Additionally, note that more variation is observed in longitudinal BMI measurements around the estimated model (Figure 4).

3.3. Data Validation Using ELISA

For the reported data above, we had a small sample size and we did not adjust for multiple testing, thus a follow-up confirmation analyses on a subset of PD biomarkers was carried out. Unfortunately, from the list of PD biomarkers identified in Table 3 above, a good validated ELISA assay that used low sera volume was available for only three PD biomarker candidates: MMP3, IGFBP2 and IGF-I. There were other ELISA assays for other candidates, but they were either not validated or required a larger volume of sera samples, which we did not have. To confirm the SomaScan® findings obtained for MMP3, IGFBP2, and IGF-1, we used a larger sample size of samples over a longer treatment duration and a wider age range of subjects. The findings from the ELISA runs are summarized in Table 4. Results for MMP-3 validated the SOMAscan® signal (same directionality; p-value for the difference in mean levels = 0.022), IGFBP-2 neared significance (p-value for difference in trajectory slopes = 0.051), while IGF-I did not validate the SOMAscan® signal. Figure 5 shows correlations between SOMAscan® data and ELISA data with the IGF-I measurements having a substantially worse Pearson correlation than MMP-3 and IGFBP-2.

4. Discussion

GCs are and continue to be the standard of care for several chronic inflammatory and auto-immune diseases including DMD [11,21,25]. In this study, and due to the interest and ongoing debates between families and clinicians regarding which GC drug is more beneficial and safe for DMD boys, we focused on a subset of Dfz- an Pred-treated DMD patients to define the differential pharmacodynamic response to these two commonly prescribed drugs using blood circulating proteins. We previously defined a set of PD biomarkers that are responsive to GC treatment in DMD [17]. In that previous study, we showed that GC affected the levels of 107 circulating serum proteins. In general, use of GCs decreased the levels of several circulating pro-inflammatory and immune response associated proteins, but caused an increase in certain proteins associated with metabolism and extracellular remodeling [17]. In this current study, we find that, in general, both Dfz and Pred engaged the PD biomarkers in a similar fashion, however, among the 107 PD biomarkers, 17 exhibited a differential longitudinal response to Dfz relative to Pred, in directionality, mean levels, or both.
A close examination of the 17 differentially altered PD biomarkers (from screening dataset A) between Dfz and Pred in term of their levels, longitudinal directionality and their potential physiological function led to the following interpretation. Among these 17 PD biomarkers, 12 were significantly reduced in the Pred-treated group relative to the Dfz-treated group after adjusting for treatment duration (the samples were age matched). These 12 decreased PD biomarkers by Pred can be classified into subgroups with the first subgroup consisting of inflammatory and immune associated proteins such as LILRB1, TNFRSF21, sRAGE, CD166, and sCD163. Previous studies have suggested that Dfz is a stronger immune suppressant than Pred [43]. In contrast, our analysis showed that for the above-mentioned five serum proteins, Pred seemed to be more immunosuppressive than Dfz. However, further studies using a larger sample size and additional cellular biomarkers are needed to claim this finding. Another subgroup that was differentially decreased by Pred relative to Dfz included bone mineralization protein GPNMB and the cell adhesion protein NCAM-L1. GPNMB, also known as osteoactvin in rats, has been shown to be implicated in bone mineralization and bone regeneration [39,40] and the larger decrease in its mean level in the Pred treated group relative to the Dfz treated group could suggest that Pred treated patients might have a higher risk of losing bone density than Dfz treated patients. This agrees with earlier studies showing that decreases in bone density was markedly spared by Dfz, as compared to Pred in both adult and children cohorts [44]. However, further studies correlating low levels of circulating GPNMB to loss of bone density in Pred treated group relative to Dfz treated group are needed to support this hypothesis. The differential decrease in the circulating levels of NCAM-L1 by Pred could be indicative of another potential side effect of Pred over Dfz. Indeed, a recent study linked low levels of circulating NCAM-L1 to risk of developing type 2 diabetes [41], which is in agreement with an earlier study showing that Pred was more diabetogenic than Dfz in a pediatric population [45]. The differential decrease of sRAGE by Pred over Dfz could also be linked to the diabetogenic effect of Pred since sRAGE has been suggested to act as a decoy that dampens the advanced glycation end product signaling (e.g., RAGE signaling) and thus increases the risk of developing diabetes [33].
However, the differential decrease in the mean levels of MCL-1, MAPK14, and CHRDL1 by Pred relative to Dfz could be considered to be evidence of possible efficacy in DMD. Indeed, MCL-1 and MAPK14 were previously reported by us [17] to be elevated in blood of untreated DMD boys relative to age matched healthy controls then decreased by GC treatment toward the levels in healthy controls. CHRDL1, however, is known to bind to BMP4 and antagonizes its function [30]. BMP4 is a member of the member of the transforming growth factor-β (TGF-β), a well-known pathway reported to be involved in DMD pathogenesis [46]. Thus, a differential decrease of CRDL1 by Pred relative to Dfz could suggest a beneficial effect but further experiments using a CRDL1 knockout animal model are needed to verify this hypothesis.
Another relevant PD biomarker that was differentially altered over time between the two drugs is IGFBP-2. It sharply decreased in its longitudinal trajectory in Dfz-treated group compared to the Pred-treated group. This observation was further validated by ELISA assay (p = 0.0507; near significance). Note that, in our ELISA confirmations, we had the smallest number of additional samples available for IGFBP-2. IGFBP-2 binds to IGF-1 and regulates growth. While IGFBP-2 was decreased by Dfz treatment, IGF-1 was relatively increased by Dfz over Pred. Unfortunately, we were unable to validate the IGF-1 SomaScan® signal using ELISA; nevertheless, an increase in IGF-1 by GC treatment might be associated with antiinflammation efficacy [31]. However, the selective longitudinal decrease of IGFBP2 by Dfz could be associated with significant growth stunting by Dfz over the Pred reported by others [20,47,48] and confirmed by us in this study. Furthermore, a gene knockout study, conducted on mouse model comparing Igfbp2-/- mice colony with Igfbp2+/+ mice colony determined the role of IGFBP2 in bone turnover and showed that Igfbp2-/- males had shorter femurs and were heavier than controls but were not insulin resistant [42]. The decrease in IGFBP2 could have dual side effects on both growth stunting and risk of bone fracture. Indeed, a recent study compared a cohort of DMD boys treated with Dfz to a cohort treated with hydrocortisone and concluded that DMD boys receiving daily dose of Dfz had a higher incidence of bone fractures with greater risk of growth stunting [47].
The second group of steroidal PD biomarkers that were differentially affected by the two drugs are those that increased following treatment with GC. FCN1, IGF-I and PDE3A were relatively increased by Pred over Dfz, while MMP3 and HJV were relatively increased by Dfz over Pred. While FCN1 mean levels were higher in the Pred treated group relative to Dfz treated group, it decreased over time with Pred treatment while remaining relatively unchanged over time following Dfz treatment. This is intriguing and requires further validation using orthogonal methods such as ELISA. Unfortunately, the available ELISA assay required the use of a larger sample volume, which was a limitation of our study.
A comparison of the relative effects of Pred and dexamethasone (Dex), another glucocorticoid, on short-term growth and bone turnover confirmed that both drugs affected short-term bone turnover and growth [20]. However, Dex may be more potent in suppressing linear growth, simulating weight gain and bone turnover compared to prednisone. Dex was more potent at depressing IGF-I levels than prednisolone. This is consistent with our SOMAscan® data, IGF-1 was significantly more depressed or decreased in Dfz treated group relative to Pred treated group. However, we were not able to confirm this finding using ELISA. The discrepancy between the SOMAscan® data and ELISA data could be associated with an epitope effect. The two techniques recognize different epitopes on the IGF-I. Furthermore, IGF-I might co-exist as free form and bound form to circulating insulin growth factor binding proteins and these might further interfere with both these affinity-based assays that relies on epitope binding. Future analyses using denatured condition and targeted mass spectrometry analysis are needed to examine the true levels of circulating IGF-I.
MMP3, also known as stromelysin-1, was dramatically increased in the blood circulation following GC treatment, as previously shown by us [17] and others in both DMD patients [24] and inflammatory bowel disease patients [25] and anti-neutrophil cytoplasmic antibody-associated vasculitis, and juvenile dermatomyositis [25]. In this study, we further show that this increase is associated more with Dfz use than Pred use. This was further confirmed by an ELISA assay on a confirmation dataset using a larger sample size with longer treatment duration. The mechanism by which GC induces levels of circulating MMP3 remains unclear. Further studies using cell culture expressing MMP3 are needed to explain the difference of action of Dfz and Pred on the expression of MMP3. Owing to the function of MMP3, i.e., degradation of extracellular matrix, this differential increase might result in more extracellular matrix remodeling in Dfz-treated DMD patients relative to Pred-treated patients. Whether this extracellular matrix remodeling is adverse or beneficial remain to be carefully examined.

5. Conclusions

Here, we investigated less invasive and objective PD biomarkers that might prove useful to monitor disease progression and response to GC therapies in DMD. We identified differences in serum levels of PD biomarkers between Pred- and Dfz-treated subjects, some of which may be associated with safety. Such blood accessible biomarkers in DMD could play an important role in clinical trials and decision making [49]. In our comparisons, we adjusted for duration of treatment, the comparison groups were age matched, and for a subset of proteins, we conducted ELISA validation testing of SOMAscan® signal. IGFB2 had a decreasing longitudinal trend associated with Dfz, at odds with the observation for Pred. This may be associated with differential growth stunting seen between the two drugs. The dramatic increase of MMP3 by Dfz relative to Pred remains to be interpreted. Further studies are needed to test the physiological significance of these Dfz and Pred differentially affected PD biomarkers. The study’s limitations include the small sample size, no adjustments for multiple testing, and that the data come from a natural history study with lots of observed variability. We have tried to overcome these limitations by performing lab validation using ELISA and using a larger sample size and longer treatment duration. For biological validation, an external cohort is needed to validate our findings. Note also that we did not investigate direct associations of biomarkers with efficacy/safety outcomes; larger sample sizes are required for such studies and are the objective of ongoing research. While preliminary, this study identified serum proteins that were altered between the Dfz and Pred groups using SOMAscan® array data and we successfully validated two biomarkers using ELISA that may be associated with adverse effects. However, further studies using larger sample sizes collected from well controlled cohorts enrolled in ongoing and future clinical trials comparing the safety and efficacy of Dfz versus Pred are needed to validate and test these differentially altered pharmacodynamic biomarkers identified herein.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4426/10/4/164/s1, Supplemental Table S1: A subset of DMD patients from the Cooperative International Neuromuscular Research Group Duchenne Natural History Study (CINRG-DNHS).

Author Contributions

Conceptualization, Y.H. and U.J.D.; Methodology, Y.H., U.J.D. and S.T.; Validation, S.T. and A.S.; Formal Analysis, S.T.; Data Curation, Y.H., U.J.D., H.G.-D. and S.T.; Writing—Original Draft Preparation, S.T.; Writing—Review & Editing, Y.H. and U.J.D.; Funding Acquisition, Y.H. The CINRG DNHS study was run by the CINRG DNHS investigators who also collected the bio-samples. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by DoD grant W81XWH-16-1-0557, NIH grant P50HD090254 and the Muscular Dystrophy Association USA (MDA353094).

Acknowledgments

The authors would like to thank the patients and their families for their participation in the study. The authors would also like to thank study team members for our participating CINRG sites: University of California Davis: C. McDonald, E. Henricson, M. Cregan, L. Johnson, J. Han, N. Joyce, A. Nicorici, D. Reddy; Sundaram Medical Foundation and Apollo Children’s Hospital, Alberta Children’s Hospital, Calgary: J. Mah, A. Chiu, T. Haig, M. Harris, M. Kornelsen, N. Rincon, K. Sanchez, L. Walker; Queen Silvia Children’s Hospital: M. Tulinius, A. Alhander, A. Ekstrom, A. Gustafsson, A. Kroksmark, U. Sterky, L. Wahlgren; Children’s National Health System: R. Leshner, N. Brody, B. Drogo, M. Leach, C. Tesi-Rocha, M. Birkmeier, B. Tadese, A. Toles and M. Thangarajh; Royal Children’s Hospital: A. Kornberg, K. Carroll, K. DeValle, R. Kennedy, V. Rodriguez, D. Villano; Hadassah Hebrew University Hospital: Y. Nevo, R. Adani, A. Bar Leve, L. Chen-Joseph, M. Daana, V. Panteleyev-Yitshak, E. Simchovitz, D. Yaffe; Instituto de Neurosciencias Fundacion Favaloro: L. Andreone, F. Bonaudo, J. Corderi, L. Levi, L. Mesa, P. Marco; Children’s Hospital of Pittsburgh of UPMC and the University of Pittsburgh: P. Clemens, H. Abdel-Hamid, R. Bendixen, C. Bise, A. Craig, K. Karnavas, C. Matthews, G. Niizawa, A. Smith, J. Weimer; Washington University: J. Anger, T. Christenson, J. Florence, R. Gadeken, P. Golumbak, B. Malkus, A. Pestronk, R. Renna, J. Schierbecker, C. Seiner, C. Wulf; Children’s Hospital of Richmond at VCU: J. Teasley, S. Blair, B. Grillo, E. Monasterio; University of Tennessee: T. Bertorini, M. Barrett-Adair, C. Benzel, K. Carter, J. Clift, B. Gatlin, R. Henegar, J. Holloway, M. Igarashi, F. Kiphut, A. Parker, A. Phillips, R. Young; Children’s Hospital of Westmead: K. North, K. Cornett, N. Gabriel, M. Harman, C. Miller, K. Rose, S. Wicks; University of Alberta: H. Kolski, L. Chen, C. Kennedy; Centro Clinico Nemo: M. Beneggi, L. Capone, A. Molteni, V. Morettini; Texas Children’s Hospital: T. Lotze, A. Gupta, A. Knight, B. Lott, R. McNeil, G. Orozco, R. Schlosser; University of Minnesota: G. Chambers, J. Day, J. Dalton, A. Erickson, M. Margolis, J. Marsh, C. Naughton; Mayo Clinic: K. Coleman-Wood, A. Hoffman, W. Korn-Petersen, N. Kuntz; University of Puerto Rico: B. Deliz, S. Espada, P. Fuste, C. Luciano, J. Torres; and the CINRG Coordinating Center: Lauren Morgenroth, M. Ahmed, A. Arrieta, N. Bartley, T. Brown-Caines, C. Carty, T. Duong, J. Feng, F. Hu, L. Hunegs, Z. Sund, W. Tang, and A. Zimmerman. The CINRG-DNHS was funded by the U.S. Department of Education/NIDRR (#H133B031118, #H133B090001); U.S. Department of Defense (#W81XWH-12-1-0417); National Institutes of Health/NIAMS (#R01AR061875); Parent Project Muscular Dystrophy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoffman, E.P.; Brown, R.H.; Kunkel, L.M. Dystrophin: The protein product of the Duchenne muscular dystrophy locus. Cell 1987, 51, 919–928. [Google Scholar] [CrossRef]
  2. Wallace, G.Q.; McNally, E.M. Mechanisms of muscle degeneration, regeneration, and repair in the muscular dystrophies. Annu. Rev. Physiol. 2009, 71, 37–57. [Google Scholar] [CrossRef] [PubMed]
  3. Ryan, N.J. Ataluren: First global approval. Drugs 2014, 74, 1709–1714. [Google Scholar] [CrossRef] [PubMed]
  4. Mendell, J.R.; Rodino-Klapac, L.R.; Sahenk, Z.; Roush, K.; Bird, L.; Lowes, L.P.; Alfano, L.; Gomez, A.M.; Lewis, S.; Kota, J.; et al. Eteplirsen for the treatment of Duchenne muscular dystrophy. Ann. Neurol. 2013, 74, 637–647. [Google Scholar] [CrossRef]
  5. Syed, Y.Y. Eteplirsen: First global approval. Drugs 2016, 76, 1699–1704. [Google Scholar] [CrossRef]
  6. Heo, Y.A. Golodirsen: First approval. Drugs 2020, 80, 329–333. [Google Scholar] [CrossRef]
  7. Frank, D.E.; Schnell, F.J.; Akana, C.; El-Husayni, S.H.; Desjardins, C.A.; Morgan, J.; Charleston, J.S.; Sardone, V.; Domingos, J.; Dickson, G.; et al. Increased dystrophin production with golodirsen in patients with Duchenne muscular dystrophy. Neurology 2020, 94, e2270–e2282. [Google Scholar] [CrossRef] [Green Version]
  8. Hoffman, E.P.; Bronson, A.; Levin, A.A.; Takeda, S.; Yokota, T.; Baudy, A.R.; Connor, E.M. Restoring dystrophin expression in Duchenne muscular dystrophy muscle: Progress in exon skipping and stop codon read through. Am. J. Pathol. 2011, 179, 12–22. [Google Scholar] [CrossRef] [Green Version]
  9. Bushby, K.; Finkel, R.; Birnkrant, D.J.; Case, L.E.; Clemens, P.R.; Cripe, L.; Kaul, A.; Kinnett, K.; McDonald, C.; Pandya, S.; et al. Diagnosis and management of Duchenne muscular dystrophy, part 1: Diagnosis, and pharmacological and psychosocial management. Lancet Neurol. 2010, 9, 77–93. [Google Scholar] [CrossRef]
  10. Griggs, R.C.; Miller, J.P.; Greenberg, C.R.; Fehlings, D.L.; Pestronk, A.; Mendell, J.R.; Moxley, R.T.; King, W.; Kissel, J.T.; Cwik, V.; et al. Efficacy and safety of deflazacort vs prednisone and placebo for Duchenne muscular dystrophy. Neurology 2016, 87, 2123–2131. [Google Scholar] [CrossRef] [Green Version]
  11. McDonald, C.M.; Henricson, E.K.; Abresch, R.T.; Duong, T.; Joyce, N.C.; Hu, F.; Clemens, P.R.; Hoffman, E.P.; Cnaan, A.; Gordish-Dressman, H.; et al. Long-term effects of glucocorticoids on function, quality of life, and survival in patients with Duchenne muscular dystrophy: A prospective cohort study. Lancet 2018, 391, 451–461. [Google Scholar] [CrossRef]
  12. Bello, L.; Gordish-Dressman, H.; Morgenroth, L.P.; Henricson, E.K.; Duong, T.; Hoffman, E.P.; Cnaan, A.; McDonald, C.M. Prednisone/prednisolone and deflazacort regimens in the CINRG Duchenne Natural History Study. Neurology 2015, 85, 1048–1055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lamb, M.M.; West, N.A.; Ouyang, L.; Yang, M.; Weitzenkamp, D.; James, K.; Ciafaloni, E.; Pandya, S.; DiGuiseppi, C.; Cunniff, C.; et al. Corticosteroid treatment and growth patterns in ambulatory males with Duchenne muscular dystrophy. J. Pediatr. 2016, 173, 207–213.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bonifati, M.D.; Ruzza, G.; Bonometto, P.; Berardinelli, A.; Gorni, K.; Orcesi, S.; Lanzi, G.; Angelini, C. A multicenter, double-blind, randomized trial of deflazacort versus prednisone in Duchenne muscular dystrophy. Muscle Nerve 2000, 23, 1344–1347. [Google Scholar] [CrossRef]
  15. Cohen, L.; Morgan, J.; Bozyk, M. Variable effects of corticosteroid treatment of serum enzyme activities in Duchenne’s muscular dystrophy. Res. Commun. Chem. Pathol. Pharmacol. 1977, 17, 529–538. [Google Scholar]
  16. FDA News Release (2017) FDA Approves Drug to Treat Duchenne Muscular Dystrophy. Available online: https://www.fda.gov/news-events/press-announcements/fda-approves-drug-treat-duchenne-muscular-dystrophy (accessed on 1 September 2020).
  17. Hathout, Y.; Liang, C.; Ogundele, M.; Xu, G.; Tawalbeh, S.M.; Dang, U.J.; Hoffman, E.P.; Gordish-Dressman, H.; Conklin, L.S.; van den Anker, J.N.; et al. Disease-specific and glucocorticoid-responsive serum biomarkers for Duchenne muscular dystrophy. Sci. Rep. 2019, 9, 12167. [Google Scholar] [CrossRef]
  18. Dang, U.J.; Ziemba, M.; Clemens, P.R.; Hathout, Y.; Conklin, L.S.; Hoffman, E.P. Serum biomarkers associated with baseline clinical severity in young steroid-naïve Duchenne muscular dystrophy boys. Hum. Mol. Genet. 2020, 29, 2481–2495. [Google Scholar] [CrossRef]
  19. Spitali, P.; Hettne, K.; Tsonaka, R.; Charrout, M.; van den Bergen, J.; Koeks, Z.; Kan, H.E.; Hooijmans, M.T.; Roos, A.; Straub, V.; et al. Tracking disease progression non-invasively in Duchenne and Becker muscular dystrophies. J. Cachexia. Sarcopenia Muscle 2018, 9, 715–726. [Google Scholar] [CrossRef]
  20. Ahmed, S.F.; Tucker, P.; Mushtaq, T.; Wallace, A.M.; Williams, D.M.; Hughes, I.A. Short-term effects on linear growth and bone turnover in children randomized to receive prednisolone or dexamethasone. Clin. Endocrinol. 2002, 57, 185–191. [Google Scholar] [CrossRef]
  21. Matthews, E.; Brassington, R.; Kuntzer, T.; Jichi, F.; Manzur, A.Y. Corticosteroids for the treatment of Duchenne muscular dystrophy. Cochrane Database Syst. Rev. 2016. [Google Scholar] [CrossRef] [Green Version]
  22. McDonald, C.M.; Sajeev, G.; Yao, Z.; McDonnell, E.; Elfring, G.; Souza, M.; Peltz, S.W.; Darras, B.T.; Shieh, P.B.; Cox, D.A.; et al. Deflazacort vs prednisone treatment for Duchenne muscular dystrophy: A meta-analysis of disease progression rates in recent multicenter clinical trials. Muscle Nerve. 2020, 61, 26–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ganapati, A.; Ravindran, R.; David, T.; Yadav, B.; Jeyaseelan, V.; Jeyaseelan, L.; Danda, D. Head to head comparison of adverse effects and efficacy between high dose deflazacort and high dose prednisolone in systemic lupus erythematosus: A prospective cohort study. Lupus 2018, 27, 890–898. [Google Scholar] [CrossRef] [PubMed]
  24. Hathout, Y.; Conklin, L.S.; Seol, H.; Gordish-Dressman, H.; Brown, K.J.; Morgenroth, L.P.; Nagaraju, K.; Heier, C.R.; Damsker, J.M.; van den Anker, J.N.; et al. Serum pharmacodynamic biomarkers for chronic corticosteroid treatment of children. Sci. Rep. 2016, 6, 31727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Conklin, L.S.; Merkel, P.A.; Pachman, L.M.; Parikh, H.; Tawalbeh, S.; Damsker, J.M.; Cuthbertson, D.D.; Morgan, G.A.; Monach, P.A.; Hathout, Y.; et al. Serum biomarkers of glucocorticoid response and safety in anti-neutrophil cytoplasmic antibody-associated vasculitis and juvenile dermatomyositis. Steroids 2018, 140, 159–166. [Google Scholar] [CrossRef]
  26. McDonald, C.M.; Henricson, E.K.; Abresch, R.T.; Han, J.J.; Escolar, D.M.; Florence, J.M.; Duong, T.; Arrieta, A.; Clemens, P.R.; Hoffman, E.P.; et al. The cooperative international neuromuscular research group duchenne natural history study-a longitudinal investigation in the era of glucocorticoid therapy: Design of protocol and the methods used. Muscle Nerve 2013, 48, 32–54. [Google Scholar] [CrossRef] [PubMed]
  27. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  28. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  29. Kuznetsova, A.; Brockhoff, P.B.; Christensen, R.H.B. lmerTest package: Tests in linear mixed effects models. J. Stat. Softw. 2017, 82, 1–26. [Google Scholar] [CrossRef] [Green Version]
  30. Lin, J.; Patel, S.R.; Cheng, X.; Cho, E.A.; Levitan, I.; Ullenbruch, M.; Phan, S.H.; Park, J.M.; Dressler, G.R. Kielin/chordin-like protein, a novel enhancer of BMP signaling, attenuates renal fibrotic disease. Nat. Med. 2005, 11, 387–393. [Google Scholar] [CrossRef]
  31. Zatorski, H.; Marynowski, M.; Fichna, J. Is insulin-like growth factor 1 (IGF-1) system an attractive target inflammatory bowel diseases? Benefits and limitation of potential therapy. Pharmacol. Rep. 2016, 68, 809–815. [Google Scholar] [CrossRef]
  32. Serrano, A.L.; Muñoz-Cánoves, P. Fibrosis development in early-onset muscular dystrophies: Mechanisms and translational implications. Semin. Cell Dev. Biol. 2017, 64, 181–190. [Google Scholar] [CrossRef]
  33. Selvin, E.; Halushka, M.K.; Rawlings, A.M.; Hoogeveen, R.C.; Ballantyne, C.M.; Coresh, J.; Astor, B.C. SRAGE and risk of diabetes, cardiovascular disease, and death. Diabetes 2013, 62, 2116–2121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Tanida, S.; Mizoshita, T.; Ozeki, K.; Katano, T.; Kataoka, H.; Kamiya, T.; Joh, T. Advances in refractory ulcerative colitis treatment: A new therapeutic target, Annexin A2. World J. Gastroenterol. 2015, 21, 8776–8786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Liu, N.; Jiang, Y.; Chung, J.Y.; Li, Y.; Yu, Z.; Kim, J.W.; Lok, J.M.; Whalen, M.J.; Wang, X. Annexin A2 deficiency exacerbates neuroinflammation and long-term neurological deficits after traumatic brain injury in mice. Int. J. Mol. Sci. 2019, 20, 6125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Defour, A.; Medikayala, S.; Van der Meulen, J.H.; Hogarth, M.W.; Holdreith, N.; Malatras, A.; Duddy, W.; Boehler, J.; Nagaraju, K.; Jaiswal, J.K. Annexin A2 links poor myofiber repair with inflammation and adipogenic replacement of the injured muscle. Hum. Mol. Genet. 2017, 26, 1979–1991. [Google Scholar] [CrossRef] [Green Version]
  37. Grønbæk, H.; Kreutzfeldt, M.; Kazankov, K.; Jessen, N.; Sandahl, T.; Hamilton-Dutoit, S.; Vilstrup, H.; Møller, H.J. Single-centre experience of the macrophage activation marker soluble (s)CD163–associations with disease activity and treatment response in patients with autoimmune hepatitis. Aliment. Pharmacol. Ther. 2016, 44, 1062–1070. [Google Scholar]
  38. Espinasse, M.-A.; Pépin, A.; Virault-Rocroy, P.; Szely, N.; Chollet-Martin, S.; Pallardy, M.; Biola-Vidamment, A. Glucocorticoid-induced leucine zipper is expressed in human neutrophils and promotes apoptosis through Mcl-1 down-regulation. J. Innate Immun. 2016, 8, 81–96. [Google Scholar] [CrossRef]
  39. Hu, X.; Zhang, P.; Xu, Z.; Chen, H.; Xie, X. GPNMB enhances bone regeneration by promoting angiogenesis and osteogenesis: Potential role for tissue engineering bone. J. Cell. Biochem. 2013, 114, 2729–2737. [Google Scholar] [CrossRef]
  40. Huang, Y.; Bai, B.; Yao, Y. Prospects of osteoactivin in tissue regeneration. Expert Opin. Ther. Targets 2016, 20, 1357–1364. [Google Scholar] [CrossRef]
  41. Belongie, K.J.; Ferrannini, E.; Johnson, K.; Andrade-Gordon, P.; Hansen, M.K.; Petrie, J.R. Identification of novel biomarkers to monitor β-cell function and enable early detection of type 2 diabetes risk. PLoS ONE 2017, 12, e0182932. [Google Scholar] [CrossRef]
  42. DeMambro, V.E.; Clemmons, D.R.; Horton, L.G.; Bouxsein, M.L.; Wood, T.L.; Beamer, W.G.; Canalis, E.; Rosen, C.J. Gender-specific changes in bone turnover and skeletal architecture in Igfbp-2-null mice. Endocrinology 2008, 149, 2051–2061. [Google Scholar] [CrossRef]
  43. Hahn, B.H.; MacDermott, R.P.; Jacobs, S.B.; Pletscher, L.S.; Beale, M.G. Immunosuppressive effects of low doses of glucocorticoids: Effects on autologous and allogeneic mixed leukocyte reactions. J. Immunol. 1980, 124, 2812–2817. [Google Scholar] [PubMed]
  44. Lippuner, K.; Casez, J.-P.; Horber, F.F.; Jaeger, P. Effects of deflazacort versus prednisone on bone mass, body composition, and lipid profile: A randomized, double blind study in kidney transplant patients. J. Clin. Endocrinol. Metab. 1998, 83, 3795–3802. [Google Scholar] [CrossRef] [PubMed]
  45. Pagano, G.; Lombardi, A.; Ferraris, G.M.; Imbimbo, B.; Cavallo Perin, P. Acute effect of prednisone and deflazacort on glucose tolerance in prediabetic subjects. Eur. J. Clin. Pharmacol. 1982, 22, 469–471. [Google Scholar] [CrossRef] [PubMed]
  46. Kemaladewi, D.U.; Pasteuning, S.; Van Der Meulen, J.W.; Van Heiningen, S.H.; Van Ommen, G.J.; Ten Dijke, P.; Aartsma-Rus, A.; ’T Hoen, P.A.C.; Hoogaars, W.M. Targeting TGF-β signaling by antisense oligonucleotide-mediated knockdown of TGF-β type I receptor. Mol. Ther.-Nucleic Acids 2014, 3, e156. [Google Scholar] [CrossRef] [PubMed]
  47. Joseph, S.; Wang, C.; Bushby, K.; Guglieri, M.; Horrocks, I.; Straub, V.; Ahmed, S.F.; Wong, S.C. Fractures and linear growth in a nationwide cohort of boys with Duchenne muscular dystrophy with and without glucocorticoid treatment: Results from the UK Northstar database. JAMA Neurol. 2019, 76, 701–709. [Google Scholar] [CrossRef] [PubMed]
  48. Marden, J.R.; Freimark, J.; Yao, Z.; Signorovitch, J.; Tian, C.; Wong, B.L. Real-world outcomes of long-term prednisone and deflazacort use in patients with Duchenne muscular dystrophy: Experience at a single, large care center. J. Comp. Eff. Res. 2020, 9, 177–189. [Google Scholar] [CrossRef]
  49. Conklin, L.S.; Damsker, J.M.; Hoffman, E.P.; Jusko, W.J.; Mavroudis, P.D.; Schwartz, B.D.; Mengle-Gaw, L.J.; Smith, E.C.; Mah, J.K.; Guglieri, M.; et al. Phase IIa trial in Duchenne muscular dystrophy shows vamorolone is a first-in-class dissociative steroidal anti-inflammatory drug. Pharmacol. Res. 2018, 136, 140–150. [Google Scholar] [CrossRef]
Figure 1. 2D structures of Prednisone and Deflazacort (top panels) and their active drugs: prednisolone and 21-desacetyl deflazacort. The similarities are clear. The red arrows point to the structural differences between the active drugs.
Figure 1. 2D structures of Prednisone and Deflazacort (top panels) and their active drugs: prednisolone and 21-desacetyl deflazacort. The similarities are clear. The red arrows point to the structural differences between the active drugs.
Jpm 10 00164 g001
Figure 2. Schematic describing the workflow for the biomarker analyses. Proteins identified as differentially affected by prednisone vs. deflazacort were identified from the SomaScan® based screening dataset. This was followed by confirmation analysis using ELISA assays of three selected biomarkers on a larger set of subjects (confirmation dataset). n: number of subjects, Ns: number of serum samples used in this study.
Figure 2. Schematic describing the workflow for the biomarker analyses. Proteins identified as differentially affected by prednisone vs. deflazacort were identified from the SomaScan® based screening dataset. This was followed by confirmation analysis using ELISA assays of three selected biomarkers on a larger set of subjects (confirmation dataset). n: number of subjects, Ns: number of serum samples used in this study.
Jpm 10 00164 g002
Figure 3. Selected examples of serum protein PD biomarkers showing difference between prednisone (Pred) and deflazacort (Dfz)-treated samples. FCN-1 has relatively higher mean RFU levels in Dfz-treated patients (p = 0.038 for mean levels; p = 0.03 for difference in trajectory slopes over time). TNFRSF21 has higher mean RFU levels in Pred- vs. Dfz-treated DMD patients (p = 0.005 mean levels; p = 0.934 for difference in trajectory slopes over time). The longitudinal trajectories of IGFBP-2 levels are different between the two drugs (p = 0.04). MMP-3 mean RFU level is elevated in DMD boys treated with Dfz, as compared to Pred (p = 0.008).
Figure 3. Selected examples of serum protein PD biomarkers showing difference between prednisone (Pred) and deflazacort (Dfz)-treated samples. FCN-1 has relatively higher mean RFU levels in Dfz-treated patients (p = 0.038 for mean levels; p = 0.03 for difference in trajectory slopes over time). TNFRSF21 has higher mean RFU levels in Pred- vs. Dfz-treated DMD patients (p = 0.005 mean levels; p = 0.934 for difference in trajectory slopes over time). The longitudinal trajectories of IGFBP-2 levels are different between the two drugs (p = 0.04). MMP-3 mean RFU level is elevated in DMD boys treated with Dfz, as compared to Pred (p = 0.008).
Jpm 10 00164 g003
Figure 4. Comparison of longitudinal trajectories between DMD patients treated with deflazacort and prednisone of BMI (p = 0.08 for difference in trajectory slopes), height (cm; p = 0.006), and weight (kg; p = 0.112).
Figure 4. Comparison of longitudinal trajectories between DMD patients treated with deflazacort and prednisone of BMI (p = 0.08 for difference in trajectory slopes), height (cm; p = 0.006), and weight (kg; p = 0.112).
Jpm 10 00164 g004
Figure 5. SOMAscan® signal confirmation using ELISA. Upper panel shows longitudinal trajectories of selected DMD serum protein PD biomarkers from ELISA assays. Lower panel shows correlation plots between SOMA and ELISA data for the samples that overlapped between the screening and confirmation data sets.
Figure 5. SOMAscan® signal confirmation using ELISA. Upper panel shows longitudinal trajectories of selected DMD serum protein PD biomarkers from ELISA assays. Lower panel shows correlation plots between SOMA and ELISA data for the samples that overlapped between the screening and confirmation data sets.
Jpm 10 00164 g005
Table 1. Characteristics of patients/samples for SOMAscan® screening dataset A. All samples were from ambulatory patients.
Table 1. Characteristics of patients/samples for SOMAscan® screening dataset A. All samples were from ambulatory patients.
TreatmentNumber of PatientsAverage Number of Visits (Min, Max)Age Range at Sample Collection (Years)Average Time between Biosample Collection (Days) (Mix, Max)Regimen
Deflazacort82.5 (2, 4)4.7–9.4465 (56, 1268)Daily on Dfz
Prednisone72.14 (2, 3)4.3–8.3582 (157, 1392)Daily on Pred
Table 2. Characteristics of patients/samples for confirmation dataset B. These serum samples were used for ELISA confirmation of MMP3, IGFBP-2 and IGF-1. All the samples were from ambulatory patients.
Table 2. Characteristics of patients/samples for confirmation dataset B. These serum samples were used for ELISA confirmation of MMP3, IGFBP-2 and IGF-1. All the samples were from ambulatory patients.
TreatmentNumber of PatientsAverage Number of Visits (Min, Max)Mean Age at Sample Collection in Years (Min-Max)Average Time between Biosample Collection (Days) (Mix, Max)
Deflazacort172.7 (2, 5)9 (4.7–15.3)1392 (370, 3058)
Prednisone122.5 (2, 3)8.5 (4.2–15.8)1685 (594, 3391)
Table 3. List of PD biomarkers that were differentially altered between Pred-treated (n = 7) and Dfz-treated (n = 8) groups (SOMAscan® screening dataset A).
Table 3. List of PD biomarkers that were differentially altered between Pred-treated (n = 7) and Dfz-treated (n = 8) groups (SOMAscan® screening dataset A).
Abbreviated Gene Name 1 (Uniprot ID)Fold Change between Untreated DMD vs. Healthy Controls with p-Value from [17]Fold-Change between DMD Subjects Treated with Pred and Dfzp-Value: Difference in Mean Levels 2p-Value: Difference in Trajectories 2Protein Function—Biological ProcessPotential Significance
LILRB1 (Q8NHL6)0.96 (0.861)−1.5 (↓) in Pred vs. Dfz0.005 0.426Immune responseSide effect (immune suppression)
TNFRSF21 (O75509)1.14 (0.454)−1.33 (↓) in Pred vs. Dfz0.005 0.934Apoptotic process, adaptive immune responseSide effect (immune suppression)
CHRDL1 (Q9BU40)1.30 (0.039)−1.2 (↓) in Pred vs. Dfz0.007 0.455Bone Morphogenetic Proteins (BMP) signaling pathwayEfficacy via action on TGF-β signaling [30]
IGF-I (P05019)0.85 (0.161)1.14 (↑) in Pred vs. Dfz0.007 0.096Promotes growthPotential efficacy associated with anti-inflammatory propriety [31]
MMP-3 (P08254)0.77 (0.488)2 (↑) in Dfz vs. Pred0.008 0.216Extracellular matrix degradationEfficacy/Side effect [17,25,32]
sRAGE/AGER (Q15109)0.60 (0.005)−1.82 (↓) in Pred vs. Dfz0.010 0.640Inflammatory (causes complications in diabetes) Side effect may be associated with diabetes risk [33]
ANXA2 (P07355)1.24 (0.246)−1.22 (↓) in Pred vs. Dfz0.011 0.054Angiogenesis, biomineral tissue development, inflammationPotential efficacy marker associated with inflammation [34,35,36]
CD166 (Q13740)0.81 (0.088)−1.25 (↓) in Pred vs. Dfz0.014 0.616Cell adhesion, adaptive immune response Side effect (immune suppression)
HJV (Q6ZVN8)0.83 (0.073)1.19 (↑) in Dfz vs. Pred0.017 0.385BMP signaling pathway, iron ion homeostasis Unknown
sCD163 (Q86VB7)0.86 (0.466)−1.30 (↓) in Pred vs. Dfz0.025 0.340InflammationEfficacy [37]
Mcl-1 (Q07820)2.31 (<0.001)−1.33 (↓) in Pred vs. Dfz0.029 0.146Apoptosis, DifferentiationInflammationEfficacy [38]
PDE3A (Q14432)1.03 (0.805)1.22 (↑) in Pred vs. Dfz0.033 0.956Cell to cell signalingUnknown
GPNMB (Q14956)0.79 (0.111)−1.32 (↓) in Pred vs. Dfz0.033 0.944Cell adhesion, bone mineralization. Side effect related to bone [39,40]
FCN1 (O00602)1.02 (0.905)1.18 (↑) in Pred vs. Dfz0.038 0.030Innate immune responseUnknown
MAPK14 (Q16539)2.10 (<0.001)−1.18 (↓) in Pred vs. Dfz0.041 0.982Potential muscle injury biomarkerEfficacy, muscle injury protein that normalized after GC treatment [17]
NCAM-L1 (P32004)0.76 (0.072)−1.28 (↓) in Pred vs. Dfz0.042 0.691Cell adhesion and differentiationSide effect, risk of developing diabetes [41]
IGFBP-2 (P18065)2.55 (<0.001)−1.16 (↓) in Dfz vs. Pred0.3290.040Growth regulation Side effect, growth stunting [42]
1 LILRB1 = Leukocyte immunoglobulin-like receptor subfamily B member 1; TNFRSF21 = Tumor necrosis factor receptor superfamily member 21; CHRDL1 = Chordin-like protein 1; IGF-I = Insulin-like growth factor I; MMP-3 = Stromelysin-1; sRAGE/AGER = Advanced glycosylation end product-specific receptor, soluble; ANXA2 = Annexin A2; CD166 = CD166 antigen; HJV = Hemojuvelin; sCD163 = Scavenger receptor cysteine-rich type 1 protein M130; Mcl-1 = Induced myeloid leukemia cell differentiation protein; PDE3A = cGMP-inhibited 3′,5′-cyclic phosphodiesterase A; GPNMB = Transmembrane glycoprotein NMB; FCN1 = Ficolin-1; MAPK14 = Mitogen-activated protein kinase 14; NCAM-L1 = Neural cell adhesion molecule L1; IGFBP-2 = Insulin-like growth factor-binding protein 2. 2 The two p-values are for the difference in mean levels and difference in longitudinal trajectories between Pred-treated and Dfz-treated subjects, respectively.
Table 4. Summary for biomarker signal confirmation using ELISA.
Table 4. Summary for biomarker signal confirmation using ELISA.
Protein Name (Uniprot ID)SOMAscan® DataELISA DataFunction
p-Value 1: Difference in Mean Levelsp-Value 1: Difference in TrajectoriesNumber of Patients/Samplesp-Value 1: Difference in Mean Levelsp-Value 1: Difference in TrajectoriesNumber of Patients/Samples
MMP-3 (P08254)0.0080.2168 Dfz, 7 Pred/35 samples0.0220.37817 Dfz, 12 Pred/76 SamplesExtracellular matrix degradation
IGFBP-2 (P18065)0.3280.048 Dfz, 7 Pred/35 samples0.7440.050710 Dfz, 10 Pred/49 SamplesRegulates growth
IGF-I (P05019)0.0070.0968 Dfz, 7 Pred/35 samples0.2460.13717 Dfz, 12 Pred/75 SamplesPromotes growth
1 Two p-values are provided for difference in mean levels and difference in longitudinal trajectory slopes between Pred-treated and Dfz-treated subjects for both SOMAscan® and ELISA® analysis, respectively.

Share and Cite

MDPI and ACS Style

Tawalbeh, S.; Samsel, A.; Gordish-Dressman, H.; Hathout, Y.; CINRG-DNHS Investigators; Dang, U.J. Comparison of Serum Pharmacodynamic Biomarkers in Prednisone-Versus Deflazacort-Treated Duchenne Muscular Dystrophy Boys. J. Pers. Med. 2020, 10, 164. https://doi.org/10.3390/jpm10040164

AMA Style

Tawalbeh S, Samsel A, Gordish-Dressman H, Hathout Y, CINRG-DNHS Investigators, Dang UJ. Comparison of Serum Pharmacodynamic Biomarkers in Prednisone-Versus Deflazacort-Treated Duchenne Muscular Dystrophy Boys. Journal of Personalized Medicine. 2020; 10(4):164. https://doi.org/10.3390/jpm10040164

Chicago/Turabian Style

Tawalbeh, Shefa, Alison Samsel, Heather Gordish-Dressman, Yetrib Hathout, CINRG-DNHS Investigators, and Utkarsh J. Dang. 2020. "Comparison of Serum Pharmacodynamic Biomarkers in Prednisone-Versus Deflazacort-Treated Duchenne Muscular Dystrophy Boys" Journal of Personalized Medicine 10, no. 4: 164. https://doi.org/10.3390/jpm10040164

APA Style

Tawalbeh, S., Samsel, A., Gordish-Dressman, H., Hathout, Y., CINRG-DNHS Investigators, & Dang, U. J. (2020). Comparison of Serum Pharmacodynamic Biomarkers in Prednisone-Versus Deflazacort-Treated Duchenne Muscular Dystrophy Boys. Journal of Personalized Medicine, 10(4), 164. https://doi.org/10.3390/jpm10040164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop