Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Oligomeric interface modulation causes misregulation of purine 5´-nucleotidase in relapsed leukemia

BMC Biology volume 14, Article number: 91 (2016) Cite this article

Abstract

Background

Relapsed acute lymphoblastic leukemia (ALL) is one of the main causes of mortality in childhood malignancies. Previous genetic studies demonstrated that chemoresistant ALL is driven by activating mutations in NT5C2, the gene encoding cytosolic 5´-nucleotidase (cN-II). However, molecular mechanisms underlying this hyperactivation are still unknown. Here, we present kinetic and structural properties of cN-II variants that represent 75 % of mutated alleles in patients who experience relapsed ALL (R367Q, R238W and L375F).

Results

Enzyme kinetics measurements revealed that the mutants are consitutively active without need for allosteric activators. This shows that hyperactivity is not caused by a direct catalytic effect but rather by misregulation of cN-II. X-ray crystallography combined with mass spectrometry-based techniques demonstrated that this misregulation is driven by structural modulation of the oligomeric interface within the cN-II homotetrameric assembly. These specific conformational changes are shared between the studied variants, despite the relatively random spatial distribution of the mutations.

Conclusions

These findings define a common molecular mechanism for cN-II hyperactivity, which provides a solid basis for targeted therapy of leukemia. Our study highlights the cN-II oligomerization interface as an attractive pharmacological target.

Background

Cytosolic purine 5´-nucleotidase (cN-II; EC 3.1.3.5) dephosphorylates purine nucleotide monophosphates with canonical activity towards inosine monophosphate (IMP) [1]. IMP represents a branch point intermediate linking major routes of purine metabolism, including de novo synthesis, the salvage pathway and oxidative degradation [2]. Therefore, cN-II serves as an essential regulatory element to maintain the balance of the nucleotide pool [3]. Unlike other nucleotidases, it is subject to complex allosteric regulation.

The activity of cN-II is regulated according to the adenylate energy charge of cells [4]. In a low energy state or anoxia, adenosine triphosphate (ATP) is depleted and free phosphate accumulates, leading to inhibition of cN-II. This hinders degradation of IMP and contributes to consequent losses of purine nucleotides via diffusion of dephosphorylated nucleosides.

The cN-II protein assembles into functional homotetramers [5], with each subunit comprising 561 amino acid residues [6]. The tetrameric assembly is formed by two identical dimers in which a dimerization site (interface A) comprises 53 intersubunit contacts. These two dimers are held together by interface B, which contains 28 interacting residues from each subunit [7].

The cN-II structure possesses the overall fold characteristic of the haloacid dehalogenase superfamily: an α/β Rossmann-like “core domain” and a smaller “cap domain.” The active site of cN-II is formed by three highly conserved regions: motif I (D52, D54, T56, and L57), motif II (T249), and motif III (K292, D351, and D359) [8]. To achieve a full catalytic activity, cN-II requires the presence of adenylate compounds, such as ATP, adenosine diphosphate or diadenosine tetraphosphate [1]. Binding of these allosteric activators induces a disorder-order transition of helix A, which spans residues 355–365. This transition places residues in the catalytic pocket into an appropriate orientation to coordinate the catalytically essential magnesium ion [9]. In addition, the cN-II sequence contains an acidic region consisting of 15 glutamates at the C terminus, which has been suggested to be involved in the protein oligomerization [10]. However, its role remains unclear due to a lack of structure-function information for this region, as X-ray crystallography studies have been feasible only with the C-terminally truncated form of cN-II (∆537–561) [7, 9].

Fluctuations in cN-II activity have been observed in various pathologies. The activity is upregulated in Lesch-Nyhan syndrome [11], an inherited disorder of purine metabolism that causes hyperuricemia and neurological symptoms. High levels of cN-II expression have also been correlated with poor prognosis in patients with hematological neoplasias and solid tumors, regardless of the therapeutic strategy [12]. In addition, cN-II activity has been implicated in decreased therapeutic efficacy of nucleoside analogs used for the treatment of cancer and viral diseases [13].

Whole-exome sequencing studies have recently revealed that somatic missense mutations in the gene encoding cN-II (NT5C2) are associated with relapsed acute lymphoblastic leukemia (ALL) [14, 15]. The reported mutations lead to biosynthesis of a hyperactive cN-II enzyme that boosts the viability of cancer cells in the presence of purine nucleoside analogs. These findings demonstrate that cN-II is involved in chemoresistance of leukemia cells and suggest that inhibition of hyperactive mutants might become an efficient strategy for increasing the success rates of the current chemotherapy approaches [16]. However, the molecular mechanism underlying the hyperactivity of the cN-II variants remains unknown.

In this study, we describe the biochemical and structural behavior of the three most common mutants: R367Q, R238W, and L375F. Together, these mutants represent approximately 75 % of cN-II variants observed in patients with relapsed ALL [14, 15, 17, 18]. Our work provides structural insight into the impaired allosteric regulation of ALL-specific cN-II mutants and forms a solid basis for development of novel antileukemic therapies.

Results

ALL-specific mutations affect regulation of cN-II

To investigate the enzymological aspects of cN-II hyperactivity, we examined kinetic properties of the most frequently occurring mutants. Kinetic parameters were determined in the presence or absence of the allosteric activator ATP. The concentration of ATP (3 mM) used in the enzyme assay reflected the physiological level of this ligand [1]. The analysis of wild-type cN-II under variable concentrations of ATP revealed that the enzyme is fully active under such conditions (Additional file 1).

Examination of the full-length wild-type cN-II revealed very low activity without ATP, as was demonstrated by the large value (≈30 mM) of the Michaelis constant (K M). Upon addition of ATP, the enzyme is stimulated, and the K M value is decreased by 12-fold while V max is not altered. The previous kinetic studies reported that the ATP binding affected either both K M and V max [1] or just V max [19]. The divergent results may be accounted to different procedures during the protein preparation as well as to a diverse biological source of the enzyme. We used human recombinant protein expressed in bacteria, while the other studies used material isolated from human placenta or calf thymus. Nevertheless, our kinetic data are consistent with a role of cN-II activity in maintaining homeostasis of purine nucleotides [20].

Analysis of the three full-length mutant enzymes revealed that the variants exhibited different Michaelis constants but similar turnover numbers (Table 1). When compared to the wild type, all studied mutants had lower K M in the absence of ATP, indicating that the hyperactive variants are catalytically efficient even in the absence of physiological activators. Moreover, the mutants exhibited an abnormal response to the presence of 3 mM ATP in the reaction mixture. The addition of ATP led to a large decrease in K M (12-fold) for the wild-type enzyme, compared to only a moderate decrease (2-fold) for the R367Q. In contrast, the R238W and L375F mutants were stimulated by ATP to a higher extent than the wild-type enzyme (23-fold and 76-fold decreases in K M, respectively). Our data suggest that the R238W and L375F mutations stimulate the catalytic activity in synergy with the action of the allosteric activator, whereas the R367Q variant adopts a fully active conformation even in the absence of ATP.

Table 1 Kinetic parameters of cN-II enzymes in the presence and absence of ATP
Full size table

In summary, our comprehensive enzyme kinetics study reveals that all hyperactive variants have increased catalytic activity, even in the absence of an allosteric effector. Our results also suggest that the mutations do not have a direct catalytic effect but rather impair the allosteric regulation of the enzyme activity.

Mutants form tetramers with altered thermostability

As a next step, we studied the effects of the activating mutations on the cN-II conformational properties. Because these mutations are distributed close to the intersubunit contact areas in the three-dimensional (3D) structure of wild-type cN-II (Fig. 1a), we investigated whether they affect conformational stability or tetramer formation. Using differential scanning fluorimetry (DSF), we observed increased thermostability of the R367Q mutant compared to the wild type, while the R238W and L375F mutants exhibited decreased melting temperature (T m) values (Fig. 1b). To further examine the thermostability of cN-II variants, we analyzed their catalytic activity under varying temperature (Additional file 2). The wild-type protein and the R367Q mutant displayed similar behavior, exhibiting the highest activity at 50 °C. Nevertheless, the R367Q mutant had a moderately increased relative residual activity in the 45–55 °C range. In contrast, the R238W and L375F mutants exhibited maximal activity at only 45 °C, showing lower catalytic efficiency at higher temperatures. These data are consistent with the DSF observations and confirm the effect of the mutations on the conformational stability of the cN-II variants. Our results demonstrate that the mutations cause an increase (R367Q) or decrease (R238W and L375F) in the structural stability of the cN-II protein.

