# Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase

Joan Giménez-Dejoz, Michal H. Kolář, Francesc Xavier Ruiz, Isidro Crespo, Alexandra Cousido-Siah, Alberto Podjarny, Oleg A. Barski, Jindřich Fanfrlík, Xavier Parés, Jaume Farrés, et al.

► To cite this version:

Joan Giménez-Dejoz, Michal H. Kolář, Francesc Xavier Ruiz, Isidro Crespo, Alexandra Cousido-Siah, et al.. Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase. PLoS ONE, 2015, 10 (7), pp.e0134506. 10.1371/journal.pone.0134506. [hal-04977179](https://hal.science/hal-04977179v1)

HAL Id: hal-04977179

https://hal.science/hal-04977179v1

Submitted on 5 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase

Joan Giménez-Dejoz¹, Michal H. Kolář²,³, Francesc X. Ruiz⁴, Isidro Crespo¹,
Alexandra Cousido-Siah⁴, Alberto Podjarny⁴, Oleg A. Barski⁵, Jindřich Fanfrlík²,
Xavier Parés¹, Jaume Farrés¹*, Sergio Porté¹*

1 Department of Biochemistry and Molecular Biology, Faculty of Biosciences, Universitat Autònoma de Barcelona, Bellaterra, Barcelona, Spain, 2 Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic, 3 Institute of Neuroscience and Medicine (INM-9) and Institute for Advanced Simulation (IAS-5), Forschungszentrum Jülich GmbH, Jülich, Germany, 4 Department of Integrative Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire-Centre de Biologie Intégrative, CNRS, INSERM, UdS, Illkirch CEDEX, France, 5 Diabetes and Obesity Center, School of Medicine, University of Louisville, Louisville, Kentucky, United States of America

* jaume.farres@uab.cat (JF); sergio.porte@gmail.com (SP)

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Giménez-Dejoz J, Kolář MH, Ruiz FX, Crespo I, Cousido-Siah A, Podjarny A, et al. (2015) Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase. PLoS ONE 10(7): e0134506.
doi:10.1371/journal.pone.0134506

**Editor:** Fernando Rodrigues-Lima, University Paris Diderot-Paris 7, FRANCE

**Received:** March 31, 2015

**Accepted:** July 9, 2015

**Published:** July 29, 2015

**Copyright:** © 2015 Giménez-Dejoz et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** JGD, IC, XP, JF and SP received support from Spanish Ministerio de Economía y Competitividad [BFU2011-24176] (http://www.mineco.gob.es) and Generalitat de Catalunya [2009 SGR 795] (http://www.gencat.cat/agaur). FXR, ACS and AP received support from Centre National de la Récherche Scientifique (http://www.cnrs.fr); Institut National de la Santé et de la Recherche Médicale (http://www.inserm.fr); Université de Strasbourg

## Abstract

Human aldo-keto reductase 1B15 (AKR1B15) is a newly discovered enzyme which shares 92% amino acid sequence identity with AKR1B10. While AKR1B10 is a well characterized enzyme with high retinaldehyde reductase activity, involved in the development of several cancer types, the enzymatic activity and physiological role of AKR1B15 are still poorly known. Here, the purified recombinant enzyme has been subjected to substrate specificity characterization, kinetic analysis and inhibitor screening, combined with structural modeling. AKR1B15 is active towards a variety of carbonyl substrates, including retinoids, with lower $k_{cat}$ and $K_m$ values than AKR1B10. In contrast to AKR1B10, which strongly prefers all-trans-retinaldehyde, AKR1B15 exhibits superior catalytic efficiency with 9-cis-retinaldehyde, the best substrate found for this enzyme. With ketone and dicarbonyl substrates, AKR1B15 also shows higher catalytic activity than AKR1B10. Several typical AKR inhibitors do not significantly affect AKR1B15 activity. Amino acid substitutions clustered in loops A and C result in a smaller, more hydrophobic and more rigid active site in AKR1B15 compared with the AKR1B10 pocket, consistent with distinct substrate specificity and narrower inhibitor selectivity for AKR1B15.

## Introduction

Members of the aldo-keto reductase (AKR) superfamily share an $(\alpha/\beta)_8$ barrel fold and are mostly monomeric (~35 kDa) NAD(P)H-dependent enzymes catalyzing the reduction of a
(https://www.unistra.fr); Région Alsace (http://www.region.alsace); Hôpital Civil de Strasbourg (http://www.chru-strasbourg.fr/Hopital-civi); Instruct (part of the European Strategy Forum of Research Infrastructures; ESFRI) (https://www.structuralbiology.eu); and French Infrastructure for Integrated Structural Biology (FRISBI) [ANR-10-INSB-05-01] (http://frisbi.eu). FXR received support from Fondation pour la Recherche Médicale [SPF20121226275] (http://www.frm.org). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

wide variety of endogenous carbonyl compounds such as carbohydrates, lipid aldehydes, prostaglandins, steroids and retinoids [1,2]. Several AKR enzymes act in phase-I drug metabolism by transforming some xenobiotics, are induced by nuclear factor erythroid 2-related factor 2 (Nrf2) under oxidative stress, and are involved in cancer chemoresistance [3,4]. In the human genome, fifteen AKR genes have been described which belong to three different gene families (AKR1, AKR6 and AKR7). The AKR1B subfamily gene cluster, located in chromosome 7q33-35, includes the AKR1B1, AKR1B10 and AKR1B15 genes. A syntenic gene cluster with four *loci* has been described in rodent species, although gene orthologs can only be unambiguously assigned for AKR1B1 [5,6].

The most studied enzyme of the AKR1B subfamily is AKR1B1 or aldose reductase, which reduces glucose to sorbitol under hyperglycemia and has been involved in the secondary complications of diabetic disease [7]. Another member, AKR1B10, is normally expressed in adrenal gland and small intestine, and induced in several types of cancer, such as non-small cell lung carcinoma and hepatoma [3]. Both enzymes have been proposed as promising oncogenic targets [8,9] and for this reason, along with the role of AKR1B1 in diabetic disease, they have been the subject of many studies in the search of selective and potent inhibitors [10–15]. Unlike other members of the subfamily, AKR1B10 is highly active in the reduction of all-*trans*-retinal-dehyde [5]. The third gene in the AKR1B cluster, AKR1B15, has been predicted in the last decade as a result of high-throughput sequencing and annotation projects, i.e. human genome. Recently, AKR1B15 has been demonstrated to be a functional gene with low expression restricted to placenta, testes and adipose tissues. The AKR1B15 gene undergoes alternative splicing giving rise to two protein isoforms, designated as AKR1B15.1 and AKR1B15.2. The former is a 316-amino acid protein encoded by AKR1B15-201 mRNA (Ensembl database) and showing 92% amino acid sequence identity with AKR1B10, whereas AKR1B15.2 (AKR1B15-001, Ensembl) has a longer N-terminus not homologous to other AKRs, and does not exhibit enzymatic activity or nucleotide binding [16]. AKR1B15.1 (henceforth referred to in this manuscript as AKR1B15 for brevity) is localized in the mitochondrial fraction and the recombinant protein was purified and characterized, showing limited *in vitro* activity with steroids and acetoacetyl-CoA [16]. Previously, AKR1B15.1 had been expressed in the insoluble fraction of mammalian cells, showing low activity with D,L-glyceraldehyde and 4-nitrobenzaldehyde [6]. Similarly to AKR1B10, the AKR1B15 gene was found to be up-regulated in the airway epithelium by smoking [17] and by exposure to sulforaphane, a known activator of the antioxidant response [18]. Interest in the AKR1B15 gene has risen lately because some allelic variants have been linked to a mitochondrial oxidative phosphorylation disease [19], serous ovarian carcinoma [20] and increased longevity [21].

With the aim of further characterizing the enzymatic function of AKR1B15, we have performed enzyme kinetics of the purified recombinant protein with retinaldehyde isomers and other typical carbonyl substrates of AKR1B10. We have also conducted a screening against potential inhibitors using compounds previously described for AKR1B1 or AKR1B10. Finally, based on the crystallographic structure of the AKR1B10 complex with NADP+ and tolrestat, we have constructed a model of the AKR1B15 active-site pocket.

# Materials and Methods

## Bacterial strains, plasmids and reagents

E. coli BL21(DE3) strain was obtained from Novagen, while plasmids pBB540 and pBB542 (containing the chaperone-coding genes grpE, clpB and dnaK, dnaJ, groESL, respectively), were a kind gift from Dr. A. de Marco [22]. The pET-28a vector containing the cDNA coding for isoform 2 of AKR1B15 (UniProt ID: C9JRZ8-2) had been described by Salabei et al. [6].
Tolrestat and sorbinil were generously provided by Prof. T.G. Flynn and Pfizer, respectively, whereas JF0064 (2,2',3,3',5,5',6,6'-octafluoro-4,4'-biphenyldiol) was obtained from Sigma-Aldrich. *Trans*-2-hexenal and 4-hydroxy-2-nonenal were commercially obtained from Cayman Chemical. All other reagents, including substrates, were purchased from Sigma-Aldrich unless otherwise indicated.

## Protein expression and purification

*E. coli* BL21(DE3) strain transformed with pET-28a/AKR1B15 was grown in 1 L of 2xYT medium in the presence of 33 µg/mL kanamycin, while *E. coli* BL21(DE3) containing pBB540, pBB542 and pET-28a/AKR1B15 was grown in 6 L of M9 minimal medium supplemented with 20% glucose as a carbon source, in the presence of 34 µg/mL chloramphenicol, 50 µg/mL spectinomycin and 33 µg/mL kanamycin. Protein expression was then induced by the addition of 1 mM IPTG (Apollo Scientific) and cells were further incubated for 4 h at 22°C. Cells were then pelleted and resuspended in ice-cold TBI buffer (150 mM NaCl, 10 mM Tris-HCl, 5 mM imidazole, pH 8.0) containing 1% (v/v) Triton X-100. In the case of the non-chaperone-expressing *E. coli* BL21(DE3) strain, the TBI buffer also contained 1% (w/v) sarkosyl. The protein was purified using a His-Trap HP nickel-charged chelating Sepharose Fast Flow (GE Healthcare) 5-mL column using an AKTA FPLC purification system. The column was washed with TBI buffer and the enzyme was eluted stepwise with 5, 60, 100 and 500 mM imidazole in TBI buffer. The enzyme fraction eluted with 100 mM imidazole was loaded onto a PD-10 column (Millipore), which removed imidazole and changed the buffer to storage buffer (200 mM potassium phosphate, pH 7.4, 5 mM EDTA, 5 mM DTT). Finally, the protein monomer was purified through gel filtration chromatography using a Superdex 75 10/300 GL column (GE Healthcare) equilibrated with the storage buffer. In the case of the protein expressed in the *E. coli* BL21 (DE3) strain, in the absence of chaperones, the TBI and storage buffers contained 0.1% (w/v) sarkosyl throughout the purification procedure. AKR1B10 and AKR1B1 were expressed and purified as described previously [23].

## Fluorimetric and spectrophotometric assays

NADPH binding was analyzed by quenching of Trp intrinsic fluorescence of 0.5 µM protein, using a Cary Eclipse (Varian) fluorimeter, in 20 mM sodium phosphate, pH 7.0, at 25°C in a final volume of 1 mL. The excitation wavelength was 290 nm and the emission wavelength was monitored at 340 nm. AKR1B10 was used as a control. The dissociation constant ($K_D$) values were calculated by nonlinear fitting of experimental data, using Grafit 5.0 (Eritacus Software), to the Morrison equation:

$$ \Delta F = \Delta F_{max} \cdot \frac{[E] + [NADPH] + K_D - \sqrt{([E] + [NADPH] + K_D)^2 - 4 \cdot [E] \cdot [NADPH]}}{2 \cdot [E]} $$

The activity towards aldehydes, with the exception of retinaldehydes, was analyzed spectrophotometrically following the decrease in the absorbance of the cofactor NADPH at 340 nm ($\epsilon_{340} = 6,220 \text{ M}^{-1}\cdot\text{cm}^{-1}$) or at 365 nm in the case of cinnamaldehyde ($\epsilon_{365} = 3,510 \text{ M}^{-1}\cdot\text{cm}^{-1}$) [24]. Activities were determined in 100 mM sodium phosphate, pH 7.0, at 25°C, using 0.2 mM NADPH in 0.2-cm path length cuvettes, with freshly prepared substrate solutions. One unit of activity is defined as the amount of enzyme required to transform 1 µmol of substrate per min at 25°C.
## HPLC enzymatic activity assay

Activity assays with retinoids were carried out using an HPLC-based methodology [25]. Briefly, retinaldehyde isomers were solubilized using glass tubes by a 10-min sonication at molar ratio 1:1 with fatty acid-free bovine serum albumin in 90 mM potassium phosphate, 40 mM potassium chloride, pH 7.4. The actual amount of solubilized retinoid was determined based on the corresponding molar absorption coefficient in aqueous solutions at the appropriate wavelength: ε<sub>400</sub> = 29,500 M<sup>-1</sup>·cm<sup>-1</sup> for all-*trans*-retinaldehyde and ε<sub>367</sub> = 26,700 M<sup>-1</sup>·cm<sup>-1</sup> for 9-*cis*-retinaldehyde [25]. For retinol isomers, which were used as standards of the reaction product, their concentration was determined in hexane using ε<sub>325</sub> = 51,770 M<sup>-1</sup>·cm<sup>-1</sup> for all-*trans*-retinol [26] and ε<sub>325</sub> = 43,765 M<sup>-1</sup>·cm<sup>-1</sup> for 9-*cis*-retinol [27]. The reactions were started by the addition of cofactor and carried out for 15 min at 37°C in a final volume of 0.5 mL. With the aim to measure the steady-state enzymatic activity, the concentration of enzyme was kept from 25- to 100-fold lower than that of the substrate for all the enzymatic assays. The reactions were stopped by the addition of 1 mL of cold methanol and after two rounds of extraction with hexane, retinoids were analyzed by HPLC as previously described [25]. All retinoid manipulations were performed under dim light.

## Determination of kinetic constants and inhibition screening

All compounds tested as inhibitors were dissolved in DMSO and assayed in a final concentration of 0.1% (v/v) DMSO using 6 mM D<sub>L</sub>-glyceraldehyde as a substrate. Kinetic constant and IC<sub>50</sub> (compound concentration that inhibits enzymatic activity by 50%) values were calculated by fitting the initial rates to the appropriate equation using Grafit 5.0 (Eritacus Software) and values were given as the mean ± standard error of three experiments. Standard error values were less than 20% of the mean values.

## Homology model and conformational ensembles

The structural model of *apo*-AKR1B15 was obtained from the AKR1B10 ternary complex crystallographic structure (PDB ID: 1ZUA), used as a template for homology modeling, by running the SCWRL program [28]. Because of the high sequence identity (92%) between the two proteins, the approach that keeps invariant the conformation of the conserved residues was adopted. The flexibility features of both AKR1B10 and AKR1B15 were studied by means of computer simulations.

The hydrogen atoms were added to the structural model of AKR1B15 as well as to the AKR1B10 crystal by the pdb2gmx tool of the GROMACS program package [29]. The *apo* as well as *holo* forms of both structures were studied as follows: the all-atom models were energy minimized employing the Amber99sb-ildn force field [30,31] for the protein, parameters of Holmberg *et al.* [32] for the cofactor, and Generalized-Born implicit solvent model with parameters of Hawkins *et al.* [33]. The minimized structures were used as the input for the tCONCOORD algorithm [34], which generates a set of independent conformations based on geometrical constraints [35]. It was designed to accurately capture the protein flexibility and the validity of the resulting conformational ensembles has been proven on a variety of proteins, including also AKR1B1 [34]. By means of tCONCOORD, we generated and analyzed the ensembles of 2500 conformations, which were subsequently used to calculate root mean square fluctuations (RMSF) of backbone atoms.

The energy minimized geometries (i.e. the starting ones for conformational sampling) were evaluated employing PROCHECK [36] and Verify 3D [37], which allow checking their stereochemical quality and calculating the percentage of conformations in favored regions obtained from the Ramachandran plots.
## AKR1B15 ternary complexes

Geometries of the ternary complexes were derived from the homology model of *holo*-AKR1B15, which was built by superimposition of the *apo* form and *holo*-AKR1B10 (PDB ID: 1ZUA), including the cofactor coordinates into *apo*-AKR1B15. Considering that all-*trans*-retinaldehyde and 9-*cis*-retinaldehyde are substrates, the distance between the oxygen of retinaldehyde and the catalytic residues (i.e. His111 and Tyr49) should be less than 3 Å for productive catalysis to occur. Under this assumption, the substrates were manually docked into the active site of *holo*-AKR1B15. Inhibitors were automatically docked using AutoDock 4.0 [38]. The inhibitor JF0064 coordinates were obtained from the PDB (PDB ID 4ICC). The target geometry was extracted from the energy minimized all-*trans*-retinaldehyde complex (see below). For docking, the target was kept rigid, while all the torsional bonds in JF0064, except for the conjugated double bonds, were free to move. The docking parameters were the same as described previously [39].

The hydrogen atoms were added to the ternary complex and the complex was energy minimized by adopting the PM6-D3H4 method [40] combined with the COSMO solvent model [41]. The residues, with at least one atom within 5 Å from any atom of the cofactor, substrate or inhibitor, were allowed to move and the rest was kept frozen but included in the Hamiltonian calculation. The PM6-D3H4 method, which has been developed to accurately describe non-covalent interactions in biomolecules, represents a well-established computational tool [42] and recently it has been used in the study of the inhibition of AKR1B1 [43]. The final pdb files of *holo*-AKR1B15 and the ternary complexes were validated with the QMEAN server [44], with the flag "ignoring the agreement terms", recommended for proteins known to have the correct fold. The volume of the active-site pocket was measured by using the POVME algorithm [45], whereas PyMOL (v.1.3; Schrödinger) was used for figure drawing.

## Results and Discussion

### Expression and purification of recombinant human AKR1B15

We attempted the expression and purification of recombinant AKR1B15 using different *E. coli* strains and procedures. In all cases, the protein appeared to be mostly present in the insoluble fraction of cell lysates. Previously, Salabei *et al*. [6] had described a successful procedure based on the use of the anionic detergent sarkosyl (sodium N-lauroylsarcosinate). In our hands and in the presence of sarkosyl, the amount of soluble protein and the final yield were acceptable (Fig 1A and 1B). However, AKR1B15 was mainly associated with higher molecular weight aggregates (Fig 1C), which were found to be enzymatically inactive. Only a small fraction of AKR1B15 was found to be in an active monomeric form. A minimal concentration of sarkosyl (0.001%, w/v) was necessary to avoid protein precipitation. Sarkosyl had been described as an enzyme inhibitor [46], and was confirmed that it inhibited AKR1B15 following a non-competitive model (data not shown). As an alternative procedure for expression, an *E. coli* BL21(DE3) strain co-expressing three chaperone systems (DnaK-DnaJ-GrpE, ClpB and GroEL-GroES) [22] was used. Under these conditions, AKR1B15 was still found to be mainly associated with the insoluble fraction, but the amount of soluble protein recovered increased significantly (Fig 1A). Peptide mass fingerprint analysis of the soluble protein fraction separated by SDS-PAGE identified the isolated band as being human AKR1B15 (data not shown). The enzyme was purified with a final yield of 1 mg protein per liter of culture, when using minimal medium. Analysis by gel filtration chromatography showed a molecular weight of 37 kDa for the purified protein, suggesting that AKR1B15 was obtained as a monomer (Fig 1C). By analyzing the purified protein on SDS-PAGE (Fig 1B, lane 2), the major band corresponded to the 37-kDa
A

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Expression and purification of recombinant human AKR1B15. (A) SDS-PAGE analysis of protein expression, showing that AKR1B15 was predominantly associated with the insoluble fraction of BL21(DE3) cell lysates. Treatment with 1% (w/v) sarkosyl (Sk) provided a much higher amount of AKR1B15 in the soluble fraction. In the case of BL21(DE3) pBB540 pBB542 cells, a protein band which is highlighted with a red oval was identified as human AKR1B15 by Peptide Mass Fingerprinting. Lanes: C, control for the soluble fraction not induced by IPTG; S, soluble fraction; and I, insoluble fraction. (B) SDS-PAGE analysis of protein purification, showing fractions eluted from the nickel affinity column chromatography using 100 mM imidazole. Lanes: 1, Protein eluted from the soluble fraction of BL21(DE3) + Sk; and 2, Protein eluted from the soluble fraction of BL21(DE3) pBB540 pBB542. (C) Elution profile from the Superdex 75 10/300 GL column chromatography. AKR1B15 purified from soluble fraction of BL21(DE3) + Sk and from soluble fraction of BL21(DE3) pBB540 pBB542 are shown in grey and black lines, respectively. Major peaks eluting at 7.9 and 11.4 mL correspond to aggregated (132 kDa) and monomer (37 kDa) protein, respectively.</figcaption></figure>

doi:10.1371/journal.pone.0134506.g001
enzyme, with some minor contaminating protein bands. Fluorescence analysis of cofactor
binding allowed us to determine the K<sub>D</sub> value for NADPH (113 ± 9 nM, Fig 2). The determined
K<sub>D</sub> value was in the same range with that of AKR1B10 (K<sub>D</sub> = 92 ± 35 nM).

Enzymatic activity and inhibition studies

It has been recently reported that AKR1B15 catalyzes the reduction of acetoacetyl-CoA and the carbonyl group at C17 position of sex steroids [16]. Additional results, obtained by using protein extracts from transfected human COS-7 cells, indicated that AKR1B15 exhibited very weak activity towards D₄,L-glyceraldehyde and 4-nitrobenzaldehyde [6], which are typical substrates of AKR enzymes. Our current analysis shows that the purified recombinant enzyme has broad substrate specificity and displays a significant enzymatic activity towards aliphatic and aromatic aldehydes and ketones (Table 1). AKR1B15 was active with D₄,L-glyceraldehyde, which was used as a substrate in the standard assay. The enzyme reduced 200 mM glucose with a much lower rate (12 mU/mg) than that exhibited by AKR1B1 (500 mU/mg), precluding the determination of kinetic constants. In contrast, AKR1B15 displayed a catalytic activity comparable or higher than that of AKR1B1 and AKR1B10 with a variety of physiological or model aldehydes and ketones of various classes. Among the compounds assayed, medium-chain (i.e. ≥ 6-carbon) aliphatic and aromatic carbonyl compounds were excellent substrates, with K<sub>m</sub> values in the low micromolar range. Importantly, AKR1B15 was also active towards retinaldehyde isomers (Table 2), which along with acrolein, trans-2-hexenal, 4-hydroxy-2-nonenal and farnesal may constitute physiological substrates for this enzyme. The 9-cis isomer of retinaldehyde was the best substrate based on the catalytic efficiency (k<sub>cat</sub>/K<sub>m</sub>) due to a very low K<sub>m</sub> value (160 nM) (Table 2), which could be determined by the use of an HPLC-based method providing higher sensitivity than the spectrophotometric assay. Regarding cofactor specificity, AKR1B15 was confirmed to be NADPH dependent (i.e. the enzymatic activity using 0.2 mM NADH was less than 5% of that with NADPH), as reported in [16], and described for most AKRs. The K<sub>m</sub> value of AKR1B15 for NADPH (5 µM) is within the same range as those of AKR1B1 and AKR1B10 (Table 1).

The comparison between AKR1B15 and the other human AKR1B enzymes reveals some relevant distinct kinetic features (Table 1). In particular, ketones and α-dicarbonyl compounds were good substrates for AKR1B15, showing higher activity and lower K<sub>m</sub> values than AKR1B10 and AKR1B1. The β-dicarbonyl compound, 2,4-pentanedione, was also reduced by AKR1B15, consistent with the reported activity with acetoacetyl-CoA, which also has two carbonyl groups in a β-position [16]. For other substrates, AKR1B15 exhibited lower K<sub>m</sub> values than AKR1B10, suggesting that AKR1B15 could play a role in aldehyde detoxification, similar to what has been suggested for other AKR1B enzymes [47]. AKR1B15 resembles AKR1B10 in having high activity with retinoids, in contrast to AKR1B1 (Table 2). A distinct feature of AKR1B15 is that it keeps similar k<sub>cat</sub> values towards all assayed substrates, including retinoids, suggesting a common rate-limiting step. In comparison, the k<sub>cat</sub> values of AKR1B1 with both retinaldehyde isomers and AKR1B10 with the 9-cis isomer are significantly lower than those for other substrates. This decrease in k<sub>cat</sub> values had been interpreted before as a change in the rate-limiting step (from cofactor dissociation to a slower step) of AKR1B1 and AKR1B10 reactions with these retinoids [23]. This seems not to be the case for AKR1B15. Finally, the higher specificity of AKR1B15 for the 9-cis isomer is unique for the human AKR1Bs, and it matches those of AKR1C3 [48], which like AKR1B15 also has bulky Phe residues at positions 299 and 304 in its active site, and chicken AKR [49].

Seven AKR inhibitors were tested against AKR1B15 using D₄,L-glyceraldehyde as substrate (Table 3). Among these inhibitors, there are five of the carboxylic acid type (tolrestat,
<figure><img src="image_5.png" /><figcaption>**Fig 2. Quenching of AKR1B15 and AKR1B10 fluorescence upon binding of NADPH.** Change of the protein fluorescence intensity (in percentage) upon addition of cofactor is shown. All proteins were used at a concentration of 0.5 µM in 20 mM sodium phosphate, pH 7.0, at 25°C. Graph symbols: AKR1B15 (diamonds), AKR1B10 (open squares).</figcaption></figure>

doi:10.1371/journal.pone.0134506.g002

epalrestat, oleanolic acid, sulindac and lithocholic acid), one of cyclic imide type (sorbinil) and one non-classical aldose reductase inhibitor or ARI (JF0064, which has been recently described by Cousido-Siah *et al* [15], Fig 3). Tolrestat, JF0064 and sulindac are potent inhibitors of AKR1B1 as well as of AKR1B10 (with IC₅₀<10 µM for both enzymes). Sorbinil and epalrestat are more selective against AKR1B1 [50,51], while oleanolic and lithocholic acid are more selective against AKR1B10 [52]. For AKR1B15, only JF0064 showed a significant inhibition (IC₅₀ = 0.034 ± 0.005 µM, Table 3), much stronger than for AKR1B1 and AKR1B10 [15]. The steroid lithocholic acid was found to be an inhibitor of AKR1B15, in agreement with the reported observation that certain steroids are good substrates for this enzyme [16].

Structural model of AKR1B15

With the aim to compare the substrate-binding site of AKR1B15 with that of AKR1B10, and upon unsuccessful attempts of protein crystallization, a 3D homology model of AKR1B15 was built (Fig 4A). The SCWRL server was chosen, as it is designed specifically for predicting side-chain conformations, provided a fixed backbone usually obtained from an experimental structure. In practical terms, this is the case of AKR1B15, given its 92% sequence identity with AKR1B10. Model quality was checked by PROCHECK analyses, indicating that most residues are in the preferred regions (94.8%), whereas residues in allowed regions and outliers are 3.5 and 1.6%, respectively. The compatibility of the atomic model (3D) was checked with its own amino acid sequence (1D) using Verify 3D analysis. This analysis shows that there is no region with negative scores, which would otherwise indicate potential problems. In addition, the analysis using the QMEAN server indicates the reliability of the model, with a QMEAN score of 0.67 out of 1, and with all the Z-scores being consistent with the good quality of the structure. Thus it may be concluded that this model is suitable for structural studies.

Model analysis

Regarding the cofactor-binding site, AKR1B15 shares with AKR1B10 all residues, except for Arg22, Met265 and His269 (Fig 5). The change of Lys22 in AKR1B10 to Arg in AKR1B15 may prevent its interaction with the pyrophosphate bridge of NADP⁺ but, due to the mobility of the Arg side chain in solution, its interaction with the cofactor cannot be excluded. Met265 and His269 keep interactions with the same groups of NADP⁺ that involve Val265 and Arg269 in
<table><caption>Table 1. Kinetic constants of AKR1B15, AKR1B10 and AKR1B1 with aldehydes and ketones.</caption><thead><tr><th rowspan="2">Substrate</th><th colspan="3">AKR1B15</th><th colspan="3">AKR1B10</th><th colspan="3">AKR1B1</th></tr><tr><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(min<sup>-1</sup>)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(mM<sup>-1</sup>.min<sup>-1</sup>)</th><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(min<sup>-1</sup>)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(mM<sup>-1</sup>.min<sup>-1</sup>)</th><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(min<sup>-1</sup>)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(mM<sup>-1</sup>.min<sup>-1</sup>)</th></tr></thead><tbody><tr><td><strong>Carbohydrate aldehydes</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>D<sub>1</sub>L-glyceraldehyde</td><td>880</td><td>10.7</td><td>12.4</td><td>6,000<sup>a</sup></td><td>35<sup>a</sup></td><td>6<sup>a</sup></td><td>50<sup>a</sup></td><td>31<sup>a</sup></td><td>660<sup>a</sup></td></tr><tr><td><strong>Aromatic aldehydes</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>pyridine-3-aldehyde</td><td>2.9</td><td>9</td><td>3,150</td><td>13<sup>b</sup></td><td>150<sup>b</sup></td><td>12,000<sup>b</sup></td><td>10<sup>b</sup></td><td>61<sup>b</sup></td><td>6,100<sup>b</sup></td></tr><tr><td>benzaldehyde</td><td>&lt;1*</td><td>12.1</td><td>&gt;12,100</td><td>9.7<sup>b</sup></td><td>104<sup>b</sup></td><td>11,000<sup>b</sup></td><td>21<sup>b</sup></td><td>91<sup>b</sup></td><td>4,300<sup>b</sup></td></tr><tr><td>cinnamaldehyde</td><td>&lt;1*</td><td>13.3</td><td>&gt;13,300</td><td>34.7</td><td>240</td><td>6,900</td><td>31</td><td>29</td><td>950</td></tr><tr><td><strong>Alkanals</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>hexanal</td><td>3.1</td><td>7.3</td><td>2,300</td><td>112<sup>b</sup></td><td>142<sup>b</sup></td><td>1,300<sup>b</sup></td><td>5<sup>b</sup></td><td>28<sup>b</sup></td><td>5,600<sup>b</sup></td></tr><tr><td><strong>Alkenals</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>acrolein</td><td>36</td><td>8.8</td><td>240</td><td>110<sup>c</sup></td><td>116<sup>c</sup></td><td>1,070<sup>c</sup></td><td>884<sup>c</sup></td><td>11<sup>c</sup></td><td>12<sup>c</sup></td></tr><tr><td>trans-2-hexenal</td><td>5</td><td>11.3</td><td>2,200</td><td>28<sup>b</sup></td><td>49<sup>b</sup></td><td>1,700<sup>b</sup></td><td>9<sup>b</sup></td><td>16<sup>b</sup></td><td>1,800<sup>b</sup></td></tr><tr><td>4-hydroxy-2-nonenal</td><td>2.2</td><td>5.2</td><td>2,500</td><td>31<sup>c</sup></td><td>121<sup>c</sup></td><td>3,900<sup>c</sup></td><td>716<sup>c</sup></td><td>50<sup>c</sup></td><td>70<sup>c</sup></td></tr><tr><td>citral</td><td>1.5</td><td>5.5</td><td>3,700</td><td>4.5</td><td>35</td><td>7,400</td><td>35</td><td>68</td><td>1,750</td></tr><tr><td>farnesal</td><td>&lt;1*</td><td>4.8</td><td>&gt;4,800</td><td>2.3<sup>b</sup></td><td>30<sup>b</sup></td><td>13,000<sup>b</sup></td><td>37<sup>b</sup></td><td>27<sup>b</sup></td><td>700<sup>b</sup></td></tr><tr><td><strong>Ketones</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>2-butanone</td><td>780</td><td>10.5</td><td>13.5</td><td></td><td>LA</td><td></td><td></td><td>LA</td><td></td></tr><tr><td>3-buten-2-one</td><td>21.3</td><td>8.2</td><td>380</td><td></td><td>LA</td><td></td><td></td><td>LA</td><td></td></tr><tr><td>3-nonen-2-one</td><td>1.7</td><td>6.76</td><td>4,000</td><td></td><td>LA</td><td></td><td></td><td>LA</td><td></td></tr><tr><td>2-cyclohexen-1-one</td><td>365</td><td>4.41</td><td>12.1</td><td></td><td>LA</td><td></td><td></td><td>LA</td><td></td></tr><tr><td><strong>α-Dicarbonyls</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>2,3-butanedione</td><td>&lt;1*</td><td>10.6</td><td>&gt;10,600</td><td>540<sup>d</sup></td><td>260<sup>d</sup></td><td>480<sup>d</sup></td><td>110<sup>d</sup></td><td>23<sup>d</sup></td><td>210<sup>d</sup></td></tr><tr><td>2,3-hexanedione</td><td>&lt;1*</td><td>9.5</td><td>&gt;9,500</td><td>51<sup>b</sup></td><td>79<sup>b</sup></td><td>1,500<sup>b</sup></td><td>17<sup>b</sup></td><td>49<sup>b</sup></td><td>2,900<sup>b</sup></td></tr><tr><td><strong>β-Dicarbonyls</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>2,4-pentanedione</td><td>40</td><td>2.2</td><td>55</td><td>58,900</td><td>8.6</td><td>0.15</td><td>8,100</td><td>16.7</td><td>2.2</td></tr><tr><td>3,5-heptanedione</td><td>1,300</td><td>5.3</td><td>3.9</td><td>&gt;50,000</td><td>-</td><td>-</td><td>12,000</td><td>26</td><td>2.2</td></tr><tr><td><strong>Cofactor</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NADPH</td><td>5.7</td><td></td><td></td><td>10</td><td></td><td></td><td>2.9<sup>e</sup></td><td></td><td></td></tr></tbody></table>

Enzymatic activity was measured spectrophotometrically. For the determination of the kinetic parameters for NADPH, 3 and 60 mM D<sub>1</sub>L-glyceraldehyde was used with AKR1B15 and AKR1B10, respectively. LA, low activity (≤ 10 mU/mg) was detected at saturating concentration of substrate for AKR1B15; ND, not determined.

*Because of very low K<sub>m</sub> value, data are only approximate.

<sup>a</sup>Data from [23].

<sup>b</sup>Data from [66].

<sup>c</sup>Data from [67].

<sup>d</sup>Data from [68].

<sup>e</sup>Data from [69].

doi:10.1371/journal.pone.0134506.t001

AKR1B10 (i.e. a hydrogen bond between the N atom of Met265 and the 2'-phosphate group, a second hydrogen bond between the Nε of His269 and the 2'-phosphate group, and a stacking interaction of His269 with the adenine ring). Interestingly, a His residue at position 269 and its interactions have also been described in the rat AKR1B14 X-ray structure [53]. The salt bridge between the side-chain of Asp217 and Lys263, acting as a safety belt in cofactor binding, and the stacking interaction of Tyr210 side chain with the nicotinamide ring are also conserved.
<table><caption>Table 2. Kinetic constants with retinaldehyde isomers.</caption><thead><tr><th>Substrate and parameter</th><th>AKR1B15</th><th>AKR1B10<sup>a</sup></th><th>AKR1B1<sup>a</sup></th></tr></thead><tbody><tr><td><strong>All-trans-retinaldehyde</strong></td><td></td><td></td><td></td></tr><tr><td>K<sub>m</sub> (µM)</td><td>1 ± 0.3</td><td>0.6 ± 0.1</td><td>1.1 ± 0.1</td></tr><tr><td>k<sub>cat</sub> (min<sup>-1</sup>)</td><td>5.4 ± 0.5</td><td>27 ± 1</td><td>0.9 ± 0.1</td></tr><tr><td>k<sub>cat</sub>/K<sub>m</sub> (mM<sup>-1</sup>·min<sup>-1</sup>)</td><td>5,300 ± 1,700</td><td>45,000 ± 7,600</td><td>1,300 ± 160</td></tr><tr><td><strong>9-cis-retinaldehyde</strong></td><td></td><td></td><td></td></tr><tr><td>K<sub>m</sub> (µM)</td><td>0.16 ± 0.03</td><td>0.7 ± 0.1</td><td>0.4 ± 0.1</td></tr><tr><td>k<sub>cat</sub> (min<sup>-1</sup>)</td><td>3.8 ± 0.2</td><td>0.9 ± 0.1</td><td>0.7 ± 0.2</td></tr><tr><td>k<sub>cat</sub>/K<sub>m</sub> (mM<sup>-1</sup>·min<sup>-1</sup>)</td><td>25,600 ± 5,300</td><td>1,300 ± 190</td><td>1,500 ± 170</td></tr></tbody></table>

Enzymatic activity was measured by using the HPLC-based method.

doi:10.1371/journal.pone.0134506.t002

This is consistent with the AKR1B15 cofactor preference for NADP(H) and the low K<sub>m</sub> value of AKR1B15 for NADPH.

On the other hand, the active site displays high divergence between AKR1B15 and AKR1B10. The residue differences between the two proteins are concentrated in loops A and C which, along with loop B, give shape to the active-site pocket (Fig 4A). Thus, the catalytic residues (Asp44, Tyr49, Lys78 and His111) and those in loop B are strictly conserved. In contrast, Ser118, Leu122, Ala131, and Gly133 (in loop A), together with Cys299, Asn300, Val301, Leu302, Gln303, Ser304, and Tyr310 (in loop C) of AKR1B10 are substituted by Thr118, Phe122, Met131, Ser133, Phe299, Asp300, Phe301, Lys302, Glu303, Phe304 and Phe310, in AKR1B15. Noteworthy, residues Phe299, Phe301, Glu303 and Phe304, along with Phe48 (Val48 in AKR1B10), participate in the AKR1B15 active-site pocket, which is thus significantly smaller (60 Å³ for AKR1B15 versus 279 Å³ for AKR1B10) and more hydrophobic, compared to the AKR1B10 pocket (Fig 4B and 4C). Some of these substitutions might not only have a consequence on the shape, volume and hydrophobicity of the active site, but also on the flexibility of the polypeptide chain, which is most relevant in the loop regions. The analysis of

<table><caption>Table 3. Inhibitory effect of different compounds on enzymatic activity.</caption><thead><tr><th rowspan="2">Inhibitor</th><th colspan="3">IC<sub>50</sub> (µM)</th></tr><tr><th>AKR1B15</th><th>AKR1B10</th><th>AKR1B1</th></tr></thead><tbody><tr><td>Tolrestat</td><td>&gt; 100</td><td>0.006</td><td>0.01<sup>a</sup></td></tr><tr><td>Sorbinil</td><td>&gt; 100</td><td>9.6<sup>b</sup></td><td>0.55<sup>b</sup></td></tr><tr><td>JF0064</td><td>0.034 ± 0.005</td><td>1.0<sup>d</sup></td><td>0.3<sup>d</sup></td></tr><tr><td>Epalrestat</td><td>&gt; 50</td><td>0.33<sup>c</sup></td><td>0.021<sup>c</sup></td></tr><tr><td>Oleanolic acid</td><td>&gt; 100</td><td>0.09<sup>b</sup></td><td>124<sup>b</sup></td></tr><tr><td>Sulindac</td><td>&gt;100</td><td>0.35<sup>c</sup></td><td>0.21<sup>c</sup></td></tr><tr><td>Lithocholic acid</td><td>16.3 ± 7.6</td><td>0.12<sup>d</sup></td><td>7.2<sup>d</sup></td></tr></tbody></table>

The enzymatic activity assay with inhibitors was performed by using D,L-glyceraldehyde as a substrate. JF0064: 2,2',3,3',5,5',6,6'-octafluoro-4,4'-biphenyldiol

doi:10.1371/journal.pone.0134506.t003

<footnote><sup>a</sup>Data from [23]</footnote>

<footnote><sup>a</sup>Data from [70].</footnote>

<footnote><sup>b</sup>Data from [52].</footnote>

<footnote><sup>c</sup>Data from [71].</footnote>

<footnote><sup>d</sup>Data from [66].</footnote>
<figure><img src="image_6.png" /><figcaption>Fig 3. Molecular structure of compound JF0064, 2,2',3,3',5,5',6,6'-octafluoro-4,4'-biphenyldiol.</figcaption></figure>

doi:10.1371/journal.pone.0134506.g003

conformational ensembles indicates that AKR1B15 would display less flexible loops A and C than AKR1B10 (Fig 4D). A detailed analysis of the AKR1B15 model shows van der Waals interactions between residues from these loops and residues from other protein regions (i.e. Trp21-Phe299_loopC, Phe48-Phe116_loopA, Phe48-Phe123_loopA, Phe48-Phe301_loopA, Trp112_loopA-Phe301_loopC, Phe116_loopA-Phe304_loopC, Phe123_loopA-Phe304_loopC, Met131_loopA-Phe304_loopC and Met131_loopA-Leu307_loopC). Such interactions result in a hydrophobic cluster (Fig 4E), which likely contributes to the lower flexibility of the AKR1B15 active site.

## Docking of substrates and the inhibitor JF0064

As it has been described above, AKR1B15 is active towards retinoids, and thus the binding mode of all-trans- and 9-cis-retinaldehyde was analyzed. The obtained models were also
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Model of AKR1B15 structure. (A) Side view of the (α/β)₈ barrel. In red, the divergent residues between AKR1B15 and AKR1B10. (B) and (C) Active-site pockets of AKR1B15 and AKR1B10, respectively. The AKR1B15-specific residues are displayed in magenta. NADP⁺ cofactor is colored in orange. The surface contour of pockets is shown in grey and orange for AKR1B15 and AKR1B10, respectively. (D) The local conformational changes in *holo* forms of AKR1B15 (red line) and AKR1B10 (blue line) derived from computer simulations, as indicated by root mean square fluctuations (RMSF) of backbone atoms. The residues of loops A, B and C are highlighted by grey background. The difference in RMSF between the two enzymes is displayed as a black line in the top. (E) AKR1B15 loops A and C indicating potential contacts between different residues are shown as sticks, which may explain the low flexibility of the protein in this region. Figures have been drawn using PyMOL.</figcaption></figure>

doi:10.1371/journal.pone.0134506.g004

analyzed with the QMEAN server and displayed similar scores (0.68/1 and 0.69/1, respectively) as the AKR1B15 holoenzyme model. The analysis showed that both substrates could be placed with their carbonyl groups at catalytic distance from the hydroxyl group of Tyr49, the Nε of His111, and the cofactor C4 atom (2.9, 2.9, and 3.2 Å, respectively). The two molecules would be positioned in a similar manner into a narrow and hydrophobic pocket, establishing contacts with Trp21, Phe48, Phe123, Trp220, and Phe301 (all-*trans*-retinaldehyde), and Phe48, Trp220, Phe299, and Phe301 (9-*cis*-retinaldehyde) (Fig 6A). A slight rearrangement of loop A and loop
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Comparison of the cofactor-binding site between the AKR1B15 model (A) and the AKR1B10 crystal structure (B). Interactions of Met265 and His269 with NADP+ in AKR1B15 are similar to those of Val265 and Arg269 in AKR1B10 (black dotted lines). His269 forms a π-stacking interaction with the adenine ring of the cofactor. The substitution of Lys22 by Arg in AKR1B15 prevents its interaction with the pyrophosphate bridge of NADP+. The salt bridge between Asp217 and Lys263 (red dotted line), acting as a safety belt in the coenzyme binding, and the π-stacking interaction of Tyr210 with the cofactor nicotinamide ring are conserved between the two AKRs. Carbon atoms of the cofactor are shown in green, whereas those of the enzyme are colored grey. Figures have been drawn using PyMOL.</figcaption></figure>

doi:10.1371/journal.pone.0134506.g005

C (Fig 6A) could allow the establishment of a hydrogen-bond or an electrostatic interaction between Lys125 and Glu303 (Fig 6B).

The particular AKR1B15 topology brings the cyclohexene ring of all-*trans* and 9-*cis*-retinal-dehyde further away from loops A and C than in AKR1B10 (Fig 6A, 6B and 6C). Thus these substrates do not interact with the residues that have been described as important for all-*trans*-retinaldehyde binding in AKR1B10, and which are in the most external part of the AKR1B10 pocket [23]. Notably, some of these residues are not interacting with the retinoid substrates in AKR1B15, i.e. Lys125 and Met131 (Ala in AKR1B10) from loop A, Glu303 (Gln in AKR1B10) and Phe304 (Ser in AKR1B10) from loop C. In AKR1B10, all-*trans*-retinaldehyde binding requires a rearrangement of Lys125, not necessary for 9-*cis*-retinaldehyde, explaining the higher *k*cat value towards the former substrate [3]. In the case of AKR1B15, Lys125 is not involved in the binding of either substrate, being consistent with similar *k*cat values. The hydrophobic pocket in the external region of the active site of the enzyme, where part of the polyene chain and the cyclohexene ring of the substrate bind, corresponds to the protein region which shows more rigidity in comparison to AKR1B10 (Fig 4D). Furthermore, this rigidity, likely due to the presence of the bulky Phe residues (Phe48, Phe299, Phe301 and Phe304), would make difficult the opening of the so-called “specificity pocket” (Fig 6B and 6C), which has been described for ARI binding in AKR1B10 and AKR1B1 and which usually accommodates the hydrophobic moiety of inhibitors [52,54]. This feature could have important functional consequences. For instance, it could explain the different substrate specificity and inhibitor selectivity of AKR1B15, since flexibility has been well established in connection with the active-site accessibility, and substrate and ligand binding [55].

The compound JF0064 is the only ARI found to significantly inhibit AKR1B15, probably due to the reduced volume of the active-site cleft and to the difficulty in opening the “specificity pocket”. In order to analyze JF0064 binding, a docking simulation along with energy
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Molecular docking of substrates or inhibitors to the active-site pocket of AKR1B15. (A) Residues implicated in binding all-trans- and 9-cis-retinaldehyde are displayed in light and dark grey sticks; while the substrates are shown in light and dark blue, respectively. The residues found in the most external part of all-trans-retinaldehyde binding channel in AKR1B10 are highlighted in yellow. The energy minimized apo-conformation is displayed in magenta cartoon. (B) and (C) Side view of the surface contour of the active-site pocket, depicted in grey and orange for AKR1B15 and AKR1B10, respectively, to show the inhibitor "specificity pocket". A thick grey curved line indicates the "specificity pocket" in AKR1B10. As it is shown, this pocket may not be opened in AKR1B15, likely due to the presence of bulky Phe residues. (D) The inhibitor JF0064 (PDB ID 4ICC) bound to AKR1B15 is displayed as sticks with C atoms in magenta, while residues interacting with the inhibitor are shown as sticks with C atoms in grey. (E) Steric hindrance preventing tolrestat (in blue) and sorbinil (in orange) from binding to the active site of AKR1B15. For this analysis, the AKR1B15 structure model was superimposed with the AKR1B10 crystallographic structures with tolrestat (PDB ID 1ZUA) and sorbinil (PDB ID 4GA8). NADP+ is colored in orange. Figures have been drawn using PyMOL.</figcaption></figure>

doi:10.1371/journal.pone.0134506.g006

minimization was performed. The model was again validated by using the QMEAN server and displayed a similar, though slightly better score than the rest of the AKR1B15 models (0.72/1). The Z-score (which is tending from negative digits to 0 when approximating to experimentally determined structures) of the AKR1B15-NADP+-JF0064 model (−0.92) is comparatively
improved with respect to the AKR1B15 holoenzyme model (–1.64). All-*trans*- and 9-*cis*-retinaldehyde complexes show intermediate values of –1.51 and –1.35, respectively.

The conformation of the AKR1B15-NADP⁺-JF0064 complex corresponding to the energy minimum is displayed in Fig 6D. Binding would occur through van der Waals interactions with a large number of hydrophobic residues, and by establishing hydrogen bonds with catalytic residues (Tyr49 and His111) and Glu303. As it has been described above for substrates, the binding of the inhibitor also induced the same rearrangement of loop A and loop C, allowing for the interaction between Lys125 and Glu303. The other tested compounds did not inhibit AKR1B15 likely because of steric hindrance (e.g. Phe301 and Phe48 may clash against sorbinil, and Phe304, along with the fact that the “specificity pocket” could not be opened, may prevent tolrestat from binding) (Fig 6C and 6E).

## Physiological significance

The kcat/Km value of AKR1B15 for 9-*cis*-retinaldehyde (25,600 mM⁻¹·min⁻¹) is the highest among all the substrates checked so far, including the steroids and 3-keto-acyl-CoAs analyzed by Weber *et al.* [16], although these authors used different conditions for enzyme purification and kinetic studies. This is reminiscent of what has been observed for other members of the AKR superfamily, such as AKR1C3, which displays a high kcat/Km value for 9-*cis*-retinaldehyde but it is also active against steroids and prostaglandins [48]. Observations of dual substrate specificity have also been made for members of the short-chain dehydrogenase/reductase superfamily [56]. Thus, it is conceivable that AKR1B15 and some of these enzymes have a multifunctional role in the pre-receptor regulation of hormonal signaling pathways. The AKR1B15 specificity for the 9-*cis* isomer may suggest a major role in the control of RAR and RXR mediated signaling, but we cannot exclude other physiological functions. Regarding the reported localization of AKR1B15 in mitochondria, there is increasing evidence that retinoid metabolism takes place in different subcellular compartments [57,58], mitochondria being one of them. Carotenoids and their aldehyde metabolites can be generated by the asymmetric cleavage of β-carotene by β-carotene-9',10'-oxygenase (BCO2), which is associated with the inner mitochondrial membrane [59]. Therefore, a putative physiological role of AKR1B15 in retinoid metabolism is compatible with its mitochondrial localization. The presence in mitochondria of other retinaldehyde reductases, such as RDH13 [57], gives further support to this notion. In addition, RDH10, an enzyme involved in retinoic acid synthesis, shifts between mitochondria associated membranes and lipid droplets during retinyl ester biosynthesis, similarly to cellular retinol-binding protein type 1 [60]. Retinol has also been pinpointed as a modulator of energy homeostasis in mitochondria by regulating oxidative phosphorylation [61]. An additional role of AKR1B15 in metabolizing lipid peroxidation products and alkenals in mitochondria cannot be ruled out [62,63].

## Conclusions

Despite its high sequence identity with AKR1B10, AKR1B15 appears to be an enzyme with a unique substrate specificity and narrower inhibitor selectivity. AKR1B15 displays distinct kinetic features with ketones, α-dicarbonyl compounds and is among the best 9-*cis*-retinaldehyde reductases within the AKR superfamily. Some of the most potent inhibitors of AKR1B1 and AKR1B10 did not inhibit AKR1B15. Amino acid substitutions clustered in residues located in loops A and C result in a smaller, more hydrophobic and more rigid active-site pocket of AKR1B15 as compared to that of AKR1B10. The structural model of AKR1B15 provides a powerful tool for the virtual screening of substrates and inhibitors for this enzyme. The distinct topology of the AKR1B15 fold should facilitate the design of more selective inhibitors, as it has
been shown for other enzyme pairs with high sequence similarity [64,65]. Finally, the finding of all-*trans*- and 9-*cis*-retinaldehyde as substrates for AKR1B15 adds further complexity to the enzymatic pathways of retinoid transformations and their cross-talk with other hormonal signaling pathways, such as that of steroids. This opens a research line to elucidate the physiological contribution of this novel human retinaldehyde reductase.

## Acknowledgments

We thank Dr. A. de Marco for kindly providing chaperone-encoding plasmids pBB540 and pBB545.

## Author Contributions

Conceived and designed the experiments: JGD XP J. Farrés SP. Performed the experiments: JGD MHK FXR ACS AP J. Fanfrlík IC SP. Analyzed the data: JGD MHK FXR ACS J. Fanfrlík IC XP J. Farrés SP. Contributed reagents/materials/analysis tools: MHK FXR OAB. Wrote the paper: JGD MHK FXR ACS J. Fanfrlík IC XP J. Farrés SP.

## References

1. Jin Y, Penning TM. Aldo-keto reductases and bioactivation/detoxication. Annu Rev Pharmacol Toxicol. 2007; 47: 263–292. PMID: 16970545

2. Barski OA, Tipparaju SM, Bhatnagar A. The aldo-keto reductase superfamily and its role in drug metabolism and detoxification. Drug Metab Rev. 2008; 40: 553–624. PMID: 18949601

3. Ruiz FX, Gallego O, Ardèvol A, Moro A, Domínguez M, Alvarez S, et al. Aldo-keto reductases from the AKR1B subfamily: retinoid specificity and control of cellular retinoic acid levels. Chem Biol Interact. 2009; 178: 171–177. doi: 10.1016/j.cbi.2008.10.027 PMID: 19014918

4. Matsunaga T, Wada Y, Endo S, Soda M, El-Kabbani O, Hara A. Aldo-Keto Reductase 1B10 and Its Role in Proliferation Capacity of Drug-Resistant Cancers. Front Pharmacol. 2012; 3: 5. doi: 10.3389/fphar.2012.00005 PMID: 22319498

5. Ruiz FX, Moro A, Gallego O, Ardèvol A, Rovira C, Petrash JM, et al. Human and rodent aldo-keto reductases from the AKR1B subfamily and their specificity with retinaldehyde. Chem Biol Interact. 2011; 191: 199–205. doi: 10.1016/j.cbi.2011.02.007 PMID: 21329680

6. Salabei JK, Li X-P, Petrash JM, Bhatnagar A, Barski OA. Functional expression of novel human and murine AKR1B genes. Chem Biol Interact. 2011; 191: 177–184. doi: 10.1016/j.cbi.2011.01.020 PMID: 21276782

7. Oates PJ. Aldose reductase, still a compelling target for diabetic neuropathy. Curr Drug Targets. 2008; 9: 14–36. PMID: 18220710

8. Ma J, Luo DX, Huang C, Shen Y, Bu Y, Markwell S, et al. AKR1B10 overexpression in breast cancer: association with tumor size, lymph node metastasis and patient survival and its potential as a novel serum marker. Int J Cancer. 2012; 131: 862–871.

9. Tammali R, Reddy ABM, Saxena A, Rychahou PG, Evers BM, Qiu S, et al. Inhibition of aldose reductase prevents colon cancer metastasis. Carcinogenesis. 2011; 32: 1259–1267. doi: 10.1093/carcin/bgr102 PMID: 21642355

10. Endo S, Matsunaga T, Kuwata K, Zhao H-T, El-Kabbani O, Kitade Y, et al. Chromene-3-carboxamide derivatives discovered from virtual screening as potent inhibitors of the tumour maker, AKR1B10. Bioorg Med Chem. 2010; 18: 2485–2490. doi: 10.1016/j.bmc.2010.02.050 PMID: 20304656

11. Chatzopoulou M, Alexiou P, Kotsampasakou E, Demopoulos VJ. Novel aldose reductase inhibitors: a patent survey (2006—present). Expert Opin Ther Pat. 2012; 22: 1303–1323. doi: 10.1517/13543776. 2012.726615 PMID: 22998509

12. Pathania S, Randhawa V, Bagler G. Prospecting for novel plant-derived molecules of Rauvolfia serpentina as inhibitors of Aldose Reductase, a potent drug target for diabetes and its complications. PLoS One. 2013; 8: e61327. doi: 10.1371/journal.pone.0061327 PMID: 23613832

13. Soda M, Hu D, Endo S, Takemura M, Li J, Wada R, et al. Design, synthesis and evaluation of caffeic acid phenethyl ester-based inhibitors targeting a selectivity pocket in the active site of human aldo-keto reductase 1B10. Eur J Med Chem. 2012; 48: 321–329. doi: 10.1016/j.ejmech.2011.12.034 PMID: 22236472
14. Zhang L, Zhang H, Zheng X, Zhao Y, Chen S, Chen Y, et al. Structural basis for the inhibition of AKR1B10 by caffeic acid phenethyl ester (CAPE). ChemMedChem. 2014; 9: 706–709. doi: 10.1002/ cmdc.201300455 PMID: 24436249

15. Cousido-Siah A, Ruiz FX, Mitschler A, Porté S, de Lera AR, Martín MJ, et al. Identification of a novel polyfluorinated compound as a lead to inhibit the human enzymes aldose reductase and AKR1B10: structure determination of both ternary complexes and implications for drug design. Acta Crystallogr D Biol Crystallogr. 2014; 70: 889–903. doi: 10.1107/S1399004713033452 PMID: 24598757

16. Weber S, Salabei JK, Möller G, Kremmer E, Bhatnagar A, Adamski J, et al. Aldo-Keto Reductase 1B15 (AKR1B15): a Mitochondrial Human Aldo-Keto Reductase with Activity towards Steroids and 3-Keto-acyl-CoA Conjugates. J Biol Chem. 2015; 209: 6531–6545

17. Hackett NR, Butler MW, Shaykhiev R, Salit J, Omberg L, Rodriguez-Flores JL, et al. RNA-Seq quantification of the human small airway epithelium transcriptome. BMC Genomics. 2012; 13: 82. doi: 10.1186/1471-2164-13-82 PMID: 22375630

18. Constantinescu S, Hecht K, Sobotzki N, Erzinger MM, Bovet C, Shay JW, et al. Transcriptomic responses of cancerous and noncancerous human colon cells to sulforaphane and selenium. Chem Res Toxicol. 2014; 27: 377–386. doi: 10.1021/tx400427t PMID: 24383545

19. Calvo SE, Compton AG, Hershman SG, Lim SC, Lieber DS, Tucker EJ, et al. Molecular diagnosis of infantile mitochondrial disease with targeted next-generation sequencing. Sci Transl Med. 2012; 4: 118ra10. doi: 10.1126/scitranslmed.3003310 PMID: 22277967

20. Reumers J, De Rijk P, Zhao H, Liekens A, Smeets D, Cleary J, et al. Optimized filtering reduces the error rate in detecting genomic variants by short-read sequencing. Nat Biotechnol. 2012; 30: 61–8.

21. Yashin AI, Wu D, Arbeev KG, Ukraintseva SV. Joint influence of small-effect genetic variants on human longevity. Aging. 2010; 2: 612–620. PMID: 20834067

22. De Marco A, Deuerling E, Mogk A, Tomoyasu T, Bukau B. Chaperone-based procedure to increase yields of soluble recombinant proteins produced in *E. coli*. BMC Biotechnol. 2007; 7: 32. PMID: 17565681

23. Gallego O, Ruiz FX, Ardèvol A, Domínguez M, Alvarez R, de Lera AR, et al. Structural basis for the high all-trans-retinaldehyde reductase activity of the tumor marker AKR1B10. Proc Natl Acad Sci USA. 2007; 104: 20764–20769. PMID: 18087047

24. Larroy C, Fernandez MR, Gonzalez E, Parés X, Biosca JA. Characterization of the Saccharomyces cerevisiae YMR318C (ADH6) gene product as a broad specificity NADPH-dependent alcohol dehydrogenase: relevance in aldehyde reduction. Biochem J. 2002; 172: 163–172.

25. Gallego O, Belyaeva OV, Porté S, Ruiz FX, Stetsenko AV, Shabrova EV, et al. Comparative functional analysis of human medium-chain dehydrogenases, short-chain dehydrogenases/reductases and aldoketo reductases with retinoids. Biochem J. 2006; 399: 101–109. PMID: 16787387

26. Kuksa V, Imanishi Y, Batten M, Palczewski K, Moise AR. Retinoid cycle in the vertebrate retina: experimental approaches and mechanisms of isomerization. Vision Res. 2003; 43: 2959–2981. PMID: 14611933

27. Belyaeva OV, Korkina OV, Stetsenko AV, Kim T, Nelson PS, Kedishvili NY. Biochemical properties of purified human retinol dehydrogenase 12 (RDH12): catalytic efficiency toward retinoids and C9 aldehydes and effects of cellular retinol-binding protein type I (CRBPI) and cellular retinaldehyde-binding protein (CRALBP) on the oxidation and reduction of retinoids. Biochemistry. 2005; 44: 7035–7047. PMID: 15865448

28. Canutescu AA, Shelenkov AA, Dunbrack RL. A graph-theory algorithm for rapid protein side-chain prediction. Protein Sci. 2003; 12: 2001–2014. PMID: 12930999

29. Pronk S, Páll S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, et al. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013; 29: 845–854. doi: 10.1093/bioinformatics/btt055 PMID: 23407358

30. Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins. 2006; 65: 712–725. PMID: 16981200

31. Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, et al. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010; 78: 1950–1958. doi: 10.1002/ prot.22711 PMID: 20408171

32. Holmberg N, Ryde U, Bülow L. Redesign of the coenzyme specificity in L-lactate dehydrogenase from bacillus stearothermophilus using site-directed mutagenesis and media engineering. Protein Eng. 1999; 12: 851–856. PMID: 10556245

33. Hawkins GD, Cramer CJ, Truhlan DG. Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium. J Phys Chem. 1996; 100: 19824–19839.
34. Seeliger D, Haas J, de Groot BL. Geometry-based sampling of conformational transitions in proteins. Structure. 2007; 15: 1482–1492. PMID: 17997973

35. De Groot BL, van Aalten DM, Scheek RM, Amadei A, Vriend G, Berendsen HJ. Prediction of protein conformational freedom from distance constraints. Proteins. 1997; 29: 240–251. PMID: 9329088

36. Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereo-chemical quality of protein structures. J Appl Crystallogr. International Union of Crystallography; 1993; 26: 283–291.

37. Eisenberg D, Lüthy R, Bowie JU. VERIFY3D: assessment of protein models with three-dimensional profiles. Methods Enzymol. 1997; 277: 396–404. PMID: 9379925

38. Goodsell DS, Morris GM, Olson AJ. Automated docking of flexible ligands: applications of AutoDock. J Mol Recognit. 9: 1–5. PMID: 8723313

39. Porté S, Ruiz FX, Giménez J, Molist I, Alvarez S, Domínguez M, et al. Aldo-keto reductases in retinoid metabolism: search for substrate specificity and inhibitor selectivity. Chem Biol Interact. 2013; 202: 186–194. doi: 10.1016/j.cbi.2012.11.014 PMID: 23220004

40. ezáč J, Fanfrlík J, Salahub D, Hobza P. Semiempirical Quantum Chemical PM6 Method Augmented by Dispersion and H-Bonding Correction Terms Reliably Describes Various Types of Noncovalent Complexes. J Chem Theory Comput. 2009; 5: 1749–1760.

41. Klamt A, Schüürmann G. COSMO: a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J Chem Soc. 1993; 2: 799–805.

42. Lepšík M, ezáč J, Kolář M, Pecina A, Hobza P, Fanfrlík J. The Semiempirical Quantum Mechanical Scoring Function for In Silico Drug Design. Chempluschem. 2013; 78: 921–931.

43. Fanfrlík J, Kolář M, Kamlar M, Humý D, Ruiz FX, Cousido-Siah A, et al. Modulation of aldose reductase inhibition by halogen bond tuning. ACS Chem Biol. 2013; 8: 2484–2492. doi: 10.1021/cb400526n PMID: 23988122

44. Benkert P, Tosatto SCE, Schomburg D. QMEAN: A comprehensive scoring function for model quality assessment. Proteins. 2008; 71: 261–277. PMID: 17932912

45. Durrant JD, de Oliveira CAF, McCammon JA. POVME: an algorithm for measuring binding-pocket volumes. J Mol Graph Model. 2011; 29: 773–776. doi: 10.1016/j.jmgm.2010.10.007 PMID: 21147010

46. Tao H, Liu W, Simmons BN, Harris HK, Cox TC, Massiah M A. Purifying natively folded proteins from inclusion bodies using sarkosyl, Triton X-100, and CHAPS. Biotechniques. 2010; 48: 61–64. doi: 10.2144/000113304 PMID: 20078429

47. Pastel E, Pointud J-C, Volat F, Martinez A, Lefrançois-Martinez A-M. Aldo-Keto Reductases 1B in Endocrinology and Metabolism. Front Pharmacol. 2012; 3: 148. doi: 10.3389/fphar.2012.00148 PMID: 22876234

48. Ruiz FX, Porté S, Gallego O, Moro A, Ardèvol A, Del Río-Espínola A, et al. Retinaldehyde is a substrate for human ald-o-keto reductases of the 1C subfamily. Biochem J. 2011; 440: 335–344. doi: 10.1042/BJ20111286 PMID: 21851338

49. Crosas B, Cederlund E, Torres D, Jornvall H, Farrés J, Parés X. A vertebrate ald-o-keto reductase active with retinoids and ethanol. J Biol Chem. 2001; 276: 19132–19140. PMID: 11278684

50. Verma M, Martin HJ, Haq W, O'Connor TR, Maser E, Balendiran GK. Inhibiting wild-type and C299S mutant AKR1B10; a homologue of aldose reductase upregulated in cancers. Eur J Pharmacol. 2008; 584: 213–221. doi: 10.1016/j.ejphar.2008.01.036 PMID: 18325492

51. Ruiz FX, Cousido-Siah A, Mitschler A, Farrés J, Parés X, Podjarny A. X-ray structure of the V301L ald-o-keto reductase 1B10 complexed with NADP+ and the potent aldose reductase inhibitor fidarestat: Implications for inhibitor binding and selectivity. Chem Biol Interact. 2013; 202: 178–185. doi: 10.1016/j.cbi.2012.12.013 PMID: 23295227

52. Zhang L, Zhang H, Zhao Y, Li Z, Chen S, Zhai J, et al. Inhibitor selectivity between ald-o-keto reductase superfamily members AKR1B10 and AKR1B1: role of Trp112 (Trp111). FEBS Lett. 2013; 587: 3681–3686. doi: 10.1016/j.febslet.2013.09.031 PMID: 24100137

53. Sundaram K, Dhagat U, Endo S, Chung R, Matsunaga T, Hara A, et al. Structure of rat aldose reductase-like protein AKR1B14 holoenzyme: Probing the role of His269 in coenzyme binding by site-directed mutagenesis. Bioorg Med Chem Lett. 2011; 21: 801–804. doi: 10.1016/j.bmcl.2010.11.086 PMID: 21168333

54. Urzhumtsev A, Tête-Favier F, Mitschler A, Barbanton J, Barth P, Urzhumtseva L, et al. A "specificity" pocket inferred from the crystal structures of the complexes of aldose reductase with the pharmaceutically important inhibitors tolrestat and sorbinil. Structure. 1997; 5: 601–612. PMID: 9195881

55. Kokkinidis M, Glykos NM, Fadouloglou VE. Protein Flexibility and Enzymatic Catalysis. Adv Protein Chem Struct Biol. 2012. pp. 181–218. doi: 10.1016/B978-0-12-398312-1.00007-X PMID: 22607756
56. Kedishvili NY. Enzymology of retinoic acid biosynthesis and degradation. J Lipid Res. 2013; 54: 1744-1760. doi: 10.1194/jlr.R037028 PMID: 23630397

57. Belyaeva OV, Korkina OV, Stetsenko AV, Kedishvili NY. Human retinol dehydrogenase 13 (RDH13) is a mitochondrial short-chain dehydrogenase/reductase with a retinaldehyde reductase activity. FEBS J. 2008; 275: 138-147. PMID: 18039331

58. Ruiz FX, Porté S, Parés X, Farrés J. Biological role of aldo-keto reductases in retinoic Acid biosynthesis and signaling. Front Pharmacol. 2012; 3: 58. doi: 10.3389/fphar.2012.00058 PMID: 22529810

59. Palczewski G, Amengual J, Hoppel CL, von Lintig J. Evidence for compartmentalization of mammalian carotenoid metabolism. FASEB J. 2014; 1-13. doi: 10.1096/fj.14-0101LTE

60. Jiang W, Napoli JL. The retinol dehydrogenase Rdh10 localizes to lipid droplets during acyl ester biosynthesis. J Biol Chem. 2013; 288: 589-597. doi: 10.1074/jbc.M112.402883 PMID: 23155051

61. Acin-Perez R, Hoyos B, Zhao F, Vinogradov V, Fischman D a, Harris R a, et al. Control of oxidative phosphorylation by vitamin A illuminates a fundamental role in mitochondrial energy homoeostasis. FASEB J. 2010; 24: 627-636. doi: 10.1096/fj.09-142281 PMID: 19812372

62. Echtay KS, Brand MD. 4-hydroxy-2-nonenal and uncoupling proteins: an approach for regulation of mitochondrial ROS production. Redox Rep. 2007; 12: 26-29. PMID: 17263904

63. Zhong H, Yin H. Role of lipid peroxidation derived 4-hydroxynonenal (4-HNE) in cancer: Focusing on mitochondria. Redox Biol. 2015; 4: 193-199. doi: 10.1016/j.redox.2014.12.011 PMID: 25598486

64. Luong C, Miller A, Barnett J, Chow J, Ramesha C, Browner MF. Flexibility of the NSAID binding site in the structure of human cyclooxygenase-2. Nat Struct Biol. 1996; 3: 927-933. PMID: 8901870

65. Kimble-Hill AC, Parajuli B, Chen CH, Mochly-Rosen D, Hurley TD. Development of selective inhibitors for aldehyde dehydrogenases based on substituted indole-2,3-diones. J Med Chem. 2014; 57: 714-722. doi: 10.1021/jm401377v PMID: 24444054

66. Endo S, Matsunaga T, Mamiya H, Ohta C, Soda M, Kitade Y, et al. Kinetic studies of AKR1B10, human aldose reductase-like protein: endogenous substrates and inhibition by steroids. Arch Biochem Bio-phys. 2009; 487: 1-9. doi: 10.1016/j.abb.2009.05.009 PMID: 19464995

67. Shen Y, Zhong L, Johnson S, Cao D. Human aldo-keto reductases 1B1 and 1B10: a comparative study on their enzyme activity toward electrophilic carbonyl compounds. Chem Biol Interact. 2011; 191: 192-198. doi: 10.1016/j.cbi.2011.02.004 PMID: 21329684

68. Calam E, Porté S, Fernández MR, Farrés J, Parés X, Biosca JA. Biocatalytic production of alpha-hydroxy ketones and vicinal diols by yeast and human aldo-keto reductases. Chem Biol Interact. 2013; 202: 195-203. doi: 10.1016/j.cbi.2012.12.006 PMID: 23295224

69. Kubiseski TJ, Flynn TG. Studies on human aldose reductase. Probing the role of arginine 268 by site-directed mutagenesis. J Biol Chem. 1995; 270: 16911-16917. PMID: 7622508

70. Ehrig T, Bohren KM, Prendergast FG, Gabbay KH. Mechanism of aldose reductase inhibition: binding of NADP+/NADPH and alrestatin-like inhibitors. Biochemistry. 1994; 33: 7157-7165. PMID: 8003482

71. Endo S, Matsunaga T, Soda M, Tajima K, Zhao H-T, El-Kabbani O, et al. Selective inhibition of the tumor marker AKR1B10 by antiinflammatory N-phenylanthranilic acids and glycyrrhetic acid. Biol Pharm Bull. 2010; 33: 886-890. PMID: 20460771