Alexandra
Males
a,
Gaetano
Speciale
b,
Spencer J.
Williams
*b and
Gideon J.
Davies
*a
aYork Structural Biology Laboratory, Department of Chemistry, The University of York, YO10 5DD York, UK. E-mail: gideon.davies@york.ac.uk
bSchool of Chemistry and Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Melbourne, Victoria 3010, Australia. E-mail: sjwill@unimelb.edu.au
First published on 13th August 2019
Enzyme transition-state mimics can act as powerful inhibitors and allow structural studies that report on the conformation of the transition-state. Here, mannoimidazole, a mimic of the transition state of mannosidase catalyzed hydrolysis of mannosides, is shown to bind in a B2,5 conformation on the Clostridium perfringens GH125 α-1,6-mannosidase, providing additional evidence of a OS2–B2,5–1S5 conformational itinerary for enzymes of this family.
GHs have been grouped into >150 families on the basis of primary sequence and ability to cleave specific glycoside substrates.4,5 As sequence defines the three-dimensional structure of proteins, family members typically adopt similar folds and possess highly conserved amino acids within their active site clefts. The chemical mechanisms through which family members perform catalysis tend to be conserved across a family and lead to substrate hydrolysis with either retention or inversion of anomeric stereochemistry.6 Inversion of stereochemistry occurs in a single step reaction when water acts as a nucleophile. Retention of stereochemistry is typically achieved through a two-step substitution mechanism involving participation by an enzymic nucleophile or by a pendant neighboring group on the substrate.6–8 Both pathways benefit from enzymatic amino acids that provide general acid catalysis to promote the departure of the anomeric groups, and general base catalysis to promote nucleophilic attack by water. In almost all cases the transition states possess oxocarbenium-ion-like character and the stability of this TS is greatest when partial double bond character is maximized between C1 and O5, which occurs when the system C5–O5–C1–C2 is planar.9 Such transition-state conformations include two boats: B2,5 and 2,5B; two half-chairs: 4H3 and 3H4 (Fig. 1a); and the related envelopes: 3E, E3, 4E and E4. As transition state mimicry is a powerful strategy for inhibitor design, efforts have been made to assign conformational itineraries, and in particular the TS conformation, for a large number of GH families.1,2
Mannosidases and mannanases catalyze the hydrolysis of mannoside linkages and these enzymes are found within 11 GH families. Enzymes within families GH2, 5, 26, 38, 76, 92, 99, 113, 125 and 134 perform catalysis through itineraries along the OS2–B2,5–1S5 region of conformational space whereas enzymes of families GH47 and 134 react through the 3S1–3H4–1C4 conformational axis (Fig. 1b). The use of mannoimidazole and analogues thereof have been instrumental in the assignment of conformational itineraries and transition state conformation (Table 1).
In the case of family GH125, the conformational itinerary remained elusive for many years. Initial X-ray structures of the foundation member CpGH125 exo-α-1,6-mannosidase (CpGH125) with the nonhydrolyzable substrate mimic methyl 1,6-α-thiomannobioside (1) (Protein Databank (PDB) ID: 3QT9) and deoxymannojirimycin (2) (PDB ID: 3QRY) revealed undistorted, ground-state 4C1 conformations (Fig. 2a).18 Instead, a computational approach involving ab initio QM/MM metadynamics was used to model the conformation of substrate, 1,6-α-mannobiose, in the enzyme active site, by in silico substitution of the sulfur in the experimentally determined X-ray structure.19 This calculation predicted that the O-glycoside favoured an OS2 conformation on-enzyme and hinted at an OS2 → B2,5‡ → 1S5 conformational itinerary, which was fully supported by QM/MM simulations of the reaction mechanism (Fig. 2b). Direct experimental support for this mechanism was obtained through determination of the Michaelis complex for the catalytically-disabled acid mutant (D220N) of CpGH125 in complex with 1,6-α-mannobiose, which revealed an OS2 conformation consistent with the computational predictions.19
Fig. 2 (a) Structures of methyl α-thiomannobioside, deoxymannojirimycin and mannoimidazole. (b) Mechanism showing distortion of mannose from an OS2 conformation to 1S5, via a B2,5 transition state, a credible conformational itinerary of GH125 α-mannosidases. (c) syn and anti trajectory of the protonator/general acid residue in respect to the carbohydrate ligand. (d) Conformational free energy landscape of mannoimidazole in complex with a variety of GH family mannosidases. 1 + 2: β-Man-1,4-MIm bound to GH26 CjMan26C and GH113 AaManA;33: MIm bound to GH2 BtMan2A;144: MIm bound to GH47 CkMan47;105 + 6: MIm bound to GH38 DmGManII;3,207: MIm bound to GH92 BtMan3990;12 * (asterisk) represent MIm bound to CpGH125 (this work). |
As the ultimate determinant of the activation barrier for catalysis, the TS is the cardinal feature of the conformational itinerary. Mannoimidazole (3) (MIm) has emerged as a powerful chemical probe for TS conformation of mannosidases owing to its flattened nature, with an sp2-hybridized anomeric centre and energetically feasible access to all of the possible stereoelectronically-compliant conformations expected for an oxocarbenium ion-like transition state (Fig. 2b).3 Indeed, for a single step reaction, determination of the conformation of the transition state allows prediction of the conformational itinerary (from substrate to product via the TS) by invoking the principle of least nuclear motion.21 This principle states that elementary reactions that involve the least change in atomic position and electronic configuration will be favoured. Given the importance of the TS conformation in a reaction coordinate, and the utility of MIm in assigning conformational itineraries, we chose to revisit family GH125 to gather additional evidence in support of the proposed conformational itinerary. Here we report an X-ray structure of mannoimidazole (3) bound to wild-type CpGH125 that provides compelling evidence in support of a ‘latitudinal’ OS2 → B2,5‡ → 1S5 conformational itinerary, complementing previous computational and structural data.
We studied the thermodynamics for binding of MIm to CpGH125 using isothermal titration calorimetry (Fig. 3a). Analysis of the calorimogram using a one site binding model allowed prediction of binding at 3.2 sites with an overall KD value of 1.6 mM. Using a sequential binding site model, individual dissociation constants can be calculated at each site yielding values of 40 μM, 105 μM, and 3.1 mM. However, whilst this is indicative of multiple binding sites, it is not clear that the individual KD values are accurately determined by the data. In order to gain evidence that binding by MIm, as shown by ITC, results in enzyme inhibition, we determined its IC50 value. At a concentration of 0.5 mM DNPMan, we determined an IC50 of 1.7 mM (Fig. 3b), consistent with the overall ITC figure of 1.6 mM; this IC50 value cannot be converted to a KI value since the substrate does not give saturation kinetics and so KM cannot be determined.
Binding stoichiometry of >1 for basic sugar-shaped heterocycles is unusual and suggests that the ligand may not solely be binding at the −1 subsite. Previous work has identified three subsites, −1, +1 and +2, in CpGH125 through a complex with 1,6-α-mannotriose (PDB ID: 5M7Y),19 and studies of a disaccharide substrate that binds in the −1/+1 subsites revealed saturation kinetics, whereas a monosaccharide substrate, 2,4-dinitrophenyl α-mannoside (DNPMan), that bound only in the −1 subsite did not.18,26 Previous poor electron density in the +2 subsite suggests that this may represent the low affinity binding site, a conclusion supported by subsequent X-ray crystallographic studies (vide infra). The two tighter affinities therefore should likely reflect binding at the −1 and +1 subsites; however, we cannot draw conclusions as to which value is assigned to each of these subsites. Nevertheless, the interesting ITC data encouraged efforts to determine the 3D structure of a complex of CpGH125 with mannoimidazole.
Preformed crystals of CpGH125 recombinantly expressed in E. coli were soaked with MIm and the structure of the complex was solved at 1.85 Å resolution (Table 2). The crystal form contains 2 monomers in the asymmetric unit. Density for MIm was observed both in the −1 and +1 subsites, but not in the +2 subsite (Fig. 4a). The MIm molecule in the −1 subsite bound in a B2,5 conformation and the MIm molecule in the +1 subsite bound in a 4H3 conformation. The occupancy of MIm in the −1 subsite was estimated to be 80% and 70% for monomer A and B, respectively; the occupancy of MIm in the +1 subsite was 100% for both monomers. These observations are consistent with tight binding of MIm at the −1/+1 subsites, and weaker binding at the +2 subsite.
Fig. 4 (a) Complex of CpGH125 with mannoimidazole in the −1 and +1 subsites. The residues of the enzyme that form hydrogen bonding interactions and close contacts with mannoimidazole (<3 Å distance) are shown in blue and the nucleophilic water molecule is depicted by a red sphere. The electron density map is a REFMAC maximum-likelihood/σA-weighted 2 Fo − Fc synthesis contoured at 0.75 electrons Å−3 and 1.20 electrons Å−3 clipped to MIm in the −1 and +1 subsite, respectively. (b) Superposition of CpGH125 in complex with MIm (monomer A) and 1,6-α-mannobiose (PDB ID: 5M7I), shown in purple. Hydrogen bonds shown in black between the catalytic residues of CpGH125 and nucleophilic water. |
MIm molecules in each site engage with a range of active site residues (Fig. 4a). Key interactions with the −1 subsite MIm include the carbonyl oxygen of Ser219 with O2, Oδ1 of Asn302 with O3, Oδ1 and Oδ2 of Asp65 with O4 and O6, respectively, and NH2 of Arg405 with O6. The general base Glu393 is 2.3 Å away from a water poised for inline attack at the ‘anomeric carbon’ of MIm. Major interactions with MIm in subsite +1 include Oδ1 of Asp195 with O2, Glu129 Oε1 with O3 and Oε2 with O4, and the terminal amine nitrogen, NZ of Lys131 and NH2 of Arg188 with O3. The catalytic amino acid, Asp220, interacts with O6 of MIm in the +1 subsite. Trp116 is stacked parallel to the imidazole of MIm, with CZ2 of Trp119 and N2 of MIm spaced 3.79 Å apart.
Overlay of the complex of MIm with wildtype CpGH125 with that of the D220N mutant of the same enzyme with 1,6-α-mannobiose (PDB ID: 5M7I) is shown in Fig. 4b. This shows that the two MIm molecules bind such that their carbohydrate rings are well-matched with the equivalent mannose rings of 1,6-α-mannobiose. Moreover, the general acid (Asp220; or variant Asn220), and general base (Glu393) residues occupy similar positions in space. This overlay highlights that the interaction of Asp220 with the +1 subsite MIm molecule, in the complex determined here, can be considered equivalent to the interaction of Asn220 with the inter-glycosidic oxygen of 1,6-α-mannobiose seen in the CpGH125-D220N complex. Both complexes include a well-defined water molecule hydrogen-bonded to Glu393 and O3 of the −1 subsite MIm that is poised to undertake a nucleophilic attack at the anomeric centre. The interaction of this water molecule with O3 assists in orientating the water molecule to perform an in-line attack on C1 (Fig. 2b).
Nerinckx generalized the concept of syn/anti protonation to consider divergence from canonical orientation in the mean plane of the ring to include orientations in a syn or anti hemisphere (Fig. 2c).27 According to this extended concept, CpGH125 can be considered a syn protonator and with the present work CpGH125 joins a growing list examples of glycosidases that can bind glycoimidazoles even without achieving lateral protonation. While binding is seen in this and related cases, presumably the potency of binding is less than what might be expected for an ideally located general acid that can achieve such an interaction, an expectation that is borne-out by much tighter binding to canonical anti-protonators such as Bt3130 and Bt3965 from Bacteroides thetaiotaomicron GH92 (Ki = 1.0 and 0.4 μM, respectively) whereas other α-mannosidases from B. thetaiotaomicron have a Ki values ranging from 40–6000 μM.12 Also, imidazole-substituted mannoimidazoles can exhibit even tighter binding, for example, the β-mannosidase BtMan2A (from family GH2) binds phenethylmannoimidazole with Ki = 57 nM compared to binding to mannoimidazole with Ki = 1400 nM.14 In this case since the inhibitor lies across both the −1 and +1 subsite, the number of interactions the inhibitor can make with the protein increases. We highlight one important anomalous case, in which a complex of a disaccharide MIm analogue with BxGH99 1,2-α-mannanase revealed an unusual 2H3/E3 conformation possibly due to steric clashes with an active site residue, Tyr252.25 In this case no KD value could be measured by ITC indicating poor affinity for the enzyme.
Data were collected on Diamond MX beamline I24 and processed using xia2 3dii.28 Integrated data were processed using the CCP4 software suite;29 the programs used were aimless,30,31 MOLREP32 with 5M7I as the reference model and cycles of REFMAC533–38 and COOT39 for refinement of the structure. CpGH125 structures figures were produced using CCP4 mg.40 Coordinate and observed structure-factor amplitudes have been deposited on the PDB with accession code 6RQK.
This journal is © The Royal Society of Chemistry 2019 |