Oligomeric interface modulation causes misregulation of purine 5´-nucleotidase in relapsed leukemia
© Hnizda et al. 2016
Received: 25 July 2016
Accepted: 23 September 2016
Published: 19 October 2016
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).
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.
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.
KeywordsNucleotidase Cancer mutations Relapsed ALL Purine metabolism Allosteric regulation
Cytosolic purine 5´-nucleotidase (cN-II; EC 126.96.36.199) dephosphorylates purine nucleotide monophosphates with canonical activity towards inosine monophosphate (IMP) . IMP represents a branch point intermediate linking major routes of purine metabolism, including de novo synthesis, the salvage pathway and oxidative degradation . Therefore, cN-II serves as an essential regulatory element to maintain the balance of the nucleotide pool . 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 . 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 , with each subunit comprising 561 amino acid residues . 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 .
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) . To achieve a full catalytic activity, cN-II requires the presence of adenylate compounds, such as ATP, adenosine diphosphate or diadenosine tetraphosphate . 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 . 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 . 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 , 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 . In addition, cN-II activity has been implicated in decreased therapeutic efficacy of nucleoside analogs used for the treatment of cancer and viral diseases .
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 . 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.
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 . 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  or just V max . 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 .
Kinetic parameters of cN-II enzymes in the presence and absence of ATP
V max [μmol.min-1mg-1]
26.4 ± 4.7
33.5 ± 9.1
1.27 ± 0.38
29.5 ± 2.1
0.88 ± 0.25
31.6 ± 1.8
2.9 ± 0.8
0.90 ± 0.14
35.3 ± 2.1
12.3 ± 3.3
31.4 ± 2.9
5.8 ± 1.7
0.88 ± 0.16
35.1 ± 3.3
6.0 ± 1.8
28.8 ± 1.3
2.8 ± 0.8
0.85 ± 0.10
32.2 ± 1.5
11.4 ± 2.1
23.4 ± 1.9
12.1 ± 3.1
0.85 ± 0.09
26.1 ± 2.1
2.1 ± 0.6
31.2 ± 2.5
0.52 ± 0.25
0.54 ± 0.15
34.8 ± 2.8
67 ± 33
36.2 ± 3.9
12.2 ± 3.9
0.78 ± 0.10
40.5 ± 3.4
3.3 ± 1.1
29.5 ± 2.3
0.16 ± 0.08
0.49 ± 0.16
33.0 ± 2.6
206 ± 105
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
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 . 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 . 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.
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.
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 , 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 . 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.
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 , 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.
Differentially modified lysine residues revealed from mass spectrometric study of protein crosslinking in the presence and absence of ATP
Relative abundance (wild type compared to R367Q)
Protein:BS2G molar ratio
Protein:BS3 molar ratio
Lys 344; Lys 359
77-23 ± 20
88-12 ± 12
47-53 ± 24
65-35 ± 1
88-12 ± 12
89-11 ± 11
62-38 ± 33
52-48 ± 43
Lys 359; Lys 361
37-63 ± 4
48-52 ± 13
20-80 ± 11
19-81 ± 5
55-45 ± 5
43-57 ± 23
53-47 ± 5
44-56 ± 3
24-76 ± 11
34-66 ± 15
17-83 ± 12
14-86 ± 14
24-76 ± 8
32-68 ± 18
30-70 ± 8
17-83 ± 8
80-20 ± 16
82-18 ± 10
94-6 ± 4
76-24 ± 13
84-16 ± 16
85-15 ± 4
93-7 ± 7
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.
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 , 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 . 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.
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 . 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 . 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) . 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 . Potent allosteric inhibitors of glutaminase  and tumor-specific IDH2  mutants are currently being tested in clinical trials .
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.
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 . 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  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 . 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  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 , 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 , and Porod volume were computed with GNOM . All used programs are included in the ATSAS software package .
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 .
The C-terminally truncated mutant proteins lacking residues 537–561 were crystallized as described for wild-type cN-II . 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) . Data processing was performed using the XDS software package . 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 . Refinement was performed with REFMAC 11.0 , including translation/libration/screwing (TLS) refinement  with 12 TLS groups at a later stage of refinement. Manual rebuilding of the model was done using Coot . The MolProbity server  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 , namely PISA , LSQKAB , and BAVERAGE. In addition, the residue interaction network was assessed with RINalyzer . Structural models were visualized using PyMOL software (DeLano Scientific). Missing loops were modeled using the ModLoop server  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 .
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 .
Crosslinked proteins were analyzed using SDS-PAGE to check that the reactions did not yield nonspecific aggregates (Additional file 14).
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.
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.
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.
The authors declare that they have no competing interests.
Open AccessThis 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.
- 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.PubMedGoogle Scholar
- Itoh R. Enzymatic properties and physiological roles of cytosolic 5 '-nucleotidase II. Curr Med Chem. 2013;20(34):4260–84.View ArticlePubMedGoogle Scholar
- Bianchi V, Spychala J. Mammalian 5'-nucleotidases. J Biol Chem. 2003;278(47):46195–8.View ArticlePubMedGoogle Scholar
- Itoh R. Regulation of cytosol 5'-nucleotidase by adenylate energy charge. Biochim Biophys Acta. 1981;659(1):31–7.View ArticlePubMedGoogle Scholar
- Itoh R. Studies on some molecular properties of cytosol 5'-nucleotidase from rat liver. Biochim Biophys Acta. 1982;716(1):110–3.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Aster JC, DeAngelo DJ. Resistance revealed in acute lymphoblastic leukemia. Nat Med. 2013;19(3):264–5.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ipata PL, Balestri F. The functional logic of cytosolic 5 '-nucleotidases. Curr Med Chem. 2013;20(34):4205–16.View ArticlePubMedGoogle Scholar
- Carlomagno T. Ligand-target interactions: what can we learn from NMR? Annu Rev Biophys Biomol Struct. 2005;34:245–66.View ArticlePubMedGoogle Scholar
- Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372(3):774–97.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Nussinov R, Ma B, Tsai CJ. Multiple conformational selection and induced fit events take place in allosteric propagation. Biophys Chem. 2014;186:22–30.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Gabizon R, Friedler A. Allosteric modulation of protein oligomerization: an emerging approach to drug design. Front Chem. 2014;2:9.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Svergun DI. Determination of the regularization parameter in indirect-transform methods using perceptual criteria. J Appl Crystallogr. 1992;25:495–503.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Kabsch W. XDS. Acta Crystallogr D. 2010;66:125–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Vagin A, Teplyakov A. Molecular replacement with MOLREP. Acta Crystallogr D. 2010;66:22–5.View ArticlePubMedGoogle Scholar
- Murshudov GN, Vagin AA, Dodson EJ. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D. 1997;53:240–55.View ArticlePubMedGoogle Scholar
- Winn MD, Isupov MN, Murshudov GN. Use of TLS parameters to model anisotropic displacements in macromolecular refinement. Acta Crystallogr D. 2001;57:122–33.View ArticlePubMedGoogle Scholar
- Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D. 2004;60:2126–32.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Fiser A, Sali A. ModLoop: automated modeling of loops in protein structures. Bioinformatics. 2003;19(18):2500–1.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Kukacka Z, Rosulek M, Strohalm M, Kavan D, Novak P. Mapping protein structural changes by quantitative cross-linking. Methods. 2015;89:112–20.View ArticlePubMedGoogle Scholar