Fig. 1
figure 1

Oligomerization and thermostability of ALL-specific cN-II mutants. a Distribution of the identified mutations within the 3D structure of cN-II (Protein Data Bank [PDB] ID 2XCX; missing loops were modeled using the ModLoop server). The mutants studied in this work (R367Q, R238W and L375F) are highlighted using a boldface font and thick lines. Each subunit is represented in a different color. b Thermostability of cN-II proteins studied by differential scanning fluorimetry. Each sample was measured twice in two independent measurements. c Size exclusion chromatography of full-length cN-II protein. Elution volumes together with estimated molecular weight are summarized in the table (e). d Blue native polyacrylamide gel electrophoresis (BN-PAGE) of the studied full-length proteins. The protein samples (5 μg) were loaded as follows: M molecular weight standard, S β-amylase from sweet potato (200 kDa), 1 wild type, 2 R367Q, 3 R238W, 4 L375F. e Table summarizing data for the cN-II variants obtained from size exclusion chromatography and small-angle X-ray scattering (SAXS). Analysis of proteins with 1 mg/ml concentration is shown. The values of Rg and MWGuinier were calculated by Guinier approximation; Dmax and Porod volume were determined using fitting of distance distribution functions

Full size image

Next, we used DSF to monitor the thermostability of cN-II proteins under varying concentrations of ATP. Our results showed that T m values of both the wild-type enzyme and the mutants increased in the presence of ATP by more than 5 °C (Additional file 3). Such a high degree of protein stabilization in an ATP concentration-dependent manner strongly suggests that the effects are caused by the binding event. To determine the binding affinity of the cN-II variants towards ATP, we tested the feasibility of isothermal calorimetry and microscale thermophoresis. However, we could not obtain dissociation constants due to a rather weak binding of the activating compounds to the cN-II proteins with apparent K d in the low millimolar range [1]. As a surrogate, we used saturation transfer difference nuclear magnetic resonance (STD NMR) spectroscopy, which is a useful technique for probing of the low-affinity ligand binding [21]. We acquired the STD spectra of the variants in the presence of variable concentrations of ATP and specifically followed the enhancement in the intensity of the ATP signal. We observed a significant increase in the STD enhancement of the ATP signals for all three studied mutants relative to that of the wild-type protein at all tested ligand concentrations (Additional file 4). These experiments clearly showed that the altered allosteric regulation of the mutants is accompanied by a higher binding affinity towards the activator.

To assess the quaternary structures of the cN-II variants (full-length as well as truncated), we used size exclusion chromatography (SEC) and blue native polyacrylamide gel electrophoresis (BN-PAGE). We observed that the wild type and all the mutants migrated as biomolecules of 250–300 kDa (Fig. 1c–e), corresponding to tetrameric assemblies. This finding was further confirmed by small-angle X-ray scattering (SAXS) measurements of the cN-II protein samples at various concentrations. This analysis was feasible only for truncated enzymes due to insufficient solubility of the full-length variants. The wild-type and mutant proteins consistently exhibited radii of gyration of ~5.0 nm and molecular weights of ~260 kDa as determined by Guinier approximation, which was in good agreement with the calculation of the P (r) function. In addition, similar Porod volume values of ~500 Å3 were obtained for all analyzed cN-II variants (Fig. 1e; representative data for the wild type are shown in Additional file 5). In summary, the results of our experiments indicate that the mutations moderately alter the protein stability but do not affect the formation of homotetramers.

Mutation-induced structural changes are found at the oligomeric interface

To gain detailed structural insight into the misregulated cN-II mutants, we determined their 3D structures using X-ray crystallography. We analyzed C-terminally truncated variants, as previously described for wild-type cN-II [7, 9]. Our kinetic analysis, together with DSF and SEC measurements, showed that the full-length and C-terminally truncated enzymes are essentially identical (Additional files 6, 7, and 8). Therefore, the truncated variants can be considered suitable models to provide relevant structural data.

Crystals of all three cN-II mutants belonged to the I222 space group with one molecule per asymmetric unit. We solved their structures by the difference Fourier technique using the isostructural wild-type structure [7]. Overall, the spatial arrangements of the mutants were almost identical to that of the wild type. The root-mean-square deviation (RMSD) values for superposition of 465 Cα atoms of individual mutants with the wild type were approximately 0.5 Å, indicating that the mutations do not induce large conformational changes in the protein. Nevertheless, detailed analysis revealed significant perturbations within the amino acid substitution microenvironment in each mutant (Fig. 2).

Fig. 2
figure 2

Local conformational changes in the cN-II mutants revealed by X-ray crystallography. Wild-type protein subunits are colored in light/dark gray, mutant subunits are colored in yellow or orange, and relevant residues are represented as sticks. Mutated residues are highlighted with boldface font. a The R367Q mutation disrupts the local hydrogen bonding network. b The R238W substitution triggers changes at the oligomeric interface. The region spanning residues 385–410 adopts a helical arrangement in the wild type, while it is unstructured in the mutant. Missing loops were modeled using the ModLoop server (highlighted in blue and magenta for wild type and mutant, respectively). c The L375F mutation alters local intersubunit contacts due to steric hindrance of the more bulky phenylalanine residue

Full size image

The R367Q mutation is located in beta-strand 14 (amino acid residues 367–371) in the proximity of helix A, which is a crucial regulatory segment of cN-II. The replacement of the positively charged arginine residue by a glutamine led to loss of residue 367 hydrogen bonding interactions with both the carbonyl group of S40 and the carboxyl group of D459. Consequently, D459 formed a compensatory polar contact with K361, a residue located in helix A, as illustrated in Fig. 2a. This interaction likely stabilizes the helix A region, which directly stimulates the catalytic activity of mutant cN-II.

The R238W mutation lies in alpha-helical region 9 (amino acid residues 231–242), which forms intersubunit contacts with a region spanning residues 385–433 of the adjacent subunit. We found structural differences between wild type and the R238W mutant mainly in residues 397–405. This segment was unstructured in the R238W mutant but adopted a helical arrangement in the wild-type protein (Fig. 2b). The helical structure was stabilized by its interactions with residues 474–480 of the adjacent subunit, and these contacts were found only in the wild-type structure. In addition, in the R238W mutant, the W238 side chain mutant was reoriented to avoid steric hindrance with F473, inducing changes within the interhelical loop at the N-terminal region (amino acid residues 3–30). Overall, the R238W mutation causes local changes, mainly at the oligomeric interface of cN-II.

The L375F mutation is located in alpha-helix 13 (amino acid residues 375–381), which interacts with a 3/10 helix (amino acid residues 485–487) in the neighboring subunit. Introduction of F375 induced a closer contact between the mutated residue and H486 of the adjacent subunit, which was accompanied by a slight displacement of the proximal L379 side chain (Fig. 2c). Consequently, this led to rearrangement of other surrounding interactions at the oligomeric interface; i.e., the NE1 atom of W382 formed a hydrogen bond with the carbonyl group of H486 in the mutant and with the D487 carbonyl in the wild type. Our data show that the L375F mutation causes local changes at intersubunit contacts.

To study distant structural effects of the mutations, we thoroughly analyzed superpositions of the hyperactive variants’ structures with the previously published structure of the wild-type enzyme apo form (PDB ID 2XCX) [9]. This revealed that major changes are consistently distributed across the oligomeric interface in the cN-II mutants (RMSD > 0.8 Å; Fig. 3). This observation indicates that local effects of the mutations propagate through the entire structure of cN-II and alter both contact areas responsible for oligomerization: interfaces A and B. The most significant changes at interface A (RMSD ≈ 1.2 Å) occurred within the regions between residues 385–433 and 470–486, which are in mutual contact and form a peripheral part of the oligomerization site (Fig. 3d). Major perturbations at interface B included an exclusive formation of antiparallel beta-sheets in the mutants, which very likely induced other changes in this region (RMSD ≈ 1.0 Å; Fig. 3c).

