Thermodynamic penalty arising from burial of a ligand polar group within a hydrophobic pocket of a protein receptor

Here, we examine the thermodynamic penalty arising from burial of a polar group in a hydrophobic pocket that forms part of the binding-site of the major urinary protein (MUP-I). X-ray crystal structures of the complexes of octanol, nonanol and 1,8 octan-diol indicate that these ligands bind with similar orientations in the binding pocket. Each complex is characterised by a bridging water molecule between the hydroxyl group of Tyr120 and the hydroxyl group of each ligand. The additional hydroxyl group of 1,8 octan-diol is thereby forced to reside in a hydrophobic pocket, and isothermal titration calorimetry experiments indicate that this is accompanied by a standard free energy penalty of +21 kJ/mol with respect to octanol and +18 kJ/mol with respect to nonanol. Consideration of the solvation thermodynamics of each ligand enables the "intrinsic" (solute-solute) interaction energy to be determined, which indicates a favourable enthalpic component and an entropic component that is small or zero. These data indicate that the thermodynamic penalty to binding derived from the unfavourable desolvation of 1,8 octan-diol is partially offset by a favourable intrinsic contribution. Quantum chemical calculations suggest that this latter contribution derives from favourable solute-solute dispersion interactions.     

Structural determinants of trypsin affinity and specificity for cationic inhibitors

The binding free energies of four inhibitors to bovine beta-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson-Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.     

Conformational selection through electrostatics: Free energy simulations of GTP and GDP binding to archaeal initiation factor 2

Archaeal Initiation Factor 2 is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and their dependence on the conformational state, molecular dynamics free energy simulations were performed. The ON state specificity was predicted to be weak, with a GTP/GDP binding free energy difference of -1 kcal/mol, favoring GTP. The OFF state specificity is larger, 4 kcal/mol, favoring GDP. The overall effects result from a competition among many interactions in several complexes. To interpret them, we use a simpler, dielectric continuum model. Several effects are robust with respect to the model details. Both nucleotides have a net negative charge, so that removing them from solvent into the binding pocket carries a desolvation penalty, which is large for the ON state, and strongly disfavors GTP binding compared to GDP. Short-range interactions between the additional GTP phosphate group and ionized sidechains in the binding pocket offset most, but not all of the desolvation penalty; more distant groups also contribute significantly, and the switch 1 loop only slightly. The desolvation penalty is lower for the more open, wetter OFF state, and the GTP/GDP difference much smaller. Short-range interactions in the binding pocket and with more distant groups again make a significant contribution. Overall, the simulations help explain how conformational selection is achieved with a single phosphate group.     

On the binding determinants of the glutamate agonist with the glutamate receptor ligand binding domain

Ionotropic glutamate receptors (GluRs) are ligand-gated membrane channel proteins found in the central neural system that mediate a fast excitatory response of neurons. In this paper, we report theoretical analysis of the ligand-protein interactions in the binding pocket of the S1S2 (ligand binding) domain of the GluR2 receptor in the closed conformation. By utilizing several theoretical methods ranging from continuum electrostatics to all-atom molecular dynamics simulations and quantum chemical calculations, we were able to characterize in detail glutamate agonist binding to the wild-type and E705D mutant proteins. A theoretical model of the protein-ligand interactions is validated via direct comparison of theoretical and Fourier transform infrared spectroscopy (FTIR) measured frequency shifts of the ligand's carboxylate group vibrations [Jayaraman et al. (2000) Biochemistry 39, 8693-8697; Cheng et al. (2002) Biochemistry 41, 1602-1608]. A detailed picture of the interactions in the binding site is inferred by analyzing contributions to vibrational frequencies produced by protein residues forming the ligand-binding pocket. The role of mobility and hydrogen-bonding network of water in the ligand-binding pocket and the contribution of protein residues exposed in the binding pocket to the binding and selectivity of the ligand are discussed. It is demonstrated that the molecular surface of the protein in the ligand-free state has mainly positive electrostatic potential attractive to the negatively charged ligand, and the potential produced by the protein in the ligand-binding pocket in the closed state is complementary to the distribution of the electrostatic potential produced by the ligand itself. Such charge complementarity ensures specificity to the unique charge distribution of the ligand.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

Selective alteration of substrate specificity by replacement of aspartic acid-189 with lysine in the binding pocket of trypsin

