Positional ordering of reacting groups contributes significantly to the efficiency of proton transfer at an antibody active site

Multiscale computation delivers organophosphorus reactivity and stereoselectivity to immunoglobulin scavengers

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.

Enhancing a de novo enzyme activity by computationally-focused ultra-low-throughput screening

Directed evolution has revolutionized protein engineering. Still, enzyme optimization by random library screening remains sluggish, in large part due to futile probing of mutations that are catalytically neutral and/or impair stability and folding. FuncLib is a novel approach which uses phylogenetic analysis and Rosetta design to rank enzyme variants with multiple mutations, on the basis of predicted stability. Here, we use it to target the active site region of a minimalist-designed, de novo Kemp eliminase. The similarity between the Michaelis complex and transition state for the enzymatic reaction makes this system particularly challenging to optimize. Yet, experimental screening of a small number of active-site variants at the top of the predicted stability ranking leads to catalytic efficiencies and turnover numbers (∼2 × 104 M-1 s-1 and ∼102 s-1) for this anthropogenic reaction that compare favorably to those of modern natural enzymes. This result illustrates the promise of FuncLib as a powerful tool with which to speed up directed evolution, even on scaffolds that were not originally evolved for those functions, by guiding screening to regions of the sequence space that encode stable and catalytically diverse enzymes. Empirical valence bond calculations reproduce the experimental activation energies for the optimized eliminases to within ∼2 kcal mol-1 and indicate that the enhanced activity is linked to better geometric preorganization of the active site. This raises the possibility of further enhancing the stability-guidance of FuncLib by computational predictions of catalytic activity, as a generalized approach for computational enzyme design.

Electric Fields and Enzyme Catalysis

What happens inside an enzyme's active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists' attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme's active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

Kemp Eliminase Activity of Ketosteroid Isomerase

Kemp eliminases represent the most successful class of computationally designed enzymes, with rate accelerations of up to 10⁹-fold relative to the rate of the same reaction in aqueous solution. Nevertheless, several other systems such as micelles, catalytic antibodies, and cavitands are known to accelerate the Kemp elimination by several orders of magnitude. We found that the naturally occurring enzyme ketosteroid isomerase (KSI) also catalyzes the Kemp elimination. Surprisingly, mutations of D38, the residue that acts as a general base for its natural substrate, produced variants that catalyze the Kemp elimination up to 7000-fold better than wild-type KSI does, and some of these variants accelerate the Kemp elimination more than the computationally designed Kemp eliminases. Analysis of the D38N general base KSI variant suggests that a different active site carboxylate residue, D99, performs the proton abstraction. Docking simulations and analysis of inhibition by active site binders suggest that the Kemp elimination takes place in the active site of KSI and that KSI uses the same catalytic strategies of the computationally designed enzymes. In agreement with prior observations, our results strengthen the conclusion that significant rate accelerations of the Kemp elimination can be achieved with very few, nonspecific interactions with the substrate if a suitable catalytic base is present in a hydrophobic environment. Computational design can fulfill these requirements, and the design of more complex and precise environments represents the next level of challenges for protein design.

Computational strategies for the design of new enzymatic functions

In this contribution, recent developments in the design of biocatalysts are reviewed with particular emphasis in the de novo strategy. Studies based on three different reactions, Kemp elimination, Diels-Alder and Retro-Aldolase, are used to illustrate different success achieved during the last years. Finally, a section is devoted to the particular case of designed metalloenzymes. As a general conclusion, the interplay between new and more sophisticated engineering protocols and computational methods, based on molecular dynamics simulations with Quantum Mechanics/Molecular Mechanics potentials and fully flexible models, seems to constitute the bed rock for present and future successful design strategies.

Quantum and Molecular Mechanical (QM/MM) Monte Carlo Techniques for Modeling Condensed-Phase Reactions