Fig. 3
figure 3

Global structural changes observed in cN-II hyperactive mutants. a Structural changes in the R367Q mutant revealed by its superposition with the wild-type cN-II (PDB ID 2XCX). The subunits of cN-II tetramer are represented as a light gray surface or dark gray ribbons. The residues with altered positions of Cα atoms are shaded according to the calculated RMSD values. The active site and mutated residues are highlighted as blue and green spheres, respectively, in one of the protein subunits. b Altered intersubunit interactions in the R367Q mutant. Cartoon representations (colored in light blue) illustrate altered intersubunit contacts in the mutant protein, i.e., loss of contact (highlighted as green spheres), formation of a compensatory interaction between subunits (red spheres), and residues with altered binding partners (orange spheres). c Detail of interface B within the wild type and the R367Q mutant. Each subunit is shown in a different color: light and dark gray (wild type) or yellow and orange (R367Q). Residues involved in altered intersubunit contacts are shown as sticks. The beta-sheet structure observed only in the mutants is highlighted. d Segments with the highest RMSD values located within interface A of the R367Q mutant. These regions are shown in cartoon representation, with subunits in different colors. The interhelical loop spanning residues 405–416 was not found in the crystal structure and was modeled using the ModLoop server for illustration. e Structural changes observed in the R238W (PDB ID 5L4Z) and L375F mutants as revealed by a superposition with the wild-type structure (surface representation) and by analysis of polar intersubunit contacts (cartoon representation). Coloring is identical as for the R367Q mutant (panels a and b). f Superposition of the structures obtained for the studied mutants with the ATP-bound wild-type cN-II structure (PDB ID 2XCW). The coloring is identical as for the structures in panels a, b, and e

Full size image

Interestingly, interface A also includes helix A (amino acid residues 355–365), which was well defined in all the mutant crystal structures. Helix A is formed only in the active state of the cN-II enzyme [9], and the stabilization of this region in the mutants is consistent with their high catalytic potency even in the absence of ATP. Furthermore, the superpositions of the mutants with the ATP-bound wild-type structure (PDB ID 2XCW) yielded reduced RMSD values in comparison to the values obtained for the wild-type structure in the absence of ATP (Fig. 3f), indicating the constitutively active state of misregulated variants.

To gain further insight into the rearrangement at the oligomeric interface of the mutants, we analyzed intersubunit contacts using Protein Interfaces, Surfaces and Assemblies (PISA) software [22]. The results confirmed differences in the pattern of intersubunit interactions at both interfaces that were accompanied by a moderate decrease in oligomerization area of the hyperactive variants (Fig. 3 and Additional file 9). In the mutants, the interface A regions were extensively rearranged and formed by unique polar interactions that were not found in the wild-type protein. These unique interactions represented nearly half of all polar intersubunit contacts in the studied cN-II variants. On the other hand, the mutants had a smaller number of hydrophobic interactions, which was reflected by a twofold lower solvation energy in interface A (Fig. 3b and Additional file 9). In contrast, the changes in interface B were located predominantly within one specific segment that was altered due to the formation of the antiparallel beta-sheet in the mutants (Fig. 3c).

In summary, the analysis of crystal structures demonstrates that all studied mutations induce topological changes almost exclusively at the oligomerization sites of the cN-II tetramer. These results show that the mutations, regardless of their position within the protein, cause essentially identical structural perturbations (Fig. 3a and e). This defines a common molecular mechanism of hyperactivation for all studied variants. X-ray crystallography revealed that the local changes induced by each point mutation propagate throughout the entire biomolecule, causing the rearrangement within the oligomeric interface of cN-II. These alterations are also accompanied by the stabilization of helix A, which had been previously identified as a major structural component in the cN-II allosteric network. Therefore, our data provide detailed mechanistic insight into how ALL-specific mutations cause hyperactivity.

Conformational changes are accompanied by altered protein dynamics

To further examine the effects of the mutations, we studied the full-length wild-type and mutant cN-II proteins in solution. Using hydrogen/deuterium exchange (HDX) mass spectrometry, we examined the microenvironment of the backbone amide protons.

Consistent with the crystallographic data, HDX showed that perturbations in the mutants are localized mainly at the oligomeric interface of the cN-II protein (Fig. 4). Moreover, all studied mutants exhibited altered HDX kinetics in essentially identical regions, supporting the assumption that the variants share a common structural mechanism underlying enzyme misregulation. We identified regions with altered exchange rates at both oligomeric interfaces. In particular, interface A exhibited faster HDX in the mutants than in the wild-type protein, suggesting a more open conformation. In contrast, interface B had slower exchange rates in the mutants, indicating a higher rigidity within this region when compared to the wild type. These observations are in agreement with the structural changes identified by X-ray crystallography.

Fig. 4
figure 4

Structural changes in the cN-II mutants captured by H/D exchange. a Perturbations of all analyzed mutants are observed mainly at the oligomeric interface. Regions with faster (red) or slower (blue) deuteration kinetics in the mutants (related to the wild-type protein) are highlighted on the cN-II tetramer. The perturbed regions are indicated only at two subunits (surface representation) of tetramer. Graphs (percentage of deuteration) show detailed kinetics of these regions (black circle wild type, blue triangle R367Q, red square R238W, green diamond L375F). b HDX provides additional insight into functionally important sites in cN-II proteins. In particular, faster exchange of active site region in the mutants and high level of deuteration for helix A in all cN-II variants are shown. In addition, HDX kinetics of structured regions that do not exhibit significant changes are shown for illustration

Full size image

Furthermore, HDX measurements provided insight into the behavior of cN-II proteins that was not captured by crystallography (Fig. 4b). Peptides spanning helix A exhibited fast HDX in the wild type as well as in the mutants, indicating a high degree of conformational flexibility for this segment (a comparison with peptides from fully structured regions of cN-II is shown in Fig. 4b). In addition, the region containing active site residue D52 (residues 47–54) exhibited increased exchange rates in each mutant. This indicates that conformational changes in hyperactive mutants are transmitted to the catalytic pocket, resulting in increased solvent accessibility at the active site. Importantly, this observation provides a direct link between the remote conformational changes and the hyperactivity of mutant enzymes.

To further analyze differences observed by HDX, we focused on the perturbed regions in crystal structures. Using RINalyzer [23], we examined residue interaction networks in the wild type and the studied mutants. This analysis revealed that most differences observed with HDX could be explained by conformational changes accompanied with altered intramolecular interactions (see comparison of HDX data with analysis of crystal structures in Additional file 10). However, a relatively small number of regions with perturbed HDX had unaltered topology in the crystal structures. This suggests that these segments had different structural flexibility that could not be captured using X-ray crystallography. Overall, HDX recorded conformational changes predominantly at the oligomeric interface of the cN-II mutants, thus supporting observations from X-ray crystallography. Moreover, HDX measurements provide complementary data about the dynamic behavior of cN-II proteins.

To further probe the flexibility of helix A, we analyzed this segment within the wild-type enzyme and the R367Q mutant, a variant locked in a catalytically active state. Because this region is positively charged and contains three lysine residues (K359, K361, and K362), we performed amine-reactive chemical crosslinking combined with mass spectrometric analysis. Crosslinking reactions of wild-type cN-II and the mutant with bis(sulfosuccinimidyl)glutarate (BS2G) and bis(sulfosuccinimidyl)suberate (BS3) provided 13 and 9 unique crosslinked and one-site modification-containing peptides, respectively (Additional file 11). Quantitative comparison of crosslinking in the wild-type and R367Q enzymes revealed several products with altered abundance (Table 2). Interestingly, the identical crosslinking pattern and similar changes between wild type and mutant were also observed in the proteins complexed with ATP. Differentially modified sites were found in the intersubunit contact area and/or close to helix A. In particular, peptides containing a K344-K359 crosslink or one-site modification of K344 were more abundant in the wild-type cN-II. This observation is consistent with the crystal structures of cN-II proteins because K344 is localized in a region with altered intersubunit contacts (Fig. 5a). Furthermore, we detected K359-K361 crosslinks and one-site modification of K359 to a greater extent in the R367Q mutant. These changes likely reflect a lower reactivity of K344 and a subsequently lower amount of K344-K359 in the mutant.