To test the role of Asp-189 which is located at the base of the substrate binding pocket in determining the specificity of trypsin toward basic substrates, this residue was replaced with a lysine residue by site-directed mutagenesis. Both rat trypsinogen and Lys-189 trypsinogen were expressed and secreted into the periplasmic space of Escherichia coli. The proteins were purified to homogeneity and activated by porcine enterokinase, and their catalytic activities were determined on natural and synthetic substrates. Lys-189 trypsin displayed no catalytic activity toward arginyl and lysyl substrates. Further, there was no compensatory change in specificity toward acidic substrates; no cleavage of aspartyl or glutamyl bonds was detected. Additional studies of substrate specificity involving gas-phase sequence analyses of digested natural substrates revealed an inherent but low chymotrypsin-like activity of trypsin. This activity was retained but modified by the Asp to Lys change at position 189. In addition to hydrolyzing phenylalanyl and tyrosyl peptide bonds, the mutant enzyme has the unique property of cleaving leucyl bonds. On the basis of computer graphic modeling studies of the Lys-189 side chain, it appears that the positively charged NH2 group is directed outside the substrate binding pocket. The resulting hydrophobic cavity may explain the altered substrate specificity of the mutant enzyme. The relatively low chymotrypsin-like activity of both recombinant enzymes may be due to distorted positioning of the scissile bond with respect to the catalytic triad rather than to the lack of sufficient interaction between the hydrophobic side chains and the substrate binding pocket of the enzyme.     

The association between a negatively charged ligand and the electronegative binding pocket of its receptor.

Many examples exist of charged amino acids that play a role in attracting or holding a charged ligand toward or inside an oppositely charged binding pocket of the protein. For example, the enzymes superoxide dismutase, triose-phosphate isomerase, and acetylcholinesterase can steer ligands toward their oppositely charged binding pockets or gorges. Interestingly, in our Brownian dynamics simulations of a phosphate-binding protein, we discovered that negatively charged phosphate (HPO(2-)(4)) could make its way into the negatively charged binding pocket. In fact, the phosphate-binding protein exhibits counterintuitive kinetics of association. That is, one would expect that the rate of association would increase on increases to the ionic strength since the interaction between the ligand, with a charge of -2, and the electronegative binding pocket would be repulsive and greater screening should reduce this repulsion and increase the rate of association. However, the opposite is seen-i.e., the rate of association decreases on increases in the ionic strength. We used Brownian dynamics techniques to compute the diffusion limited association rate constants between the negatively charged phosphate ligand and several open forms of PBP (wild-type and several mutants based on an x-ray structure of open-form PBP, mutant T141D). With the appropriate choices of reaction criteria and molecular parameters, the ligand was able to diffuse into the binding pocket. A number of residues influence binding of the ligand within the pocket via hydrogen bonds or salt bridges. Arg135 partially neutralizes the charges on the HPO(2-)(4) ligand in the binding pocket, allowing it to enter. It is also found that the positive electrostatic patches above and below the binding entrance of PBP contribute the major attractive forces that direct the ligand toward the surface of the protein near the binding site.     

A model for the C5a receptor and for its interaction with the ligand [corrected]

A model of the C5a receptor was built based on the assumption that the seven membrane-spanning helices of known inward/outward direction are in an arrangement roughly similar to that in bacteriorhodopsin. Guidelines for the positioning of the helices were cysteine pairing, 'ridges into grooves' interdigitation of side chains and aromatic cluster formation. The chain segments protruding from the membrane are too short for folding into an independent ectodomain. The only longer segment (179-202) is tied down in its centre onto the membrane by a disulphide bridge and, thereby, made into two short loops as well. Ideas of the interaction of the C5a receptor with its ligand were derived mainly from the search for accommodation of the functionally essential arginine residues 40 and 74 of C5a. Asp82 is the only charged residue in a pocket approximately 20 A below the receptor surface and is conserved in the rhodopsin superfamily. It commends itself for binding Arg74 which is the tip of the flexible C-terminal chain of C5a, and rules out Arg40 in the structurally well-defined part of the molecule. The latter may bind to Glu180 at the bottom of a more shallow pocket which happens to resemble the substrate-binding site of trypsin.     

Prediction of Protein-Protein Interaction Sites Using Electrostatic Desolvation Profiles

