Molecular basis of rutin inhibition of protein disulfide isomerase (PDI) by combined in silico and experimental methods

Protein disulfide isomerase (PDI) is a founding member of the thiol isomerase family, and is recently found to play critical roles in thrombus formation. The development of effective PDI inhibitors is of great significance, and attracts strong interest. We previously showed that rutin bound directly to PDI and inhibited PDI activities, leading to the suppression of platelet aggregation and fibrin generation in a mouse model. A close analog of rutin, isoquercetin, is currently in advanced phase clinical trials. However, the molecular interaction between rutin and PDI is unknown and is difficult to study by X-ray crystallography due to the weak interaction. Here, we generated a molecular model of PDI:rutin complex by molecular docking and thorough molecular dynamics (MD) simulations. We then validated the complex model through a number of different experimental methods. We mutated the key residues predicted by the model and analyzed the mutants by an optimized isothermal titration calorimetry (ITC) method and a functional assay (insulin reduction assay). The results consistently showed that the PDI residues H354, L355 and E359 are important in the binding of rutin. These residues are next to the canonical major substrate binding site of the b′ domain, and were not conserved across the members of thiol isomerases, explaining the specificity of rutin for PDI among vascular thiol isomerases. Furthermore, the inhibitory activities of three rutin analogues were evaluated using an insulin reduction assay. The results supported that the second sugar ring at the side chain of rutin was not necessary for the binding to PDI. Together, this work provides the structural basis for the inhibitory mechanism of rutin to PDI, and offers a promising strategy for the design of new generation inhibitors with higher binding affinity to PDI for therapeutic applications.