Table 2 Differentially modified lysine residues revealed from mass spectrometric study of protein crosslinking in the presence and absence of ATP
Full size table
Fig. 5
figure 5

High degree of flexibility for helix A in hyperactive cN-II mutants. a Crosslinked lysine residues located around helix A as identified by mass spectrometry. Crosslinks K344-K359 and K359-K361 were detected even in R367Q, despite the constrained position of residues within the helical region. Solid or dashed arrows indicate lysine crosslinks with higher or lower quantity, respectively, in the particular variant. The R367Q-specific polar intersubunit contact between N117 and K344 is highlighted as a red dashed line. The loop spanning residues 355–365 was modeled in the wild type using the ModLoop server due to missing electron density in this region. b Alternative conformations of the residues distributed within or around the helix A of the cN-II mutants as observed by X-ray crystallography (highlighted as orange or magenta sticks)

Full size image

We identified a K359-K361 crosslink at the region forming helix A. However, we would expect the structurally constrained positions of lysine residues within this helical segment to prevent crosslinking (Fig. 5a). The K359-K361 product was present in the wild-type and R367Q enzymes, indicating that the helix A moves dynamically in both proteins. X-ray crystallography showed that this region is unstructured in the wild type but well-ordered in the mutant. Nevertheless, we identified several residues within helix A or in its vicinity in alternative conformations in the crystals of all mutant variants (Fig. 5b). Overall, these data show that helix A fluctuates between the ordered and disordered state, even in the active conformation of cN-II (i.e., wild type in the presence of ATP and the hyperactive variants) with a properly ordered helical segment in crystal structure. This indicates that the cN-II activity is regulated by a conformational selection or by the shifts in population of the pre-existing conformers, which are well-recognized mechanisms for protein allostery [24, 25].

In summary, mass spectrometry-based approaches provided further insight into the behavior of the studied proteins in solution. These results are consistent with those of X-ray crystallography, identifying intersubunit contact areas including helix A as key regions for regulation of cN-II activity. Mass spectrometric analyses also provided additional structural information, in particular regarding mutation-induced changes within the active site and the conformational flexibility of helix A.

Discussion

Whole-exome sequencing studies have shown that activating mutations in NT5C2 confer chemoresistance in ALL [14, 15]. To elucidate the molecular mechanism underlying the hyperactivity of cN-II mutants, we scrutinized the kinetic and structural behavior of three representative variants (R367Q, R238W, and L375F). We found that the ALL-specific variants are catalytically efficient even in the absence of an allosteric activator, demonstrating that the mutations cause misregulation of cN-II activity.

Using X-ray crystallography combined with HDX and protein crosslinking, we demonstrated that this misregulation is caused by conformational changes at the oligomeric interface of the protein. Our data showed that cN-II activation is associated with an altered pattern of intersubunit interactions accompanied by differences in dynamic behavior at the oligomeric interface. These changes also include stabilization of helix A, a segment that plays a crucial role in cN-II activation.

The role of intersubunit contacts in enzyme regulation was further demonstrated by targeted analysis of the available structural data (Additional file 12). In particular, the activator binding site and helix A are localized within interface A, and the access channel to the active site is formed by interface B. Moreover, superposition of the available crystal structures of wild-type cN-II (PDB IDs 2XCX and 2XCW) reveals specific changes in residues located at the oligomeric interface depending on the presence or absence of allosteric activators. In addition, the distribution of misregulation-causing mutations within the cN-II structure strongly suggests that the oligomerization sites play an essential role in allosteric regulation (Fig. 1a). A structure of wild-type cN-II in complex with noncompetitive inhibitors [26], which induce major perturbations at the regions forming intersubunit contacts, also demonstrates the importance of the oligomeric interface for cN-II regulation. Overall, we brought together several lines of evidence indicating that the oligomeric interface is involved in regulation of cN-II.

Helix A is a crucial region for cN-II allosteric regulation. Importantly, our study provided novel insight into its behavior. Structural comparison of wild-type cN-II and the hyperactive mutants revealed that helix A is well ordered only in an active conformation of cN-II. However, mass spectrometry-based techniques showed that helix A retains a certain degree of flexibility even in the active state of the enzyme. This observation was further supported by crystal structures of hyperactive mutants that revealed alternative conformations of several residues within or around helix A (Fig. 5). However, the techniques applied could not fully capture differences in the rapid order-disorder transition between the wild type and mutants. One may speculate that an active form of cN-II contains a transiently ordered helix A that is further stabilized in the crystallized proteins.

Here, we thoroughly characterized the most common ALL-specific cN-II mutants. The studied variants are representative of the structural topology of the majority of all described mutations [14, 15, 17, 18]. The mutations are localized in two regions: either close to helix A (represented by R367Q) or directly at the intersubunit contact area (represented by R238W and L375F; see Additional file 12). The positioning of the mutated residues within the protein structure corresponded to their kinetic properties and altered thermostability. In particular, the R367Q variant adopted a fully active conformation that exhibited higher thermostability than the wild-type enzyme. On the other hand, the R238W and L375F mutations stimulated catalytic activity in synergy with ATP binding and caused a decrease in protein thermostability. Based on these observations, the activating mutations can be divided into two major groups. One group comprises mutations that directly stabilize the neighborhood of helix A (R39Q, K359Q, S360P, and R367Q). This mode of action was previously suggested for K359Q, while other mutants lacked structural cues explaining hyperactivity [14]. The second group includes oligomeric interface mutations that alter intersubunit contacts (R39Q, R238W, R238G, R238L, R291W, K359Q, L375F, D407A, D407Y, S408R, P414S, and S445F). The R39Q and K359Q mutations belong to both groups, as they may stabilize helix A while affecting intersubunit contacts. These findings imply that the cN-II structure contains two hotspot regions where point mutations can cause enzyme hyperactivity.

Note that the relapsed ALL-specific NT5C2 mutations are heterozygous and complement the wild-type allele. Consequently, cN-II tetramers can be formed by two different subunits. It is tempting to speculate that mutations induce hyperactivity even in the wild-type subunit via allosteric transmission through the oligomeric interface. Exploration of this phenomenon is beyond the scope of the present study, but should be investigated in the future.

Our enzyme kinetics analysis showed that cN-II mutants are not properly regulated and are catalytically effective even in the absence of activators. It is tempting to speculate that hyperactive mutants circumvent physiological feedback inhibition of the enzyme in order to persistently inactivate nucleoside analogs in chemoresistant leukemic cells. As a consequence, the elevated cN-II activity affects nucleotide metabolism homeostasis. We propose that misregulation of purine nucleotidase may represent a common mechanism underlying chemoresistant leukemia driven by mutations in NT5C2 (Fig. 6).

Fig. 6
figure 6

Mechanism underlying chemoresistant ALL driven by mutations in NT5C2. Mutations located either in the vicinity of helix A or at the oligomerization sites induce misregulation of cN-II activity through conformational modulation of the oligomeric interface, which boosts the viability of tumor cells during chemotherapy

Full size image

Interestingly, misregulation of mutated PRPS1, an enzyme responsible for the first step of de novo purine synthesis, was recently confirmed as another cause for chemoresistance in relapsed ALL [27]. This strongly supports the notion that misregulation of purine metabolism affects the outcome of chemotherapy drugs targeting nucleotide metabolism. However, the list of proteins that play a role in drug metabolism remains incomplete. For example, MutT homolog 2 (MTH2) was recently implicated in the metabolism of purine nucleoside analogs [28, 29]. Indeed, the R139C variant of MTH2 causes intolerance to mercaptopurine during the treatment of ALL or Crohn’s disease [30, 31]. These findings jointly demonstrate that drug metabolism and its regulation, including the effects of genetic variations, need to be further investigated to improve the efficacy of currently used therapeutic strategies.