Protein-protein complex formation involves removal of water from the interface region. Surface regions with a small free energy penalty for water removal or desolvation may correspond to preferred interaction sites. A method to calculate the electrostatic free energy of placing a neutral low-dielectric probe at various protein surface positions has been designed and applied to characterize putative interaction sites. Based on solutions of the finite-difference Poisson equation, this method also includes long-range electrostatic contributions and the protein solvent boundary shape in contrast to accessible-surface-area-based solvation energies. Calculations on a large set of proteins indicate that in many cases (>90%), the known binding site overlaps with one of the six regions of lowest electrostatic desolvation penalty (overlap with the lowest desolvation region for 48% of proteins). Since the onset of electrostatic desolvation occurs even before direct protein-protein contact formation, it may help guide proteins toward the binding region in the final stage of complex formation. It is interesting that the probe desolvation properties associated with residue types were found to depend to some degree on whether the residue was outside of or part of a binding site. The probe desolvation penalty was on average smaller if the residue was part of a binding site compared to other surface locations. Applications to several antigen-antibody complexes demonstrated that the approach might be useful not only to predict protein interaction sites in general but to map potential antigenic epitopes on protein surfaces.     

Mutational analysis of the primary substrate specificity pocket of complement factor B. Asp(226) is a major structural determinant for p(1)-Arg binding

Factor B is a serine protease, which despite its trypsin-like specificity has Asn instead of the typical Asp at the bottom of the S(1) pocket (position 189, chymotrypsinogen numbering). Asp residues are present at positions 187 and 226 and either one could conceivably provide the negative charge for binding the P(1)-Arg of the substrate. Determination of the crystal structure of the factor B serine protease domain has revealed that the side chain of Asp(226) is within the S(1) pocket, whereas Asp(187) is located outside the pocket. To investigate the possible role of these atypical structural features in substrate binding and catalysis, we constructed a panel of mutants of these residues. Replacement of Asp(187) caused moderate (50-60%) decrease in hemolytic activity, compared with wild type factor B, whereas replacement of Asn(189) resulted in more profound reductions (71-95%). Substitutions at these two positions did not significantly affect assembly of the alternative pathway C3 convertase. In contrast, elimination of the negative charge from Asp(226) completely abrogated hemolytic activity and also affected formation of the C3 convertase. Kinetic analyses of the hydrolysis of a P(1)-Arg containing thioester by selected mutants confirmed that residue Asp(226) is a primary structural determinant for P(1)-Arg binding and catalysis.     

Efficient electrostatic solvation model for protein-fragment docking

A method is presented for the fast evaluation of the binding energy of a protein-small molecule complex with electrostatic solvation. It makes use of a fast preprocessing step based on the assumption that the main contribution to electrostatic desolvation upon ligand binding originates from the displacement of the first shell of water molecules. For a rigid protein, the precomputation of the energy contributions on a set of grids allows the estimation of the energy in solution of about 300 protein-fragment binding modes per second on a personal computer. The docking procedure is applied to five rigid binding sites whose size ranges from 17 residues to a whole protein of 107 amino acids. Using a library of 70 mainly rigid molecules, known micromolar inhibitors or close analogs are docked and prioritized correctly. The docking based rank-ordering of the library requires about 5 h and is proposed as a complementary approach to structure-activity relationships by nuclear magnetic resonance. Proteins 2001;42:256-268.     

Computation of electrostatic complements to proteins: a case of charge stabilized binding.

Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.     

Molecular Determinants of Snurportin 1 Ligand Affinity and Structural Response upon Binding

The transport of large biomolecules such as proteins and RNA across nuclear pore complexes is a field of strong interest and research. Although the basic mechanisms are fairly well understood, the details of the underlying intermolecular interaction within these transport complexes are still unclear. The recognition dynamics and energetics of cargo binding to the transport receptor are not yet resolved. Here, the binding of dimethylated RNA-caps to snurportin 1 is studied by molecular-dynamics simulations. The simulations reveal a strong structural response of the protein upon RNA-cap release. In particular, major rearrangements occur in regions already intrinsically flexible in the holo structure. Additionally, the difference in free energy of binding to snurportin 1 between the two methylation states of the RNA-cap, responsible for the directionality of the transport is quantified. In particular, desolvation of the ligand is revealed as the key-step in binding to snurportin 1. These findings suggest that the binding of m₃G-capped RNA is mainly driven by the enhanced water entropy gain of the solvation shell.     

Probing molecular docking in a charged model binding site

