Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment.

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid-base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

Distinct immune stimulatory effects of anti-human VISTA antibodies are determined by Fc-receptor interaction.

VISTA (PD-1H) is an immune regulatory molecule considered part of the next wave of immuno-oncology targets. VISTA is an immunoglobulin (Ig) superfamily cell surface molecule mainly expressed on myeloid cells, and to some extent on NK cells and T cells. In previous preclinical studies, some VISTA-targeting antibodies provided immune inhibitory signals, while other antibodies triggered immune stimulatory signals. Importantly, for therapeutic antibodies, the isotype backbone can have a strong impact on antibody function. To elucidate the mode of action of immune stimulatory anti-VISTA antibodies, we studied three different anti-human VISTA antibody clones, each on three different IgG isotypes currently used for therapeutic antibodies: unaltered IgG1 (IgG1-WT), IgG1-KO (IgG1-LL234,235AA-variant with reduced Fc-effector function), and IgG4-Pro (IgG4- S228P-variant with stabilized hinge region). Antibody functionality was analysed in mixed leukocyte reaction (MLR) of human peripheral blood mononuclear cells (PBMCs), as a model system for ongoing immune reactions, on unstimulated human PBMCs, as a model system for a resting immune system, and also on acute myeloid leukemia (AML) patient samples to evaluate anti-VISTA antibody effects on primary tumor material. The functions of three anti-human VISTA antibodies were determined by their IgG isotype backbones. An MLR of healthy donor PBMCs was effectively augmented by anti-VISTA-IgG4-Pro and anti-VISTA-IgG1-WT antibodies, as indicated by increased levels of cytokines, T cell activation markers and T cell proliferation. However, in a culture of unstimulated PBMCs of single healthy donors, only anti-VISTA-IgG1-WT antibodies increased the activation marker HLA-DR on resting myeloid cells, and chemokine levels. Interestingly, interactions with different Fc-receptors were required for these effects, namely CD64 for augmentation of MLR, and CD16 for activation of resting myeloid cells. Furthermore, anti-VISTA-IgG1-KO antibodies had nearly no impact in any model system. Similarly, in AML patient samples, anti-VISTA-antibody on IgG4-Pro backbone, but not on IgG1-KO backbone, increased interactions, as a novel readout of activity, between immune cells and CD34+ AML cancer cells. In conclusion, the immune stimulatory effects of antagonistic anti-VISTA antibodies are defined by the antibody isotype and interaction with different Fc-gamma-receptors, highlighting the importance of understanding these interactions when designing immune stimulatory antibody therapeutics for immuno-oncology applications.

Dead-End Elimination with a Polarizable Force Field Repacks PCNA Structures.

A balance of van der Waals, electrostatic, and hydrophobic forces drive the folding and packing of protein side chains. Although such interactions between residues are often approximated as being pairwise additive, in reality, higher-order many-body contributions that depend on environment drive hydrophobic collapse and cooperative electrostatics. Beginning from dead-end elimination, we derive the first algorithm, to our knowledge, capable of deterministic global repacking of side chains compatible with many-body energy functions. The approach is applied to seven PCNA x-ray crystallographic data sets with resolutions 2.5-3.8 Å (mean 3.0 Å) using an open-source software. While PDB_REDO models average an Rfree value of 29.5% and MOLPROBITY score of 2.71 Å (77th percentile), dead-end elimination with the polarizable AMOEBA force field lowered Rfree by 2.8-26.7% and improved mean MOLPROBITY score to atomic resolution at 1.25 Å (100th percentile). For structural biology applications that depend on side-chain repacking, including x-ray refinement, homology modeling, and protein design, the accuracy limitations of pairwise additivity can now be eliminated via polarizable or quantum mechanical potentials.

Disassembly of all SNARE complexes by N-ethylmaleimide-sensitive factor (NSF) is initiated by a conserved 1:1 interaction between α-soluble NSF attachment protein (SNAP) and SNARE complex.