In summary, our data provide comprehensive insight into the allosteric regulation of cN-II, which could be used for the design of novel antileukemic drugs. Our study implies that the oligomeric interface in the cN-II protein represents a potential target for pharmacological intervention. Importantly, the intersubunit contact area contains effector sites that bind physiological ligands [7, 9, 32]. Recently, low-affinity noncompetitive inhibitors have been developed to target this site in wild-type cN-II [26]. In addition, oligomerization sites have been successfully targeted in various screening campaigns against medicinally relevant proteins (e.g., p53, tumor necrosis factor, and G-protein coupled receptors) [33]. Moreover, allosteric sites of metabolic enzymes, such as pyruvate kinase, isocitrate dehydrogenase (IDH), and glutaminase, have been proposed as important targets for specific cancer therapies [34]. Potent allosteric inhibitors of glutaminase [35] and tumor-specific IDH2 [36] mutants are currently being tested in clinical trials [37].

Conclusion

Our study clearly demonstrates that ALL-specific mutations cause misregulation of cN-II due to allosteric modulation of the oligomeric interface, as illustrated in Fig. 6. This defines a common molecular mechanism underlying chemoresistance driven by hyperactive purine nucleotidase. Allosteric inhibition of cN-II mutants may become an effective therapeutic approach for treatment of ALL.

Methods

Materials

Except where stated otherwise, all compounds were purchased from Sigma-Aldrich.

Cloning and plasmid preparation

All plasmids were derived from pET28 expressing C-terminally truncated human cN-II (Δ 536–561) with an N-terminal hexahistidine tag (Source Bioscience, UK) as previously described [7]. Site-directed mutagenesis was performed using the QuikChange Site-Directed Mutagenesis Protocol (Agilent Technologies). To prepare plasmids producing full-length proteins, the region encoding the C-terminal tail was prepared synthetically (Eurofins Scientific) and ligated into the constructs encoding the truncated variants.

Protein expression and purification

Protein samples were prepared according to previously described procedures [7] with several modifications. Fractions collected from nickel affinity chromatography that contained cN-II protein were further purified using ion-exchange columns. Full-length proteins were dialyzed against 50 mM Bis-Tris propane (pH 7.0) containing 10 % glycerol and 2 mM tris(2-carboxyethyl)phosphine (TCEP; buffer A) and subsequently loaded onto a MonoQ column (Amersham Biosciences) and eluted with a gradient of NaCl (0–1 M) in buffer A over 30 column volumes. The C-terminally truncated variants were dialyzed against 50 mM sodium phosphate (pH 6.7) containing 10 % glycerol and 2 mM TCEP (buffer B), loaded onto a MonoS column (Amersham Biosciences), and eluted with a gradient of NaCl (0–1 M) in buffer B over 30 column volumes. Fractions containing cN-II were collected and dialyzed against 50 mM sodium phosphate (pH 7.4) containing 100 mM NaCl, 10 % glycerol, and 2 mM TCEP and stored at –80 °C.

Proteins were purified to homogeneity as judged by silver-stained SDS-PAGE. The identity of the purified cN-II proteins was validated by mass spectrometry-based peptide mass fingerprinting.

Assay of cN-II activity

The catalytic activity of cN-II was assessed with inosine monophosphate (IMP) as a substrate according to a previously described procedure [1]. The nucleotidase activity (0.1 μM cN-II subunit) was measured in the presence or absence of 3 mM MgATP, with the IMP concentrations ranging between 10 μM to 60 mM, in 60 mM imidazole buffer (pH 7.4) containing 150 mM KCl, 1 mM MgCl2, and 0.5 mM dithiothreitol (DTT). The reaction was performed at 37 °C for 1–20 minutes and then chilled on ice and quenched with 0.15 M ethylenediaminetetraacetic acid. The rate of enzymatic conversion was analyzed using reversed-phase high-performance liquid chromatography. Components in the reaction mixture were separated on a Gemini 5 μm C18 110 Å column (Phenomenex) by isocratic elution with 0.1 M triethylammonium bicarbonate, and signals were obtained using UV detection at λ = 262 nm. The enzyme activity was evaluated by integrating of the peaks with mobility corresponding to the standards of substrate and product. Nonlinear data fitting was performed using OriginPro 8.5 (OriginLab). The linearity of the assay was tested in the presence and absence of ATP during a reaction time of 20 minutes as illustrated in Additional file 1. Kinetic properties of the wild-type cN-II, including a response upon ATP binding, were in agreement with a previous study on human recombinant enzymes [10] and with reports deposited in the BRaunschweig ENzyme DAtabase (BRENDA).

Differential scanning fluorimetry

Proteins (0.1 mg/ml) were dissolved in 50 mM sodium phosphate (pH 7.4) containing 100 mM NaCl, 2 mM TCEP and 8× SYPRO Orange. To analyze the effect of ATP on protein thermostability, the proteins were incubated with the ligand for 10 minutes at room temperature prior to addition of the SYPRO Orange dye. Using a LightCycler 480 (Roche) real-time detection system, the proteins were heated from 20 °C to 90 °C in increments of 0.5 °C and with 1-minute hold intervals. The signal was monitored by fluorescence detection (excitation wavelength of 470 nm and emission wavelength of 570 nm). The melting temperatures (T m) of proteins were determined as minima from first derivative curves.

Saturation transfer difference NMR

All data for cN-II proteins were collected at 25 °C using a Bruker Avance III™ HD 850 MHz spectrometer equipped with a 5 mm CPTCI 1H/13C/15N/D Z-GRD cryoprobe. In order to compare the relative ATP affinities towards the cN-II variants, 1H STD NMR spectra were acquired using 350-μl samples of 1.5 μM cN-II proteins dissolved in 50 mM sodium phosphate buffer (pH 7.4) containing 100 mM NaCl, 2 mM TCEP, and 5 % D2O/95 % H2O. The frequency used for non-selective irradiation of the protein signals was 666 Hz (0.78 ppm), using a 50-ms shaped pulse Eburp2.1000 at a power of 40 dB. The spectra were acquired with a free induction decay (FID) resolution of 0.83 Hz, and the typical experimental time was 20 minutes.

Size exclusion chromatography

Proteins were dissolved in 50 mM sodium phosphate (pH 7.4) containing 100 mM NaCl and 2 mM TCEP and analyzed using a Superdex 200 10/300 GL column. The column was calibrated with a high molecular weight protein standards kit (Amersham Biosciences). For analysis, 200 μl of cN-II protein samples (0.5 mg/ml) were injected.

Small-angle X-ray scattering

Truncated proteins of varying concentrations (from 0.25 to 2 mg/ml) dissolved in 50 mM sodium phosphate (pH 7.4) containing 100 mM NaCl, 2 mM TCEP, and 5 % glycerol were analyzed at the BM29 beamline at the European Synchrotron Radiation Facility (ESRF) (Grenoble) with a Pilatus 1 M detector. The measurements were performed at a sample-detector distance of 2.864 m and a wavelength of 0.992 Å. Data were processed with PRIMUS [38], and the scattering curves were analyzed using Guinier approximation to calculate the forward scattering I (0) and radius of gyration (R g). A solution of β-amylase from sweet potato (200 kDa) was used as a reference sample. Distance distribution functions P (r) , maximum particle dimensions D max [39], and Porod volume were computed with GNOM [40]. All used programs are included in the ATSAS software package [41].

Blue native polyacrylamide gel electrophoresis (BN-PAGE)

BN-PAGE was performed using 4–16 % polyacrylamide Bis-Tris precast gels (Invitrogen) at room temperature according to the manufacturer’s protocol with high molecular weight native marker kit (Amersham Biosciences) and β-amylase from sweet potato as protein standards. A 5-μg aliquot of protein was loaded in each lane as described previously [42].

X-ray crystallography

