Abstract
Significance
Developed quantum mechanics/molecular mechanics (QM/MM) supercomputation is used to identify maturation of immunoglobulin (Ig) reactivity through optimization of amino acids in the catalytic core, focused on enhanced nucleophilic attack on phosphorus. A computationally optimized variant,
l
-Leu47Lys, delivers 340-fold reactivity enhancement for paraoxon relative to wild-type Ig-paraoxonase resulting from substrate-induced ionization of the Tyr37 nucleophile. Variant antibody stereoselectivity toward a P-chiral phenylphosphonate toxin results from starkly contrasting solvent-assisted or general base catalytic mechanisms relative to the parent Ig-paraoxonase. These mechanisms are defined by crystallographic and QM/MM analyses as resulting from changed substrate orientation and participating amino acid conformations within the catalytic site.
Quantum mechanics/molecular mechanics (QM/MM) maturation of an immunoglobulin (Ig) powered by supercomputation delivers novel functionality to this catalytic template and facilitates artificial evolution of biocatalysts. We here employ density functional theory-based (DFT-b) tight binding and funnel metadynamics to advance our earlier QM/MM maturation of A17 Ig-paraoxonase (WTIgP) as a reactibody for organophosphorus toxins. It enables regulation of biocatalytic activity for tyrosine nucleophilic attack on phosphorus. The single amino acid substitution
l
-Leu47Lys results in 340-fold enhanced reactivity for paraoxon. The computed ground-state complex shows substrate-induced ionization of the nucleophilic
l
-Tyr37, now H-bonded to
l
-Lys47, resulting from repositioning of
l
-Lys47. Multiple antibody structural homologs, selected by phenylphosphonate covalent capture, show contrasting enantioselectivities for a P-chiral phenylphosphonate toxin. That is defined by crystallographic analysis of phenylphosphonylated reaction products for antibodies A5 and WTIgP. DFT-b analysis using QM regions based on these structures identifies transition states for the favored and disfavored reactions with surprising results. This stereoselection analysis is extended by funnel metadynamics to a range of WTIgP variants whose predicted stereoselectivity is endorsed by experimental analysis. The algorithms used here offer prospects for tailored design of highly evolved, genetically encoded organophosphorus scavengers and for broader functionalities of members of the Ig superfamily, including cell surface-exposed receptors.
