Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2016; 7: 11248.
Published online 2016 April 14. doi:  10.1038/ncomms11248
PMCID: PMC4834639

Steroid binding to Autotaxin links bile salts and lysophosphatidic acid signalling


Autotaxin (ATX) generates the lipid mediator lysophosphatidic acid (LPA). ATX-LPA signalling is involved in multiple biological and pathophysiological processes, including vasculogenesis, fibrosis, cholestatic pruritus and tumour progression. ATX has a tripartite active site, combining a hydrophilic groove, a hydrophobic lipid-binding pocket and a tunnel of unclear function. We present crystal structures of rat ATX bound to 7α-hydroxycholesterol and the bile salt tauroursodeoxycholate (TUDCA), showing how the tunnel selectively binds steroids. A structure of ATX simultaneously harbouring TUDCA in the tunnel and LPA in the pocket, together with kinetic analysis, reveals that bile salts act as partial non-competitive inhibitors of ATX, thereby attenuating LPA receptor activation. This unexpected interplay between ATX-LPA signalling and select steroids, notably natural bile salts, provides a molecular basis for the emerging association of ATX with disorders associated with increased circulating levels of bile salts. Furthermore, our findings suggest potential clinical implications in the use of steroid drugs.

Autotaxin (ATX) is a secreted lysophospholipase D (lysoPLD) that produces the bioactive lipid lysophosphatidic acid (LPA) from extracellular lysophospholipids, predominantly lysophosphatidylcholine (LPC)1,2. ATX is a unique member of the ectonucleotide pyrophosphatase/phosphodiesterase (ENPP) family of enzymes that hydrolyse phosphodiester bonds in various substrates, including nucleoside triphosphates, lysophospholipids and choline phosphate esters3.

The bioactive product of ATX, LPA, acts on six distinct G protein-coupled receptors (LPA1–6) that activate multiple signalling pathways4. The biological outcome of ATX-LPA signalling is remarkably diverse, depending on LPA receptor expression patterns and tissue context, and includes the stimulation of cell migration, proliferation and survival4,5,6. In pathophysiology, ATX-LPA signalling has been implicated in tumour progression, pulmonary fibrosis, neuropathic pain, cardiovascular disease and cholestatic pruritus, making the ATX-LPA signalling axis an attractive therapeutic target7.

Crystal structures of ATX8,9 (ENPP2) revealed a catalytic domain with a bimetallic active site adjacent to a catalytic threonine. Substrate binding takes place in a shallow hydrophilic groove that accommodates nucleotides as well as the glycerol moiety of lysophospholipids, and in a deep hydrophobic pocket binding the lysophospholipid acyl chain. In addition, the catalytic domain together with the first of two somatomedin-β (SMB) domains, form a tunnel adjacent to the active site (Fig. 1a,b).

Figure 1
The tunnel of ATX.

The function of this unique tunnel, which is absent in other ENPP family members like ENPP1 (refs 10, 11), remains enigmatic. One hypothesis based on structural and mutagenesis data suggests that the tunnel serves to deliver LPA to its cognate G protein-coupled receptors (GPCRs)9. Recent structural analysis shows that access to the LPA1-binding pocket occurs from the extracellular space12, in contrast to the proposed access route of the related lipid mediator sphingosine 1-phosphate (S1P) to its receptor S1P1 (ref. 13). ATX isoforms have affinity for cell-surface integrins8 and heparan sulphate proteoglycans14. In addition, kinetic analysis revealed that the half-time of the ATX:LPA complex is several minutes15, allowing ATX:LPA to diffuse and possibly recognize cell surface receptors, bringing the LPA product closer to the cell membrane. Collectively, these observations place ATX at the centre of an intrinsic network of interactions that could explain its extensive roles in diverse physiological and disease processes. Until now, no physiological regulators of ATX activity have been identified.

Our on-going studies of ATX structure and catalytic mechanisms reveals that recombinant ATX contained steroids bind to the tunnel. We discover that binding of at least one class of steroids, bile salts, modulates ATX activity. Our findings reveal an unexpected novel interplay between LPA and bile salts, and possibly other steroids, and offer testable hypotheses on the function of the ATX-LPA signalling axis in various physiological contexts.


ATX binds steroids in its tunnel

During the course of our structural investigation of the ATX catalytic mechanism, we obtained a 1.6-Å resolution structure of ATX (Table 1). Following structure refinement, some residual density in the tunnel became apparent (Fig. 2a,b). The density was compatible with the four-ring system bearing two axially oriented methyl groups, characteristic of steroid moieties. Indeed, modelling of cholesterol in the density yielded a nearly perfect fit (Fig. 2c). To explain remaining difference density, adjacent to the C7 position of the B ring, we modelled 7α-hydroxycholesterol (7HCS, Fig. 2d,e).

Figure 2
7α-OH steroids are bound in the tunnel of ATX.
Table 1
Crystallographic data.

Notably, when we determined the first ATX structure8 (2XR9) at 2.05 Å resolution, we have noted that ‘residual electron density occupies this tunnel, probably indicating the presence of a molecule that we were unable to model (possibly a lipid acquired from the growth medium)'. When we re-examined the maps corresponding to that structure, with the hindsight that the residual electron density could be 7HCS, it was indeed trivial to fit a 7HCS moiety to that residual density in the tunnel (Supplementary Fig. 1). This indicated that 7HCS has always been present in our protein preparations, and it is certainly a bound ligand in the 2XR9 ATX structure, which we were unable to identify at the time. Indeed, the presence of cholesterol (~15%), hydroxysterols (~80%) and 7HCS specifically (~75% of total hydroxysterols and ~60% of the total) at near stoichiometric levels to ATX, in several preparations, was confirmed by mass spectrometry (Table 2 and Supplementary Fig. 2). The mass spectrometry data also suggest that the predominant 7HCS species is the α epimer. The 7HCS structure and the revised 2XR9 structure strongly suggest that 7HCS is most likely the only bound molecule, albeit the exact nature of the aliphatic chain could not be uniquely identified.