Introduction
Protein disulde isomerase (PDI; EC 5.3.4.1) is a member of the thioredoxin superfamily of redox proteins, 1 which contains more than 20 different members. 2 PDI normally exists and functions in endoplasmic reticulum (ER), 3 but under certain circumstances can escape from ER to other locations, e.g. plasma membrane. 4 PDI can be either an oxidoreductase or a molecular chaperone. 5 PDI facilitates the folding of nascent proteins or misfolded proteins to ensure proper disulde bond formation through its oxidase and reductase activities. 6 PDI is composed by multiple domains, including a short signal peptide at the amino terminal, four thioredoxin-like domains (namely a, b, b 0 , and a 0 ), one linker between b 0 and a 0 domains (x linker), a short acidic C-terminal extension and the ER retention signal sequence (KDEL) at the carboxyl terminal (Fig. 1A). 7 The a and a 0 domains of PDI are considered to be the catalytic domains containing the conserved active CXXC motifs. 8 The b and b 0 domains of PDI are thought to function as substrate recognition and binding sites. 9 PDI is involved in a number of physiological and pathological processes. 10 It has been demonstrated that PDI is overexpressed in diseases with ER stress. 11,12 PDI is a promising target for cancer therapy. 13 PDI has also been shown to play a vital role in thrombus formation. 14 PDI is rapidly secreted from platelets and endothelial cells and is critical for platelet accumulation and brin generation at sites of injury in a mouse thrombosis model, demonstrating that PDI is also a feasible target for antithrombotic treatment. 15 There has been great interest in developing PDI inhibitors in the past few decades. A few types of PDI inhibitors were reported previously, which include antibodies, peptides, and small organic molecules. [16][17][18][19] Many of the small molecular inhibitors covalently bind to the catalytic cysteine residues within the active sites and are irreversible inhibitors of PDI. It is also a challenge to develop specic small-molecule inhibitors of PDI as the majority of the inhibitors act broadly on other thiol isomerases. Rutin (Fig. 1E) was found to be a specic and reversible inhibitor of PDI. 20 Rutin selectively bound and inhibited PDI activities in a dose-dependent way with an IC 50 of at micromolar range. 20,21 Although the interaction between rutin and PDI is weaker than other small-molecule inhibitors, for example, juniferdin that inhibits PDI at 156 nM, 22 the in vivo inhibitory effect of rutin on PDI was well documented in a number of studies. Inhibition of PDI by rutin blocked both Fig. 1 Overviews of PDI proteins and small-molecule inhibitors mentioned in this paper. (A) Domain organization diagram of PDI and MBP-fused PDI proteins. PDI contains four thioredoxin-like domains, labeled as a, a 0 , b and b 0 from the amino terminal to the carboxyl terminal. There are additional x linker between b 0 and a 0 domains and short C-terminal extension after a 0 domain. In addition, a and a 0 domains have a conserved sequence, termed CXXC motif. And normally, there are signal peptide sequence at the N-terminus and KDEL at the C-terminus of PDI. For MBPfused PDI, MBP was fused at the N-terminus of PDI. We constructed b 0 x and b 0 xa 0 recombinant proteins of human and murine PDI, namely MBP-PDI-b 0 x, MBP-PDI-b 0 xa 0 , MBP-mPDI-b 0 x and MBP-mPDI-b 0 xa 0 , respectively. MBP-fused PDI proteins were isolated by affinity chromatography with the His-tag purification resin and further purified by gel filtration chromatography. (B) Cartoon representation of the crystal structure of human PDI (abb 0 xa 0 , PDB ID: 4EL1). The a, b, b 0 and a 0 domains were colored red, cyan, green and magenta, respectively. The x linker was displayed in blue. The four cysteine residues of the two active sites (CXXC motifs) were represented by spheres. (C) The reduced SDS-PAGE analysis of the purified MBP-fused PDI proteins. The molecular weight of MBP-PDI-b 0 x, MBP-PDI-b 0 xa 0 , MBP-mPDI-b 0 x and MBP-mPDI-b 0 xa 0 are 59.4 kDa, 71.8 kDa, 59.2 kDa and 71.4 kDa, respectively. (D) Analysis of the purified wt MBP-PDI-b 0 xa 0 by gel filtration chromatography on Superdex75. (E) and (F), chemical structure of rutin (E) and isoquercetin (F). Rutin and isoquercetin has typical structure of flavonoid with two benzene rings (A and B rings) and a pyranic ring (C ring). Isoquercetin is glycosylated at position 3, while rutin core is attached with a rhamnose and a glucose at position 3.
platelet aggregation and brin generation in vivo at a dose of 0.5 mg kg À1 in mouse models. Furthermore, the inhibitory activities of rutin on PDI can be rapidly reversed in vivo by the infusion of PDI b 0 x fragment in mouse model. 21 Recently, isoquercetin, a close analogue of rutin with only one glucose moiety in its structure (Fig. 1F), was tested on human through oral administration (1 gram) and was demonstrated to reduce plasma thrombin generation activity by about half through inhibition of PDI. 23 Even though rutin appears to be a potential antithrombotic agent, it still has low potency and low bioavailability. In addition, the detailed molecular interaction between PDI and rutin is still largely unknown, and is difficult to study by X-ray crystallography, partly due to the weak nature of this interaction. We previously localized the binding site of rutin to the b 0 xa 0 fragment of PDI. 21 In this study, we used a combined molecular simulation and molecular biology method to tackle this problem by rst generating a molecular model of rutin binding to PDI through thorough computational study, and validating the model using rigorous experimental assays. The ndings obtained demonstrated the feasibility of this approach, and provided a structural template to design better and more potent inhibitors of PDI.

Materials
Mouse P4HB Gene ORF cDNA clone in cloning vector purchased from SinoBiological Inc. Constructed plasmid vector pMCSG9 (containing His-tag, TEV site and MBP) and Escherichia coli (E. coli) BL21 (DE3) cells were from our laboratory. Rutin was purchased from Sangon Biotech (Shanghai) Co., Ltd. Kaempferitrin was purchased from Chengdu Micxy Chemical Co., Ltd. Tiliroside was purchased from Shanghai Yuanye Bio-Technology Co., Ltd. 2 0 -O-Galloylhyperin was purchased from Chengdu Herbpurify Co., Ltd. Insulin was purchased from Shanghai Yeasen Bio-Technology Co., Ltd. We analysed the purity of these small-molecule chemicals by HPLC (Fig. S1 †). All these chemical reagents have high purity for scientic study.

