Abstract Duchenne muscular dystrophy (DMD) is caused by loss of dystrophin in muscle, and while all patients share the primary gene and biochemical defect, there is considerable patient–patient variability in clinical symptoms. We sought to develop multivariate models of serum protein biomarkers that explained observed variation, using functional outcome measures as proxies for severity. Serum samples from 39 steroid-naïve DMD boys 4 to <7 years enrolled into a clinical trial of vamorolone were studied (NCT02760264). Four assessments of gross motor function were carried out for each participant over a 6-week interval, and their mean was used as response for biomarker models. Weighted correlation network analysis was used for unsupervised clustering of 1305 proteins quantified using SOMAscan® aptamer profiling to define highly representative and connected proteins. Multivariate models of biomarkers were obtained for time to stand performance (strength phenotype; 17 proteins) and 6 min walk performance (endurance phenotype; 17 proteins) including some shared proteins. Identified proteins were tested with associations of mRNA expression with histological severity of muscle from dystrophinopathy patients (n = 28) and normal controls (n = 6). Strong associations predictive of both clinical and histological severity were found for ERBB4 (reductions in both blood and muscle with increasing severity), SOD1 (reductions in muscle and increases in blood with increasing severity) and CNTF (decreased levels in blood and muscle with increasing severity). We show that performance of DMD boys was effectively modeled with serum proteins, proximal strength associated with growth and remodeling pathways and muscle endurance centered on TGFβ and fibrosis pathways in muscle.