Two competing passageways for the exit of acetic acid and acetate ion in
Torpedo californica (Tc)acetylcholinesterase (ACHE) were studied by examining free energies of passage through two potentialtra
jectories using the umbrella sampling technique as implemented in CHARMM. The coordinates for migrationwere defined as the distance from Ser200 O
![](/images/gifchars/gamma.gif)
, one through the 20-Å long active-site gorge ending with Trp279and a 14-Å long route ending at Arg244. The free energies were calculated in successive windows 0.5 Åwide for 40-90 ps. The potential of mean force (PMF) was calculated along the coordinate for migration.The PMF for the migration of acetic acid decreases by ~8 kcal/mol after 8-Å travel through the main gorge.The PMF profile for acetate ion migration falls to a 6 kcal/mol lower value than for acetic acid migration inthe main gorge. The free energy barrier for the migration of acetate ion is 1.5 kcal/mol due to a constrictionformed by Tyr121, Phe290, Phe330, and Phe331 in the main gorge. The interaction between acetic acid/acetate ion and the OH group of Tyr121 appears to guide product release through the main gorge. Acetate ionremains hydrogen-bonded to Tyr121 until it approaches Trp279 when it is expelled into bulk water. Aceticacid encounters a 6 kcal/mol barrier through the alternate pathway, while the PMF for acetate ion drops ~27kcal/mol when it approaches Arg244. This is the lowest energy path. Full molecular dynamics simulations,free of restraint for 170 ps, result in the migration of acetate ion through the short channel but not through themain gorge. The results indicate that if the nascent acetic acid should ionize within 3.5 Å from Ser200 O
![](/images/gifchars/gamma.gif)
,it would be more likely to exit via the alternate channel than through the main gorge.