Molecular dynamics was used to simulate the transition state for the first chemical reaction step (TS1) of cocaine hydrolysis catalyzed by human butyrylcholinesterase (BChE) and its mutants. The simulated results demonstrate that the overall hydrogen bonding between the carbonyl oxygen of (؊)-cocaine benzoyl ester and the oxyanion hole of BChE in the TS1 structure for (؊)-cocaine hydrolysis catalyzed by A199S͞S287G͞A328W͞Y332G BChE should be significantly stronger than that in the TS1 structure for (؊)-cocaine hydrolysis catalyzed by the WT BChE and other simulated BChE mutants. Thus, the transition-state simulations predict that A199S͞ S287G͞A328W͞Y332G mutant of BChE should have a significantly lower energy barrier for the reaction process and, therefore, a significantly higher catalytic efficiency for (؊)-cocaine hydrolysis. The theoretical prediction has been confirmed by wet experimental tests showing an Ϸ(456 ؎ 41)-fold improved catalytic efficiency of A199S͞S287G͞A328W͞Y332G BChE against (؊)-cocaine. This is a unique study to design an enzyme mutant based on transitionstate simulation. The designed BChE mutant has the highest catalytic efficiency against cocaine of all of the reported BChE mutants, demonstrating that the unique design approach based on transition-state simulation is promising for rational enzyme redesign and drug discovery. molecular dynamics ͉ rational design ͉ transition-state stabilization ͉ cocaine ͉ enzyme-substrate binding C ocaine is recognized as the most reinforcing of all drugs of abuse (1-3). The disastrous medical and social consequences of cocaine addiction have made the development of an effective pharmacological treatment a high priority (4-6). However, cocaine mediates its reinforcing and toxic effects by blocking neurotransmitter reuptake, and the classical pharmacodynamic approach has failed to yield small-molecule receptor antagonists because of the difficulties inherent in blocking a blocker (1-5). An alternative to receptor-based approaches is to interfere with the delivery of cocaine to its receptors or accelerate its metabolism in the body (5,(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17). An ideal molecule for this purpose should be a potent enzyme catalyzing the hydrolysis of cocaine into biologically inactive metabolites. The dominant pathway for cocaine metabolism in primates is butyrylcholinesterase (BChE)-catalyzed hydrolysis at the benzoyl ester group (Fig. 3, which is published as supporting information on the PNAS web site), and the metabolites are all biologically inactive (5, 18). Clearly, BChE-catalyzed hydrolysis of cocaine at the benzoyl ester is the metabolic pathway most suitable for amplification. However, the catalytic activity of this plasma enzyme is Ϸ1,000-fold lower against the naturally occurring (Ϫ)-cocaine than that against the biologically inactive (ϩ)-cocaine enantiomer (19)(20)(21)(22). (ϩ)-cocaine can be cleared from plasma in seconds, before partitioning into the CNS, whereas (Ϫ)-cocaine has a plasma half-life of Ϸ45-90 min, long enough for ...