Vesicle trafficking in eukaryotic cells is facilitated by SNARE-mediated membrane fusion. The ATPase NSF (N-ethylmaleimide-sensitive factor) and the adaptor protein α-SNAP (soluble NSF attachment protein) disassemble all SNARE complexes formed throughout different pathways, but the effect of SNARE sequence and domain variation on the poorly understood disassembly mechanism is unknown. By measuring SNARE-stimulated ATP hydrolysis rates, Michaelis-Menten constants for disassembly, and SNAP-SNARE binding constants for four different ternary SNARE complexes and one binary complex, we found a conserved mechanism, not influenced by N-terminal SNARE domains. α-SNAP and the ternary SNARE complex form a 1:1 complex as revealed by multiangle light scattering. We propose a model of NSF-mediated disassembly in which the reaction is initiated by a 1:1 interaction between α-SNAP and the ternary SNARE complex, followed by NSF binding. Subsequent additional α-SNAP binding events may occur as part of a processive disassembly mechanism.

Processive ATP-driven substrate disassembly by the N-ethylmaleimide-sensitive factor (NSF) molecular machine.

SNARE proteins promote membrane fusion by forming a four-stranded parallel helical bundle that brings the membranes into close proximity. Post-fusion, the complex is disassembled by an AAA+ ATPase called N-ethylmaleimide-sensitive factor (NSF). We present evidence that NSF uses a processive unwinding mechanism to disassemble SNARE proteins. Using a real-time disassembly assay based on fluorescence dequenching, we correlate NSF-driven disassembly rates with the SNARE-activated ATPase activity of NSF. Neuronal SNAREs activate the ATPase rate of NSF by ∼26-fold. One SNARE complex takes an average of ∼5 s to disassemble in a process that consumes ∼50 ATP. Investigations of substrate requirements show that NSF is capable of disassembling a truncated SNARE substrate consisting of only the core SNARE domain, but not an unrelated four-stranded coiled-coil. NSF can also disassemble an engineered double-length SNARE complex, suggesting a processive unwinding mechanism. We further investigated processivity using single-turnover experiments, which show that SNAREs can be unwound in a single encounter with NSF. We propose a processive helicase-like mechanism for NSF in which ∼1 residue is unwound for every hydrolyzed ATP molecule.

Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site.

Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol-Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations.

Noncanonical E2 recruitment by the autophagy E1 revealed by Atg7-Atg3 and Atg7-Atg10 structures.

Core functions of autophagy are mediated by ubiquitin-like protein (UBL) cascades, in which a homodimeric E1 enzyme, Atg7, directs the UBLs Atg8 and Atg12 to their respective E2 enzymes, Atg3 and Atg10. Crystallographic and mutational analyses of yeast (Atg7-Atg3)(2) and (Atg7-Atg10)(2) complexes reveal noncanonical, multisite E1-E2 recognition in autophagy. Atg7's unique N-terminal domain recruits distinctive elements from the Atg3 and Atg10 'backsides'. This, along with E1 and E2 conformational variability, allows presentation of 'frontside' Atg3 and Atg10 active sites to the catalytic cysteine in the C-terminal domain from the opposite Atg7 protomer in the homodimer. Despite different modes of binding, the data suggest that common principles underlie conjugation in both noncanonical and canonical UBL cascades, whereby flexibly tethered E1 domains recruit E2s through surfaces remote from their active sites to juxtapose the E1 and E2 catalytic cysteines.

Quantitative, directional measurement of electric field heterogeneity in the active site of ketosteroid isomerase.

