<p> We show that a combination of different theoretical methods is a viable approach to calculate and explain the relative binding affinities of inhibitors of the human butyrylcholinesterase enzyme. We probe structural properties of the enzyme-inhibitor complex in the presence of dialkyl phenyl phosphates and derivatives that include changes to the aromatic group and alkane-to-cholinyl substitutions that help these inhibitors mimic physiological substrates. Monte Carlo docking allowed for the identification of three regions within the active site of the enzyme where substituents of the phosphate group could be structurally stabilized. Computational clustering was used to identify distinct binding modes and their relative stabilities. Molecular dynamics suggest an essential asparagine residue not previously characterized as strongly influencing inhibitor strength which may serve as a crucial component in catalytic and inhibitory activity. This study provides a framework for suggesting future inhibitors that we expect will be effective at sub-micromolar concentrations. </p><p>