Starch is a branched polymer of glucose with two components, both of which have (1 → 4)-α linear links and (1 → 6)-α branch points: amylopectin, of high molecular weight with many short branches, and amylose, of lower molecular weight and only a few long-chain branches. Granule-bound starch synthase I (GBSSI) is one of the main enzymes controlling amylose synthesis and chain-length distribution. As production of different GBSSI mutants is time-consuming and laborious, molecular dynamics (MD) simulations are used here to predict the binding of different GBSSI mutants to a representative amylose fragment. The simulations were atomistic, with explicit solvent and docking, a method successfully used to understand the binding of wild-type GBSSI to amylose fragments. The binding of GBSSI to G5 (a pentasaccharide amylose fragment) is combined with free-energy calculations employing a thermodynamic integration method to predict the effects of mutations on enzyme activity. Ten GBSSI mutants with different enzyme activities were analyzed to find the structural and energy changes among different single amino-acid mutants and their possible relationship to starch characteristics. Comparing the structural changes and the relative binding free energy of G5 to the wild type GBSSI and GBSSI mutants, it was found that mutants with negative binding energy (lower than -2.0 kcal/mol) are more likely to have higher enzyme activity and amylose content compared to the wild type. This theoretical paper used simulations and robust free energy calculations to interpret in planta data with potential predictions as to what mutants might be generated to give desired properties. This study can be used to help develop grains with improved functional properties.