Expression and purication of recombinant proteins
MBP fusion proteins expressed in E. coli strain BL21 (DE3) aer they were induced by isopropyl-b-D-thiogalactoside (IPTG). All recombinant proteins were captured by affinity chromatography with His-tag purication resin (Chelating Sepharose Fast Flow from GE Healthcare), followed by gel ltration chromatography through Superdex75 HR 10/300 column (GE Healthcare) equilibrated with protein buffer (pH 7.4) (20 mM Tris containing 150 mM NaCl, 1 mM ethylenediaminetetraacetic acid (EDTA), 1 mM b-mercaptoethanol (b-ME) and 5% glycerin). The puried recombinant proteins were stored at À80 C.

Fluorescence-based direct binding assay
The relevant recombinant proteins were incubated with rutin in black 384-well plate at room temperature for 30 min before measurements. The concentrations of proteins were always kept at 18 mM and rutin was 55 mM. The total assay volume was lled to 100 mL using the protein buffer (pH 7.4). The uorescence emission spectrum of every mixture were measured with the excitation at 430 nm and the emission at 460 nm at room temperature on a BioTek Synergy microplate reader with a sensitivity of 100 or 120 and measurements of per data with 100 times.

Insulin reduction assay
The insulin reduction assay was performed to assess the enzymatic activity of PDI and its MBP-fused recombinant proteins. PDI can reduce insulin in the presence of dithiothreitol (DTT) along with increase in turbidity of the reaction mixture. We thus monitored spectrophotometrically the aggregation of reduced insulin through the increase of turbidity. The reductase activity was measured in 100 mL in presence of 600 mM insulin, 1.5 mM PDI and recombinant PDI proteins, 1 mM DTT, and with protein buffer (pH 7.4), at room temperature over 100 min by absorption at 650 nm. For inhibition assay, 100 mM of rutin, kaempferitrin, tiliroside and 2 0 -O-galloylhyperin or control buffer was added to the reaction system.

Isothermal titration calorimetry (ITC)
Microcalorimetry titration of rutin to PDI recombinant and mutant proteins used a MicroCal ITC200 microcalorimeter (GE Healthcare). PDI recombinant proteins, which concentration ranges from 25 mM to 50 mM, were added into the cell (about 300 mL). Rutin, with concentration ranging from 500 mM to 1 mM, was injected from syringe to titrate the protein in cell (almost 40 mL). The experiment was performed at 25 C with a equilibration time of 10 min, the reference power was set to 6.25 mcal s À1 . The process including 25 injections, each of 1, 2 or 2.5 mL rutin, in a duration of 3 s, and the interval time between two injections was 180 s. Data was analyzed with MicroCal Origin (MicroCal).

Molecular docking
The molecular docking calculations of rutin onto PDI were performed with AutoDock 4.2 package 24 using the Lamarckian genetic algorithm. 25 The initial 3D structures of PDI was obtained from the PDB database (PDB ID: 4EL1, 26 resolution at 2.88 A). The geometry of rutin was optimized using GAMESS 27 with the B3LYP functional 28 and 6-31G(*) basis set. 29 The grid map, centered at PDI was calculated using 126 Â 126 Â 126 points with 0.375 A spacing. Two-thousand docking runs were carried out. The binding energy and cluster analyses were conducted using the default parameters implemented in Auto-DockTools 1.5. 24 The best docking result was dened by the model that has the best binding affinity within the dominant cluster.

