Antibodies (Abs) are a crucial component of the immune system and are often used as diagnostic and therapeutic agents. The need for high-affinity and high-specificity antibodies in research and medicine is driving the development of computational tools for accelerating antibody design and discovery. We report a diverse set of antibody binding data with accompanying structures that can be used to evaluate methods for modeling antibody interactions. Our Antibody-Bind (AB-Bind) database includes 1101 mutants with experimentally determined changes in binding free energies (ΔΔG) across 32 complexes. Using the AB-Bind data set, we evaluated the performance of protein scoring potentials in their ability to predict changes in binding free energies upon mutagenesis. Numerical correlations between computed and observed ΔΔG values were low (r = 0.16-0.45), but the potentials exhibited predictive power for classifying variants as improved vs weakened binders. Performance was evaluated using the area under the curve (AUC) for receiver operator characteristic (ROC) curves; the highest AUC values for 527 mutants with |ΔΔG| > 1.0 kcal/mol were 0.81, 0.87, and 0.88 using STATIUM, FoldX, and Discovery Studio scoring potentials, respectively. Some methods could also enrich for variants with improved binding affinity; FoldX and Discovery Studio were able to correctly rank 42% and 30%, respectively, of the 80 most improved binders (those with ΔΔG < -1.0 kcal/mol) in the top 5% of the database. This modest predictive performance has value but demonstrates the continuing need to develop and improve protein energy functions for affinity prediction.

The flavin-mediated enzymatic oxidation of a CN bond in amino acids can occur through hydride transfer, carbanion, or polar nucleophilic mechanisms. Previous results with D-arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) using multiple deuterium kinetic isotope effects (KIEs) and computational studies established preferred binding of the substrate protonated on the α-amino group, with cleavages of the NH and CH bonds occurring in asynchronous fashion, consistent with the three possible mechanisms. The hydroxyl groups of Y53 and Y249 are ≤4 Å from the imino and carboxylate groups of the reaction product iminoarginine, suggesting participation in binding and catalysis. In this study, we have investigated the reductive half-reactions of the Y53F and Y249F variants of PaDADH using substrate and solvent deuterium KIEs, solvent viscosity and pH effects, and quantum mechanical/molecular mechanical computational approaches to gain insights into the catalytic roles of the tyrosines and evaluate whether their mutations affect the transition state for substrate oxidation. Both Y53F and Y249F enzymes oxidized D-arginine with steady-state kinetic parameters similar to those of the wild-type enzyme. Rate constants for flavin reduction (k(red)) with D-leucine, a slow substrate amenable to rapid kinetics, were 3-fold smaller than the wild-type value with similar pKa values for an unprotonated group of ∼10.0. Similar pKa values were observed for (app)Kd in the variant and wild-type enzymes. However, cleavage of the substrate NH and CH bonds in the enzyme variants occurred in synchronous fashion, as suggested by multiple deuterium KIEs on k(red). These data can be reconciled with a hydride transfer mechanism, but not with carbanion and polar nucleophilic mechanisms.

Pseudomonas aeruginosa D‐arginine dehydrogenase (PaDADH) is a flavin dependent enzyme that can be classified as a member of D‐amino acid oxidase, sarcosine oxidase, and related enzyme families. It catalyzes the oxidation of D‐amino acids except D‐aspartate, D‐glutamate and glycine, to imino acids. Mechanistic details of C‐N bond oxidations of D‐amino acids by PaDADH are still elusive. Multiple mechanistic possibilities such as carbanion, hydride transfer and polar nucleophilic mechanisms could be considered. Previous work on PaDADH showed that NH and CH bond cleavages during amine oxidation are asynchronous with no further details available on the chemical mechanism. Here, we mutated two conserved active site residues Y53 and Y249 to phenylalanine in an attempt to reveal the mechanistic details of amine oxidation via pH profiles and deuterium isotope effects on rapid reaction kinetics. To further minimize mechanistic possibilities, QM/MM methodology was used to probe the relative local electrophilicity of key FAD sites within PaDADH active site. Our results are consistent with D‐amino acid oxidation occurring through the direct transfer of a hydride ion from the substrate to FAD without stabilization of reaction intermediates or FAD‐derived adducts, in which the active site tyrosines do not have a direct role in catalysis in PaDADH.Grant Funding Source: NSF MCB‐1121695 (G.G.) and a GSU Molecular Basis for Disease Fellowship (S.G.).

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nerve system that hydrolyzes acetylcholine (ACh) and terminates synaptic signals by reducing the effective concentration of ACh in the synaptic clefts. Organophosphate compounds irreversibly inhibit AChEs, leading to irreparable damage to nerve cells. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the covalent inhibition mechanism between AChE and the nerve toxin soman and determined its free energy profile for the first time. Our results indicate that phosphonylation of the catalytic serine by soman employs an addition-elimination mechanism, which is highly associative and stepwise: in the initial addition step, which is also rate-limiting, His440 acts as a general base to facilitate the nucleophilic attack of Ser200 on the soman's phosphorus atom to form a trigonal bipyrimidal pentacovalent intermediate; in the subsequent elimination step, Try121 of the catalytic gorge stabilizes the leaving fluorine atom prior to its dissociation from the active site. Together with our previous characterization of the aging mechanism of soman inhibited AChE, our simulations have revealed detailed molecular mechanistic insights into the damaging function of the nerve agent soman.

