key: cord-0055214-ugocr73g
authors: Tanawattanasuntorn, Tanotnon; Thongpanchang, Tienthong; Rungrotmongkol, Thanyada; Hanpaibool, Chonnikan; Graidist, Potchanapond; Tipmanee, Varomyalin
title: (−)-Kusunokinin as a Potential Aldose Reductase Inhibitor: Equivalency Observed via AKR1B1 Dynamics Simulation
date: 2020-12-21
journal: ACS Omega
DOI: 10.1021/acsomega.0c05102
sha: f196a6db38382e3725d5f76872aa14014fd5d269
doc_id: 55214
cord_uid: ugocr73g

[Image: see text] (−)-Kusunokinin performed its anticancer potency through CFS1R and AKT pathways. Its ambiguous binding target has, however, hindered the next development phase. Our study thus applied molecular docking and molecular dynamics simulation to predict the protein target from the pathways. Among various candidates, aldo-keto reductase family 1 member B1 (AKR1B1) was finally identified as a (−)-kusunokinin receptor. The predicted binding affinity of (−)-kusunokinin was better than the selected aldose reductase inhibitors (ARIs) and substrates. The compound also had no significant effect on AKR1B1 conformation. An intriguing AKR1B1 efficacy, with respect to the known inhibitors (epalrestat, zenarestat, and minalrestat) and substrates (UVI2008 and prostaglandin H(2)), as well as a similar interactive insight of the enzyme pocket, pinpointed an ARI equivalence of (−)-kusunokinin. An aromatic ring and a γ-butyrolactone ring shared a role with structural counterparts in known inhibitors. The modeling explained that the aromatic constituent contributed to π–π attraction with Trp111. In addition, the γ-butyrolactone ring bound the catalytic His110 using hydrogen bonds, which could lead to enzymatic inhibition as a consequence of substrate competitiveness. Our computer-based findings suggested that the potential of (−)-kusunokinin could be furthered by in vitro and/or in vivo experiments to consolidate (−)-kusunokinin as a new AKR1B1 antagonist in the future.

Chemotherapy is commonly used in the treatment of cancer. Chemotherapeutic agent development is ongoing, so as to increase treatment efficacy and reduce adverse effects. 1 The major source of synthetic chemotherapeutic agents is active plant-derived chemicals. Various types of natural products have been successfully developed into practical chemotherapeutic drugs. 2 Lignans, a dimeric phenylpropanoid from natural products, provide a wide range of biological functions. In particular, flaxseed lignans have been revealed to have anticancer effects by multiple targets of the hallmarks of cancer. 3 Well-known cases include lignan-based anticancer drugs, namely, etoposide and other podophyllotoxin semisynthetic derivatives. 4, 5 Another reported lignan, arctigenin and its glycoside, arctiin, hindered the cancer cell growth. 6, 7 Arctigenin prevented the osteoarthritis progression via targeting of the PI3K/AKT pathway. 8 Arctiin possessed antiangiogenic effects on human myeloma cells via inhibition of STAT3 phosphorylation. 9 In addition, arctigenin and arctiin also exhibited anti-inflammatory effects through the inhibition of reactive oxygen species (ROS) production and COX-2. 10 (−)-Kusunokinin, one of the reported lignan compounds, exhibited anticancer activity on various human cancer cell lines, including breast cancer, colon cancer, and cholangiocarcinoma. 11 This compound inhibited tumor growth with no side effects and exhibited the normal complete blood count parameters and clinical chemistry of the renal and liver function test in the breast cancer rat model. 12 The anticancer mechanism of (−)-kusunokinin could occur through the partial binding with CSF1R, resulting in the suppression of AKT and its downstream proteins, including c-myc, cyclinB1, and cyclinD1. 12, 13 Computational modeling of (−)-kusunokinin was carried out to investigate its protein target based on a similar case of bursehernin, 13, 14 a structural isomer of (−)-kusunokinin. Colony stimulating factor 1 receptor (CSF1R) was reported as a possible target of these compounds, in which the stabilized π−π interaction with aromatic amino acids played a crucial role. This observation suggested that the selectivity of (−)-kusunokinin depended on the π−π stacking with aromatic amino acids at the binding pocket, such as tyrosine, phenylalanine, or tryptophan. In addition, the (−)-kusunokinin structure could perform interactions via various hydrogenbond acceptor counterparts, namely, an oxygen atom(s) of the methoxy group (−OCH 3 ) or/and an oxygen atom at the ketone group on the γ-butyrolactone ring. 13, 14 As we had already noticed some clues in the interaction and specificity of (−)-kusunokinin, we thus speculated that another protein containing an aromatic-rich domain could behave with (−)-kusunokinin in the same manner. Moreover, due to the partial binding with CSF1R, the downstream molecules of the CSF1R and AKT signaling pathways were selected for screening for the possibility of (−)-kusunokinin target proteins in this study. Herein, we used the molecular docking technique to predict an additional target protein of (−)-kusunokinin. The time-dependent conformational behavior of the (−)-kusunokinin−target complex was elucidated using molecular dynamics simulation. Then, the binding manner of (−)-kusunokinin was also evaluated and compared with the reported inhibitors. Finally, our present work provided more information regarding the target protein of (−)-kusunokinin with aldo-keto reductase family 1 member B1 (AKR1B1) and the downstream molecules of CSF1R and AKT. (−)-Kusunokinin bound this protein at the anion binding pocket. These findings could become constructive for furthering the (−)-kusunokinin development phase as a specifically targeted anticancer inhibitor.