The C-terminally truncated mutant proteins lacking residues 537–561 were crystallized as described for wild-type cN-II [7]. The proteins (8 mg/ml) in 100 mM sodium citrate (pH 6.0) containing 2 mM TCEP were crystallized using the hanging drop vapor diffusion technique. The precipitant for each mutant was as follows: R367Q, 100 mM MES/imidazole (pH 6.5) containing 200 mM NaCl, 30 % glycerol, and 10 % PEG 4 k; R238W, 100 mM MOPS/HEPES-Na (pH 7.5) containing 100 mM NaCl, 15 % glycerol, and 10 % PEG 4 k; L375F, 200 mM MES/imidazol (pH 6.5) containing 100 mM NaCl, 30 % glycerol, and 10 % PEG 4 k. The crystals grew in 1.5-μl drops of 2:1 protein:precipitant solution for 1 week at room temperature and subsequently were transferred into a drop of perfluoropolyether for 20 s and then flash-frozen in liquid nitrogen. Diffraction data were collected on BL14.1 operated by the Joint Berlin MX Laboratory at the BESSY II electron storage ring (Berlin-Adlershof, Germany) [43]. Data processing was performed using the XDS software package [44]. The 3D structure of the mutants was solved by molecular replacement with the structure of wild-type enzyme (PDB ID 2J2C) as a template using the MOLREP program [45]. Refinement was performed with REFMAC 11.0 [46], including translation/libration/screwing (TLS) refinement [47] with 12 TLS groups at a later stage of refinement. Manual rebuilding of the model was done using Coot [48]. The MolProbity server [49] was used for assessment of the final model quality. Structural statistics are shown in Additional file 13. The protein structures were analyzed using software included in the CCP4 package [50], namely PISA [22], LSQKAB [51], and BAVERAGE. In addition, the residue interaction network was assessed with RINalyzer [23]. Structural models were visualized using PyMOL software (DeLano Scientific). Missing loops were modeled using the ModLoop server [52] to illustrate the possible position of these regions in the protein structure. These are not part of the models deposited in the PDB under accession codes: 5K7Y (R367Q), 5L4Z (R238W), and 5L50 (L375F).

Hydrogen/deuterium exchange mass spectrometry

Full-length proteins (1 mg/ml) were dissolved in 20 mM sodium carbonate (pH 7.4) containing 100 mM potassium chloride, 1 mM magnesium chloride, and 0.5 mM TCEP. Proteins were diluted to a final concentration of 0.2 mg/ml with the same buffer prepared in D2O (pD 7.4) to initiate HDX that was subsequently quenched at varying time points by adding the same volume of 0.25 M glycine-HCl buffer (pH 2.3) and immediately freezing the samples on liquid nitrogen. Liquid chromatography-mass spectrometry (LC-MS) analysis, including on-line pepsin digestion, data acquisition, processing, and evaluation, was performed according to a previously described procedure [53].

Chemical crosslinking studied by mass spectrometry

Full-length proteins (0.2 mg/ml) dissolved in 50 mM triethylammonium bicarbonate (pH 7.4) containing 50 mM NaCl, 1 mM MgCl2, and 0.5 mM DTT were incubated with the crosslinking agents BS2G and BS3 in protein:crosslinker ratios of 1:20 and 1:50, respectively. Reactions were carried out at room temperature for 1 h and then quenched with 0.15 mM ethanolamine. Crosslinked proteins were reduced with DTT, alkylated with iodoacetamide, and digested using trypsin. To compare the crosslinking pattern of wild type and R367Q, isotopically labeled crosslinkers were used, and a quantitative comparison of light/heavy crosslinked peptides was performed by LC-MS as described elsewhere [54].

Crosslinked proteins were analyzed using SDS-PAGE to check that the reactions did not yield nonspecific aggregates (Additional file 14).