Molecular dynamics (MD) simulation
The model of PDI in complex with rutin obtained from docking calculations was inserted into box of explicit water and 150 mM NaCl (corresponding to the salt concentration used in our binding assays). The protonation states of residues were assigned according to the corresponding pK a values calculated by using the H++ webserver. 30 Sodium ions were added to counterbalance the charge of the complex. The system contained 86 132 atoms in total. It was underwent MD simulations using the GROMACS 4.6.5 code 31 on in-house GPU clusters. The force elds of PDI and rutin were set to AMBER ff99SB-ILDN force eld 32-34 and general AMBER force eld (GAFF), 35 respectively. TheÅqvist potential 36 and TIP3P model 37 were used for the ions and for the water molecules, respectively. All bond lengths were constrained by LINCS algorithm. 38 Periodic boundary conditions were applied. Electrostatic interactions were calculated using the Particle Mesh-Ewald (PME) method, 39 and van der Waals and Coulomb interactions were truncated at 10 A. The system underwent 1000 steps of steepest-descent energy minimization with 1000 kJ mol À1 A À2 harmonic position restraints on the protein, followed by 2500 steps of steepest-descent and 2500 steps of conjugate-gradient minimization without restraints. The system was then gradually heated from 0 K up to 298 K in 20 steps of 2 ns. Aer that, two parallel productive MD simulations of 6000 ns and 4000 ns long with different microscopic initial conditions carried out in the NPT ensemble, respectively. The most representative structure was identied by the cluster analysis 31 over the equilibrated trajectories (aer 1000 ns). To assess the convergence of the simulated trajectory, we considered both the time evolution of backbone atoms' root-mean-square deviation (RMSD, see Fig. S2 †) of PDI and their projections on the top essential dynamical spaces obtained from a standard covariance analysis. Following Hess's criterion, 40 these projections were next compared with those expected for a random reference. The observed negligible overlap (i.e. cosine content close to 0, data not shown) conrms a posteriori adequate sampling of the complex conformations around the equilibrium position. Hydrogen bond interactions were dened to be present if the atomic distance between the acceptor and donor atoms is below 3.5 A and the angle among the hydrogen-donor-acceptor atoms is below 30 degree. Hydrophobic interactions were dened to be present if the center-of-mass distance between side chains are smaller than 4.5 A. 41

Rutin bound and inhibited MBP-fused recombinant human and mouse PDI proteins
We rst tried X-ray crystallography to dene the binding of PDI with rutin. With this original structural biology goal, we generated recombinant proteins of PDI and its fragments fused with maltose binding protein (MBP). MBP is a 43 kDa protein commonly used as fusion tag to increase the crystallizability of fused protein. In addition, MBP tag augments water solubility of target proteins and facilitates expression of recombinant proteins. 42 For human PDI, we constructed MBP-PDI-b 0 x and MBP-PDI-b 0 xa 0 recombinant proteins. All proteins were expressed in E. coli and puried by Ni-NTA chromatography followed by gel ltration chromatography ( Fig. 1C and D).
We performed additional experiments to study if these MBPfused recombinant proteins can bind to rutin and have reductase functions. We previously developed a direct binding assay of rutin to PDI. 21 Rutin emits uorescence when excited at 430 nm alone and has a maximal emission peak at 550 nm, while human PDI has little uorescence excited at 430 nm alone. With addition of PDI into rutin, the uorescence intensity of the ligand signicantly increased and the wavelength of peak emission appeared blue-shied to 530 nm ( Fig. 2A). Our results for MBP-fusion proteins showed that rutin indeed bound to MBP-PDI-b 0 x and MBP-PDI-b 0 xa 0 fusion proteins (Fig. 2B). The IC 50 of these MBP-fusion proteins are similar with that of the full-length human PDI (Fig. S3 †). As a negative control to rule out the possibility of MBP interaction with rutin, we measured the uorescence of rutin in the presence of three other different MBP-fusion proteins, and found no change of uorescence at all (Fig. S4 †). We also found the reductase activities of MBP-PDI-b 0 xa 0 through insulin reduction assays (Fig. 3C). Taken together, these results demonstrated that MBP fusion did not perturb PDI function and rutin binding, which is not surprising as MBP is a compact protein and usually does not affect the activity of target proteins.
Next, we studied the molecular interaction of PDI from murine origin (mPDI) with rutin. Although the effect of rutin was demonstrated in a murine thrombosis model, 14,20 it is unknown whether rutin can inhibit mPDI directly. We rst generated MBP-mPDI-b 0 x and MBP-mPDI-b 0 xa 0 recombinant proteins and puried them to high homogeneity (Fig. 1C). We then carried out direct binding assays, and found that rutin bound to mPDI (Fig. 2C). The binding affinities of murine proteins (IC 50 of about 30 mM) were similar to those of human counterparts (Fig. S3 †). The reductase activity of murine MBP-mPDI-b 0 xa 0 was also comparable to that of the human counterparts ( Fig. 3C and E). Thus, mPDI possessed similar characteristics of human PDI in interacting with rutin.
Furthermore, we constructed a number of MBP-PDI-b 0 xa 0 proteins with site-directed mutations on K328, H354, L355 and E359 (see below for the design of these mutants). These mutants were also expressed and puried by the methods described above (Fig. S5 †).