A recent review (Acc. Chem. Res. 2010, 43:142-151) examined our use and development of a combined quantum and molecular mechanical (QM/MM) technique for modelling organic and enzymatic reactions. Advances included the PDDG/PM3 semiempirical QM (SQM) method, computation of multi-dimensional potentials of mean force (PMF), incorporation of on-the-fly QM in Monte Carlo simulations, and a polynomial quadrature method for rapidly treating proton-transfer reactions. The current article serves as a follow up on our progress. Highlights include new reactions, alternative SQM methods, a polarizable OPLS force field, and novel solvent environments, e.g., "on water" and room temperature ionic liquids. The methodology is strikingly accurate across a wide range of condensed-phase and antibody-catalyzed reactions including substitution, decarboxylation, elimination, isomerization, and pericyclic classes. Comparisons are made to systems treated with continuum-based solvents and ab initio or density functional theory (DFT) methods. Overall, the QM/MM methodology provides detailed characterization of reaction paths, proper configurational sampling, several advantages over implicit solvent models, and a reasonable computational cost.

Use of anion-aromatic interactions to position the general base in the ketosteroid isomerase active site

Although the cation-pi pair, formed between a side chain or substrate cation and the negative electrostatic potential of a pi system on the face of an aromatic ring, has been widely discussed and has been shown to be important in protein structure and protein-ligand interactions, there has been little discussion of the potential structural and functional importance in proteins of the related anion-aromatic pair (i.e., interaction of a negatively charged group with the positive electrostatic potential on the ring edge of an aromatic group). We posited, based on prior structural information, that anion-aromatic interactions between the anionic Asp general base and Phe54 and Phe116 might be used instead of a hydrogen-bond network to position the general base in the active site of ketosteroid isomerase from Comamonas testosteroni as there are no neighboring hydrogen-bonding groups. We have tested the role of the Phe residues using site-directed mutagenesis, double-mutant cycles, and high-resolution X-ray crystallography. These results indicate a catalytic role of these Phe residues. Extensive analysis of the Protein Data Bank provides strong support for a catalytic role of these and other Phe residues in providing anion-aromatic interactions that position anionic general bases within enzyme active sites. Our results further reveal a potential selective advantage of Phe in certain situations, relative to more traditional hydrogen-bonding groups, because it can simultaneously aid in the binding of hydrophobic substrates and positioning of a neighboring general base.

SABER: a computational method for identifying active sites for new reactions

A software suite, SABER (Selection of Active/Binding sites for Enzyme Redesign), has been developed for the analysis of atomic geometries in protein structures, using a geometric hashing algorithm (Barker and Thornton, Bioinformatics 2003;19:1644-1649). SABER is used to explore the Protein Data Bank (PDB) to locate proteins with a specific 3D arrangement of catalytic groups to identify active sites that might be redesigned to catalyze new reactions. As a proof-of-principle test, SABER was used to identify enzymes that have the same catalytic group arrangement present in o-succinyl benzoate synthase (OSBS). Among the highest-scoring scaffolds identified by the SABER search for enzymes with the same catalytic group arrangement as OSBS were L-Ala D/L-Glu epimerase (AEE) and muconate lactonizing enzyme II (MLE), both of which have been redesigned to become effective OSBS catalysts, demonstrated by experiments. Next, we used SABER to search for naturally existing active sites in the PDB with catalytic groups similar to those present in the designed Kemp elimination enzyme KE07. From over 2000 geometric matches to the KE07 active site, SABER identified 23 matches that corresponded to residues from known active sites. The best of these matches, with a 0.28 Å catalytic atom RMSD to KE07, was then redesigned to be compatible with the Kemp elimination using RosettaDesign. We also used SABER to search for potential Kemp eliminases using a theozyme predicted to provide a greater rate acceleration than the active site of KE07, and used Rosetta to create a design based on the proteins identified.

Evaluation of the energetics of the concerted acid-base mechanism in enzymatic catalysis: the case of ketosteroid isomerase

Structures of enzymes invariably reveal the proximity of acidic and basic residues to reactive sites on the substrate, so it is natural and common to suggest that enzymes employ concerted mechanisms to catalyze their difficult reactions. Ketosteroid isomerase (KSI) has served as a paradigm of enzymatic proton transfer chemistry, and its catalytic effect has previously been attributed to concerted proton transfer. We employ a specific inhibitor that contains an IR probe that reports directly and quantitatively on the ionization state of the ligand when bound in the active site of KSI. Measurement of the fractional ionization provides a missing link in a thermodynamic cycle that can discriminate the free energy advantage of a concerted versus nonconcerted mechanism. It is found that the maximum thermodynamic advantage that KSI could capture from a concerted mechanism (ΔΔG° = 0.5 kcal mol(-1)) is quite small.

Oxidoreductase activity of multifunctional monoclonal antibody B13-DE1

