Summary Substantial improvements in enzyme activity demand multiple mutations at spatially proximal positions in the active site. Such mutations, however, often exhibit unpredictable epistatic (non-additive) effects on activity. Here, we describe FuncLib - an automated method for designing multipoint mutations at enzyme active sites using phylogenetic analysis and Rosetta design calculations. We applied FuncLib to two unrelated enzymes, a phosphotriesterase and an acetyl-CoA synthetase. All designs were active and most showed activity profiles that significantly differed from the wild type and from one another. Several dozen designs with only 3-6 active-site mutations exhibited 10-4,000-fold higher efficiencies with a range of alternative substrates, including the hydrolysis of the toxic organophosphate nerve agents soman and cyclosarin and the synthesis of butyryl-CoA. FuncLib is implemented as a web-server (http://FuncLib.weizmann.ac.il); it circumvents iterative, high-throughput screens and opens the way to design highly efficient and diverse catalytic repertoires.

Organophosphate nerve agents are extremely lethal compounds. Rapid in vivo organophosphate clearance requires bioscavenging enzymes with catalytic efficiencies of >10(7) (M(-1) min(-1)). Although serum paraoxonase (PON1) is a leading candidate for such a treatment, it hydrolyzes the toxic S(p) isomers of G-agents with very slow rates. We improved PON1's catalytic efficiency by combining random and targeted mutagenesis with high-throughput screening using fluorogenic analogs in emulsion compartments. We thereby enhanced PON1's activity toward the coumarin analog of S(p)-cyclosarin by ∼10(5)-fold. We also developed a direct screen for protection of acetylcholinesterase from inactivation by nerve agents and used it to isolate variants that degrade the toxic isomer of the coumarin analog and cyclosarin itself with k(cat)/K(M) ∼ 10(7) M(-1) min(-1). We then demonstrated the in vivo prophylactic activity of an evolved variant. These evolved variants and the newly developed screens provide the basis for engineering PON1 for prophylaxis against other G-type agents.

The ability to redesign enzymes to catalyze non-cognate chemical transformations would have wide-ranging applications. We developed a computational method for repurposing the reactivity of active site functional groups of metalloenzymes to catalyze new reactions. Using this method, we engineered a zinc-containing murine adenosine deaminase to catalyze the hydrolysis of a model organophosphate with a catalytic efficiency kcat/Km ~104 M−1s−1 after directed evolution. In the high-resolution crystal structure of the enzyme, all but one of the designed residues adopt the designed conformation. The designed enzyme efficiently catalyzed the hydrolysis of the RP-isomer of a coumarinyl analog of the nerve agent cyclosarin, and showed striking substrate selectivity for coumarinyl leaving-groups. Computational redesign of native enzyme active sites complements directed evolution methods and offers a general approach for exploring their untapped catalytic potential for new reactivities.

VX and its Russian (RVX) and Chinese (CVX) analogues rapidly inactivate acetylcholinesterase and are the most toxic stockpile nerve agents. These organophosphates have a thiol leaving group with a choline-like moiety and are hydrolyzed very slowly by natural enzymes. We used an integrated computational and experimental approach to increase Brevundimonas diminuta phosphotriesterase's (PTE) detoxification rate of V-agents by 5000-fold. Computational models were built of the complex between PTE and V-agents. On the basis of these models, the active site was redesigned to be complementary in shape to VX and RVX and to include favorable electrostatic interactions with their choline-like leaving group. Small libraries based on designed sequences were constructed. The libraries were screened by a direct assay for V-agent detoxification, as our initial studies showed that colorimetric surrogates fail to report the detoxification rates of the actual agents. The experimental results were fed back to improve the computational models. Overall, five rounds of iterating between experiment and model refinement led to variants that hydrolyze the toxic SP isomers of all three V-agents with kcat/KM values of up to 5 × 10(6) M(-1) min(-1) and also efficiently detoxify G-agents. These new catalysts provide the basis for broad spectrum nerve agent detoxification.

A preferred strategy for preventing nerve agents intoxication is catalytic scavenging by enzymes that hydrolyze them before they reach their targets. Using directed evolution, we simultaneously enhanced the activity of a previously described serum paraoxonase 1 (PON1) variant for hydrolysis of the toxic S(P) isomers of the most threatening G-type nerve agents. The evolved variants show ≤340-fold increased rates and catalytic efficiencies of 0.2-5 × 10(7) M(-1) min(-1). Our selection for prevention of acetylcholinesterase inhibition also resulted in the complete reversion of PON1's stereospecificity, from an enantiomeric ratio (E) < 6.3 × 10(-4) in favor of the R(P) isomer of a cyclosarin analog in wild-type PON1, to E > 2,500 for the S(P) isomer in an evolved variant. Given their ability to hydrolyze G-agents, these evolved variants may serve as broad-range G-agent prophylactics.