A molecular model of PDI:rutin complex revealed by in silico simulations
We originally wanted to crystalize the MBP-fused PDI or mPDI or their fragments in complex with rutin in order to understand the molecular interaction of rutin with PDI. Unfortunately, we did not obtain any diffracting crystals despite thousands of crystallization trials, presumably due to the weak interaction between PDI and rutin.
We then resorted to in silico methods to generate a molecular model through molecular docking calculations and molecular dynamics (MD) simulations based on the available crystal structure of PDI (PDB ID: 4EL1 26 ). Both the current studies using MBP-fused PDI and previous direct binding assays 21 showed that the b 0 domain of PDI is the major binding site of rutin. Our docking and MD studies also revealed that rutin indeed bound to PDI mainly at the b 0 domain, but also contacted with the x linker (Fig. 4A). In our structural model, the quercetin moiety of rutin formed stable hydrogen bonds with the backbone oxygen atom of E303 from the helix a3 of b 0 domain and with the side-chains of H354 and E359 from x linker (Fig. 4B, quercetin site). The stabilities of such intermolecular interactions were supported by the analysis of the occupancies of formation and average distances of hydrogen bonds across MD simulations (Table S1 †), as well as by the mutagenesis studies (see following section). Moreover, one additional stable hydrogen bond was observed between the hydroxyl group at position 2 of the glucose of rutin and the backbone oxygen atom of L355 ( Fig. 4B and Table S1 †).
In contrast to the stable interactions established between PDI and the quercetin and glucose moieties of rutin, the putative hydrogen bonds contributed by the rhamnose moiety of rutin and PDI residues (K328, E330 and K352, named rhamnose site) were non-optimal and transient across the MD simulations, as indicated by the relatively low occupancies of formation of hydrogen bonds and large standard deviations of the average distances between the donor and acceptor atoms (Table S1 †). To further assess the effect of the rhamnose moiety of rutin on binding to PDI, we performed in silico-based ADME analysis on both rutin and isoquercetin (Table S2 †). The results showed that the omission of this ring may not apparently change the ADME properties of rutin. Therefore, the stability analysis of the putative interactions between PDI and rutin suggested that quercetin glycosylated at 3 position with glucose may be the major contributor to the binding affinity, while the rhamnose moiety may not contribute signicantly.

