A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the binding specificity of acetylcholinesterase (AChE) with its natural substrate acetylcholine (ACh), a family of substrate analogues, and choline. Calculated docking energies are well correlated to experimental k(cat)/K(M) values, as well as to experimental binding affinities of a related series of TMTFA inhibitors. The "esteratic" and "anionic" subsites are found to act together to achieve substrate binding specificity. We find that the presence of ACh in the active site of AChE not only stabilizes the setup of the catalytic triad but also tightens both subsites to achieve better binding. The docking energy gained from this induced fit is 0.7 kcal/mol for ACh. For the binding of the substrate tailgroup to the anionic subsite, both the size and the positive charge of the tailgroup are important. The removal of the positive charge leads to a weaker binding of 1 kcal/mol loss in docking energy. Substituting each tail methyl group with hydrogen results in both an incremental loss in docking energy and also a decrease in the percentage of structures docked in the active site correctly set up for catalysis.