2.1. Protein Structure Preparation. Structures of 41 proteins associated with the CSF1R and AKT pathways were retrieved from the Protein Data Bank (https://www.rcsb.org/) in a PDB format file. All atoms were selected based on the highest occupancy. Crystallographic waters were removed and polar hydrogen atoms were added using AutoDockTools 1.5.6 (ADT). 15 The protein structures were written into a PDBQT format file.

2.2. Ligand Preparation. The structure files of (−)-kusunokinin and the other ligands used were obtained from the PubChem database. (−)-Kusunokinin was taken from a structure with CID 384876. The crystal-bound (native) ligands were directly separated from the crystallography structure using the Visual Molecular Dynamics (VMD) package. 16 All ligands used in this study are summarized in Table S1 . The structure was then converted into a PDB file using the Biovia Discovery Studio 2019 Client. 17 Hydrogen atoms were added to all ligands. Gasteiger charges and default atom parameters were automatically assigned from ADT. Finally, the structure was written into the PDBQT file format using ADT.

2.3. Molecular Docking Parameters. Molecular docking studies were performed using the AutoDock4 program. During the process, the protein structure was set as a rigid molecule with a flexible ligand. The docking site on the protein target was defined by a grid box of 120 × 120 × 120 cubic angstrom (Å 3 ), with the grid spacing at the center of the protein structure. Other parameters followed the default values in ADT. Fifty genetic algorithm (GA) runs with a population size of 200 were executed for conformational sampling. The conformation with the lowest binding energy was considered as the best pose. The interaction of the best docked pose was compared with the native ligands in the crystallographic structure.

2.4. Molecular Docking Study and Binding Energy Calculation. The reproducibility of the AutoDock4 program was validated via a self-redocking experiment of the native reference ligand. The all-atom root-mean-square deviation (RMSD) cutoff of less than 2 Å between docked and native crystal structures was set as the criterion for a successful procedure. The AutoDock4 program score was recorded from the binding free energy between the ligand and the protein. 18, 19 2.5. Molecular Dynamics Simulation. The AKR1B1 structure selected was the PDB structure code 1US0 (a resolution of 0.66 Å). All water and other solvents were removed. Only AKR1B1 with the cofactor dihydro-nicotinamide-adenine-dinucleotide phosphate (NADPH) was retained. All polar hydrogen in the best ligand docked pose was eliminated. AMBER16 force field was exploited to model the drug−AKR1B1 complex, similar to previous reports of AKR1B1 molecular dynamics study. 20−23 The hydrogen atom was instead replaced via the leap module implemented in the AMBER16 package. 24 An ionizable side chain of the amino acid was considered at pH 7. The side chain in lysine or arginine contained +1e charge, while the glutamate and aspartate side chains contained −1e charge. The imidazole ring in histidine was neutral. No disulfide was observed in the AKR1B1 structure. All ligand molecules, namely, (−)-kusunokinin, AKR1B1 inhibitors, and a substrate (PGH 2 ), were parameterized using the Gaussian16 package, 25 based on calculations of restrained electrostatic atomic potential (RESP). NADPH parameters were adopted from a previous report. 26 The structure of ligand-free AKR1B1 or AKR1B1 with the best docked conformer was solvated by TIP3P water with a distance of 14 Å from the surface, leading to approximately 12000 water molecules. Sodium ion (Na + ) was used to neutralize the system, and 18 NaCl pairs were included to yield a 0.1 M NaCl solution. The aforementioned AKR1B1−ligand complex was used as an initial coordinate for the molecular dynamics process. The AKR1B1 solution was first equilibrated in a canonical (NVT) ensemble using Langevin dynamics at 310 K (37°C). The harmonic restraint potential was applied to the positions of AKR1B1, NADPH, and the ligand, with a force constant of 200, 100, 50, 20, and 10 kcal/mol-Å 2 , respectively. Each restrained NVT ensemble lasted for 200 ps using a time step of 1 fs. The system was switched into an NPT ensemble via the Berendsen algorithm at 310 K and 1.013 bar.

The molecular dynamics (MD) simulation was run for 150 ns with a time step of 2 fs. An MD trajectory with 2500 equidistant snapshots was obtained. The last 60-ns-MD conformation (1500 snapshots) was taken for conformation and binding analysis. Root-mean-square deviation (RMSD) calculation and structure visualization were conducted using the VMD package. Root-mean-square fluctuation (RMSF) was computed, and the pattern of distance geometry was acquired using the cpptraj module.

2.6. Binding Free Energy Evaluation. Average binding free energy (ΔG) was calculated using the molecular mechanics/generalized Born-surface area (MM/GBSA) 27−30 method from MD trajectories. Average molecular mechanics energy was obtained from the AMBER functional formula. 18 Table  S2 . (−)-Kusunokinin exhibited the best predicted binding potential toward aldo-keto reductase family 1 member B1 (AKR1B1), with a binding energy of −11.11 kcal/mol. The other candidate was the dual specificity mitogen-activated protein kinase kinase2 (MEK2), with a binding energy of −10.14 kcal/mol.

In spite of the closeness in the predicted binding energy between AKR1B1 and MEK2, AKR1B1 piqued our interest not only for its lower binding energy but also due to some supporting data regarding its lignan-based inhibitor. AKR1B1, an NADPH-dependent enzyme, is a target for diabetic complication treatment through the regulation of the polyol metabolic glucose pathway. 32 Recently, AKR1B1 was found in various human cancers, 33 and it could be the novel target for cancer aggressiveness reduction. 34−36 Arctigenin and arctiin exhibited potency in the aldose reductase (AKR1B1) 

http://pubs.acs.org/journal/acsodf Article inhibition. 37−39 Since the lignans share a feature for the aromatic ring structure, the π−π stacking with the aromatic binding site of AKR1B1 would then be expected to be like CSF1R. 13, 14 3.2. Preliminary Atomistic View of (−)-Kusunokinin on AKR1B1 Binding. (−)-Kusunokinin was suitably inserted in the AKR1B1 active pocket cleft, which was surrounded by several hydrophobic and aromatic residues (Figure 2 ). The site was categorized into three apparent pockets: an anion binding pocket (Tyr48, His110, and Trp111), a small hydrophobic pocket (Trp20, Val47, and Tyr48), and a large hydrophobic specific pocket (Trp79, Trp111, Phe115, Phe122, Trp219, Cys298, Leu300, Cys303, and Tyr309). These active pockets resembled previously reported AKR1B1 binding with βaminophenylpropanoic acid derivatives. 23 Later, a molecular docking study of (−)-kusunokinin against AKR1B1 was performed and compared with the other available AKR1B1 substrates, 40,41 commercial carboxylic acid aldose reductase inhibitors (ARIs), commercial hydantoin ARIs, and some suggested potential ARIs, respectively. 42−44 The result is summarized in Table 1 . Interestingly, the binding energy of (−)-kusunokinin (−11.11 kcal/mol) was in a similar range to the reported ARIs, from −6.92 to −11.54 kcal/mol. Moreover, (−)-kusunokinin provided a better predicted energy than substrates (PGH 2 and retinol), substrate analogs (UVI2008 and PGA1), and some reported potential ARIs. (−)-Kusunokinin was likewise inserted in the active pocket with a similar position and orientation to the well-known ARIs (Figure 3 ).

Trp20 and Trp111 sustained ARI conformation by forming π−π stacking with the aromatic rings of ARIs. Furthermore, Tyr48 and His110 also formed hydrogen bonds with the AKR1B1-bound ligand. These interactions, found in most ARIs, are concordant with previous reports that Trp20, Tyr48, His110, and Trp111 are present in the anion binding pocket. 32 The results gave a hint that the aromatic ring-containing compounds could be the potential ARIs due to π−π stacking formation at the active pocket.

Since molecular docking was performed on the basis of rigid proteins with a flexible ligand, the study of molecular dynamics simulation was later carried out to investigate the timedependent characteristics of AKR1B1 and its bound compounds.

3.3. Molecular Dynamics Trajectory Analysis. Molecular dynamics simulation was performed to reflect the more realistic dynamic behaviors of the drug-free AKR1B1 or drug-AKR1B1 complex. All of the MD trajectories of all systems became stable after 90 ns (Figure 4a ). The flexibility pattern of ligand−AKR1B1 complexes was similar mostly throughout the structure. A distinctive exception was observed in the AKR1B1 residue 210 to 230 in PGH 2 -AKR1B1 and UVI2008-AKR1B1, which was distinctively more flexible (Figure 4b) .

The next point was to investigate whether the bound ligand can affect the AKR1B1 conformation. The distance between a center of geometry in each amino acid and an origin point was plotted to give the pattern of distance. The reference structure was ligand-free AKR1B1 MD simulation. The same distance pattern indicated a similar conformation of the assigned amino acid. However, a different distance pattern pinpointed a different conformation from the AKR1B1. Herein, similar patterns among the seven AKR1B1 simulations were observed ( Figure 5 ). The conclusive finding was that (−)-kusunokinin and other ARIs did not affect the AKR1B1 conformation.

AKR1B1 binding affinity of (−)-kusunokinin was considered from an average MM/GBSA binding free energy with respect to other experimentally known ARIs, which were zenarestat, minalrestat, and epalrestat ( Table 2) . MM/GBSA calculations revealed that the binding free energy of (−)-kusunokinin (−42.23 kcal/mol) was comparable to other selected ARIs, ranging from −29.63 to −38.96 kcal/mol, and suggested that (−)-kusunokinin would bind with AKR1B1 like other inhibitors.

An oxygen atom at the ketone group on the γ-butyrolactone ring of (−)-kusunokinin interacted with His110 via hydrogen bonding. Meanwhile, the aromatic ring on 1,2-dimethoxybenzene was stabilized by π−π stacking with Trp111 throughout the MD progression (Figure 6a,b) . The γ-butyrolactone ring acted equivalently to the carboxylic acid group of epalrestat and zenarestat to access into a key catalytic His110 in an AKR1B1 anion binding site. The hydantoin ring of minalrestat, however, formed temporary hydrogen bonds with His110. This ACS Omega http://pubs.acs.org/journal/acsodf Article observed finding is concordant with the reported role of His110 with carboxylic acid or hydantoin ARIs. 45 Furthermore, Trp111 reportedly acted as a key residue in ARI specificity. 32, 46 Our study showed that an aromatic ring of (−)-kusunokinin was stabilized by π−π stacking with Trp111 ( Figure 6a ), similar to selected ARIs (Figure 6b−d) and UVI2008 (an AKR1B1 substrate analog) (Figure 6e ). Although no aromatic ring was found in the PGH 2 structure, Trp111 was also responsible in PGH 2 binding, using π−alkyl interaction with the PGH 2 hydrocarbon chain (Figure 6f) . Besides, Tyr48 could temporarily form hydrogen bonds with (−)-kusunokinin and ARIs.

3.4. Comparative Behavior of (−)-Kusunokinin with ARIs and AKR1B1 Substrate. AKR1B1 requires NADPH cofactor to catalyze the various substrates, including prostaglandin H 2 (PGH 2 ). PGH 2 is converted into PGF 2α and relays the signals in the NF-κB pathway to regulate the epithelial−mesenchymal transition (EMT) process on various cancers. The catalytic imidazole ring in His110 acts as a proton donor for the dioxabicycloheptane ring of PGH 2 to produce PGF 2α , 4040 as shown in Figure 7a .

The molecular dynamics simulations showed that the PGH 2 dioxabicycloheptane ring made contact with a transferable proton of NADPH and His110, as expected. The stable hydrogen bond formation with a catalytic His110 due to (−)-kusunokinin, as well as epalrestat, directly prevented PGH 2 access into the catalytic site (Figure 7b ). However, no stable hydrogen bond with His110 was noticed in the other ARIs. Instead, aromatic residue Trp111 was found to be the principal stabilizer of (−)-kusunokinin, carboxylic and hydantoin ARIs. Mutual interaction with Trp111 would be a rather major contributor to AKR1B1 inhibition because the ligand hindered PGH 2 access into the catalytic site, no matter whether or not the compound could form a hydrogen bond with His110 ( Figure 8) .

The binding behaviors of (−)-kusunokinin were closer to carboxylic ARI than hydantoin ARIs. In Figure 8 , (−)-kusunokinin could bind AKR1B1 via its carbonyl oxygen in the γbutyrolactone ring and aromatic ring (Figure 8a ). The γbutyrolactone ring resembled a carboxylic group carboxylic ARI. Carboxylic ARIs have exerted carboxylate groups into anion binding pockets, 47 including His110, similar to the carbonyl oxygen of (−)-kusunokinin (Figure 8a,b) . This structural information could also be supported by ARIreported lignans (arctigenin and arctiin) 37−39 because their core structures contain the γ-butyrolactone ring and aromatic entity (Figure 9 ). Although many ARIs were synthesized, several ARIs were withdrawn in consequence of the severe adverse effects of hepatotoxicity 48−50 and hypersensitivity. 51, 52 Only commercial epalrestat is used in Asian countries for diabetic neuropathy treatment 53 and metastatic triple-negative breast cancer in clinical trials. 54 However, its poor efficiency and adverse effects remain a problem in other countries. 55 Hence, the search for new, efficient ARIs with low toxicity is important.

Related to (−)-kusunokinin action in CSF1R and AKT pathways, AKR1B1 was found as a tentative (−)-kusunokinin inhibitors. An aromatic ring acted as an AKR1B1 binding counterpart through π−π stacking with Trp11. An obstruction in AKR1B1 function was also expected because of the hydrogen bond formation of the γ-butyrolactone ring with a crucial catalytic His110. Due to persuasive in silico 

From chemotherapy to biological therapy: A review of novel concepts to reduce the side effects of systemic cancer treatment

Plants as a source of anti-cancer agents

Flaxseed lignans as important dietary polyphenols for cancer prevention and treatment: Chemistry, pharmacokinetics, and molecular targets

Podophyllotoxin: distribution, sources, applications and new cytotoxic derivatives

Podophyllotoxin: a novel potential natural anticancer agent

Arctigenin inhibits prostate tumor growth in high-fat diet fed mice through dual actions on adipose tissue and tumor

Arctiin induces cell growth inhibition through the down-regulation of cyclin D1 expression

Arctigenin prevents the progression of osteoarthritis by targeting PI3K/Akt/NF-κB axis: In vitro and in vivo studies

Arctiin is a pharmacological inhibitor STAT3 phosphorylation at tyrosine 705 residue and potentiates bortezomib-induced apoptotic and antiangiogenic effects in human multiple myeloma cells

Overview of the anti-inflammatory effects, pharmacokinetic properties and clinical efficacies of arctgenin and arctiin from Arctium lappa L

Anticancer activity of synthetic (±)-kusunokinin and its derivative (±)-bursehernin on human cancer cell lines

-)-kusunokinin inhibits breast cancer in N-nitrosomethylureainduced mammary tumor rats

