ADP-glucose pyrophosphorylase (AGPase) is a key allosteric enzyme in plant starch biosynthesis. Plant AGPase is a heterotetrameric enzyme that consists of large (LS) and small subunits (SS), which are encoded by two different genes. Computational and experimental studies have revealed that the heterotetrameric assembly of AGPase is thermodynamically weak. Modeling studies followed by the mutagenesis of the LS of the potato AGPase identified a heterotetramer-deficient mutant, LSR88A. To enhance heterotetrameric assembly, LSR88A cDNA was subjected to error-prone PCR, and second-site revertants were identified according to their ability to restore glycogen accumulation, as assessed with iodine staining. Selected mutations were introduced into the wild-type (WT) LS and co-expressed with the WT SS in Escherichia coli glgC(-). The biochemical characterization of revertants revealed that (LSSSWT)-S-I90V, (LSSSWT)-S-Y378C and (LSSSWT)-S-D410G mutants displayed enhanced heterotetrameric assembly with the WT SS. Among these mutants, (LSSSWT)-S-Y378C AGPase displayed increased heat stability compared with the WT enzyme. Kinetic characterization of the mutants indicated that the (LSSSWT)-S-I90V and (LSSSWT)-S-Y378C AGPases have comparable allosteric and kinetic properties. However, the (LSSSWT)-S-D410G mutant exhibited altered allosteric properties of being less responsive and more sensitive to 3-phosphoglyceric acid activation and inorganic phosphate inhibition. This study not only enhances our understanding of the interaction between the SS and the LS of AGPase but also enables protein engineering to obtain enhanced assembled heat-stable variants of AGPase, which can be used for the improvement of plant yields.