Understanding the electrostatic forces and features within highly heterogeneous, anisotropic, and chemically complex enzyme active sites and their connection to biological catalysis remains a longstanding challenge, in part due to the paucity of incisive experimental probes of electrostatic properties within proteins. To quantitatively assess the landscape of electrostatic fields at discrete locations and orientations within an enzyme active site, we have incorporated site-specific thiocyanate vibrational probes into multiple positions within bacterial ketosteroid isomerase. A battery of X-ray crystallographic, vibrational Stark spectroscopy, and NMR studies revealed electrostatic field heterogeneity of 8 MV/cm between active site probe locations and widely differing sensitivities of discrete probes to common electrostatic perturbations from mutation, ligand binding, and pH changes. Electrostatic calculations based on active site ionization states assigned by literature precedent and computational pK(a) prediction were unable to quantitatively account for the observed vibrational band shifts. However, electrostatic models of the D40N mutant gave qualitative agreement with the observed vibrational effects when an unusual ionization of an active site tyrosine with a pK(a) near 7 was included. UV-absorbance and (13)C NMR experiments confirmed the presence of a tyrosinate in the active site, in agreement with electrostatic models. This work provides the most direct measure of the heterogeneous and anisotropic nature of the electrostatic environment within an enzyme active site, and these measurements provide incisive benchmarks for further developing accurate computational models and a foundation for future tests of electrostatics in enzymatic catalysis.

High-resolution analysis of Zn(2+) coordination in the alkaline phosphatase superfamily by EXAFS and x-ray crystallography.

Comparisons among evolutionarily related enzymes offer opportunities to reveal how structural differences produce different catalytic activities. Two structurally related enzymes, Escherichia coli alkaline phosphatase (AP) and Xanthomonas axonopodis nucleotide pyrophosphatase/phosphodiesterase (NPP), have nearly identical binuclear Zn(2+) catalytic centers but show tremendous differential specificity for hydrolysis of phosphate monoesters or phosphate diesters. To determine if there are differences in Zn(2+) coordination in the two enzymes that might contribute to catalytic specificity, we analyzed both x-ray absorption spectroscopic and x-ray crystallographic data. We report a 1.29-Å crystal structure of AP with bound phosphate, allowing evaluation of interactions at the AP metal site with high resolution. To make systematic comparisons between AP and NPP, we measured zinc extended x-ray absorption fine structure for AP and NPP in the free-enzyme forms, with AMP and inorganic phosphate ground-state analogs and with vanadate transition-state analogs. These studies yielded average zinc-ligand distances in AP and NPP free-enzyme forms and ground-state analog forms that were identical within error, suggesting little difference in metal ion coordination among these forms. Upon binding of vanadate to both enzymes, small increases in average metal-ligand distances were observed, consistent with an increased coordination number. Slightly longer increases were observed in NPP relative to AP, which could arise from subtle rearrangements of the active site or differences in the geometry of the bound vanadyl species. Overall, the results suggest that the binuclear Zn(2+) catalytic site remains very similar between AP and NPP during the course of a reaction cycle.

Reintroducing electrostatics into macromolecular crystallographic refinement: application to neutron crystallography and DNA hydration.

Most current crystallographic structure refinements augment the diffraction data with a priori information consisting of bond, angle, dihedral, planarity restraints, and atomic repulsion based on the Pauli exclusion principle. Yet, electrostatics and van der Waals attraction are physical forces that provide additional a priori information. Here, we assess the inclusion of electrostatics for the force field used for all-atom (including hydrogen) joint neutron/X-ray refinement. Two DNA and a protein crystal structure were refined against joint neutron/X-ray diffraction data sets using force fields without electrostatics or with electrostatics. Hydrogen-bond orientation/geometry favors the inclusion of electrostatics. Refinement of Z-DNA with electrostatics leads to a hypothesis for the entropic stabilization of Z-DNA that may partly explain the thermodynamics of converting the B form of DNA to its Z form. Thus, inclusion of electrostatics assists joint neutron/X-ray refinements, especially for placing and orienting hydrogen atoms.

Polarizable Atomic Multipole X-Ray Refinement: Particle Mesh Ewald Electrostatics for Macromolecular Crystals.