Validation of PDI:rutin complex model by site-directed mutagenesis
To validate the binding site of rutin on PDI predicted by in silico simulations, we constructed a list of mutants of human MBP-PDI-b 0 xa 0 on the basis of the structural model (Fig. 4): K328E, H354A, L355R, E359A, and E359K. Among them, H354A and E359A mutations were intended to verify the quercetin binding site on PDI. Other mutations, including the introduction of charges to hydrophobic residue or charge-reversal mutation, drastically changed amino acid properties to induce stronger effect of the mutations.
We rst used the uorescence-based binding assays to measure the direct binding of rutin to the mutants. When added to rutin, mutant proteins like K328E and H354A showed similar uorescence proles. The uorescence intensities of K328E and H354A mutants were also similar to that of the wildtype (wt) MBP-PDI-b 0 xa 0 protein (Fig. 2D). These results showed that our uorescence-based direct binding assay was not able to distinguish the differences introduced by the mutations, reecting the low sensitivity of this direct binding assay.
Next, we resorted to a more sensitive method, isothermal titration calorimetry (ITC), to measure the direct binding. Due to the weak interaction, this method requires large amounts of high-purity protein samples (about 240 mL at 20-50 mM) to be titrated by the highest concentration of rutin (its water solubility of about 200 mM), so that the release of heat as a result of complex formation reached plateau phase rapidly and was enough to be measured. A problem we encountered is that the starting phase of the S-shape curve of heat measurement cannot always be resolved using the typical titration program of ITC, thus affecting the accuracy of dissociation constant (K d ) measurement from the S-shape ITC curve. This was due to the use of equal volume of small molecule titrant ($2.5 mL) used in the typical ITC titration. Therefore, we established an optimized ITC titration scheme by using variable drop sizes, which started with small drops (1 mL at rst ten times of injection) and ended with large drops (2 mL). We thus obtained S-shape ITC curves with well-dened upper and lower platforms (Fig. 5 and S6 †), which allowed the accurate measurements of K d values.
At rst, we carried out this improved ITC method to measure the K d of rutin to wt MBP-PDI-b 0 xa 0 . In two independent ITC experiments, the K d was found to be 2.39 mM and 2.31 mM, respectively (Fig. 5). This result validated that our improved ITC procedure is quite sensitive and demonstrated that the ITC results were reproducible. We then applied this ITC method for rutin binding to wt MBP-PDI-b 0 xa 0 and its mutants (Table 1 and Fig. S6 †). The results showed that the stoichiometric ratio for rutin binding to wt MBP-PDI-b 0 xa 0 protein was about 1 : 1, which was consistent with our previously reported measurements of the full-length human PDI 21 and further conrmed that the binding site of rutin was indeed within the b 0 xa 0 domains of PDI. In general, the effect of mutation on rutin binding to PDI was modest (

Measurements of reductase activities of the site-directed mutants further validated the PDI:rutin complex model
We further used the insulin reduction assay to measure the reductase activities of these site-directed MBP-PDI-b 0 xa 0 mutants. This assay has typically good reproducibility as shown by a control experiment with three repetitions (Fig. S7 †). We found the PDI mutants maintained reductase activities in the presence of dithiothreitol (DTT). The H354A mutation lowered the reductase activity of wt MBP-PDI-b 0 xa 0 (Fig. 3G) while other mutants had similar reductase activities with that of wt MBP-PDI-b 0 xa 0 (Fig. 3C, F and H-J). Addition of rutin inhibited the reductase activities of both wt MBP-PDI-b 0 xa 0 and the mutants ( Fig. 3C and F-J). As control, we showed that rutin and other other rutin analogues (see section below) did not perturb the insulin reductase assay (Fig. S8 †). We quantied the inhibition of rutin on these MBP-PDI-b 0 xa 0 mutants by measuring the increase of turbidity (OD 650 at 100 min) upon addition of rutin (Fig. 3K). We found that, compared with the wt MBP-PDI-b 0 xa 0 , the inhibitory efficiency of rutin on H354A, E359A and E359K decreased about 20% while L355R showed modest decrease. Together, these results are in agreement with our in silico simulations and ITC assays, and supported that the residues H354 and E359 are important in the binding of PDI to rutin.

Inhibitory activities of rutin analogues against PDI
We identied a series of rutin analogues and purchased three of them including kaempferitrin, tiliroside and 2 0 -O-galloylhyperin where the second sugar ring is different from rutin. Their inhibitory activities against the full-length human PDI and PDI fragment (wt MBP-PDI-b 0 xa 0 ) were also measured by insulin reduction assays. The reductase activities of both PDI and MBP-PDI-b 0 xa 0 were inhibited by these analogues as observed in the assays of rutin (Fig. 6). The inhibitory activity of kaempferitrin to the full-length human PDI was smaller than that of rutin to PDI. The results indicated that the second sugar ring at position 3 of rutin did not contribute to the binding to PDI. This was in agreement with the results obtained in the previous study of the structure-activity relationship between avonols and their inhibitory potency against PDI, where the IC 50 values of the inhibitors with or without the second sugar ring at 3 position were similar. 20 This was also consistent with our in silico simulations and ITC assays, both of which suggested that the residue K328 involved in the interaction with the second sugar ring at position 3 of rutin did not function in the binding of rutin to PDI.

Protein exibility of PDI
PDI is a dynamic protein with two redox-dependent conformations clearly identied by crystal structural studies of either human 26,43 or yeast PDI. 44 In the crystal structure of the oxidized fungal PDI b 0 xa 0 bound to its substrate, a-synuclein, the protein took an open conformation with an exposed major substrate binding site in b 0 domain bound to the substrate. At the reduced condition, there was no direct binding of the substrate, a-synuclein, presumably due to the closed conformation. 45 Recent MD studies conrmed such redox-regulation of PDI conformations, but revealed a much larger degree of exibility. 46,47 Smallangle X-ray scattering (SAXS) experiments, 48,49 including ours, 21 were also carried out to study the solution conformation of PDI. All these results suggest a hypothesis of redox-regulated substrate binding of PDI where PDI is in an open conformation in the oxidized state (in the presence of GSSG) to expose its hydrophobic substrate binding pocket, and becomes compact in the reduced form to expel the substrate away from the substrate binding site. In our current study, we used the reducing condition in all assays because the insulin reduction assay was carried out in the presence of DTT, and the recombinant protein contains two free cysteines and needs reducing condition to maintain its activities. We additionally did an experiment to study if the redox state will affect the binding site of rutin. We found that the oxidative condition (GSSG) gave a higher uorescence signal compared to the reducing condition (DTT) in the direct binding assay of rutin to wild type PDI (Fig. S9 †). For PDI H354 mutant, such redox-induced uorescence enhancement has a magnitude identical to wt PDI (Fig. S9 †). This indicated that the binding sites of rutin on PDI most likely not affected by the redox state of PDI. This is understandable because the binding site locates far away from the free cysteine residues, which are most like one key sensing elements for the redox state (GSSG/ DTT).
There are only limited numbers of crystal structures of the full-length human PDI, presumably reecting its dynamic nature and the difficulty of PDI crystallization. The resolved crystal structures showed that the four domains of human PDI arranged in a horseshoe shape, wherein the a and a 0 domains faced each other at the open end of the horseshoe while b and b 0 domains formed its closed end (Fig. 1B). The structure of PDI in oxidized and reduced states also showed that the b 0 xa 0 region of PDI was mainly responsible for the redox-regulation of its conformations. 26 A tryptophan residue of the x linker (20 amino acids long) bound to the substrate binding site, 50 demonstrating one example of conformational exibility induced by the x linker.
We narrowed down the rutin binding site to the b 0 xa 0 fragment of PDI, and used this fragment and its fusion with MBP for the co-crystallization with rutin in order to reduce conformational entropy penalty and increase the crystallizability. However, we still were unable to obtain diffracting crystals of rutin in complex with MBP-fused recombinant PDI fragments despite extensive efforts, likely due to the weak nature of the complex. Thus, we resorted to thorough computational methods to predict the binding model of rutin bound to PDI. In light of relatively low binding affinity of rutin to PDI, we use two carefully designed assays to validate the binding model. Both assays consistently showed that the H354A mutation lowered the binding affinity of PDI fragment to rutin and the activity of PDI fragment. We did not combine the single mutants together for further validation. A better approach will be nuclear magnetic resonance (NMR) for further validation of the PDI binding residues of rutin, but such method requires major effort.

Structural basis of PDI interacting with rutin
All four domains of PDI participate in the substrate recognition, 51-53 and the non-catalytic domains, especially the b 0 domain, are the major substrate binding site. 54 The PDI b 0 domain contains a characteristic thioredoxin-like fold containing ve bsheet strands and four a-helices. The major substrate binding site was formed between the a1 and a3 helices and constituted of conserved amino acids (F240, F249, I301, F304, F305, L320 and M324) (Fig. 4B). This binding site is largely hydrophobic in nature, and appears not to be the potential binding site of rutin because of the amphipathic nature of rutin structure.
We previously demonstrated that the b 0 x domain of PDI is the primary binding site of rutin. 21 Our current in silico study and site-directed mutagenesis showed that the binding site was uniquely located at the interface between b 0 and x linkers. The residues H354, L355 and E359 from the b 0 x domain were key for rutin binding to PDI (Table 1 and Fig. 4B). Notably, these key residues were not conserved across the extracellular thiol isomerases, indicating a candidate interface to target PDI specically (Fig. 4C). This was consistent with previous experimental data that rutin did not inhibit other thiol isomerases, including ERp5, ERp57, ERp72 or thioredoxin. 20 Our model also showed that the second sugar ring at 3 position of rutin was not important for binding to PDI. Three analogues of rutin which mainly differ on the second sugar ring at 3 position showed similar inhibitory activities against the full-length human PDI and wt MBP-PDI-b 0 xa 0 (Fig. 6). In agreement with this, mutation on its binding residue (K328E) had no effect on the binding of PDI (Table 1 and Fig. 3F).
The structural model of rutin binding of PDI may drive a new strategy for structural optimization of rutin-based PDI inhibitors by increasing their binding affinities to PDI. For example, it is valuable to perform thorough scaffold hopping on the quercetin and glucose moieties of rutin. Structure-based pharmacophore strategies can be also used to identify or design novel moieties targeting PDI rhamnose site with enhanced stabilities.

Potential mechanism of rutin inhibition on the redox activities of PDI
Rutin binds and inhibits PDI selectively and reversibly. According to our current PDI:rutin model, the binding sites of rutin on PDI were primarily located at the interface of b 0 domain and x linker of PDI. Rutin bound near the major substrate binding site of PDI, but did not directly bind to the major substrate binding site (Fig. 4). However, rutin was not too far away from the opening of substrate binding site (4.3 A to residue F304). Thus, rutin may perturb the conformation of the major substrate binding site, and inhibit substrate binding to PDI. On the other hand, rutin binding at the x linker may hinder the exibility of PDI, which is critical for the catalytic activity of PDI. In addition, our MD simulations also showed that rutin induced a closed conformation of PDI, blocking the accessibility of the substrate binding site (Fig. S10 †). Collectively, it appears that it is the allosteric effect of rutin that is mainly responsible for inhibiting the redox activities of PDI.

Conclusion
Rutin demonstrated effective anticoagulant activity in human by inhibiting PDI activity. However, the rutin does not have good bioavailability and high dose is required. In fact, rutin inhibits PDI activity with an IC 50 at micromolar range. Identi-cation of molecular interaction between rutin and PDI is critical for further improvement of pharmacological prole of rutin. We developed a strategy to study such weak interaction combining computational and experimental methods. The results demonstrate the feasibility of the current strategy. In addition, rutin shows a unique binding site right next to the canonical major substrate binding site located at the b 0 domain. Thus, this work provides the structural basis for the inhibitory mechanism of rutin to PDI, and will be useful for further development of potent inhibitors with higher binding affinity to PDI for therapeutic applications.

Conflicts of interest
There are no conicts to declare.