A model binding site was used to investigate charge-charge interactions in molecular docking. This simple site, a small (180A(3)) engineered cavity in cyctochrome c peroxidase (CCP), is negatively charged and completely buried from solvent, allowing us to explore the balance between electrostatic energy and ligand desolvation energy in a system where many of the common approximations in docking do not apply. A database with about 5300 molecules was docked into this cavity. Retrospective testing with known ligands and decoys showed that overall the balance between electrostatic interaction and desolvation energy was captured. More interesting were prospective docking scre"ens that looked for novel ligands, especially those that might reveal problems with the docking and energy methods. Based on screens of the 5300 compound database, both high-scoring and low-scoring molecules were acquired and tested for binding. Out of 16 new, high-scoring compounds tested, 15 were observed to bind. All of these were small heterocyclic cations. Binding constants were measured for a few of these, they ranged between 20microM and 60microM. Crystal structures were determined for ten of these ligands in complex with the protein. The observed ligand geometry corresponded closely to that predicted by docking. Several low-scoring alkyl amino cations were also tested and found to bind. The low docking score of these molecules owed to the relatively high charge density of the charged amino group and the corresponding high desolvation penalty. When the complex structures of those ligands were determined, a bound water molecule was observed interacting with the amino group and a backbone carbonyl group of the cavity. This water molecule mitigates the desolvation penalty and improves the interaction energy relative to that of the "naked" site used in the docking screen. Finally, six low-scoring neutral molecules were also tested, with a view to looking for false negative predictions. Whereas most of these did not bind, two did (phenol and 3-fluorocatechol). Crystal structures for these two ligands in complex with the cavity site suggest reasons for their binding. That these neutral molecules do, in fact bind, contradicts previous results in this site and, along with the alkyl amines, provides instructive false negatives that help identify weaknesses in our scoring functions. Several improvements of these are considered.     

Contributions of water transfer energy to protein‐ligand association and dissociation barriers: Watermap analysis of a series of p38α MAP kinase inhibitors

In our previous work, we proposed that desolvation and resolvation of the binding sites of proteins can serve as the slowest steps during ligand association and dissociation, respectively, and tested this hypothesis on two protein‐ligand systems with known binding kinetics behavior. In the present work, we test this hypothesis on another kinetically‐determined protein‐ligand system—that of p38α and eight Type II BIRB 796 inhibitor analogs. The k_on values among the inhibitor analogs are narrowly distributed (10⁴ ≤ k_on ≤ 10⁵ M^?1 s^?1), suggesting a common rate‐determining step, whereas the k_off values are widely distributed (10^?1 ≤ k_off ≤ 10^?6 s^?1), suggesting a spectrum of rate‐determining steps. We calculated the solvation properties of the DFG‐out protein conformation using an explicit solvent molecular dynamics simulation and thermodynamic analysis method implemented in WaterMap to predict the enthalpic and entropic costs of water transfer to and from bulk solvent incurred upon association and dissociation of each inhibitor. The results suggest that the rate‐determining step for association consists of the transfer of a common set of enthalpically favorable solvating water molecules from the binding site to bulk solvent. The rate‐determining step for inhibitor dissociation consists of the transfer of water from bulk solvent to specific binding site positions that are unfavorably solvated in the apo protein, and evacuated during ligand association. Different sets of unfavorable solvation are evacuated by each ligand, and the observed dissociation barriers are qualitatively consistent with the calculated solvation free energies of those sets.     

Crystal structure of the kringle 2 domain of tissue plasminogen activator at 2.4-A resolution.

The crystal structure of the kringle 2 domain of tissue plasminogen activator was determined and refined at a resolution of 2.43 A. The overall fold of the molecule is similar to that of prothrombin kringle 1 and plasminogen kringle 4; however, there are differences in the lysine binding pocket, and two looping regions, which include insertions in kringle 2, take on very different conformations. Based on a comparison of the overall structural homology between kringle 2 and kringle 4, a new sequence alignment for kringle domains is proposed that results in a division of kringle domains into two groups, consistent with their proposed evolutionary relation. The crystal structure shows a strong interaction between a lysine residue of one molecule and the lysine/fibrin binding pocket of a noncrystallographically related neighbor. This interaction represents a good model of a bound protein ligand and is the first such ligand that has been observed in a kringle binding pocket. The structure shows an intricate network of interactions both among the binding pocket residues and between binding pocket residues and the lysine ligand. A lysine side chain is identified as the positively charged group positioned to interact with the carboxylate of lysine and lysine analogue ligands. In addition, a chloride ion is located in the kringle-kringle interface and contributes to the observed interaction between kringle molecules.     