Refinement of macromolecular models from X-ray crystallography experiments benefits from prior chemical knowledge at all resolutions. As the quality of the prior chemical knowledge from quantum or classical molecular physics improves, in principle so will resulting structural models. Due to limitations in computer performance and electrostatic algorithms, commonly used macromolecules X-ray crystallography refinement protocols have had limited support for rigorous molecular physics in the past. For example, electrostatics is often neglected in favor of nonbonded interactions based on a purely repulsive van der Waals potential. In this work we present advanced algorithms for desktop workstations that open the door to X-ray refinement of even the most challenging macromolecular data sets using state-of-the-art classical molecular physics. First we describe theory for particle mesh Ewald (PME) summation that consistently handles the symmetry of all 230 space groups, replicates of the unit cell such that the minimum image convention can be used with a real space cutoff of any size and the combination of space group symmetry with replicates. An implementation of symmetry accelerated PME for the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) force field is presented. Relative to a single CPU core performing calculations on a P1 unit cell, our AMOEBA engine called Force Field X (FFX) accelerates energy evaluations by more than a factor of 24 on an 8-core workstation with a Tesla GPU coprocessor for 30 structures that contain 240 000 atoms on average in the unit cell. The benefit of AMOEBA electrostatics evaluated with PME for macromolecular X-ray crystallography refinement is demonstrated via rerefinement of 10 crystallographic data sets that range in resolution from 1.7 to 4.5 Å. Beginning from structures obtained by local optimization without electrostatics, further optimization using AMOEBA with PME electrostatics improved agreement of the model with the data (Rfree was lowered by 0.5%), improved geometric features such as favorable (ϕ, ψ) backbone conformations, and lowered the average potential energy per residue by over 10 kcal/mol. Furthermore, the MolProbity structure validation tool indicates that the geometry of these rerefined structures is consistent with X-ray crystallographic data collected up to 2.2 Å, which is 0.9 Å better than the actual mean quality (3.1 Å). We conclude that polarizable AMOEBA-assisted X-ray refinement offers advantages to methods that neglect electrostatics and is now efficient enough for routine use.

Polarizable atomic multipole x-ray refinement: hydration geometry and application to macromolecules.

We recently developed a polarizable atomic multipole refinement method assisted by the AMOEBA force field for macromolecular crystallography. Compared to standard refinement procedures, the method uses a more rigorous treatment of x-ray scattering and electrostatics that can significantly improve the resultant information contained in an atomic model. We applied this method to high-resolution lysozyme and trypsin data sets, and validated its utility for precisely describing biomolecular electron density, as indicated by a 0.4-0.6% decrease in the R- and R(free)-values, and a corresponding decrease in the relative energy of 0.4-0.8 Kcal/mol/residue. The re-refinements illustrate the ability of force-field electrostatics to orient water networks and catalytically relevant hydrogens, which can be used to make predictions regarding active site function, activity, and protein-ligand interaction energies. Re-refinement of a DNA crystal structure generates the zigzag spine pattern of hydrogen bonding in the minor groove without manual intervention. The polarizable atomic multipole electrostatics model implemented in the AMOEBA force field is applicable and informative for crystal structures solved at any resolution.

Dissecting the paradoxical effects of hydrogen bond mutations in the ketosteroid isomerase oxyanion hole.

The catalytic importance of enzyme active-site interactions is frequently assessed by mutating specific residues and measuring the resulting rate reductions. This approach has been used in bacterial ketosteroid isomerase to probe the energetic importance of active-site hydrogen bonds donated to the dienolate reaction intermediate. The conservative Tyr16Phe mutation impairs catalysis by 10(5)-fold, far larger than the effects of hydrogen bond mutations in other enzymes. However, the less-conservative Tyr16Ser mutation, which also perturbs the Tyr16 hydrogen bond, results in a less-severe 10(2)-fold rate reduction. To understand the paradoxical effects of these mutations and clarify the energetic importance of the Tyr16 hydrogen bond, we have determined the 1.6-A resolution x-ray structure of the intermediate analogue, equilenin, bound to the Tyr16Ser mutant and measured the rate effects of mutating Tyr16 to Ser, Thr, Ala, and Gly. The nearly identical 200-fold rate reductions of these mutations, together with the 6.4-A distance observed between the Ser16 hydroxyl and equilenin oxygens in the x-ray structure, strongly suggest that the more moderate rate effect of this mutant is not due to maintenance of a hydrogen bond from Ser at position 16. These results, additional spectroscopic observations, and prior structural studies suggest that the Tyr16Phe mutation results in unfavorable interactions with the dienolate intermediate beyond loss of a hydrogen bond, thereby exaggerating the apparent energetic benefit of the Tyr16 hydrogen bond relative to the solution reaction. These results underscore the complex energetics of hydrogen bonding interactions and site-directed mutagenesis experiments.