The monoclonal antibody B13-DE1 binds fluorescein, several fluorescein derivatives, and three peptide mimotopes. Our results revealed that this antibody also catalyzed the redox reaction of resazurin to resorufin, which are both structurally related to fluorescein. By using sodium sulfite as a reducing agent, the antibody B13-DE1 lowered the activation energy of this reaction. The Michaelis-Menten constant and turnover number of the catalyzed reaction were determined as 4.2 µmol/l and 0.0056 s(-1) , respectively. Because the results showed that fluorescein inhibited the catalytic activity of the antibody, we assume that the antigen-binding site and the catalytic active site are identical.

Role of water in the multifaceted catalytic antibody 4B2 for allylic isomerization and Kemp elimination reactions

Specificity toward a single reaction is a well-known characteristic of catalytic antibodies. However, contrary to convention, catalytic antibody 4B2 possesses the ability to efficiently catalyze two unrelated reactions: a Kemp elimination and an allylic isomerization of a beta,gamma-unsaturated ketone. To elucidate how this multifaceted antibody operates, mixed quantum and molecular mechanics calculations coupled to Monte Carlo simulations were carried out. The antibody was determined to derive its adaptability for the mechanistically different reactions through the rearrangement of water molecules in the active site into advantageous geometric orientations for enhanced electrostatic stabilization. In the case of the Kemp elimination, a general base, Glu L34, carried out the proton abstraction from the isoxazole ring of 5-nitro-benzisoxazole while water molecules delivered specific stabilization at the transition state. The role of water was found to be more pronounced in the allylic isomerization because the solvent actively participated in the stepwise mechanism. A rate-limiting abstraction of the alpha-proton from the beta,gamma-unsaturated ketone via Glu L34 led to the formation of a neutral dienol intermediate, which was rapidly reprotonated at the gamma-position via a solvent hydronium ion. Preferential channeling of H(3)O(+) in the active site ensured a stereoselective proton exchange from the alpha- to the gamma-position, in good agreement with deuterium exchange NMR and HPLC experiments. Ideas for improved water-mediated catalytic antibody designs are presented. In a technical advancement, improvements to a recent polynomial fitting and integration technique utilizing free energy perturbation theory delivered greater accuracy and speed gains.

Advances in quantum and molecular mechanical (QM/MM) simulations for organic and enzymatic reactions

Application of combined quantum and molecular mechanical (QM/MM) methods focuses on predicting activation barriers and the structures of stationary points for organic and enzymatic reactions. Characterization of the factors that stabilize transition structures in solution and in enzyme active sites provides a basis for design and optimization of catalysts. Continued technological advances allowed for expansion from prototypical cases to mechanistic studies featuring detailed enzyme and condensed-phase environments with full integration of the QM calculations and configurational sampling. This required improved algorithms featuring fast QM methods, advances in computing changes in free energies including free-energy perturbation (FEP) calculations, and enhanced configurational sampling. In particular, the present Account highlights development of the PDDG/PM3 semi-empirical QM method, computation of multi-dimensional potentials of mean force (PMF), incorporation of on-the-fly QM in Monte Carlo (MC) simulations, and a polynomial quadrature method for efficient modeling of proton-transfer reactions. The utility of this QM/MM/MC/FEP methodology is illustrated for a variety of organic reactions including substitution, decarboxylation, elimination, and pericyclic reactions. A comparison to experimental kinetic results on medium effects has verified the accuracy of the QM/MM approach in the full range of solvents from hydrocarbons to water to ionic liquids. Corresponding results from ab initio and density functional theory (DFT) methods with continuum-based treatments of solvation reveal deficiencies, particularly for protic solvents. Also summarized in this Account are three specific QM/MM applications to biomolecular systems: (1) a recent study that clarified the mechanism for the reaction of 2-pyrone derivatives catalyzed by macrophomate synthase as a tandem Michael-aldol sequence rather than a Diels-Alder reaction, (2) elucidation of the mechanism of action of fatty acid amide hydrolase (FAAH), an unusual Ser-Ser-Lys proteolytic enzyme, and (3) the construction of enzymes for Kemp elimination of 5-nitrobenzisoxazole that highlights the utility of QM/MM in the design of artificial enzymes.

An aspartate and a water molecule mediate efficient acid-base catalysis in a tailored antibody pocket

