Using polarizable molecular mechanics, a recent study [de Courcy et al. J. Am. Chem. Soc., 2010, 132, 3312] has compared the relative energy balances of five competing inhibitors of the FAK kinase. It showed that the inclusion of structural water molecules was indispensable for an ordering consistent with the experimental one. This approach is now extended to compare the binding affinities of four active site ligands to the Type I Zn-metalloenzyme phosphomannose isomerase (PMI) from Candida albicans. The first three ones are the PMI substrate β-D-mannopyranose 6-phosphate (β-M6P) and two isomers, α-D-mannopyranose 6-phosphate (α-M6P) and β-D-glucopyranose 6-phosphate (β-G6P). They have a dianionic 6-phosphate substituent and differ by the relative configuration of the two carbon atoms C1 and C2 of the pyranose ring. The fourth ligand, namely 6-deoxy-6-dicarboxymethyl-β-D-mannopyranose (β-6DCM), is a substrate analogue that has the β-M6P phosphate replaced by the nonhydrolyzable phosphate surrogate malonate. In the energy-minimized structures of all four complexes, one of the ligand hydroxyl groups binds Zn(II) through a water molecule, and the dianionic moiety binds simultaneously to Arg304 and Lys310 at the entrance of the cavity. Comparative energy-balances were performed in which solvation of the complexes and desolvation of PMI and of the ligands are computed using the Langlet-Claverie continuum reaction field procedure. They resulted into a more favorable balance in favor of β-M6P than α-M6P and β-G6P, consistent with the experimental results that show β-M6P to act as a PMI substrate, while α-M6P and β-G6P are inactive or at best weak inhibitors. However, these energy balances indicated the malonate ligand β-6DCM to have a much lesser favorable relative complexation energy than the substrate β-M6P, while it has an experimental 10-fold higher affinity than it on Type I PMI from Saccharomyces cerevisiae. The energy calculations were validated by comparison with parallel ab initio quantum chemistry on model binding sites extracted from the energy-minimized PMI-inhibitor complexes. We sought to improve the models upon including explicit water molecules solvating the dianionic moieties in their ionic bonds with the Arg304 and Lys310 side-chains. Energy-minimization resulted in the formation of three networks of structured waters. The first water of each network binds to one of the three accessible anionic oxygens. The networks extend to PMI residues (Asp17, Glu48, Asp300) remote from the ligand binding site. The final comparative energy balances also took into account ligand desolvation in a box of 64 waters. They now resulted into a large preference in favor of β-6DCM over β-M6P. The means to further augment the present model upon including entropy effects and sampling were discussed. Nevertheless a clear-cut conclusion emerging from this as well as our previous study on FAK kinase is that both polarization and charge-transfer contributions are critical elements of the energy balances.

Type I phosphomannose isomerases (PMIs) are zinc-dependent metalloenzymes involved in the reversible isomerization of D-mannose 6-phosphate (M6P) and D-fructose 6-phosphate (F6P). 5-Phospho-D-arabinonohydroxamic acid (5PAH), an inhibitor endowed with nanomolar affinity for yeast (Type I) and Pseudomonas aeruginosa (Type II) PMIs (Roux et al., Biochemistry 2004; 43:2926-2934), strongly inhibits human (Type I) PMI (for which we report an improved expression and purification procedure), as well as Escherichia coli (Type I) PMI. Its K(i) value of 41 nM for human PMI is the lowest value ever reported for an inhibitor of PMI. 5-Phospho-D-arabinonhydrazide, a neutral analogue of the reaction intermediate 1,2-cis-enediol, is about 15 times less efficient at inhibiting both enzymes, in accord with the anionic nature of the postulated high-energy reaction intermediate. Using the polarizable molecular mechanics, sum of interactions between fragments ab initio computed (SIBFA) procedure, computed structures of the complexes between Candida albicans (Type I) PMI and the cyclic substrate β-D-mannopyranose 6-phosphate (β-M6P) and between the enzyme and the high-energy intermediate analogue inhibitor 5PAH are reported. Their analysis allows us to identify clearly the nature of each individual active site amino acid and to formulate a hypothesis for the overall mechanism of the reaction catalyzed by Type I PMIs, that is, the ring-opening and isomerization steps, respectively. Following enzyme-catalyzed ring-opening of β-M6P by zinc-coordinated water and Gln111 ligands, Lys136 is identified as the probable catalytic base involved in proton transfer between the two carbon atoms C1 and C2 of the substrate D-mannose 6-phosphate.