Polarizable atomic multipole X-ray refinement: application to peptide crystals.

Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (FFT) of Cartesian Gaussian multipoles. Relative to direct summation, the FFT approach can give a speedup of more than an order of magnitude for aspherical refinement of ultrahigh-resolution data sets. Use of a sublattice formalism makes the method highly parallelizable. Application of the Cartesian Gaussian multipole scattering model to a series of four peptide crystals using multipole coefficients from the AMOEBA force field demonstrates that AMOEBA systematically underestimates electron density at bond centers. For the trigonal and tetrahedral bonding geometries common in organic chemistry, an atomic multipole expansion through hexadecapole order is required to explain bond electron density. Alternatively, the addition of interatomic scattering (IAS) sites to the AMOEBA-based density captured bonding effects with fewer parameters. For a series of four peptide crystals, the AMOEBA-IAS model lowered R(free) by 20-40% relative to the original spherically symmetric scattering model.

A potent peptidomimetic inhibitor of botulinum neurotoxin serotype A has a very different conformation than SNAP-25 substrate.

Botulinum neurotoxin serotype A is the most lethal of all known toxins. Here, we report the crystal structure, along with SAR data, of the zinc metalloprotease domain of BoNT/A bound to a potent peptidomimetic inhibitor (K(i)=41 nM) that resembles the local sequence of the SNAP-25 substrate. Surprisingly, the inhibitor adopts a helical conformation around the cleavage site, in contrast to the extended conformation of the native substrate. The backbone of the inhibitor's P1 residue displaces the putative catalytic water molecule and concomitantly interacts with the "proton shuttle" E224. This mechanism of inhibition is aided by residue contacts in the conserved S1' pocket of the substrate binding cleft and by the induction of new hydrophobic pockets, which are not present in the apo form, especially for the P2' residue of the inhibitor. Our inhibitor is specific for BoNT/A as it does not inhibit other BoNT serotypes or thermolysin.

Comparative enzymology in the alkaline phosphatase superfamily to determine the catalytic role of an active-site metal ion.

Mechanistic models for biochemical systems are frequently proposed from structural data. Site-directed mutagenesis can be used to test the importance of proposed functional sites, but these data do not necessarily indicate how these sites contribute to function. In this study, we applied an alternative approach to the catalytic mechanism of alkaline phosphatase (AP), a widely studied prototypical bimetallo enzyme. A third metal ion site in AP has been suggested to provide general base catalysis, but comparison of AP with an evolutionarily related enzyme casts doubt on this model. Removal of this metal site from AP has large differential effects on reactions of cognate and promiscuous substrates, and the results are inconsistent with general base catalysis. Instead, these and additional results suggest that the third metal ion stabilizes the transferred phosphoryl group in the transition state. These results establish a new mechanistic model for this prototypical bimetallo enzyme and demonstrate the power of a comparative approach for probing biochemical function.

Arginine coordination in enzymatic phosphoryl transfer: evaluation of the effect of Arg166 mutations in Escherichia coli alkaline phosphatase.

Arginine residues are commonly found in the active sites of enzymes catalyzing phosphoryl transfer reactions. Numerous site-directed mutagenesis experiments establish the importance of these residues for efficient catalysis, but their role in catalysis is not clear. To examine the role of arginine residues in the phosphoryl transfer reaction, we have measured the consequences of mutations to arginine 166 in Escherichia coli alkaline phosphatase on hydrolysis of ethyl phosphate, on individual reaction steps in the hydrolysis of the covalent enzyme-phosphoryl intermediate, and on thio substitution effects. The results show that the role of the arginine side chain extends beyond its positive charge, as the Arg166Lys mutant is as compromised in activity as Arg166Ser. Through measurement of individual reaction steps, we construct a free energy profile for the hydrolysis of the enzyme-phosphate intermediate. This analysis indicates that the arginine side chain strengthens binding by approximately 3 kcal/mol and provides an additional 1-2 kcal/mol stabilization of the chemical transition state. A 2.1 A X-ray diffraction structure of Arg166Ser AP is presented, which shows little difference in enzyme structure compared to the wild-type enzyme but shows a significant reorientation of the bound phosphate. Altogether, these results support a model in which the arginine contributes to catalysis through binding interactions and through additional transition state stabilization that may arise from complementarity of the guanidinum group to the geometry of the trigonal bipyramidal transition state.