Design of catalysts featuring multiple functional groups is a desirable, yet formidable goal. Antibody 13G5, which accelerates the cleavage of unactivated benzisoxazoles, is one of few artificial enzymes that harness an acid and a base to achieve efficient proton transfer. X-ray structures of the Fab-hapten complexes of wild-type 13G5 and active-site variants now afford detailed insights into its mechanism. The parent antibody preorganizes Asp(H35) and Glu(L34) to abstract a proton from substrate and to orient a water molecule for leaving group stabilization, respectively. Remodeling the environment of the hydrogen bond donor with a compensatory network of ordered waters, as seen in the Glu(L34) to alanine mutant, leads to an impressive 10(9)-fold rate acceleration over the nonenzymatic reaction with acetate, illustrating the utility of buried water molecules in bifunctional catalysis. Generalization of these design principles may aid in creation of catalysts for other important chemical transformations.

Origin of the activity drop with the E50D variant of catalytic antibody 34E4 for Kemp elimination

In enzymes, multiple structural effects cooperatively lead to the high catalytic activity, while individually these effects can be small. The design of artificial enzymes requires the understanding and ability to manipulate such subtle effects. The 34E4 catalytic antibody, catalyzing Kemp elimination of 5-nitrobenzisoxazole, and its Glu50Asp (E50D) variant are the subject of the present investigation. This removal of only a methylene group yields an approximately 30-fold reduction in the rate for the catalyzed Kemp elimination. Here, the aim is to understand this difference in the catalytic performance. The mechanism of Kemp elimination catalyzed by 34E4 and the E50D mutant is elucidated using QM/MM Monte Carlo simulations and free energy perturbation theory. In both proteins, the reaction is shown to follow a single-step, concerted mechanism. In the mutant, the activation barrier rises by 2.4 kcal/mol, which corresponds to a 62-fold rate deceleration, which is in good agreement with the experimental data. The positions and functionality of the residues in the active site are monitored throughout the reaction. It is concluded that the looser contact with the base, shorter base-Asn58 contact, less favorable pi-stacking with Trp91 in the transition state of the reaction, and different solvation pattern all contribute to the reduction of the reaction rate in the E50D variant of 34E4.

Testing geometrical discrimination within an enzyme active site: constrained hydrogen bonding in the ketosteroid isomerase oxyanion hole

Enzymes are classically proposed to accelerate reactions by binding substrates within active-site environments that are structurally preorganized to optimize binding interactions with reaction transition states rather than ground states. This is a remarkably formidable task considering the limited 0.1-1 A scale of most substrate rearrangements. The flexibility of active-site functional groups along the coordinate of substrate rearrangement, the distance scale on which enzymes can distinguish structural rearrangement, and the energetic significance of discrimination on that scale remain open questions that are fundamental to a basic physical understanding of enzyme active sites and catalysis. We bring together 1.2-1.5 A resolution X-ray crystallography, (1)H and (19)F NMR spectroscopy, quantum mechanical calculations, and transition-state analogue binding measurements to test the distance scale on which noncovalent forces can constrain the structural relaxation or translation of side chains and ligands along a specific coordinate and the energetic consequences of such geometric constraints within the active site of bacterial ketosteroid isomerase (KSI). Our results strongly suggest that packing and binding interactions within the KSI active site can constrain local side-chain reorientation and prevent hydrogen bond shortening by 0.1 A or less. Further, this constraint has substantial energetic effects on ligand binding and stabilization of negative charge within the oxyanion hole. These results provide evidence that subtle geometric effects, indistinguishable in most X-ray crystallographic structures, can have significant energetic consequences and highlight the importance of using synergistic experimental approaches to dissect enzyme function.

Conformational isomerism can limit antibody catalysis

Ligand binding to enzymes and antibodies is often accompanied by protein conformational changes. Although such structural adjustments may be conducive to enzyme catalysis, much less is known about their effect on reactions promoted by engineered catalytic antibodies. Crystallographic and pre-steady state kinetic analyses of antibody 34E4, which efficiently promotes the conversion of benzisoxazoles to salicylonitriles, show that the resting catalyst adopts two interconverting active-site conformations, only one of which is competent to bind substrate. In the predominant isomer, the indole side chain of Trp(L91) occupies the binding site and blocks ligand access. Slow conformational isomerization of this residue, on the same time scale as catalytic turnover, creates a deep and narrow binding site that can accommodate substrate and promote proton transfer using Glu(H50) as a carboxylate base. Although 34E4 is among the best catalysts for the deprotonation of benzisoxazoles, its efficiency appears to be significantly limited by this conformational plasticity of its active site. Future efforts to improve this antibody might profitably focus on stabilizing the active conformation of the catalyst. Analogous strategies may also be relevant to other engineered proteins that are limited by an unfavorable conformational pre-equilibrium.

