Conformational changes in HIV-1 proteinase: effect of protonation of the active center on conformation of HIV-1 proteinase in water

At weak acidic pH, where HIV-1 proteinase is most stable and active, its catalytic Asp 25/25' dyad shares one proton. At a physiological pH the dyad is deprotonated, however, 2 ns molecular dynamics simulations of the HIV-1 protease with monoprotonated and deprotonated Asp25/25' dyad is performed, in order to investigate the influence of Asp25/25' protonation state on the proteinase dynamics. For net charge neutralization the 4 Cl- ions were included. In case of deprotonated active site the significant tertiary structure deviation of HIV-1 PR structure from crystal structure is observed, while in the monoprotonated one the tertiary structure fluctuates near starting structure. Possible mechanism of the influence of the Asp25/25' protonation state on proteinase dynamics is discussed.     

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Insights into the structure and dynamics of the dinuclear zinc beta-lactamase site from Bacteroides fragilis

Herein, we report quantum chemical calculations and molecular dynamics (MD) simulations of the dinuclear form of the Bacteroides fragilis zinc beta-lactamase. We studied four different configurations which differ in the protonation state of the Asp103 residue and in the presence or absence of a Zn1-OH-Zn2 bridge. The flexibility of the Zn1-OH-Zn2 bridge was studied by means of quantum mechanical (QM) calculations on cluster models while the relative stabilities of the different configurations were estimated from QM linear scaling calculations on the enzyme. Contacts between important residues (Cys104, Asp69, Lys185, etc.), the solvation of the zinc ions, and the conformation of the active site beta-hairpin loop were characterized by the MD analyses. The influence of the buried sodium ion close to the Zn2 position was investigated by carrying out a secondary simulation where the sodium ion was replaced with an internal water molecule. The comparative structural analyses among the different MD trajectories augmented with energetic calculations have demonstrated that the B. fragilis protein efficiently binds the internal Na(+) ion observed crystallographically. Moreover, we found that when Asp103 is unprotonated, a rigid Zn1-OH-Zn2 bridge results, while for neutral Asp103, a fluctuating Zn1-Zn2 distance was possible via the breaking and formation of the Zn1-OH-Zn2 bridge. The mechanistic implications of these observations are discussed in detail.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Litchi chinensis-derived terpenoid as anti-HIV-1 protease agent: structural design from molecular dynamics simulations

The molecular structures of the binding between human immunodeficiency virus-1 protease (HIV-1PR) and various inhibitors including existing extensive natural products extracts have been investigated for anti-HIV drug development. In this study, the binding of HIV-1PR and a terpenoid from Litchi chinensis extracts (3-oxotrirucalla-7,24-dien-21-oic acid) was investigated in order to clarify the inhibition effectiveness of this compound. Molecular dynamics (MD) simulations of HIV-1PR complex with 3-oxotrirucalla-7,24-dien-21-oic acid were performed including water molecules. The MD simulation results indicated the formation of hydrogen bonds between the oxygen atoms of the inhibitor and the catalytic aspartates, which are commonly found in inhibitors–protease complexes. On the other hand, no hydrogen bonding of this particular inhibitor to the flap region was found. In addition, the radial distribution function of water oxygens around the catalytic carboxylate nitrogens of Asp29 and Asp30 suggests that at least one or two water molecules are in the active site region whereas direct interaction of the inhibitor was found for catalytic carboxylate oxygen of Asp25. The results of this simulation, in comparison with the structures of other HIV-PR inhibitor complexes, could lead to a better understanding of the activity of 3-oxotrirucalla-7,24-dien-21-oic acid.     

Determinants of conductance of a bacterial voltage-gated sodium channel