Table 2
Mass spectrometry analysis of ATX.

The steroid moiety is oriented away from the active site with the 3-OH group facing the solvent on the opposing surface of the molecule. The steroid rings are well packed, with the conserved Trp254, Trp260 and Phe274 likely participating in π–π interactions. The 7-OH group is clearly only in a single anomeric configuration (α) as suggested by the mass spectrometry data; the preference is explained by a hydrogen bond to the Νζ atom of Trp260. The selectivity towards the 7-OH substituted steroid is notable, as a Protein Data Bank (PDB) search suggests this is the only example of specific hydroxysterol recognition by a carrier protein.

Bile salts bind ATX and inhibit LPA production

As sterols were not present during ATX purification, while ATX was selective for a 7α-substituted steroid presumably acquired from the cells and/or cell culture medium, we went on to investigate how distinct steroids may modulate ATX activity (Fig. 3). Hydroxysterols like 7HCS or 7α,25-dihydroxycholesterol, as well as various other steroids (testosterone, dexamethasone, prednisolone), had no effect on ATX-mediated LPC hydrolysis (Fig. 3b and Supplementary Fig. 3). Strikingly, however, the conjugated bile salt tauroursodeoxycholate (TUDCA) inhibited ATX activity with an apparent IC50 of 10 μM (Fig. 3c). Structure–activity relationship studies using bile salts with different substitutions of the steroid rings and different conjugates of the acidic tail, revealed that: (i) bile salts with the C7 hydroxyl in the α-configuration (taurochenodeoxycholate) or in the β-configuration (TUDCA) inhibit ATX activity with similar potency (Fig. 3 c,d); (ii) the presence of a hydroxyl in C12 of the C ring, either in the presence of the C7 hydroxyl (that is, taurocholate) or in the absence of the C7 hydroxyl (that is, taurodeoxycholate) abrogate the inhibitory activity of TUDCA (Fig. 3e,f); and (iii) the conjugation of the acidic tail did not significantly change potency, as shown by using glycochenodeoxycholate (which is conjugated to glycine instead of taurine) or ursodeoxycholate (lacking conjugation; Fig. 3g,h).

Figure 3
Activity of ATX in the presence of selected steroids.

To understand the selectivity of bile salts on ATX activity, we co-crystallized ATX with TUDCA and determined the structure to 2.0 Å resolution (Table 1). Following molecular replacement and model adjustment, there was a clear difference electron density in the tunnel, markedly different from the density for 7HCS (Fig. 4a,b). Owing to the reduction of the C5–C6 double bond of the steroid B-ring, the four-ring system in bile salts is no longer planar, but adopts an L-shaped conformation (Fig. 4b). In comparison with the 7HCS placement, the steroid moiety of TUDCA has ‘slid' inside the tunnel by about 2.5 Å towards the active site (Fig. 4c,d). This movement ‘aligns' the D-ring better with Phe274, likely resulting to more favourable π–π interactions. The Νζ of Trp260 now forms a hydrogen bond with the 3-OH group of the A-ring, instead of the 7-OH that has moved away. In addition, the conserved His251 flips towards the tunnel exit, acting as a ‘lid'. Overall, the ring placement indicates a stronger interaction of TUDCA compared with 7HCS in the tunnel. Given that the tunnel may serve as an LPA exit route9, this observation might provide a mechanism for the inhibitory activity of TUDCA. Unlike the steroid moiety, the acidic tail is not very well resolved in the electron density. Further refinement allowed modelling of three district conformations of the TUDCA acidic tail (Supplementary Fig. 4a,b).

Figure 4
The bile salt TUDCA binds the tunnel of ATX.

Strikingly, structural alignment of the TUDCA structure with the 7HCS structure or the native structure (2XR9), showed that the SMB domains move ‘outwards' to accommodate TUDCA in the tunnel (Supplementary Fig. 4c); a movement of the ATX SMB domains relative to the catalytic core has not been observed before. Thus, TUDCA binding induces a conformation change, which may have implications for SMB domain interactions and binding site dynamics.

The observed multiple conformations of the acidic tail, contrasting with the well-defined binding of the bile salt steroid moiety, agrees well with the structure–activity data, and prompted us to examine further the inhibition mechanism of bile salts towards ATX.

Bile salts are non-competitive inhibitors of LPC hydrolysis

Testing the inhibitory effect of different TUDCA concentrations against increasing concentrations of LPC, showed a reduction in Vmax with no significant change in the kM of the reaction (Fig. 5a); consistently, kinetic modelling confirmed that a competitive model has a probability of <0.01% to be correct (see the Methods for details and equations). As TUDCA is a partial inhibitor of ATX (Fig. 3b), we used a partial mixed inhibition model to explain the kinetic data (see the Methods for details). Kinetic modelling suggested that while TUDCA has an appreciable [var kappa]I of 9±3 μM, even at saturating concentrations ~40% of ATX activity towards LPC is retained. Crucially, mathematical modelling indicated that TUDCA is a noncompetitive inhibitor, implying that LPC or LPA can co-exist in a ternary complex with ATX.

Figure 5
TUDCA acts as a non-competitive inhibitor of LPC hydrolysis.

To test this hypothesis, we crystallized a complex of ATX, TUDCA and LPA and we determined a structure of this ternary complex to 1.8 Å resolution (Table 1). In this complex, TUDCA is well defined in density (Supplementary Fig. 5) and the steroid moiety lies in the tunnel (Fig. 5b–d), but the conjugated tail is in a single well-defined conformation pointing away from the active site (Fig. 5e). LPA(18:1) binds in a similar pose as shown earlier for mouse ATX9 (3NKP); whereas the glycerol moiety is well defined in density, the aliphatic tail is less well defined (Supplementary Fig. 5) and likely exists in multiple conformations. This structure reveals the molecular basis of TUDCA action as a non-competitive inhibitor, and confirms that tunnel blockade by the bile salt, rather than active site de-stabilization, underlies its inhibitory activity.