Inhibition of CSF1R and AKT by (±)-kusunokinin hinders breast cancer cell proliferation

Structure-guided cancer blockade between bioactive bursehernin and proteins: Molecular docking and molecular dynamics study

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

VMD: Visual molecular dynamics

Atomistic insight and modeled elucidation of conessine towards Pseudomonas aeruginos efflux pump

A Semiempirical Free Energy Force Field with Charge-Based Desolvation

Structural basis for the high all-trans-retinaldehyde reductase activity of the tumor marker AKR1B10

Peŕez-Sala, D. Molecular Interactions and Implications of Aldose Reductase Inhibition by PGA1 and Clinically Used Prostaglandins

Exploring the Interactional Details between Aldose Reductase (AKR1B1) and 3-Mercapto-5H-1,2,4-triazino[5,6-b]indole-5-acetic Acid through Molecular Dynamics Simulations

Discovery of new selective human aldose reductase inhibitors through virtual screening multiple binding pocket conformations

Gaussian 16, revision C.01; Gaussian, Inc., 2016. (26) Ryde, U. Molecular dynamic simulations of alcohol dehydrogenase with varying coordination number of the catalytic zinc ion

Calculating structures and free energies of complex molecules: Combining molecular mechanics and continuum models

Free energy calculations: applications to chemical and biochemical phenomena

Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations

Free energy analysis of the conformational preferences of A and B forms of DNA in solution

Why are lopinavir and ritonavir effective against the newly emerged Coronavirus 2019? Atomistic insights into the inhibitory mechanisms

Aldose Reductase Inhibitors as Potential Therapeutic Drugs of Diabetic Complications

Expression of the aldo-ketoreductases AKR1B1 and AKR1B10 in human cancers

Polyol pathway links glucose metabolism to the aggressiveness of cancer cells

AKR1B1 promotes basal-like breast cancer progression by a positive feedback loop that activates the EMT program

Inhibition of aldose reductase prevents growth factor-induced G1-S phase transition through the AKT/phosphoinositide 3-kinase/E2F-1 pathway in human colon cancer cells

Effects of lignans extracted from Eucommia ulmoides and aldose reductase inhibitor epalrestat on hypertensive vascular remodeling

Inhibitory effect of total lignan from Fructus arctii on aldose reductase

Bioacitve Constituents from Chinese Natural Medicines. XV. Inhibitory Effect on Aldose Reductase and Structures of Saussureosides A and B from Saussurea medusa

Catalytic mechanism of the primary human prostaglandin F2α synthase, aldo-keto reductase 1B1-Prostaglandin D2 synthase activity in the absence of NADP(H)

Human aldose reductase and human small intestine aldose reductase are efficient retinal reductases: Consequences for retinoid metabolism

Inhibition of human aldose reductase-like protein (AKR1B10) by α-and γ-mangostins, major components of pericarps of mangosteen

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

Selective inhibition of the tumor marker aldo-keto reductase family member 1B10 by oleanolic acid

Molecular modeling studies of the binding modes of aldose reductase inhibitors at the active site of human aldose reductase

Potency of inhibitors depends upon the accessibility of their aromatic rings within the hydrophobic specificity pocket: a novel avenue for future aldose reductase inhibitor design

Longterm effects of tolrestat on symptomatic diabetic sensory polyneuropathy

Characterization of a novel aldose reductase inhibitor, FR74366, and its effects on diabetic cataract and neuropathy in the rat

phthalazineacetic acid (zopolrestat) and congeners

spiro hydantoin aldose reductase inhibitor

Hydantoin bioisosteres. in vivo active spiro hydroxy acetic acid aldose reductase inhibitors

Stratified analyses for selecting appropriate target patients with diabetic peripheral neuropathy for long-term treatment with an aldose reductase inhibitor, epalrestat

New insights into the implication of epigenetic alterations in the EMT of triple negative breast cancer

Aldose reductase inhibitors in the treatment of diabetic peripheral neuropathy: a review