Through molecular dynamics (MD) and free energy simulations in electric fields, we examine the factors influencing conductance of bacterial voltage-gated sodium channel Na_vMs. The channel utilizes four glutamic acid residues in the selectivity filter (SF). Previously, we have shown, through constant pH and free energy calculations of pKa values, that fully deprotonated, singly protonated, and doubly protonated states are all feasible at physiological pH, depending on how many ions are bound in the SF. With 173 MD simulations of 450 or 500 ns and additional free energy simulations, we determine that the conductance is highest for the deprotonated state and decreases with each additional proton bound. We also determine that the pKa value of the four glutamic residues for the transition between deprotonated and singly protonated states is close to the physiological pH and that there is a small voltage dependence. The pKa value and conductance trends are in agreement with experimental work on bacterial Na_v channels, which show a decrease in maximal conductance with lowering of pH, with pKa in the physiological range. We examine binding sites for Na⁺ in the SF, compare with previous work, and note a dependence on starting structures. We find that narrowing of the gate backbone to values lower than the crystal structure's backbone radius reduces the conductance, whereas increasing the gate radius further does not affect the conductance. Simulations with some amount of negatively charged lipids as opposed to purely neutral lipids increases the conductance, as do simulations at higher voltages.     

Active site binding modes of curcumin in HIV-1 protease and integrase

Structure models for the interaction of curcumin with HIV-1 integrase (IN) and protease (PR) were investigated using computational docking. Curcumin was found to bind preferentially in similar ways to the active sites of both IN and PR. For IN, the binding site is formed by residues Asp64, His67, Thr66, Glu92, Thr93, Asp116, Ser119, Asn120, and Lys159. Docked curcumin contacts the catalytic residues adjacent to Asp116 and Asp64, and near the divalent metal (Mg2+). In the PR docking, the curcumin structure fitted well to the active site, interacting with residues Asp25, Asp29, Asp30, Gly27', Asp29', and Asp30'. The results suggest that o-hydroxyl and/or keto-enol structures are important for both IN and PR inhibitory actions. The symmetrical structure of curcumin seems to play an important role for binding to the PR protein, whereas the keto-enol and only one side of the terminal o-hydroxyl showed tight binding to the IN active site.     

Human immunodeficiency virus-1 protease. 2. Use of pH rate studies and solvent kinetic isotope effects to elucidate details of chemical mechanism

The pH dependence of the peptidolytic reaction of recombinant human immunodeficiency virus type 1 protease has been examined over a pH range of 3-7 for four oligopeptide substrates and two competitive inhibitors. The pK values obtained from the pKis vs pH profiles for the unprotonated and protonated active-site aspartyl groups, Asp-25 and Asp-25', in the monoprotonated enzyme form were 3.1 and 5.2, respectively. Profiles of log V/K vs pH for all four substrates were "bell-shaped" in which the pK values for the unprotonated and protonated aspartyl residues were 3.4-3.7 and 5.5-6.5, respectively. Profiles of log V vs pH for these substrates were "wave-shaped" in which V was shifted to a constant lower value upon protonation of a residue of pK = 4.2-5.2. These results indicate that substrates bind only to a form of HIV-1 protease in which one of the two catalytic aspartyl residues is protonated. Solvent kinetic isotope effects were measured over a pH (D) range of 3-7 for two oligopeptide substrates, Ac-Arg-Ala-Ser-Gln-Asn-Tyr-Pro-Val-Val-NH2 and Ac-Ser-Gln-Asn-Tyr-Pro-Val-Val-NH2. The pH-independent value for DV/K was 1.0 for both substrates, and DV = 1.5-1.7 and 2.2-3.2 at low and high pH (D), respectively. The attentuation of both V and DV at low pH (D) is consistent with a change in rate-limiting step from a chemical one at high pH (D) to one in which a product release step or an enzyme isomerization step becomes partly rate-limiting at low pH (D). Proton inventory data is in accord with the concerted transfer of two protons in the transition state of a rate-limiting chemical step in which the enzyme-bound amide hydrate adduct collapses to form the carboxylic acid and amine products.     

Cooperative fluctuations of unliganded and substrate-bound HIV-1 protease: a structure-based analysis on a variety of conformations from crystallography and molecular dynamics simulations