These observations reinforce the hypothesis that the tunnel serves as a product exit route9; partial inhibition by bile salts further suggests that the tunnel is an alternate but not obligatory exit route for LPA release. It is important to note that, although TUDCA acts as an inhibitor over a wide variety of LPC substrates with different aliphatic tails (Supplementary Fig. 6a–d), TUDCA did not inhibit the hydrolysis of nucleotide substrates (Supplementary Fig. 6e). This further supports the hypothesis that TUDCA acts by blocking LPA exit through the tunnel. Moreover, TUDCA completely inhibited hydrolysis of the fluorescent lysolipid analogue FS3 (ref. 16; KI=1.4 μM, Supplementary Fig. 6f) in a competitive manner (see the Methods for details). Although it is unknown how ATX binds FS3, it seems likely that one of the two FS3 acyl chains with the attached fluorophore is accommodated in the ATX tunnel, making TUDCA-mediated inhibition for FS3 competitive in nature. We conclude that the tunnel is an allosteric site that binds bile salts to reduce LPA production rates, most likely by slowing down product release.

Bile salts attenuate ATX function ex vivo and in situ

To validate that bile salts act as allosteric inhibitors in a physiological setting, we examined their effect on modulating LPA production in human serum. Consistent with our in vitro findings, TUDCA inhibits endogenous ATX activity in human serum ex vivo (Fig. 6a) with an apparent IC50 of ~30 μM. Also in serum, TUDCA behaves as a partial inhibitor, leaving ~40% residual activity. Next, we examined the ability of bile salts to modulate ATX-LPA signalling, by examining their effect on ATX-dependent LPA receptor activation in cultured cells (in situ), using agonist-induced LPA1 receptor internalization as readout17,18. As shown in Fig. 6b,c, at low LPC substrate levels, TUDCA reduces ATX-mediated LPAR1 endocytosis in HeLa cells, consistent with reduced and/or delayed LPA production and delivery.

Figure 6
Physiological effects of bile salts on human serum and LPA receptor activation.


Our finding that bile salts such as TUDCA and glycochenodeoxycholate act as allosteric inhibitors of ATX lysophospholipase D activity, suggests that ATX may interact with steroids in vivo. Under pathological conditions, notably cholestasis, bile salt concentrations can reach >100 μM in the blood stream19. Thus, the observed binding constants of bile salts to ATX (~10 μM) suggest that bile salts are likely to modulate ATX activity under natural conditions. As the lifetime of the ATX:LPA complex is relatively long (tens of seconds15) allowing the complex to diffuse over a distance of ~65 μm before dissociating, bile salt-mediated stabilization of the ATX:LPA complex could enable further spreading of the LPA signal and allow LPA-loaded ATX to interaction with integrin receptors8 or heparan proteoglycans14 in a more productive manner.

Modulation of ATX activity by bile salts is an intriguing finding given that cholestatic itch is accompanied by a significant increase in serum ATX activity19. In addition, intradermal injection of LPA in mice leads to a transient but significant increase in scratch behaviour, strongly suggesting that LPA produced by ATX initiates or potentiates itch19. As cholestasis is characterized by accumulation of bile salts in the circulation, this is at first sight contradictory to our findings. However, delaying LPA release by prolonging the lifetime of the ATX:LPA complex, might have stimulatory effects in select cells, such as sensory neurons involved in itch perception.

In the context of cholestatic disease, it is noteworthy that UDCA is often used for the treatment of women with cholestasis of pregnancy (ICP), which develops in ~1% of pregnant women during the third trimester of pregnancy, to relieve itch20. As ATX activity in the serum of these patients is significantly increased and serves as a diagnostic marker21, our results suggest that the antipruritogenic effect of UDCA may be attributed, at least in part, to inhibition of plasma ATX activity. This notion is also supported by our finding that TUDCA inhibits serum ATX ex vivo (Fig. 6a).

In addition, UDCA is the only FDA-approved drug for the treatment of primary biliary cholangitis, a disorder that leads to progressive cholestasis and is often associated with chronic itch22. Moreover, TUDCA has been shown to alleviate symptoms of arthritis in rats23; interestingly, ATX has been implicated in the pathogenesis of modelled arthritis24. Anti-inflammatory effects of UDCA and its conjugated derivatives have also been recognized in other non-hepatic disorders in which ATX could play a role, such as neurite growth impairment25, colitis26 and neovascularization27. Furthermore, UDCA has been explored as a drug to alleviate symptoms of several neurodegenerative diseases including ALS, Parkinson's disease and Huntington's28; although molecular mechanisms that implicate an apoptotic function and ER stress relief have been suggested29,30, the mechanism of action of UDCA and its derivatives remains an active field of research. The hypothesis that the therapeutic effect of ursodeoxycholates is partially achieved also by unwittingly interfering with ATX function and LPA signalling warrants further investigation.

Various steroidal and non-steroidal agonists of the bile salt receptors TGR5 and FXR are being investigated as candidate drugs for the treatment of chronic liver disease, hepatocellular cancer and extrahepatic inflammatory and metabolic diseases31. More specifically developments have targeted obesity, type 2 diabetes, hypertriglyceridaemia and atherosclerosis32. In addition, obetichocolic acid is being evaluated in clinical trials as a treatment for non-cirrhotic, non-alcoholic steatohepatitis33. The interaction of ATX with bile salts, even regardless of their ability to modulate LPA signalling, might be important for delivering such drugs to their cognate receptors, much as already proposed for ATX and LPA7,9.

Oxysterols, which clearly bind ATX but did not appear to change ATX activity in vitro, using LPC(18:1) as substrate, have emerged as signalling molecules linked to breast cancer pathophysiology34. We cannot exclude the possibility that certain oxysterols (or other steroids) might be affecting the processing of other physiologically relevant unsaturated LPC species with longer acyl chains (that are not available for testing).