Rab and Arl GTPase family members cooperate in the localization of the golgin GCC185.

GCC185 is a large coiled-coil protein at the trans Golgi network that is required for receipt of transport vesicles inbound from late endosomes and for anchoring noncentrosomal microtubules that emanate from the Golgi. Here, we demonstrate that recruitment of GCC185 to the Golgi is mediated by two Golgi-localized small GTPases of the Rab and Arl families. GCC185 binds Rab6, and mutation of residues needed for Rab binding abolishes Golgi localization. The crystal structure of Rab6 bound to the GCC185 Rab-binding domain reveals that Rab6 recognizes a two-fold symmetric surface on a coiled coil immediately adjacent to a C-terminal GRIP domain. Unexpectedly, Rab6 binding promotes association of Arl1 with the GRIP domain. We present a structure-derived model for dual GTPase membrane attachment that highlights the potential ability of Rab GTPases to reach binding partners at a significant distance from the membrane via their unstructured and membrane-anchored, hypervariable domains.

Crystal structures of the luciferase and green fluorescent protein from Renilla reniformis.

Due to its ability to emit light, the luciferase from Renilla reniformis (RLuc) is widely employed in molecular biology as a reporter gene in cell culture experiments and small animal imaging. To accomplish this bioluminescence, the 37-kDa enzyme catalyzes the degradation of its substrate coelenterazine in the presence of molecular oxygen, resulting in the product coelenteramide, carbon dioxide, and the desired photon of light. We successfully crystallized a stabilized variant of this important protein (RLuc8) and herein present the first structures for any coelenterazine-using luciferase. These structures are based on high-resolution data measured to 1.4 A and demonstrate a classic alpha/beta-hydrolase fold. We also present data of a coelenteramide-bound luciferase and reason that this structure represents a secondary conformational form following shift of the product out of the primary active site. During the course of this work, the structure of the luciferase's accessory green fluorescent protein (RrGFP) was also determined and shown to be highly similar to that of Aequorea victoria GFP.

Structural and functional comparisons of nucleotide pyrophosphatase/phosphodiesterase and alkaline phosphatase: implications for mechanism and evolution.

The rapid expansion of the amount of genomic and structural data has provided many examples of enzymes with evolutionarily related active sites that catalyze different reactions. Functional comparisons of these active sites can provide insight into the origins of the enormous catalytic proficiency of enzymes and the evolutionary changes that can lead to different enzyme activities. The alkaline phosphatase (AP) superfamily is an ideal system to use in making such comparisons given the extensive data available on both nonenzymatic and enzymatic phosphoryl transfer reactions. Some superfamily members, such as AP itself, preferentially hydrolyze phosphate monoesters, whereas others, such as nucleotide pyrophosphatase/phosphodiesterase (NPP), preferentially hydrolyze phosphate diesters. We have measured rate constants for NPP-catalyzed hydrolysis of phosphate diesters and monoesters. NPP preferentially catalyzes diester hydrolysis by factors of 10(2)-10(6), depending on the identity of the diester substrate. To identify features of the NPP active site that could lead to preferential phosphate diester hydrolysis, we have determined the structure of NPP in the absence of ligands and in complexes with vanadate and AMP. Comparisons to existing structures of AP reveal bimetallo cores that are structurally indistinguishable, but there are several distinct structural features outside of the conserved bimetallo site. The structural and functional data together suggest that some of these distinct functional groups provide specific substrate binding interactions, whereas others tune the properties of the bimetallo active site itself to discriminate between phosphate diester and monoester substrates.