References

  1. Spychala J, Madrid-Marina V, Fox IH. High Km soluble 5'-nucleotidase from human placenta. Properties and allosteric regulation by IMP and ATP. J Biol Chem. 1988;263(35):18759–65.

    CAS  PubMed  Google Scholar 

  2. Itoh R. Enzymatic properties and physiological roles of cytosolic 5 '-nucleotidase II. Curr Med Chem. 2013;20(34):4260–84.

    Article  CAS  PubMed  Google Scholar 

  3. Bianchi V, Spychala J. Mammalian 5'-nucleotidases. J Biol Chem. 2003;278(47):46195–8.

    Article  CAS  PubMed  Google Scholar 

  4. Itoh R. Regulation of cytosol 5'-nucleotidase by adenylate energy charge. Biochim Biophys Acta. 1981;659(1):31–7.

    Article  CAS  PubMed  Google Scholar 

  5. Itoh R. Studies on some molecular properties of cytosol 5'-nucleotidase from rat liver. Biochim Biophys Acta. 1982;716(1):110–3.

    Article  CAS  PubMed  Google Scholar 

  6. Oka J, Matsumoto A, Hosokawa Y, Inoue S. Molecular cloning of human cytosolic purine 5'-nucleotidase. Biochem Biophys Res Commun. 1994;205(1):917–22.

    Article  CAS  PubMed  Google Scholar 

  7. Wallden K, Stenmark P, Nyman T, Flodin S, Graslund S, Loppnau P, Bianchi V, Nordlund P. Crystal structure of human cytosolic 5'-nucleotidase II: insights into allosteric regulation and substrate recognition. J Biol Chem. 2007;282(24):17828–36.

    Article  CAS  PubMed  Google Scholar 

  8. Allegrini S, Scaloni A, Careddu MG, Cuccu G, D'Ambrosio C, Pesi R, Camici M, Ferrara L, Tozzi MG. Mechanistic studies on bovine cytosolic 5'-nucleotidase II, an enzyme belonging to the HAD superfamily. Eur J Biochem. 2004;271(23-24):4881–91.

    Article  CAS  PubMed  Google Scholar 

  9. Wallden K, Nordlund P. Structural basis for the allosteric regulation and substrate recognition of human cytosolic 5'-nucleotidase II. J Mol Biol. 2011;408(4):684–96.

    Article  CAS  PubMed  Google Scholar 

  10. Spychala J, Chen V, Oka J, Mitchell BS. ATP and phosphate reciprocally affect subunit association of human recombinant High Km 5'-nucleotidase. Role for the C-terminal polyglutamic acid tract in subunit association and catalytic activity. Eur J Biochem. 1999;259(3):851–8.

    Article  CAS  PubMed  Google Scholar 

  11. Pesi R, Micheli V, Jacomelli G, Peruzzi L, Camici M, Garcia-Gil M, Allegrini S, Tozzi MG. Cytosolic 5'-nucleotidase hyperactivity in erythrocytes of Lesch-Nyhan syndrome patients. Neuroreport. 2000;11(9):1827–31.

    Article  CAS  PubMed  Google Scholar 

  12. Galmarini CM, Graham K, Thomas X, Calvo F, Rousselot P, El Jafaari A, Cros E, Mackey JR, Dumontet C. Expression of high Km 5'-nucleotidase in leukemic blasts is an independent prognostic factor in adults with acute myeloid leukemia. Blood. 2001;98(6):1922–6.

    Article  CAS  PubMed  Google Scholar 

  13. Mazzon C, Rampazzo C, Scaini MC, Gallinaro L, Karlsson A, Meier C, Balzarini J, Reichard P, Bianchi V. Cytosolic and mitochondrial deoxyribonucleotidases: activity with substrate analogs, inhibitors and implications for therapy. Biochem Pharmacol. 2003;66(3):471–9.

    Article  CAS  PubMed  Google Scholar 

  14. Tzoneva G, Perez-Garcia A, Carpenter Z, Khiabanian H, Tosello V, Allegretta M, Paietta E, Racevskis J, Rowe JM, Tallman MS, et al. Activating mutations in the NT5C2 nucleotidase gene drive chemotherapy resistance in relapsed ALL. Nat Med. 2013;19(3):368–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Meyer JA, Wang J, Hogan LE, Yang JJ, Dandekar S, Patel JP, Tang Z, Zumbo P, Li S, Zavadil J, et al. Relapse-specific mutations in NT5C2 in childhood acute lymphoblastic leukemia. Nat Genet. 2013;45(3):290–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Aster JC, DeAngelo DJ. Resistance revealed in acute lymphoblastic leukemia. Nat Med. 2013;19(3):264–5.

    Article  CAS  PubMed  Google Scholar 

  17. Ma X, Edmonson M, Yergeau D, Muzny DM, Hampton OA, Rusch M, Song G, Easton J, Harvey RC, Wheeler DA, et al. Rise and fall of subclones from diagnosis to relapse in pediatric B-acute lymphoblastic leukaemia. Nat Commun. 2015;6:6604.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Kunz JB, Rausch T, Bandapalli OR, Eilers J, Pechanska P, Schuessele S, Assenov Y, Stutz AM, Kirschner-Schwabe R, Hof J, et al. Pediatric T-lymphoblastic leukemia evolves into relapse by clonal selection, acquisition of mutations and promoter hypomethylation. Haematologica. 2015;100:1442–50.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Pesi R, Turriani M, Allegrini S, Scolozzi C, Camici M, Ipata PL, Tozzi MG. The bifunctional cytosolic 5'-nucleotidase: regulation of the phosphotransferase and nucleotidase activities. Arch Biochem Biophys. 1994;312(1):75–80.

    Article  CAS  PubMed  Google Scholar 

  20. Ipata PL, Balestri F. The functional logic of cytosolic 5 '-nucleotidases. Curr Med Chem. 2013;20(34):4205–16.

    Article  CAS  PubMed  Google Scholar 

  21. Carlomagno T. Ligand-target interactions: what can we learn from NMR? Annu Rev Biophys Biomol Struct. 2005;34:245–66.

    Article  CAS  PubMed  Google Scholar 

  22. Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372(3):774–97.

    Article  CAS  PubMed  Google Scholar 

  23. Doncheva NT, Assenov Y, Domingues FS, Albrecht M. Topological analysis and interactive visualization of biological networks and protein structures. Nat Protoc. 2012;7(4):670–85.

    Article  CAS  PubMed  Google Scholar 

  24. del Sol A, Tsai CJ, Ma B, Nussinov R. The origin of allosteric functional modulation: multiple pre-existing pathways. Structure. 2009;17(8):1042–50.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Nussinov R, Ma B, Tsai CJ. Multiple conformational selection and induced fit events take place in allosteric propagation. Biophys Chem. 2014;186:22–30.

    Article  CAS  PubMed  Google Scholar 

  26. Marton Z, Guillon R, Krimm I, Preeti, Rahimova R, Egron D, Jordheim LP, Aghajari N, Dumontet C, Perigaud C, et al. Identification of noncompetitive inhibitors of cytosolic 5 '-nucleotidase II using a fragment-based approach. J Med Chem. 2015;58(24):9680–96.

    Article  CAS  PubMed  Google Scholar 

  27. Li BS, Li H, Bai Y, Kirschner-Schwabe R, Yang JJ, Chen Y, Lu G, Tzoneva G, Ma XT, Wu TM, et al. Negative feedback-defective PRPS1 mutants drive thiopurine resistance in relapsed childhood ALL. Nat Med. 2015;21(6):563–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Carter M, Jemth AS, Hagenkort A, Page BDG, Gustafsson R, Griese JJ, Gad H, Valerie NCK, Desroses M, Bostrom J, et al. Crystal structure, biochemical and cellular activities demonstrate separate functions of MTH1 and MTH2. Nat Commun. 2015;6:7871.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Moriyama T, Nishii R, Perez-Andreu V, Yang WJ, Klussmann FA, Zhao XJ, Lin TN, Hoshitsuki K, Nersting J, Kihira K, et al. NUDT15 polymorphisms alter thiopurine metabolism and hematopoietic toxicity. Nat Genet. 2016;48(4):367–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yang SK, Hong M, Baek J, Choi H, Zhao WT, Jung YS, Haritunians T, Ye BD, Kim KJ, Park SH, et al. A common missense variant in NUDT15 confers susceptibility to thiopurine-induced leukopenia. Nat Genet. 2014;46(9):1017–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yang JJ, Landier W, Yang WJ, Liu CC, Hageman L, Cheng C, Pei DQ, Chen YJ, Crews KR, Kornegay N, et al. Inherited NUDT15 variant is a genetic determinant of mercaptopurine intolerance in children with acute lymphoblastic leukemia. J Clin Oncol. 2015;33(11):1235–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Pesi R, Allegrini S, Careddu MG, Filoni DN, Camici M, Tozzi MG. Active and regulatory sites of cytosolic 5 '-nucleotidase. Febs J. 2010;277(23):4863–72.

    Article  CAS  PubMed  Google Scholar 

  33. Gabizon R, Friedler A. Allosteric modulation of protein oligomerization: an emerging approach to drug design. Front Chem. 2014;2:9.

    Article  PubMed  PubMed Central  Google Scholar 

  34. DeLaBarre B, Hurov J, Cianchetta G, Murray S, Dang L. Action at a distance: allostery and the development of drugs to target cancer cell metabolism. Chem Biol. 2014;21(9):1143–61.

    Article  CAS  PubMed  Google Scholar 

  35. DeLaBarre B, Gross S, Fang C, Gao Y, Jha A, Jiang F, Song JJ, Wei WT, Hurov JB. Full-length human glutaminase in complex with an allosteric inhibitor. Biochemistry. 2011;50(50):10764–70.

    Article  CAS  PubMed  Google Scholar 

  36. Wang F, Travins J, DeLaBarre B, Penard-Lacronique V, Schalm S, Hansen E, Straley K, Kernytsky A, Liu W, Gliser C, et al. Targeted inhibition of mutant IDH2 in leukemia cells induces cellular differentiation. Science. 2013;340(6132):622–6.

    Article  CAS  PubMed  Google Scholar 

  37. Gross MI, Demo SD, Dennison JB, Chen L, Chernov-Rogan T, Goyal B, Janes JR, Laidig GJ, Lewis ER, Li J, et al. Antitumor activity of the glutaminase inhibitor CB-839 in triple-negative breast cancer. Mol Cancer Ther. 2014;13(4):890–901.

    Article  CAS  PubMed  Google Scholar 

  38. Konarev PV, Volkov VV, Sokolova AV, Koch MHJ, Svergun DI. PRIMUS: a Windows PC-based system for small-angle scattering data analysis. J Appl Crystallogr. 2003;36:1277–82.

    Article  CAS  Google Scholar 

  39. Konig S, Svergun D, Koch MHJ, Hubner G, Schellenberger A. Synchrotron radiation solution X-ray scattering study of the pH dependence of the quaternary structure of yeast pyruvate decarboxylase. Biochemistry. 1992;31(37):8726–31.

    Article  CAS  PubMed  Google Scholar 

  40. Svergun DI. Determination of the regularization parameter in indirect-transform methods using perceptual criteria. J Appl Crystallogr. 1992;25:495–503.

    Article  Google Scholar 

  41. Petoukhov MV, Franke D, Shkumatov AV, Tria G, Kikhney AG, Gajda M, Gorba C, Mertens HDT, Konarev PV, Svergun DI. New developments in the ATSAS program package for small-angle scattering data analysis. J Appl Crystallogr. 2012;45:342–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Filoni DN, Pesi R, Allegrini S, Camici M, Tozzi MG. A native electrophoretic technique to study oligomerization and activity of cytosolic 5'-nucleotidase II. Anal Bioanal Chem. 2013;405(27):8951–4.

    Article  CAS  PubMed  Google Scholar 

  43. Mueller U, Darowski N, Fuchs MR, Forster R, Hellmig M, Paithankar KS, Puhringer S, Steffien M, Zocher G, Weiss MS. Facilities for macromolecular crystallography at the Helmholtz-Zentrum Berlin. J Synchrotron Radiat. 2012;19:442–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kabsch W. XDS. Acta Crystallogr D. 2010;66:125–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Vagin A, Teplyakov A. Molecular replacement with MOLREP. Acta Crystallogr D. 2010;66:22–5.

    Article  CAS  PubMed  Google Scholar 

  46. Murshudov GN, Vagin AA, Dodson EJ. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D. 1997;53:240–55.

    Article  CAS  PubMed  Google Scholar 

  47. Winn MD, Isupov MN, Murshudov GN. Use of TLS parameters to model anisotropic displacements in macromolecular refinement. Acta Crystallogr D. 2001;57:122–33.

    Article  CAS  PubMed  Google Scholar 

  48. Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D. 2004;60:2126–32.

    Article  PubMed  Google Scholar 

  49. Lovell SC, Davis IW, Adrendall WB, de Bakker PIW, Word JM, Prisant MG, Richardson JS, Richardson DC. Structure validation by C alpha geometry: phi, psi and C beta deviation. Proteins. 2003;50(3):437–50.

    Article  CAS  PubMed  Google Scholar 

  50. Winn MD, Ballard CC, Cowtan KD, Dodson EJ, Emsley P, Evans PR, Keegan RM, Krissinel EB, Leslie AGW, McCoy A, et al. Overview of the CCP4 suite and current developments. Acta Crystallogr D. 2011;67:235–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Krissinel E, Henrick K. Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensions. Acta Crystallogr D. 2004;60:2256–68.

    Article  CAS  PubMed  Google Scholar 

  52. Fiser A, Sali A. ModLoop: automated modeling of loops in protein structures. Bioinformatics. 2003;19(18):2500–1.

    Article  CAS  PubMed  Google Scholar 

  53. Trcka F, Durech M, Man P, Hernychova L, Muller P, Vojtesek B. The assembly and intermolecular properties of the Hsp70-Tomm34-Hsp90 molecular chaperone complex. J Biol Chem. 2014;289(14):9887–901.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Kukacka Z, Rosulek M, Strohalm M, Kavan D, Novak P. Mapping protein structural changes by quantitative cross-linking. Methods. 2015;89:112–20.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Katsiaryna Tratsiak, MSc., for technical help; Dr. Zdeněk Kukačka and Alan Kádek, MSc., for helpful assistance during mass spectrometric measurements; and Dr. Ivan Rosenberg and Dr. Petr Těšina for useful advice and comments during the entire project, including manuscript preparation.