LPA binds to albumin. How the protein-bound LPA pool, the bio-active LPA pool (which can also be the ATX-bound LPA) and the LPC substrate pool are correlated is poorly understood. Recently, Hla and co-workers have shown that altering the balance of apolipoprotein M and free S1P influences lymphopoiesis and neuroinflammation35. By analogy LPA, and steroids might affect the balance between plasma protein and ATX-bound pools of LPA, selectively promoting or inhibiting LPA signalling.

Finally, steroids specific for the ATX tunnel may offer an alternative strategy for developing ATX inhibitors as partial inhibitors might prove of importance. In addition, blocking the ATX tunnel by molecules that slow down LPA release, might even promote LPA signalling, by allowing more efficient delivery of LPA to its receptors.

In conclusion, our study reveals a previously unknown interplay between ATX-mediated LPA production and bile salts acting as allosteric modulators of ATX activity, with potential physiological and clinical implications.


ATX production and crystallization

ATX from rat (rATX) and human (hATX) proteins were produced in HEK 293 Flp-In cells as described previously36. The high-resolution crystals that contained 7HCS, rATX (3.5 mg ml−1) was first pre-incubated with 23.5 mg ml−1 heparin sulphate; 1 μl of this solution was then mixed with 1 μl of reservoir solution containing 20% (w/v) PEG 3350, 0.1 M NH4I and 0.3 M NaSCN. After 1 week, the first crystals formed, to finally obtain a crystal of 0.15 × 0.02 × 0.02 mm3. This unusually large crystal diffracted to high resolution; although we did not investigate in detail why, these were some of our best-diffracting crystals at the time. We believe that we obtained them through the combination of a good protein preparation (the quality of protein preps often vary), optimizing little practical details of ATX crystallization (as this is done manually there is variation), and reducing the pH close to 6.0 due to the acidic heparin (which we measured at the end in the reservoir solution). For co-crystallization of rATX with TUDCA, TUDCA was first mixed in 10 M excess with rATX (3.5 mg ml−1) and incubated for 30 min. For co-crystallization of rATX with TUDCA and LPA(18:1) (Avanti Polar Lipids), rATX (3.5 mg ml−1) was first mixed with a 3 M excess of LPA(18:1) and 10 M excess of TUDCA, and incubated for 30 min. In both cases, the best crystals were obtained at room temperature (293 K) in a 24-well optimization screen (18–20% PEG 3350, 0.1–0.4 M NaSCN and 0.1–0.4 M NH4I) by mixing 1 μl of protein:inhibitor solution and 1 μl of reservoir solution. All crystals were vitrified in cryoprotectant, which consisted of reservoir solution with the addition of 20% (v/v) glycerol.

Crystallographic data and methods

Structural studies were performed with rATX. The X-ray diffraction data for the complexes were collected at ESRF on beamlines ID29 (7HCS, TUDCA:LPA) and ID23 (TUDCA), respectively, at 100 K. The crystallographic diffraction data were recorded on a PILATUS detector to resolution of 1.60 Å for the 7HCS complex, 2.20 Å for the TUDCA complex and 1.80 Å for the TUDCA:LPA complex. All data were processed and integrated with XDS37 and scaled with AIMLESS38. The structures were determined by molecular replacement using PHASER39 with the structure of ATX (PDB entry 2XR9 (ref. 8)) as the search model. Model building and subsequent refinement were performed iteratively with COOT40, REFMAC5 (ref. 41) and PDB_REDO42.

The 7HCS structure was first extensively refined, before we attempted modelling of cholesterol and then 7HCS, as described earlier. For the TUDCA structure, TUDCA was initially modelled in a single conformation in one of the two ATX molecules in the asymmetric unit; after refinement, a second molecule was modelled in the ATX molecule in the second asymmetric unit. After extensive refinement, the TUDCA moiety in the first asymmetric unit was modelled in three different conformations to account for residual difference density. Following that, some neighbouring residues were also modelled in multiple conformations. In the TUDCA:LPA complex, the TUDCA and LPA moieties were modelled directly in the difference density map after initial refinement. Structure validation was carried out by MolProbity, and shows that 97% of all residues where within the Ramachandran plot favoured regions for all three structures, while one residue was an outlier; the clash score and the MolProbity43 score for all three structures belong to the 100th percentile. Additional crystallographic parameters and model quality indicators are summarized in Table 1. The structure models and experimental diffraction data were deposited at the PDB under codes 5DLT for 7HCS, 5DLV for the TUDCA complex and 5DLW for the TUDCA–LPA complex.

Structural alignment of the TUDCA structure with the 7HCS structure or the 2XR9 structure performed by RAPIDO44, suggested that an optimal alignment should include two rigid bodies: the first consisting of most of the SMB domains (residues 56–100 and 110–120) and a loop of the phosphodiesterase domain (275–284), and the second of the rest of the protein. Both rigid bodies are similar between the two structures (root-mean-square (RMS) deviation 0.5 Å). However, superposing on the second, larger, rigid body shows that the SMB domains have an RMS deviation of 2.3 Å to each other, showing that they moved ‘outwards' to accommodate TUDCA in the tunnel.

Mass spectrometry analysis

Recombinant hATX preparations or buffer controls were extracted using acidified organic solvents as reported previously for lipid extractions from cells and tissues45 with inclusion of deuterium-labelled cholesterol of oxysterol internal standards as detailed below. Samples were evaporated to dryness and reconstituted in methanol for mass spectrometry analysis.