Computational enzyme design is an emerging field that has yielded promising success stories, but where numerous challenges remain. Accurate methods to rapidly evaluate possible enzyme design variants could provide significant value when combined with experimental efforts by reducing the number of variants needed to be synthesized and speeding the time to reach the desired endpoint of the design. To that end, extending our computational methods to model the fundamental physical-chemical principles that regulate activity in a protocol that is automated and accessible to a broad population of enzyme design researchers is essential. Here, we apply a physics-based implicit solvent MM-GBSA scoring approach to enzyme design and benchmark the computational predictions against experimentally determined activities. Specifically, we evaluate the ability of MM-GBSA to predict changes in affinity for a steroid binder protein, catalytic turnover for a Kemp eliminase, and catalytic activity for α-Gliadin peptidase variants. Using the enzyme design framework developed here, we accurately rank the most experimentally active enzyme variants, suggesting that this approach could provide enrichment of active variants in real-world enzyme design applications.

Enzyme design is an important area of ongoing research with a broad range of applications in protein therapeutics, biocatalysis, bioengineering, and other biomedical areas; however, significant challenges exist in the design of enzymes to catalyze specific reactions of interest. Here, we develop a computational protocol using an approach that combines molecular dynamics, docking, and MM-GBSA scoring to predict the catalytic activity of enzyme variants. Our primary focuses are to understand the molecular basis of substrate recognition and binding in an S-stereoselective ω-aminotransferase (ω-AT), which naturally catalyzes the transamination of pyruvate into alanine, and to predict mutations that enhance the catalytic efficiency of the enzyme. The conversion of (R)-ethyl 5-methyl-3-oxooctanoate to (3S,5R)-ethyl 3-amino-5-methyloctanoate in the context of several ω-AT mutants was evaluated using the computational protocol developed in this work. We correctly identify the mutations that yield the greatest improvements in enzyme activity (20-60-fold improvement over wild type) and confirm that the computationally predicted structure of a highly active mutant reproduces key structural aspects of the variant, including side chain conformational changes, as determined by X-ray crystallography. Overall, the protocol developed here yields encouraging results and suggests that computational approaches can aid in the redesign of enzymes with improved catalytic efficiency.

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nervous system that hydrolyzes neurotransmitter acetylcholine (ACh) and terminates synaptic signals. The catalytic serine of AChE can be phosphonylated by soman, one of the most potent nerve agents, and subsequently undergo an aging reaction. This phosphonylation and aging process leads to irreversible AChE inhibition, results in accumulation of excess ACh at the synaptic clefts, and causes neuromuscular paralysis. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the aging mechanism of soman phosphonylated AChE and determined its free energy profile. This aging reaction starts with the scission of the O2-Cα bond, which is followed by methyl migration, and results in a tertiary carbenium intermediate. At the transition state, the scissile O2-Cα bond is already cleaved with an average O-C distance of 3.2 ± 0.3 Å and the migrating methyl group is shared between Cα and Cβ carbons with C-C distances of 1.9 ± 0.1 and 1.8 ± 0.1 Å, respectively. The negatively charged phosphonate group is stabilized by a salt bridge with the imidazole ring of the catalytic histidine. A major product of aging, 2,3-dimethyl-2-butanol can be formed swiftly by the reaction of a water molecule. Our characterized mechanism and simulation results provide new detailed insights into this important biochemical process.