Funding

This work was supported by grant 15-06582S from the Czech Science Foundation and in part by the Ministry of Education of the Czech Republic (program “NAVRAT” LK11205, program “NPU I” LO1304, and program “InterBioMed” LO1302). Institutional support was provided by projects RVO 61388963 and 68378050 of the Academy of Sciences of the Czech Republic. PN and PM were supported by the European Regional Development Fund (CZ.1.05/1.1.00/02.0109) and the Ministry of Education of the Czech Republic (LD15089).

We acknowledge beamline access at MX14.1 of the BESSY (HZB Berlin, Germany) and BM29 of the ESRF (Grenoble, France) to collect crystal diffraction and SAXS data, respectively. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

The structures described here were deposited in the Protein Data Base [PDB:5K7Y, 5L4Z, and 5L50].

All supporting data used for preparation of this paper are available in the additional files.

Authors’ contributions

AH designed the entire study, performed and analyzed experiments, interpreted data, and wrote the manuscript. JS, PP, and PR performed and analyzed X-ray crystallography. MF prepared all plasmids. MS performed and analyzed the biochemical experiments. LV performed and analyzed the NMR experiments. PM and PN designed, performed, and analyzed the mass spectrometry experiments. VV designed and analyzed the NMR experiments. PR, MF, and VV revised all versions of the manuscript. All authors read, revised, and accepted the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

  1. Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, Prague 6, 166 10, Czech Republic

    Aleš Hnízda, Jana Škerlová, Petr Pachl, Martina Šinalová, Lukáš Vrzal, Pavlína Řezáčová & Václav Veverka

  2. Institute of Molecular Genetics, Academy of Sciences of the Czech Republic, Videnska 1083, Prague 4, 142 20, Czech Republic

    Milan Fábry & Pavlína Řezáčová

  3. Institute of Microbiology, Academy of Sciences of the Czech Republic, Videnska 1083, Prague 4, 142 20, Czech Republic

    Petr Man & Petr Novák

Authors
  1. Aleš Hnízda
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jana Škerlová
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Milan Fábry
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Petr Pachl
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Martina Šinalová
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Lukáš Vrzal
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Petr Man
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Petr Novák
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Pavlína Řezáčová
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Václav Veverka
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Aleš Hnízda or Václav Veverka.

Additional files

Additional file 1:

Validation of the enzyme assay. a The cN-II activity under varying concentrations of ATP. The analysis was performed using 1 mM IMP as substrate. The points represent mean values with standard deviation from four measurements. b The test of linearity for catalyzed reaction for cN-II in the presence and absence of ATP. The analysis was performed with 5 mM IMP as substrate. (TIF 421 kb)

Additional file 2:

Temperature dependence of cN-II catalytic activity. The activity values are relative to the values obtained at 40 °C with 1 mM IMP as substrate. The points represent mean values from two independent measurements. (TIF 452 kb)

Additional file 3:

Thermostability of cN-II proteins in the presence of ATP studied by differential scanning fluorimetry. a Representative curves of the wild-type and mutant proteins in the presence or absence of 3 mM ATP. b The T m values of the cN-II proteins under varying concentrations of ATP. Each point represents the mean value from two measurements. (TIF 1192 kb)

Additional file 4:

1H STD NMR spectra of ATP at various concentrations in the presence of cN-II variants. The ATP reference 1H NMR spectrum (standard atom numbering) is shown at the bottom. (TIF 1286 kb)

Additional file 5:

Representative data from SAXS measurements: analysis of wild-type cN-II (1 mg/ml). a Processed scattering curve. b Guinier analysis. c Distance distribution analysis including calculation of the P (r) function. (TIF 836 kb)

Additional file 6:

Kinetic parameters of the C-terminally truncated enzymes in the presence and absence of ATP. The n value refers to the calculated Hill coefficient. (DOCX 15 kb)

Additional file 7:

Differential scanning fluorimetry of the C-terminally truncated cN-II variants in the presence or absence of 3 mM ATP. Truncated proteins are stabilized upon ATP binding, as described for full-length variants. Mean values and standard deviations from duplicates of two independent measurements are listed. (DOCX 14 kb)

Additional file 8:

Size exclusion chromatography of the C-terminally truncated cN-II proteins that were studied by X-ray crystallography. Each variant formed tetramers as reported for full-length proteins. The column was calibrated using a high molecular weight protein standard kit (Amersham Biosciences). (TIF 284 kb)

Additional file 9:

Oligomeric interface of cN-II crystals analyzed by PISA. NHB and NSB represent the number of identified hydrogen bonds and salt bridges, respectively. ∆i G is the energy of solvation. (DOCX 15 kb)

Additional file 10:

Regions with perturbed deuteration kinetics analyzed in crystal structures using RINalyzer. The ‘+’ and ‘-’ symbols refer to the presence or absence of structural changes, respectively. (DOCX 15 kb)

Additional file 11:

Cross-linking of the wild-type enzyme and R367Q mutant: modified residues identified and quantified using mass spectrometry. Relative abundance refers to a quantitative comparison between products obtained from wild-type and R367Q. The type of crosslinking (i.e., intermolecular, intramolecular or hanging) was inferred from the mass spectrometric data or from the crystal structure. Each number represents the mean value and standard deviation from three experiments. (DOCX 18 kb)

Additional file 12:

Structural analysis demonstrates an important role for the oligomeric interface in cN-II allosteric regulation. a Position of helix A (red dots) and the ATP binding site (orange dots) within the cN-II tetramer. b The access channel towards the active site is formed by interface B. Each subunit is highlighted in a different color; the active site residue D52 is shown in red. c Superposition of wild type in its free state and in complex with ATP (PDB IDs 2XCX and 2XCW). Differences in positions of Cα atoms are highlighted in the structure of the apo form (PDB ID 2XCX). The active site residue (D52) is shown as blue spheres. d Location of ALL-specific mutations reveals two hotspot regions. Each subunit is depicted in a different color, and helix A is shown as an orange cylinder. Mutated residues are highlighted as red spheres within one subunit shown in ribbon representation (wheat colored). (TIF 10513 kb)

Additional file 13:

Crystal parameters, data collection and refinement statistics. (DOCX 20 kb)

Additional file 14:

Representative SDS-PAGE of the wild-type cN-II protein upon crosslinking reaction with bis(sulfosuccinimidyl)glutarate (BS2G) and bis(sulfo-succinimidyl)suberate (BS3). M and N refer to molecular weight marker and non-modified sample, respectively. The protein:crosslinker molar ratio is indicated at the top of each lane. The arrows indicate the number of subunits contributing to crosslinking products. (TIF 1253 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hnízda, A., Škerlová, J., Fábry, M. et al. Oligomeric interface modulation causes misregulation of purine 5´-nucleotidase in relapsed leukemia. BMC Biol 14, 91 (2016). https://doi.org/10.1186/s12915-016-0313-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12915-016-0313-y

Keywords

BMC Biology

ISSN: 1741-7007

Contact us