Cholesterol and hydroxysterols were analysed using an AB Sciex 4000-Qtrap hybrid linear ion trap triple quadrupole mass spectrometer in multiple reaction monitoring mode connected to a Shimadzu HPLC. Cholesterol was analysed in positive multiple reaction monitoring mode using a Machery Nagel Nucleodur C8 gravity column, 5 μm, 125 × 2 mm2 with a flow rate of 0.5 ml min−1. Solvent A consists of 75/25 of methanol/water with 0.5% formic acid and 0.1% 5 mM ammonium formate and solvent B consists of 99/1 of methanol/water with 0.5% formic acid and 0.1% 5 mM ammonium formate. The cholesterol measurements were conducted exactly as reported previously with quantification accomplished by reference to a deuterated internal standard45. Oxysterols were measured after dimethylglycine (DMG) derivatization using an Agilent Zorbax Eclipse XDB-C18 column, 5 μm, 4.6 × 150 mm2 with a flow rate of 0.5 ml min−1. Solvent A was 85% methanol with 5 mM ammonium acetate and solvent B consisted of 80:20:2 of chloroform/water/50 mM ammonium acetate. In brief, we obtained authentic standards for 25-, 24-, 27, 7-α and 4-β OH cholesterol and deuterated 25-OH cholesterol and prepared DMG derivatives exactly as previously described46. We obtained product ion spectra from the doubly and singly charged DMG derivatives of these compounds that were essentially identical to those reported previously with some differences in product ion abundances that are likely attributable to differences between our mass spectrometer system and that used by Jiang et al.46. The DMG derivatives of these authentic standards were used in conjunction with known information reported46 to develop mass spectrometry methods that measure product ions that are common to these derivatives and ions that are diagnostic of the particular oxysterols that were then used in selected reaction monitoring mode HPLC-coupled tandem mass spectrometry assays using chromatography methods that could separate these standards. Quantification of oxysterols was accomplished by reference to the deuterated 25-OH cholesterol internal standard as explained below.

Cholesterol was detected in all ATX preparations analysed. Extracted ion chromatograms for two precursor product ion pairs, which were used to detect and quantify the deuterated cholesterol internal standard and ATX-associated cholesterol, are shown in Supplementary Fig. 2, while Table 2 shows representative data from one ATX preparation.

As reported by Jiang et al.46, bis-DMG derivatives of hydroxylated cholesterol generate doubly charged molecular ions (m/z=287.3), which fragment to produce protonated DMG (m/z=104.1) and production of singly charged mono-DMG derivative ions (m/z=367.2). Further fragmentation/rearrangements of the m/z=367.2 ion and its products through processes that are described in detail by Jiang et al.46 generates ions of m/z=232.2 and m/z=255.2 that permit discrimination of 7-OH cholesterol from other oxysterols. Accordingly, we established a selected ion-monitoring mode HPLC LC MS/MS assay that monitored these ion pairs. Supplementary Fig. 2 shows extracted ion chromatograms for these four ion pairs identifying 7-OH cholesterol as the major oxysterols associated with the ATX preparation we analysed. Jiang et al.46 also reported relative differences in signal between the 387.2 and 255.2 ions that discriminate the 7-α and 7-β OH cholesterol epimers. The relative intensities of these ions in the data shown in the figure suggest that the major species is the 7-α epimer, as also suggested by the structural data. Quantification was accomplished by comparing the integrated peak area for the 387.2/104.1 ion pair to that of the corresponding ion pair from the deuterated 25-OH cholesterol internal standard. That ATX preparation analysed contained 1.48 pmol of 7-OH cholesterol. Analysis of extracted ion chromatograms for the precursor and product ions that are common to the bis DMG derivatives of hydroxylated cholesterols with comparison to the retention times of suggested that the preparation contained significantly lower levels of 24-, 25- and 27-OH cholesterols (0.47 pmol).

Biochemical assays and modelling of kinetic data

All biochemical studies were performed with hATX. ATX lysoPLD activity was measured by choline release from LPC47. 20 nM ATX (prepared from HEK 293 Flp-In cells) was incubated with 150 μM LPC(18:1) in a final volume of 100 μl buffer containing 50 mM Tris, pH 7.4, 0.01% Triton X-100, 50 mM CaCl2, 1 U ml−1 choline oxidase, 2U ml−1 horseradish peroxidase, 2 mM homovanilic acid (HVA). The relative amount of released choline was measured by HVA fluorescence in a 96-well plate (Corning). Fluorescent intensity was determined at λex/λem=320/450 nm every 30 s for 90 min with a Fluorostar plate reader (BMG Labtech). Data analysis was performed using GraphPad Prism.

Nucleotide phosphodiesterase activity was measured using pNP-TMP (Sigma-Aldrich) as substrate in 50 mM Tris (pH 7.8), 140 mM NaCl, 5 mM KCl, 1 mM CaCl2, 1 mM MgCl2, 1 mg ml−1 fatty acid-free BSA. After addition of recombinant ATX, the amount of liberated para-nitrophenolate (pNP) was determined by the absorbance at 405 nm. The relative amount of released pNP was determined by the absorbance at 405 nM every 30 s for 90 min with a Fluorostar plate reader (BMG Labtech). Data analysis was performed using GraphPad Prism.

ATX activity against the FS3 substrate16 was measured by incubating 20 nM ATX with the indicated range of FS3 in a final volume of 50 μl in 50 mM Tris buffer (pH 7.4) containing 0.2 mg ml−1 fatty acid-free BSA. Fluorescent intensity was determined at λex/λem=485/520 nm every 30 s for 90 min with a Fluorostar plate reader (BMG Labtech).

The experiments for determining IC50 for different bile salt inhibitors were performed by adding the named bile salt and performing serial threefold dilutions; for each bile salt concentration, the velocity of the reaction was monitored as a function of time in the lysoPLD assay as described above. The linear part of the reaction velocity was visually estimated (typically between 10 min and 1 h) and the observed fluorescence signal (described above) as a function of time was modelled for all inhibitor concentrations simultaneously using the formula:

An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m1.jpg

where Ft is the observed fluorescence signal relative to the choline product at each time point (t), F0 is the background fluorescence signal at the start of each measurement, vmax and vmin are the fitted values for the minimum and maximum relative velocity (common for all inhibitor concentrations) and cinh is the corresponding inhibitor concentration for each of the experiments.

At the same time, we have also derived an observed velocity (V) by linear regression of the fluorescence signal for each inhibitor concentration and modelled that as a function of inhibitor concentration using the formula:

An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m2.jpg

for cross-checking the result of equation (1) and for creating the inhibition graphs presented in various figures.

For initial comparison between competitive, uncompetitive and non-competitive inhibition by TUDCA, using the LPC and the FS3 assays, we measured the relative velocity for increasing LPC or FS3 substrate concentrations, for no-inhibitor and for three different inhibitor concentrations. The relative velocity for each data point was measured by linear regression as described above, and a Michaelis–Menten curve was created for each inhibitor concentration. We then used the following equations to describe each mode of inhibition:

An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m3.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m4.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m5.jpg

where V is the observed velocity and cLPC is the corresponding LPC concentration for each data point, and ci is the inhibitor concentration for each curve; KM and Vmax are the Michaelis–Menten constants and the estimated maximum velocity in a relative scale, and Ki is the inhibition constant. Note that in equation (5) Ki denotes the product of Ki (which is very high as the inhibitor actually does not bind alone) and a constant that defines the apparent Ki of the ternary enzyme–substrate–inhibitor complex.

The preferred model was estimated using pairwise comparisons and we chose the best-fitting model utilizing Akaike's Informative Criteria to select the model most likely to have generated the data. For the LPC hydrolysis data, the competitive model had <0.01% probability to be correct when compared with either the noncompetitive or the uncompetitive models, whereas the noncompetitive model showed a moderate preference (75%) over the uncompetitive model (25%). For the FS3 data, the competitive model had a probability of more than 99.99% to be correct compared with other models.

For understanding the inhibition of LPC hydrolysis by TUDCA in further detail, we used a partial mixed inhibition model defined by the equation

An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m6.jpg

where in addition to previous definitions that can be derived from equations (4–6), we define Part as the partiality of the inhibition and α as an estimated constant that determines mechanism: if α=1, the formula becomes identical to equation (4) and the inhibition mode is noncompetitive; if a is very large, then the model approaches a competitive model equation (3); if a is very small but greater than zero, the model approaches an uncompetitive model.

For more robust fitting, we also did an experiment where we varied the inhibitor concentration for a fixed amount of substrate, which we modelled using an equation similar to equation (2):

An external file that holds a picture, illustration, etc.
Object name is ncomms11248-m7.jpg

equations (6) and (7) were fit simultaneously, each in the data set(s) they describe (with a relative scale factor applied to (7) to account for experiment differences), with the parameters that are determined in both experiment (KM, Vmax, Ki, Part) constrained to have identical values between experiments. This procedure allowed a more robust fit to determine all parameters from all available experiments. The value of α refined to 2.3, suggesting a noncompetitive inhibition model, where TUDCA binds to ATX both in the absence and in the presence of LPC (or LPA), albeit with a very modest twofold preference for the substrate-free protein.

All modelling was performed using the regression analysis tools within the Graphpad/Prism software. The values of all derived parameters indicate mean values as obtained by nonlinear regression analysis using the formulas listed above; triplicate data were used as separate observations during the fit; ±values indicate the standard error of the mean.

Effects of bile salts in the ex vivo activity of ATX

An ATX choline release activity assay (see above for details) was done with three different sera from normal pregnant women (informed consent was obtained from all human participants; study protocol approved by the Medical Ethical Committee of the Amsterdam Medical Center). We chose those because ATX activity is relatively high without increased serum levels of bile salts21. In place of purified ATX, 20 μl of serum were added in a total assay volume of 150 μl), aiming to dilute the serum as little as possible. Before starting the assay, the serum was pre-incubated with each indicated TUDCA concentration (1:2 dilution) for 5 min. The indicated TUDCA concentration is the final concentration in the assay.

LPAR internalization assay

HeLa cells (obtained from American Type Culture Collection) were grown in DMEM supplemented with 10% fetal bovine serum at 37 °C under 5% CO2. For LPA1 internalization studies, HeLa cells were transient transfected with GFP-LPA1 plasmid using X-tremeGene 9 (Roche) reagent according to the manufacturer's instructions. At 10 h after transfection, the cells were rinsed with SF-DMEM and incubated in the same medium for 16–24 h before further treatments. Cells were treated with 0.1 μM LPC and 5 nM ATX or with 0.1 μM LPC and 5 nM ATX and 500 μM of TUDCA for 15 min. Coverslips were fixed with 4% paraformaldehyde for 10 min, washed twice with PBS and mounted using Immu-Mount. Slides were examined on a Leica TCS-SP5 confocal microscope (× 62 objective) and analysed using ImageJ software.

Cell culture and reagents

HEK 293 Flp-In (purchased from Invitrogen) and HeLa cells were maintained in DMEM (Gibco, Life Technologies) supplemented with 10% FCS (HyClone, Thermo Scientific) 2 mM glutamine, 100 U ml−1 penicillin and 100 μg ml−1 streptomycin (Gibco, Life Technologies). Choline oxidase, HVA, horseradish peroxidase fatty acid-free BSA, Dexamethason, PEG3350, NaSCN, NH4I, pNP-TMP (p-nitrophenyl thymidine 5′-monophosphate), taurochenodeoxycholic acid, taurodeoxycholic acid, glycochenodeoxycholic acid, ursodeoxycholic acid and LPA were purchased from Sigma-Aldrich. LPC and 7α,25-dihydroxycholesterol from Avanti polar lipids. FS3 (LysoPLD/ATX substrate) from Echelon and TUDCA from Calbiochem.

Additional information

Accession code: The crystal structures have been deposited in the Protein Data Bank with accession codes: 5DLT, 5DLV, 5DLW.

How to cite this article: Keune, W.-J. et al. Steroid binding to Autotaxin links bile salts and lysophosphatidic acid signalling. Nat. Commun. 7:11248 doi: 10.1038/ncomms11248 (2016).

Supplementary Material

Supplementary Information:

Supplementary Figures 1-6


We thank Alexander Fish for help with biophysical experiments and data interpretation. Diffraction experiments were performed on beamlines ID23 and ID29 at the European Synchrotron Radiation Facility (ESRF), Grenoble, France. We are grateful to Eleonore von Castelmur for data collection and initial analysis.


