The interaction of a series of indole-2-carboxamide compounds with human liver glycogen phosphorylase a (HLGPa) have been studied employing molecular docking and 3D-QSAR approaches. The Lamarckian Genetic Algorithm (LGA) of AutoDock 3.0 was employed to locate the binding orientations and conformations of the inhibitors interacting with HLGPa. The binding models were demonstrated in the aspects of inhibitor's conformation, subsite interaction, and hydrogen bonding. The very similar binding conformations of these inhibitors show that they interact with HLGPa in a very similar way. Good correlations between the calculated interaction free energies and experimental inhibitory activities suggest that the binding conformations of these inhibitors are reasonable. The structural and energetic differences in inhibitory potencies of indole-2-carboxamide compounds were reasonably explored. Using the binding conformations of indole-2-carboxamides, consistent and highly predictive 3D-QSAR models were developed by CoMFA and CoMSIA analyses. The q2 values are 0.697 and 0.622 for CoMFA and CoMSIA models, respectively. The predictive ability of these models was validated by four compounds that were not included in the training set. Mapping these models back to the topology of the active site of HLGPa leads to a better understanding of the vital indole-2-carboxamide-HLGPa interactions. Structure-based investigations and the final 3D-QSAR results provide clear guidelines and accurate activity predictions for novel inhibitor design.