The dynamics of HIV-1 protease, both in unliganded and substrate-bound forms have been analyzed by using an analytical method, Gaussian network model (GNM). The method is applied to different conformations accessible to the protein backbone in the native state, observed in crystal structures and snapshots from fully atomistic molecular dynamics (MD) simulation trajectories. The modes of motion obtained from GNM on different conformations of HIV-1 protease are conserved throughout the MD simulations. The flaps and 40's loop of the unliganded HIV-1 protease structure are identified as the most mobile regions. However, in the liganded structure these flaps lose mobility, and terminal regions of the monomers become more flexible. Analysis of the fast modes shows that residues important for stability are in the same regions of all the structures examined. Among these, Gly86 appears to be a key residue for stability. The contribution of residues in the active site region and flaps to the stability is more pronounced in the substrate-bound form than in the unliganded form. The convergence of modes in GNM to similar regions of HIV-1 protease, regardless of the conformation of the protein, supports the robustness of GNM as a potentially useful and predictive tool.     

Molecular dynamics simulations of the free and inhibitor-bound cruzain systems in aqueous solvent: insights on the inhibition mechanism in acidic pH

The major cysteine protease of Trypanosoma cruzi, cruzain (CRZ), has been described as a therapeutic target for Chagas’ disease, which affects millions of people worldwide. Thus, a series of CRZ inhibitors has been studied, including a new competitive inhibitor, Nequimed176 (NEQ176). Nevertheless, the structural and dynamic basis for CRZ inhibition remains unclear. Hoping to contribute to this ever-growing understanding of timescale dynamics in the CRZ inhibition mechanism, we have performed the first study using 100 ns of molecular dynamics (MD) simulations of two CRZ systems in an aqueous solvent under pH 5.5: CRZ in the apo form (ligand free) and CRZ complexed to NEQ176. According to the MD simulations, the enzyme adopts an open conformation in the apo form and a closed conformation in the NEQ176–CRZ complex. We also suggest that this closed conformation is related to the hydrogen-bonding interactions between NEQ176 and CRZ, which occurs through key residues, mainly Gly66, Met68, Asn69, and Leu160. In addition, the cross-correlation analysis shows evidence of the correlated motions among Ala110–Asp140, Leu160–Gly189, and Glu190–Gly215 subdomains, as well as, the movements related to Ala1–Thr59 and Asp60–Pro90 regions seem to be crucial for CRZ activity.     

New algorithm for calculation of N-H peptide bond order parameter by the method of molecular dynamic modeling

A new algorithm was proposed for calculation of N-H bond order parameter from molecular dynamics (MD) trajectories. The MD simulation of the HIV-1 protease (3.4.23.16) with monoprotonated active centre was performed for the algorithm verification. It has been shown that the protease in aqueous solution at pH 3.5-6.5 adopts a set of conformations, which are intermediate between the "open" and "closed" ones. The N-H bond order parameters calculated from the MD trajectory are in agreement with experimental NMR data.     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

The critical role of the 185-189-loop in the factor Xa interaction with Na+ and factor Va in the prothrombinase complex

The S1 site (Asp(189)) of factor Xa (fXa) is located on a loop (residues 185-189) that contains three solvent-exposed charged residues (Asp(185), Lys(186), and Glu(188)) below the active-site pocket of the protease. To investigate the role of these residues in the catalytic function of fXa, we expressed three mutants of the protease in which the charges of these residues were neutralized by their substitutions with Ala (D185A, K186A, and E188A). Kinetic studies revealed that E188A has a normal catalytic activity toward small synthetic and natural substrates and inhibitors of fXa; however, the same activities were slightly ( approximately 2-fold) and dramatically ( approximately 20-50-fold) impaired for the D185A and K186A mutants, respectively. Further studies revealed that the affinity of D185A and K186A for interaction with Na(+) has also been altered, with a modest impairment ( approximately 2-fold) for the former and a dramatic impairment for the latter mutant. Both prothrombinase and direct binding studies indicated that K186A also has an approximately 6-fold impaired affinity for factor Va. Interestingly, a saturating concentration of factor Va restored the catalytic defect of K186A in reactions with prothrombin and the recombinant tick anticoagulant peptide that is known to interact with the Na(+) loop of fXa, but not with other substrates. These results suggest that factor Va interacts with 185-189-loop for fXa, which is energetically linked to the Na(+)-binding site of the protease.     

Computational study of IAG-nucleoside hydrolase: determination of the preferred ground state conformation and the role of active site residues