Interpreting trends in the binding of cyclic ureas to HIV-1 protease

The design of new HIV protease inhibitors requires an improved understanding of the physical basis of inhibitor/protein binding. Here, the binding affinities of seven aliphatic cyclic ureas to HIV-1 protease are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed U-shaped trend of binding free energy as a function of aliphatic chain length. Interestingly, the decrease in affinity for the longest chains is attributable primarily to the energy cost of partly desolvating charged aspartic and arginine groups at the mouths of the active site. Even aliphatic chains too short to contact these charged groups directly are subject to considerable desolvation penalties. We are not aware of other systems where binding affinity trends have been attributed to long-ranged electrostatic desolvation of ionized groups. A generalized Born/surface area solvation model yields a much smaller change in desolvation energy with chain length and, therefore, does not reproduce the experimental binding affinity trends. This result suggests that the generalized Born model should be used with caution for complex, partly desolvated systems like protein binding sites. We also find that changing the assumed protonation state of the active site aspartyl dyad significantly affects the computed binding affinity trends. The protonation state of the aspartyl dyad in the presence of cyclic ureas is discussed in light of the observation that the monoprotonated state reproduces the experimental results best.     

Thrombin inhibitor design: X-ray and solution studies provide a novel P1 determinant

The crystal structures of proflavin and 6-fluorotryptamine thrombin have been completed showing binding of both ligands at the active site S1 pocket. The structure of proflavin:thrombin was confirmatory, while the structure of 6-fluorotryptamine indicated a novel binding mode at the thrombin active site. Furthermore, speculation that the sodium atom identified in an extended solvent channel beneath the S1 pocket may play a role in binding of these ligands was investigated by direct proflavin titrations as well as chromogenic activity measurements as a function of sodium concentration at constant ionic strength. These results suggested a linkage between the sodium site and the S1 pocket. This observation could be due to a simple ionic interaction between Asp189 and the sodium ion or a more complicated structural rearrangement of the thrombin S1 pocket. Finally, the unique binding mode of 6-fluorotryptamine provides ideas toward the design of a neutrally charged thrombin inhibitor.     

Flexible docking of DNA fragments and actinocin derivatives

We have applied molecular docking methods to systems containing nucleic acids as targets and biologically active substances as ligands. The complexes of DNA fragments and actinocin derivatives with different lengths of aminoalkyl side chains were obtained by molecular docking. It was observed that actinocin derivatives could form energetically favourable complexes with DNA both as intercalators and minor groove binders. It was shown that small changes in the binding energy (～1?kcal/mol) could result in complexes with substantially different structure. The complexes of actinocin derivatives and DNA fragments were stabilized by hydrogen bonding upon intercalation and minor groove binding. It was found that the change of solvent-accessible surface area upon binding of the actinocin derivative to DNA linear increased with the growth of methylene groups' number in ligand side chains. The solvation energy change upon binding of actinocin derivatives to DNA calculated by the WSAS method was favourable in the case of small uncharged ligands and unfavourable for positively charged ligands.     

Probing ligand recognition in the decarboxylase antibody 21D8: implications for the catalytic mechanism.

Antibody 21D8, which was elicited with a naphthalene-1,5-disulfonate hapten, catalyzes the medium-sensitive decarboxylation of 5-nitro-3-carboxybenzisoxazole. Structural studies on the hapten-antibody complex show that the active site contains two anion binding pockets separated by a hydrophobic region. To gain further insight into the ligand binding and catalytic mechanism of 21D8, six site-directed mutants were prepared, four for investigating the role of each of the two hapten sulfonate binding sites and two for examining packing interactions between bound ligands and the binding pocket. With the exception of an Arg(L46)Met substitution in the more deeply buried sulfonate binding pocket, modification of the active site resulted in reductions in catalytic efficiency (k(cat)/k(uncat)), ranging between 3- and 23-fold. Importantly, and contrary to predictions based on computational docking experiments, the differential effects of the individual mutations on the K(m), K(TS), and K(product) parameters suggest that only substrate binding modes which place the carboxylate group in the more solvent-exposed sulfonate binding site are catalytically relevant. Such an orientation would permit a potentially significant interaction between the developing oxyanion in the transition state and the side chain of Arg(L96). Incomplete desolvation of the carboxylate in this orientation may also help explain the modest efficiency of 21D8 compared to the most accelerating aprotic dipolar organic solvents.