Quantitative structure-property relationship (QSPR) analysis is a promising approach to correlate structural features with properties of glass materials that lack long-range order and usually have complex structures. By using carefully chosen descriptors based on structural models generated from molecular dynamics (MD) simulations, correlations with properties and insights on glass behaviors can be obtained. Zirconia can significantly alter glass properties including chemical durability, even in a small amount, and hence plays an important role in vitrification of nuclear waste where long-term chemical durability is desired. In this study, borosilicate glasses with the composition of xZrO2-(61 - x)SiO2-17B2O3-18Na2O-4CaO with x = 0, 1, 2, 4, 6, and 8 were simulated using classical MD simulations with the recently developed composition-dependent potentials. Short-range (e.g., bond distance and coordination numbers) and medium-range (e.g., Qn distribution, network connectivity, and ring-size distribution) structural features altered by ZrO2 were obtained and analyzed. The use of a descriptor ( Fnet descriptor) that combines short-range structural characteristics, from MD simulations, and the cation-oxygen single bond strength was found to provide excellent linear correlations with the density and initial dissolution rate of these glasses. The results show that by combining MD simulations and QSPR analysis the composition and structural effect on the properties of complex multicomponent glasses can be elucidated, thus suggesting that this is a promising approach for future glass research and new composition design.