Novel enzymes through design and evolution

The generation of enzymes with new catalytic activities remains a major challenge. So far, several different strategies have been developed to tackle this problem, including site-directed mutagenesis, random mutagenesis (directed evolution), antibody catalysis, computational redesign, and de novo methods. Using these techniques, a broad array of novel enzymes has been created (aldolases, decarboxylases, dehydratases, isomerases, oxidases, reductases, and others), although their low efficiencies (10 to 100 M(-1) s(-l)) compared to those of the best natural enzymes (10(6) to 10(8) M(-1) s(-1)) remains a significant concern. Whereas rational design might be the most promising and versatile approach to generating new activities, directed evolution seems to be the best way to optimize the catalytic properties of novel enzymes. Indeed, impressive successes in enzyme engineering have resulted from a combination of rational and random design.

Survey of the year 2005 commercial optical biosensor literature

We identified 1113 articles (103 reviews, 1010 primary research articles) published in 2005 that describe experiments performed using commercially available optical biosensors. While this number of publications is impressive, we find that the quality of the biosensor work in these articles is often pretty poor. It is a little disappointing that there appears to be only a small set of researchers who know how to properly perform, analyze, and present biosensor data. To help focus the field, we spotlight work published by 10 research groups that exemplify the quality of data one should expect to see from a biosensor experiment. Also, in an effort to raise awareness of the common problems in the biosensor field, we provide side-by-side examples of good and bad data sets from the 2005 literature.

Tunable pK of Amino Acid Residues at the Air−Water Interface Gives an L-zyme (Langmuir Enzyme)

Various amino acid-carrying amphiphiles were synthesized, and the pK values of the attached amino acid residues were investigated at the air-water interface and in aqueous vesicles using pi-A isotherm measurements, (1)H NMR titration, and IR spectroscopy in reflection-adsorption mode. The epsilon-amino group of the Lys residue embedded at the air-water interface displays a significant pK shift (4 or 5 unit) compared with that observed in bulk water, while the pK shift in aqueous vesicles was not prominent (ca. 1 unit). Moreover, pK values of the amino acids at the air-water interface can be tuned simply by control of the subphase ionic strength as well as by molecular design of the amphiphiles. A simple equation based on the dominant contribution by the electrostatic energy to the pK shift reproduces well the surface pressure difference between protonated and unprotonated species, suggesting a reduction in the apparent dielectric constant at the air-water interface. Hydrolysis of a p-nitrophenyl ester derivative was used as a model reaction to demonstrate the use of the Lys-functionalized monolayer. Efficient hydrolysis was observed, even at neutral pH, after tuning of pK for the Lys residue in the monolayer, which is a similar case to that occurring in biological catalysis.

Structural origins of efficient proton abstraction from carbon by a catalytic antibody

Antibody 34E4 catalyzes the conversion of benzisoxazoles to salicylonitriles with high rates and multiple turnovers. The crystal structure of its complex with the benzimidazolium hapten at 2.5-angstroms resolution shows that a combination of hydrogen bonding, pi stacking, and van der Waals interactions is exploited to position both the base, Glu(H50), and the substrate for efficient proton transfer. Suboptimal placement of the catalytic carboxylate, as observed in the 2.8-angstroms structure of the Glu(H50)Asp variant, results in substantially reduced catalytic efficiency. In addition to imposing high positional order on the transition state, the antibody pocket provides a highly structured microenvironment for the reaction in which the carboxylate base is activated through partial desolvation, and the highly polarizable transition state is stabilized by dispersion interactions with the aromatic residue Trp(L91) and solvation of the leaving group oxygen by external water. The enzyme-like efficiency of general base catalysis in this system directly reflects the original hapten design, in which a charged guanidinium moiety was strategically used to elicit an accurately positioned functional group in an appropriate reaction environment and suggests that even larger catalytic effects may be achievable by extending this approach to the induction of acid-base pairs capable of bifunctional catalysis.