The mechanism of action of inosine-adenosine-guanosine nucleoside hydrolase (IAG-NH) has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Special attention has been given to the role of leaving group pocket residues and conformation of the bound substrate at the active site of IAG-NH. We also describe the positioning of the residues of an important flexible loop at the active site, which was previously unobservable by X-ray crystallography due to high B-factors. Five MD simulations have been performed with the Enzyme x Substrate complexes: Enzyme x anti-Adenosine with Asp40-COOH [E(40H) x Ade(a)], Enzyme x anti-Adenosine with Asp40-COO- [E(40) x Ade(a)], Enzyme x syn-Adenosine with Asp40-COOH [E(40H) x Ade(s)], Enzyme x syn-Adenosine with Asp40-COO- [E(40) x Ade(s)], and Enzyme x anti-Inosine with Asp40-COO- [E(40) x Ino(a)]. Overall, the structures generated from the MD simulation of E(40H) x Ade(s) preserve the catalytically important hydrogen bonds as well as electrostatic and hydrophobic interactions to provide a plausible catalytic structure. When deprotonated Asp40 (Asp4-COO-) is present, the active site is open to water solvent which interferes with the base stacking between Trp83 and nucleobase. A calculation using Poisson-Boltzmann equation module supports that Asp40 indeed has an elevated pK(app). Solvent accessible surface area (SASA) calculations on all the five MD structures shows that systems with protonated Asp40, namely, E(40H) x Ade(a) and E(40H) x Ade(s), have zero SASA. It is found that a water molecule is hydrogen-bonded to the N7 of the nucleobase and is probably the essential general acid to protonate the departing nucleobase anion. The N7-bonded water is in turn hydrogen-bonded to waters in a channel, held in place by the residues of the flexible loop, Tyr257, His247, and Cys245. Using normal-mode analysis with elastic network model, we find that the flexible loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Population density analysis for determining the protonation state of the catalytic dyad in BACE1-tertiary carbinamine-based inhibitor complex

BACE1 is an aspartyl protease with a very relevant role in medicinal chemistry related to Alzheimer Disease since it has demonstrated to be a promising therapeutic target for inhibition and possible control for the progress of the peptide accumulation characteristic of this pathology. The enzymatic activity of this protein is given by the aspartic dyad, Asp93 and Asp289, which can adopt several protonation states depending on the chemical nature of its inhibitors, this is, monoprotonated, diprotonated and di-deprotonated states. In the present study, the analysis of the population density, for a series of protein-inhibitor molecular dynamics simulations, was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of tertiary carbinamine (TC) transition state analog inhibitors. The results revealed that the monoprotonated Asp289i state, in which the Asp93 and Asp289 residue side chains are deprotonated and protonated on the inner oxygen, respectively, is the most preferred in the presence of TC family inhibitors. This result was obtained after evaluating, for all 9 possible protonation state configurations, the individual and combined population densities of a set of parameters sensitive to protonation state of the Aspartic dyad, using an X-ray experimental BACE1/TC crystallographic structure as reference. This case study demonstrates again the usefulness of the concept of population density as a quantitative tool to establish the most stable system settings, among all possible, by measuring the level of occurrence of simultaneous events obtained from a sampling over time. These results will help to clear the phenomena related to the TCs inhibitory pathway, as well as assist in the design of better TC inhibitors against Alzheimer’s protease.     

Active-site mobility in human immunodeficiency virus, type 1, protease as demonstrated by crystal structure of A28S mutant.

