This study investigates the effectiveness of hundreds of texture features extracted from voxel-based dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) parametric maps for early prediction of breast cancer response to neoadjuvant chemotherapy (NAC). In total, 38 patients with breast cancer underwent DCE-MRI before (baseline) and after the first of the 6-8 NAC cycles. Quantitative pharmacokinetic (PK) parameters and semiquantitative metrics were estimated from DCE-MRI time-course data. The residual cancer burden (RCB) index value was computed based on pathological analysis of surgical specimens after NAC completion. In total, 1043 texture features were extracted from each of the 13 parametric maps of quantitative PK or semiquantitative metric, and their capabilities for early prediction of RCB were examined by correlating feature changes between the 2 MRI studies with RCB. There were 1069 pairs of feature-map combinations that showed effectiveness for response prediction with 4 correlation coefficients >0.7. The 3-dimensional gray-level cooccurrence matrix was the most effective feature extraction method for therapy response prediction, and, in general, the statistical features describing texture heterogeneity were the most effective features. Quantitative PK parameters, particularly those estimated with the shutter-speed model, were more likely to generate effective features for prediction response compared with the semiquantitative metrics. The best feature-map pair could predict pathologic complete response with 100% sensitivity and 100% specificity using our cohort. In conclusion, breast tumor heterogeneity in microvasculature as measured by texture features of voxel-based DCE-MRI parametric maps could be a useful biomarker for early prediction of NAC response.