Author contributions W.-J.K. was involved in performing and/or performed biochemical experiments and all structural and cell-based experiments, and prepared figures; J.H. determined the structure that lead to the identification of bound 7a-oxycholesterol and was involved in pilot biochemical experiments; A.K. made the initial discovery for bile acid inhibition of ATX, and R.B. and D.T. demonstrated inhibitory characteristics of various bile salt species; T.H. performed cell-based experiments and crystallization experiments; R.P.J. helped in identifying the nature of the (at the time) unknown ligand in the electron density maps and finalized refinement of all crystal structures; M.S. performed and A.J.M. supervised the cholesterol/oxycholesterol quantification experiments; E.M.-R. performed cell-based experiments; W.H.M. consulted in all stages of this research and co-supervised cell-based studies; R.P.O.E. supervised research that independently discovered the role of bile salts in inhibiting ATX; A.P. supervised research that independently demonstrated binding of steroids and supervised all structural studies, performed mathematical modelling of the kinetic data, participated in structure determination and refinement, and prepared figures; R.P.O.E. and A.P. jointly supervised the bile salt biochemical experiments after joining the lab efforts; A.P. wrote the paper with assistance from W.-J.K., W.H.M. and R.P.O.E. R.P.O.E. and A.P. contributed equally to this publication.


  • Umezu-Goto M. et al. . Autotaxin has lysophospholipase D activity leading to tumor cell growth and motility by lysophosphatidic acid production. J. Cell Biol. 158, 227–233 (2002). [PMC free article] [PubMed]
  • Tokumura A. et al. . Identification of human plasma lysophospholipase D, a lysophosphatidic acid-producing enzyme, as autotaxin, a multifunctional phosphodiesterase. J. Biol. Chem. 277, 39436–39442 (2002). [PubMed]
  • Bollen M., Gijsbers R., Ceulemans H., Stalmans W. & Stefan C. Nucleotide pyrophosphatases/phosphodiesterases on the move. Crit. Rev. Biochem. Mol. Biol. 35, 393–432 (2000). [PubMed]
  • Yung Y. C., Stoddard N. C. & Chun J. LPA receptor signaling: pharmacology, physiology, and pathophysiology. J. Lipid Res. 55, 1192–1214 (2014). [PMC free article] [PubMed]
  • van Meeteren L. A. & Moolenaar W. H. Regulation and biological activities of the autotaxin-LPA axis. Prog. Lipid Res. 46, 145–160 (2007). [PubMed]
  • Choi J. W. et al. . LPA receptors: subtypes and biological actions. Annu. Rev. Pharmacol. Toxicol. 50, 157–186 (2010). [PubMed]
  • Perrakis A. & Moolenaar W. H. Autotaxin: structure-function and signaling. J. Lipid Res. 55, 1010–1018 (2014). [PMC free article] [PubMed]
  • Hausmann J. et al. . Structural basis of substrate discrimination and integrin binding by autotaxin. Nat. Struct. Mol. Biol. 18, 198–204 (2011). [PMC free article] [PubMed]
  • Nishimasu H. et al. . Crystal structure of autotaxin and insight into GPCR activation by lipid mediators. Nat. Struct. Mol. Biol. 18, 205–212 (2011). [PubMed]
  • Jansen S. et al. . Structure of NPP1, an ectonucleotide pyrophosphatase/phosphodiesterase involved in tissue calcification. Struct. Folding Design 20, 1948–1959 (2012). [PubMed]
  • Kato K. et al. . Crystal structure of Enpp1, an extracellular glycoprotein involved in bone mineralization and insulin signaling. Proc. Natl Acad. Sci. USA 109, 16876–16881 (2012). [PubMed]
  • Chrencik J. E. et al. . Crystal structure of antagonist bound human lysophosphatidic acid receptor 1. Cell 161, 1633–1643 (2015). [PMC free article] [PubMed]
  • Hanson M. A. et al. . Crystal structure of a lipid G protein-coupled receptor. Science 335, 851–855 (2012). [PMC free article] [PubMed]
  • Houben A. J. S. et al. . The polybasic insertion in autotaxin α confers specific binding to heparin and cell surface heparan sulfate proteoglycans. J. Biol. Chem. 288, 510–519 (2013). [PMC free article] [PubMed]
  • Saunders L. P. et al. . Kinetic analysis of Autotaxin reveals substrate-specific catalytic pathways and a mechanism for lysophosphatidic acid distribution. J. Biol. Chem. 286, 30141 (2011). [PMC free article] [PubMed]
  • Ferguson C. G. et al. . Fluorogenic phospholipid substrate to detect lysophospholipase D/autotaxin activity. Org. Lett. 8, 2023–2026 (2006). [PMC free article] [PubMed]
  • Lee C.-W., Rivera R., Gardell S., Dubin A. E. & Chun J. GPR92 as a new G12/13- and Gq-coupled lysophosphatidic acid receptor that increases cAMP, LPA5. J. Biol. Chem. 281, 23589–23597 (2006). [PubMed]
  • Murph M. M., Scaccia L. A., Volpicelli L. A. & Radhakrishna H. Agonist-induced endocytosis of lysophosphatidic acid-coupled LPA1/EDG-2 receptors via a dynamin2- and Rab5-dependent pathway. J. Cell Sci. 116, 1969–1980 (2003). [PubMed]
  • Kremer A. E. et al. . Lysophosphatidic acid is a potential mediator of cholestatic pruritus. Gastroenterology 139, 1008–1018 1018.e1 (2010). [PubMed]
  • Bacq Y. et al. . Efficacy of ursodeoxycholic acid in treating intrahepatic cholestasis of pregnancy: a meta-analysis. Gastroenterology 143, 1492–1501 (2012). [PubMed]
  • Kremer A. E. et al. . Autotaxin activity has a high accuracy to diagnose intrahepatic cholestasis of pregnancy. J. Hepatol. 62, 897–904 (2015). [PubMed]
  • Hohenester S., Oude Elferink R. P. J. & Beuers U. Primary biliary cirrhosis. Semin. Immunopathol. 31, 283–307 (2009). [PMC free article] [PubMed]
  • Liu M. et al. . Effects of taurochenodeoxycholic acid on adjuvant arthritis in rats. Int. Immunopharmacol. 11, 2150–2158 (2011). [PubMed]
  • Nikitopoulou I. et al. . Autotaxin expression from synovial fibroblasts is essential for the pathogenesis of modeled arthritis. J. Exp. Med. 209, 925–933 (2012). [PMC free article] [PubMed]
  • Silva S. L. et al. . Neuritic growth impairment and cell death by unconjugated bilirubin is mediated by NO and glutamate, modulated by microglia, and prevented by glycoursodeoxycholic acid and interleukin-10. Neuropharmacology 62, 2398–2408 (2012). [PubMed]
  • Martínez-Moya P. et al. . Dose-dependent antiinflammatory effect of ursodeoxycholic acid in experimental colitis. Int. Immunopharmacol. 15, 372–380 (2013). [PubMed]
  • Cho J. G. et al. . Tauroursodeoxycholic acid, a bile acid, promotes blood vessel repair by recruiting vasculogenic progenitor cells. Stem Cells 33, 792–805 (2015). [PubMed]
  • Vang S., Longley K., Steer C. J. & Low W. C. The unexpected uses of urso- and tauroursodeoxycholic acid in the treatment of non-liver diseases. Glob. Adv. Health Med. 3, 58–69 (2014). [PMC free article] [PubMed]
  • Benz C. et al. . Effect of tauroursodeoxycholic acid on bile-acid-induced apoptosis and cytolysis in rat hepatocytes. J. Hepatol. 28, 99–106 (1998). [PubMed]
  • Ozcan U. et al. . Chemical chaperones reduce ER stress and restore glucose homeostasis in a mouse model of type 2 diabetes. Science 313, 1137–1140 (2006). [PMC free article] [PubMed]
  • Schaap F. G., Trauner M. & Jansen P. L. M. Bile acid receptors as targets for drug development. Nat. Rev. Gastroenterol. Hepatol. 11, 55–67 (2014). [PubMed]
  • Thomas C., Pellicciari R., Pruzanski M., Auwerx J. & Schoonjans K. Targeting bile-acid signalling for metabolic diseases. Nat. Rev. Drug Discov. 7, 678–693 (2008). [PubMed]
  • Neuschwander-Tetri B. A. et al. . Farnesoid × nuclear receptor ligand obeticholic acid for non-cirrhotic, non-alcoholic steatohepatitis (FLINT): a multicentre, randomised, placebo-controlled trial. Lancet 385, 956–965 (2015). [PMC free article] [PubMed]
  • Nelson E. R., Wardell S. E., Jasper J. S. & Park S. 27-Hydroxycholesterol links hypercholesterolemia and breast cancer pathophysiology. Science 342, 1094–1098 (2013). [PMC free article] [PubMed]
  • Blaho V. A. et al. . HDL-bound sphingosine-1-phosphate restrains lymphopoiesis and neuroinflammation. Nature 523, 342–346 (2015). [PMC free article] [PubMed]
  • Hausmann J. et al. . Mammalian cell expression, purification, crystallization and microcrystal data collection of autotaxin/ENPP2, a secreted mammalian glycoprotein. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 66, 1130–1135 (2010). [PMC free article] [PubMed]
  • Kabsch W. XDS. Acta Crystallogr. D Biol. Crystallogr. 66, 125–132 (2010). [PMC free article] [PubMed]
  • Evans P. R. An introduction to data reduction: space-group determination, scaling and intensity statistics. Acta Crystallogr. D Biol. Crystallogr. 67, 282–292 (2011). [PMC free article] [PubMed]
  • McCoy A. Solving structures of protein complexes by molecular replacement with Phaser. Acta Crystallogr. D Biol. Crystallogr. 63, 32–41 (2007). [PubMed]
  • Emsley P., Lohkamp B., Scott W. G. & Cowtan K. Features and development of Coot. Acta Crystallogr. D Biol. Crystallogr. 66, 486–501 (2010). [PMC free article] [PubMed]
  • Murshudov G. N. et al. . REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr. D Biol. Crystallogr. 67, 355–367 (2011). [PMC free article] [PubMed]
  • Joosten R. P., Long F., Murshudov G. N. & Perrakis A. The PDB_REDO server for macromolecular structure model optimization. IUCrJ 1, 213–220 (2014). [PMC free article] [PubMed]
  • Chen V. B. et al. . MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr. D Biol. Crystallogr. 66, 12–21 (2010). [PMC free article] [PubMed]
  • Mosca R. & Schneider T. R. RAPIDO: a web server for the alignment of protein structures in the presence of conformational changes. Nucleic Acids Res. 36, W42–W46 (2008). [PMC free article] [PubMed]
  • Onono F. et al. . Efficient use of exogenous isoprenols for protein isoprenylation by MDA-MB-231 cells is regulated independently of the mevalonate pathway. J. Biol. Chem. 288, 27444–27455 (2013). [PMC free article] [PubMed]
  • Jiang X., Ory D. S. & Han X. Characterization of oxysterols by electrospray ionization tandem mass spectrometry after one-step derivatization with dimethylglycine. Rapid Commun. Mass Spectrom. 21, 141–152 (2007). [PMC free article] [PubMed]
  • Albers H. M. H. G. et al. . Boronic acid-based inhibitor of autotaxin reveals rapid turnover of LPA in the circulation. Proc. Natl Acad. Sci. USA 107, 7257–7262 (2010). [PubMed]
  • Potterton L. et al. . Developments in the CCP4 molecular-graphics project. Acta Crystallogr. D Biol. Crystallogr. 60, 2288–2294 (2004). [PubMed]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group