The mutation Ala28 to serine in human immunodeficiency virus, type 1, (HIV-1) protease introduces putative hydrogen bonds to each active-site carboxyl group. These hydrogen bonds are ubiquitous in pepsin-like eukaryotic aspartic proteases. In order to understand the significance of this difference between HIV-1 protease and homologous, eukaryotic aspartic proteases, we solved the three-dimensional structure of A28S mutant HIV-1 protease in complex with a peptidic inhibitor U-89360E. The structure has been determined to 2.0 A resolution with an R factor of 0.194. Comparison of the mutant enzyme structure with that of the wild-type HIV-1 protease bound to the same inhibitor (Hong L, Treharne A, Hartsuck JA, Foundling S, Tang J, 1996, Biochemistry 35:10627-10633) revealed double occupancy for the Ser28 hydroxyl group, which forms a hydrogen bond either to one of the oxygen atoms of the active-site carboxyl or to the carbonyl oxygen of Asp30. We also observed marked changes in orientation of the Asp25 catalytic carboxyl groups, presumably caused by the new hydrogen bonds. These observations suggest that catalytic aspartyl groups of HIV-1 protease have significant conformational flexibility unseen in eukaryotic aspartic proteases. This difference may provide an explanation for some unique catalytic properties of HIV-1 protease.     

Mechanism of aspartyl proteinase action. VII. Noncovalent complexes of HIV-1 aspartyl proteinase with substrate and substrate-like inhibitors

A computer model of a noncovalent complex of HIV-1 aspartyl protease with substrate-like inhibitor JG-365 was a priori constructed by using the approaches of theoretical conformational analysis and molecular mechanics. The root mean square deviation of the calculated conformation of the inhibitor from the X-ray diffraction analysis data was 0.87 A. These results enabled the a priori calculation of the structure of noncovalent complex of HIV-1 protease with a hexapeptide fragment of its native specific substrate Ser-Gln-Asn-Tyr-Pro-Ile-Val. The only possible orientation of the cleavable peptide bond in this and the nucleophilic water molecule relative to the catalytically active Asp residues of the enzyme (Asp25 and Asp125) was found that provides for the chemical transformation of the substrate to a tetrahedral intermediate. An action mechanism of enzymes of this class was proposed on the basis of the analysis of calculated distances. We showed that neither steric distortion of the cleavable bond nor the formation of unfavorable contacts in molecules of the enzymes and their substrates accompany the optimum orientation of substrate molecules at the active sites of HIV-1 aspartyl proteases and rhizopuspepsin.     

Mechanism of action of aspartic proteinases. V. Conformational characteristics of fragments of substrate-binding sites in rhizopuspepsin and HIV-1 proteinase

The conformational states of side chains of catalytic Asp residues in active sites of HIV-1 protease and rhizopuspepsin in the potential field of free enzymes were studied by using theoretical conformational analysis. Structural factors that stabilize the conformation of these residues in free enzymes were revealed. Methods of molecular mechanics were used to estimate the stabilization energy of the Met46-Phe53 labile fragments of HIV-1 protease in the potential field of their nearest surrounding amino acid residues for the conformations characteristic of the free protein and similar to that of the protein in enzyme-inhibitor complexes. In solution, the conformational state of the fragments of the free enzyme was concluded to be similar to that observed in the enzyme complex with the ligand and different from that determined by X-ray diffraction analysis. This difference was ascribed to the effect of crystal packing.     

Molecular dynamics simulations of a branched tetradecasaccharide substrate in the active site of a xyloglucan endo-transglycosylase

Molecular dynamics simulations of the tetradecasaccharide XXXGXXXG in complex with the hybrid aspen xyloglucan endo-transglycosylase PttXET16-34 have been performed and analysed with respect to structure, dynamics, flexibility and ligand interactions. Notably, the charge state of the so-called ‘helper residue’ aspartate 87 (Asp87), which lies between the catalytic nucleophile [glutamate 85 (Glu85)] and general acid/base (Glu89) residues on the same beta strand, had a significant effect on PttXET16-34 active site structure. When Asp87 was deprotonated, electrostatic repulsion forced the nucleophile away from C1 of the sugar ring in subsite ? 1 and the proton–donating ability of Glu89 was also weakened due to the formation of a hydrogen bond with Asp87, whereas the protonation of Asp87 resulted in the formation of a hydrogen bond with the catalytic nucleophile and correct positioning of the catalytic machinery. The results suggest that catalysis in glycoside hydrolase family 16, and by extension clan GH-B enzymes, is optimal when the catalytic nucleophile is deprotonated for nucleophilic attack on the substrate, whereas the ‘helper residue’ and general acid/base residue are both in their conjugate acid forms to align the nucleophile and deliver a proton to the departing sugar, respectively.