Structural similarities and functional differences clarify evolutionary relationships between tRNA healing enzymes and the myelin enzyme CNPase
BMC Biochemistry volume 18, Article number: 7 (2017)
Eukaryotic tRNA splicing is an essential process in the transformation of a primary tRNA transcript into a mature functional tRNA molecule. 5′-phosphate ligation involves two steps: a healing reaction catalyzed by polynucleotide kinase (PNK) in association with cyclic phosphodiesterase (CPDase), and a sealing reaction catalyzed by an RNA ligase. The enzymes that catalyze tRNA healing in yeast and higher eukaryotes are homologous to the members of the 2H phosphoesterase superfamily, in particular to the vertebrate myelin enzyme 2′,3′-cyclic nucleotide 3′-phosphodiesterase (CNPase).
We employed different biophysical and biochemical methods to elucidate the overall structural and functional features of the tRNA healing enzymes yeast Trl1 PNK/CPDase and lancelet PNK/CPDase and compared them with vertebrate CNPase. The yeast and the lancelet enzymes have cyclic phosphodiesterase and polynucleotide kinase activity, while vertebrate CNPase lacks PNK activity. In addition, we also show that the healing enzymes are structurally similar to the vertebrate CNPase by applying synchrotron radiation circular dichroism spectroscopy and small-angle X-ray scattering.
We provide a structural analysis of the tRNA healing enzyme PNK and CPDase domains together. Our results support evolution of vertebrate CNPase from tRNA healing enzymes with a loss of function at its N-terminal PNK-like domain.
Eukaryotic tRNA processing is an essential process, by which newly synthesized immature pre-tRNA matures into functional tRNA. The processing of tRNA begins with endonucleolytic cleavage of the pre-tRNA into an intron and two tRNA halves that are healed and sealed by tRNA splicing. The 5′-phosphate ligation pathway, used by yeast and plants, involves the 5′-phosphate of the 3′-tRNA half as the junction phosphate of the new phosphodiester linkage [1, 2]. 5′-phosphate ligation requires three enzymatic activities: a cyclic phosphodiesterase (CPDase), a polynucleotide kinase (PNK), and a ligase [1, 3,4,5]. In contrast to yeast and plants, animal cells employ two different ligation pathways. The 3′-phosphate ligation pathway utilizes the 3′-phosphate of the 5′-tRNA half as the junction phosphate; this pathway was first detected in HeLa cell extracts . The second pathway used by animal cells is a yeast-type 5′-phosphate ligation pathway, which has been detected in HeLa cell extracts and in the lancelet, Branchiostoma floridae [7, 8]. tRNA healing enzymes appear to be distant homologues of the 2H phosphoesterase superfamily, which is defined by the presence of two conserved H-x-T/S-x (x is a hydrophobic residue) motifs [8,9,10].
The central kinase module of yeast (Saccharomyces cerevisiae) Trl1 and the N-terminal PNK domain of lancelet PNK/CPDase show resemblance to the bacteriophage T4 PNK and belong to the P-loop phosphotransferase superfamily. They contain the signature Walker A motif (P-loop), which is an NTP-binding site in most NTP-dependent phosphotransferases [8, 11,12,13,14,15]. Initial functional analysis of the kinase module of yeast Trl1, using single alanine mutations in the P-loop, revealed that the P-loop motif (401GCGKT405) is a determinant of NTP binding . GTP is the preferred in vivo physiological substrate, and yeast Trl1 contains a single NTP binding site . GTP dependence of the yeast Trl1 kinase domain has also been verified in vitro .
The CPDase domain of both yeast Trl1 and lancelet PNK/CPDase bears resemblance to phosphoesterases of the 2H superfamily [8,9,10]. Although 5′-P RNA ligases of fungi, plants, and animals share essential mechanistic features and key residues required for their PNK and CPDase activities, their overall sequence similarity is low. The presence of the P-loop or Walker A motif (G-x-G-K-T/S, x is any residue) in the N-terminal domain and the two H-x-T/S-x motifs in the C-terminal domain of vertebrate 2′,3′-cyclic nucleotide 3′-phosphodiesterase (CNPase) indicates that the PNK/CPDase domains of yeast Trl1 (Sc PNK/CPDase) and the lancelet (Bf) PNK/CPDase are homologues of vertebrate CNPase [8, 18].
In most vertebrates, CNPase is abundantly expressed in the myelin sheath, a multilayered proteolipid membrane system. CNPase constitutes 4% of total myelin protein in the central nervous system (CNS) and 0.4% in the peripheral nervous system (PNS) . CNPase-deficient mice develop progressive motor deficits and die prematurely due to diffuse brain axonal swelling and neurodegeneration . However, it is not known whether the late onset of neurodegeneration is caused by the absence of the CNPase protein as a structural component of the myelin sheath or the absence of CNPase enzymatic activity. CNPase has been implicated in neurodegenerative disorders, such as multiple sclerosis and schizophrenia [21, 22]. A recent study showed that CNPase inhibits the assembly of infectious particles of several primate lentiviruses, including the human immunodeficiency viruses HIV-1 and HIV-2, by binding to the structural protein Gag .
The enzymatic activity of CNPase, the hydrolysis of 2′,3′-cyclic nucleotides to 2′-nucleotides, was detected in brain tissue in the 1960s . Atomic structures of the C-terminal catalytic domain of human, rat, and mouse CNPase have illuminated the catalytic role of the H-x-(T/S)-x motifs and other functionally important groups in the active site [18, 25,26,27,28,29]. The structure and function of the N-terminal domain of CNPase remain poorly characterized. The expression and purification of the N-terminal domain is difficult compared to the C-terminal domain, limiting the availability and characterization of the molecule . Although the N-terminal domain of CNPase contains the NTP-binding P-loop motif, and a study by Stingo et al. claims that nucleotide binding occurs in the N-terminal domain, it has not been experimentally proven, whether the interaction and hydrolysis take place exclusively in the N-terminal domain, as only full-length CNPase was used in the experiments . The N-terminal domain of CNPase binds RNA, mediates dimerization, and interacts with calmodulin (CaM) in a calcium-dependent manner [27, 32]. CNPase has recently been found to bind microfilaments and act as an antagonist for myelin basic protein in myelin membrane compaction , suggesting a structural role – possibly independent of enzymatic CNPase activity.
Although the structure of the C-terminal phosphodiesterase domain of vertebrate CNPase has been determined, no high-resolution structures of either full-length CNPase or its N-terminal domain are available. The functional properties of vertebrate CNPase, such as interactions with membranes, cytoskeletal proteins, and RNA have been studied; however, the link between its function in the myelin sheath and its role in RNA binding is not clear. Hence, structural and functional characterization of PNK/CPDase would benefit, not only the field of tRNA splicing, but also studies on other members of the 2H phosphoesterase family. Since BfPNK/CPDase could be involved in a new class of 5′-P RNA ligation, characterization of this enzyme might improve our understanding of the growing complexity of animal tRNA splicing. Despite genetic and biochemical analysis of PNK and CPDase domains of yeast Trl1, no structural information is available for any domain of the enzyme. Therefore, detailed biophysical and biochemical analysis of ScTrl1 PNK/CPDase and BfPNK/CPDase is warranted. We generated multiple expression constructs encoding different segments of the PNK/CPDase domains of yeast Trl1, as well as full-length lancelet PNK/CPDase. An array of biophysical and biochemical methods was employed to characterize the proteins, including assays for activity, thermal stability, and folding of the PNK/CPDase domains.
Sequence and ligation independent cloning
The Protein Crystallographic Construct Design (ProteinCCD) metaserver was used to choose fragments from different domains of yeast and lancelet healing enzymes based on information acquired from various prediction servers for secondary structure, disorder, coiled coils, transmembrane segments, conserved domains, and domain linkers . Multiple expression constructs encoding the PNK and CPDase domains of yeast tRNA ligase and full-length lancelet PNK/CPDase were assembled by sequence and ligation independent cloning (SLIC). The primers used are listed in (Additional file 1). Briefly, the DNA fragments of interest were amplified by polymerase chain reaction (PCR) from pET28a and pET20b plasmids harbouring ScTrl1 and BfPNK/CPDase, respectively . The amplified products and linearized pET-his3C-LIC-amp vector (a kind gift from the Netherlands Cancer Institute) were separately treated with T4 DNA polymerase. The insert and vector were annealed, and the resulting plasmid was used to transform E. coli NEB5α cells (New England Biolabs, Germany). Transformed colonies were screened by colony PCR for recombinant plasmids that were purified, verified by sequencing, and used for protein expression.
Protein expression and purification
Large-scale expression of ScTrl1 PNK/CPDase, ScTrl1 CPDase, and BfPNK/CPDase was performed in E. coli Rosetta(DE3) cells (Novagen, Germany) cultured in LB medium containing 100 μg/ml ampicillin and 34 μg/ml chloramphenicol. After reaching an OD600 of 0.5–0.6, expression was induced with 0.15 mM isopropyl β-D-1 thiogalactopyranoside (IPTG) for 16 h at +18 °C. The cells were harvested by centrifugation and resuspended in lysis buffer (50 mM HEPES, pH 7.5, 200 mM NaCl, 10 mM β-mercaptoethanol (β-ME), 20 μg/ml DNase, Ambion RNase cocktail containing RNase A and RNase T1(Life Technologies, Germany), 5 mM MgCl2, and cOmplete mini EDTA-free protease inhibitors (Roche, Germany)). The suspension was incubated for 20 min at +4 °C on a tube rotator. The cells were disrupted by sonication, and debris was removed by centrifugation at 35,000 g for 30 min at +4 °C. The supernatant was applied to a gravity-flow Ni-NTA column, pre-equilibrated with lysis buffer. The columns were rotated horizontally for 1 h to ensure binding of the protein to the matrix. The column was washed with lysis buffer containing 50 mM imidazole, and bound protein was eluted with lysis buffer including 500 mM imidazole. The eluted fractions were studied by SDS-PAGE, and the fractions containing the protein of expected size were pooled and dialyzed against the lysis buffer (without imidazole). The N-terminal hexahistidine tag was cleaved using recombinant 3C protease at +4 °C overnight. The cleaved protein was further purified by Ni-NTA affinity chromatography and dialyzed against the lysis buffer without imidazole. The dialyzed proteins were concentrated and applied either to a HiLoad 16/60 Superdex 200 preparative grade column (column volume: 120 ml; injection volume: 1 ml) or to a Superdex 75 10/300 GL analytical grade column (column volume: 24 ml; injection volume: 100 μl), pre-equilibrated with 50 mM HEPES (pH 7.5), 200 mM NaCl, 10 mM β-ME. Equilibration, injection, and elution were all carried out at a flow rate of 1 ml/min. The protein-containing peaks were analyzed by SDS-PAGE, and the fractions containing the proteins of interest were pooled, concentrated, flash-frozen, and stored at −80 °C. The identity of the purified proteins was verified by tryptic peptide mapping using mass spectrometry at the Biocenter Oulu Proteomics Core Facility.
Full-length mouse CNPase (MmCNPase, residues 20–420 of isoform 2), the N-terminal PNK domain (MmCNP_N, residues 20–180), and the C-terminal catalytic CPDase domain (MmCNP_C, residues 179–398) were purified using Ni-NTA chromatography, followed by His-tag cleavage using TEV protease, a second Ni-NTA step, and size exclusion chromatography (SEC), essentially as described before .
Synchrotron radiation circular dichroism spectroscopy
The folding of ScTrl1 PNK/CPDase, BfPNK/CPDase, and ScTrl1 CPDase was studied by SRCD spectroscopy on the UV-CD12 beamline at the ANKA Synchrotron. The proteins purified by SEC were dialyzed into a buffer compatible with SRCD analysis (50 mM potassium phosphate, pH 6.5). The range of protein concentrations used was between 3.5 and 5 mg/ml. Three spectra were collected at a scan rate of 14 nm/min, in a demountable 13-μm CaF2 cuvette at a wavelength range between 260 and 170 nm at 1 nm intervals. A spectrum from the corresponding buffer was measured for all the samples and used for background correction. The data were processed with CDtool  and analyzed using Dichroweb .
PNK activity assay
The PNK activity of the N-terminal PNK domain of ScTrl1 PNK/CPDase and BfPNK/CPDase was assayed using the synthetic oligonucleotides ribo A20 (A20) and deoxy A20 (dA20), and a mixture of the two (A20 + dA20). T4 PNK and ScTrl1 CPDase were used as positive and negative controls, respectively. The reaction mixture contained 70 mM Tris, pH 7.5, 10 mM MgCl2, and 5 mM DTT, 2 μg of enzyme, and 100 pmol A20 or dA20 or both. 55 pmol of radioactively labelled [γ32P] ATP were added to each reaction mixture, and the volume of the mixtures was adjusted to 40 μl with diethylpyrocarbonate-treated water. The samples were incubated at +37 °C, and 8-μl aliquots were removed at different time points (1, 15, 30, and 60 min). Each sample was quenched with 8 μl of 2X urea sample loading buffer (Invitrogen Novex) and heated at +90 °C for 4 min. The samples were analyzed by electrophoresis using a 15% Mini-PROTEAN TBE-Urea gel (Bio-Rad), at 150 V for 85 min, and visualized using PhosphorImager analysis.
CPDase activity assay
Kinetic measurements were carried out as previously described [30, 37]. In this assay, CPDase hydrolyzes the phosphodiester bond in β-nicotinamide adenine dinucleotide 2′,3′-cyclic monophosphate (2′,3′-cNADP+), and the resulting nicotinamide adenine dinucleotide phosphate (NADP+) is reduced to NADPH by glucose-6-phosphate dehydrogenase and used to transform glucose-6-phosphate to 6-phosphoglucanolactone. The quantity of NADPH formed during the coupled enzymatic reaction is a direct measure of CPDase activity .
The assay mixture contained 100 mM MES, pH 6.0, 3 mM MgCl2, 5 mM glucose-6-phosphate, and 0.6 U glucose-6-phosphate dehydrogenase. 500 ng of the purified enzyme was mixed with varying concentrations of the substrate, 2′,3′-cNADP+ (0, 0.05, 0.1, 0.2, 0.5, and 1.0 mM), and NADPH production was measured spectrophotometrically at 340 nm, using the absorption coefficient 6.22 cm−1mM−1. The measurements were carried out in triplicate at +25 °C using 96-well flat-bottom transparent plates with the TECAN Infinite M200 fluorescence spectrophotometer and iControl software.
Thermal stability assay
Fluorescence-based thermal shift assays were carried out in triplicate, essentially as described . Measurements were performed with a CFX96 RealTime PCR system (BioRad), using a 96-well thin-wall PCR plate. 1–5 μg of protein and 1–2 X SYPRO Orange dye were included in a total volume of 25 μl per condition. The plate was sealed with optical-quality adhesive film (BioRad) and heated from +25 °C to +99 °C in 0.5 °C increments. A total of 80 different conditions, varying in pH, salt type and concentration, and the presence of ligands, were screened. The fluorescence of SYPRO Orange was measured using excitation at 490 nm and emission at 575 nm. The apparent thermal melting points (Tm) were determined for each construct under all 80 conditions.
Small-angle X-ray scattering
SAXS measurements were carried out on the EMBL beamline X33 at DESY, Hamburg (Germany), the I911-4 SAXS beamline of MAX-Lab in Lund (Sweden), and the P12 BioSAXS beamline at PETRA-III, DESY, Hamburg. The ScTrl1 PNK/CPDase and BfPNK/CPDase eluted from SEC both as a dimer and as a monomer. These fractions were collected separately, concentrated, and used for SAXS measurements; only the data from the monomeric fraction were analyzed in detail, as dimerization was likely an artifact of non-specific disulphide formation. Sample concentrations were 1–10 mg/ml. Monomeric BSA was measured first as a MW standard. Solvent scattering from the corresponding buffer was measured identically before and after each sample, and the average background scattering was subtracted. The data were analyzed with the ATSAS  suite, as described [41, 42]. The data were processed using PRIMUS . Distance distributions were calculated using GNOM , and ab initio bead modeling was done with DAMMIF . DAMAVER  was used for model averaging. GASBOR  was used for building chain-like ab initio models, and MONSA  was used to assemble a 2-phase model of ScTrl1 PNK/CPDase, employing data from both ScTrl1 PNK/CPDase and ScTrl1 CPDase. The models were compared in PyMOL.
Results and discussion
Protein expression and purification
A total of nine expression clones including eight for different segments of yeast Trl1 and one for lancelet PNK/CPDase were prepared (Fig. 1). All constructs contain an N-terminal hexahistidine tag followed by a 3C protease cleavage site and the corresponding PNK/CPDase domain. The expression level of the different constructs was screened in three E. coli cell lines [BL21(DE3), BL21(DE3) CodonPlus RIPL, and Rosetta(DE3)] under different expression conditions, including various inducer concentrations and post-induction growth temperatures. The constructs encoding different regions of the Sc PNK domain did not show any expression in the three cell lines tested; the remaining five constructs could be expressed. The same pattern was observed previously for the yeast ScTrl1 PNK domain and the N-terminal domain of mouse CNPase [16, 30]. The solubility of the proteins expressed in Rosetta (DE3) was remarkably high. Screening growth conditions post-induction revealed longer incubation at lower temperatures to be the optimal expression strategy. For large-scale production of ScTrl1 PNK/CPDase and BfPNK/CPDase, up to 8 l of culture were used to obtain ~5 mg of protein, whereas 2 l of culture was sufficient to achieve similar quantities of the ScCPDase domains.
The expressed proteins were initially purified by Ni-NTA affinity chromatography and SEC (Additional file 2 & Figs. 2 and 3). ScTrl1 PNK/CPDase and BfPNK/CPDase eluted as double peaks with elution volumes corresponding to dimeric and monomeric forms. The fractions from each of the two peaks were collected separately, concentrated, and analyzed by SDS PAGE. Both peaks contained pure protein. The vertebrate homologue, mouse CNPase, exhibits similar oligomeric behaviour, and in this case, the N-terminal RNA-binding domain mediated dimerization . The addition of either 10 mM β-ME or 5 mM DTT in the lysis and purification buffers inhibited dimerization of ScTrl1 PNK/CPDase and BfPNK/CPDase, indicating that the dimers are linked by a disulphide. Most likely the physiologically relevant form of each protein is monomeric.
The elution profiles of the CPDase constructs contained a single symmetric peak (Fig. 3b, c), and no oligomerization or RNA binding was detected. The absence of dimerization and RNA binding in the CPDase domain suggest that these processes require the N-terminal PNK domain. Similarly, the N-terminal PNK-like domain of mouse CNPase is involved not only in dimerization, but also in RNA binding . The purified proteins were easily degraded upon storage at +4 °C. Hence, immediate flash-freezing and storage of the purified proteins in the freezer was essential. Since the stability of ScTrl1 CPDase was found to be comparatively better than the other two CPDase protein variants, it was used in further biophysical and biochemical experiments.
PNK/CPDase enzymes interact with Eschericia coli nucleic acids
During the purification of ScTrl1 PNK/CPDase and BfPNK/CPDase, a large peak was always present right after the void volume of the SEC column (Fig. 4a). Since the absorbance at 260 nm was very high for this peak, co-purification of a nucleic acid, possibly RNA, from the expression host, was suspected. To identify the nucleic acid bound to the protein, the fractions eluted around the void volume in SEC were treated with DNase and RNase A. Agarose gels of nuclease-treated samples show that the bright smear in the larger elution fraction and in the same fraction treated with DNase was not found in the presence of RNase A (Fig. 4b). This finding confirms that the co-purified nucleic acid impurity contained RNA. In addition, DNA was present in the BfPNK/CPDase sample. When a combination of RNase A and RNase T1 was added in the lysis buffer during subsequent purification, the A260 peak height remarkably decreased, additionally proving that the contaminant was RNA.
The behaviour of PNK/CPDase proteins in these assays was highly similar to the uracil-DNA degrading factor and mouse CNPase [27, 49]. CNPase has also been shown to interact with RNA in vitro and co-purify with poly(A)+ RNA; the catalytic domain of CNPase has been shown to be sufficient for binding with single stranded RNA homopolymers [50, 51]. However, a pulldown assay with mouse CNPase and poly(A)-sepharose indicated that the N-terminal domain binds RNA more efficiently than the C-terminal domain . The absence of the peak closely following the void volume during SEC purification of the yeast CPDase domain indicates the N-terminal PNK domain is required for RNA binding, at least in the case of E. coli RNA.
Bacterially expressed PNK/CPDase and CPDase proteins are folded
The folding state of the purified proteins was analyzed by SRCD spectroscopy. Visual inspection of the SRCD spectra indicates that the PNK/CPDase and CPDase proteins were folded, containing a mixture of α helix and β strand. SRCD spectra from ScTrl1 PNK/CPDase and BfPNK/CPDase show a similar shape (Fig. 5a), indicating the presence of similar amounts of regular secondary structure. The spectra also indicate that the N-terminal truncation of ScTrl1 CPDase is less folded than ScTrl1 CPDase (Fig. 5b). The latter is indicative of sub-optimal folding of the truncated construct.
Ligand binding improves thermal stability
Since the stability of the monomers of both ScTrl1 PNK/CPDase and BfPNK/CPDase during and after purification was initially low, the identification of a stabilizing buffer condition was necessary. A thermal stability assay was used to identify suitable buffer conditions that would offer increased stability of the protein. The melting curves and corresponding Tm values of ScTrl1 PNK/CPDase, ScTrl1 CPDase, and BfPNK/CPDase indicate that all three proteins are most stable around pH 7.5. The combination of 50 mM HEPES, pH 7.5, and 150 mM NaCl was identified as an optimal stabilizing buffer.
To analyze the effect of ligands on the thermal stability of the proteins, two known substrates of CPDase, 2′,3′-cCMP and 2′,3′-cNADP+, and the product NADP+ were included in the screen. The Tm values indicate that 2′,3′-cNADP+ improves the thermal stability of all three proteins tested (Fig. 6, Table 1); however, 2′,3′-cCMP and NADP+ do not influence the stability of any of them. The highest melting points obtained for ScTrl1 PNK/CPDase and BfPNK/CPDase were +63 °C and +48 °C, respectively, both in the presence of 2′,3′-cNADP+. ScTrl1 CPDase was found to be the most stable among the proteins tested, with a melting temperature of +72 °C in a buffer containing 50 mM HEPES, pH 7.5, 150 mM NaCl, and 1 mM 2′,3′-cNADP+. Mouse CNPase is more stable at pH 5.5 with higher salt than at neutral pH with lower salt, which stabilizes the yeast and lancelet proteins . The thermal stability of ScTrl1 CPDase is also higher than that of the mouse CNPase catalytic domain and its mutated variants [28, 29].
PNK/CPDase enzymes possess polynucleotide kinase activity
The assay to test the PNK activity of the purified ScTrl1 PNK/CPDase and BfPNK/CPDase involves transfer of the 32P-labelled terminal phosphate from [γ32P] ATP to a 20-mer 5′-OH synthetic RNA (A20) or DNA oligonucleotide (dA20). To study the substrate preferences of the enzymes on RNA and DNA oligonucleotide substrates, the transfer of the radiolabelled phosphate to a mixture of A20 and dA20 was also analyzed. The time-dependent transfer of the phosphate from [γ32P] ATP to the 5′-OH groups of the synthetic oligonucleotides indicates that both ScTrl1 PNK/CPDase and BfPNK/CPDase possess PNK activity (Fig. 7a-b). ScTrl1 CPDase did not catalyze phosphoryl transfer (Fig. 7c). The results indicate that ScTrl1 PNK/CPDase prefers A20 (Fig. 7a) and BfPNK/CPDase strongly prefers dA20 (Fig. 7b).
BfPNK/CPDase was previously shown to prefer DNA over RNA, whereas another lancelet PNK, Bf Clp1, was reported to act exclusively on RNA . This finding suggested a possible role for BfPNK/CPDase in DNA repair. The fact that we saw DNA co-purifying with BfPNK/CPDase and that the enzyme preferred the DNA substrate in the PNK assay support this hypothesis. The decrease in size of the labelled dA20 oligonucleotide during the assay could be related to an unknown activity of BfPNK/CPDase towards DNA. The homologue of yeast and lancelet PNK/CPDase, mammalian CNPase, was found to rescue a yeast strain with an inactivating mutation in the CPDase domain of yeast tRNA ligase, but not a strain with a mutated kinase domain . Thus, it appears that the N-terminal PNK-like domain of vertebrate CNPase is not functionally a PNK. This is supported by our results from the PNK activity assays. MmCNPase, as well as its two domains separately, were inactive, while the homologous yeast and lancelet PNK/CPDase proteins were active (Fig. 7a, b, Additional file 3). This finding also supports the hypothesis that vertebrates might employ the 5′-P RNA ligation pathway only as an alternate to 3′-P RNA ligation and, thus, require only the exclusive polynucleotide kinase (Clp1) for the healing reaction, and not the N-terminal PNK-like domain of CNPase.
PNK/CPDase and CPDase enzymes possess cyclic nucleotide phosphodiesterase activity
The CPDase activity of ScTrl1 PNK/CPDase, BfPNK/CPDase, and ScTrl1 CPDase was assayed using 2′,3′-cNADP+ as substrate in a coupled enzyme assay [37, 38]. The results indicate that all tested constructs are active in the phosphodiesterase reaction (Table 2) and that the activity is comparable to that of the different constructs of mouse and rat CNPase. Thus, CPDase activity is a common denominator for 2H enzymes across kingdoms of life.
PNK/CPDase enzymes are elongated in solution
SAXS measurements were carried out to determine low-resolution structural models of the proteins and to analyze their oligomeric state in solution. Although SAXS data were collected for all the five expressed constructs, only three of them, the ScTrl1 PNK/CPDase, the ScTrl1 CPDase, and the BfPNK/CPDase, yielded good-quality SAXS data, and only these were used for further processing and modeling. The best datasets were collected between 3 and 4 mg/ml, with a strong scattering signal and without significant interparticle effects, as shown by the linearity of the Guinier plot (Fig. 8a, b).
The monomeric ScTrl1 PNK/CPDase and BfPNK/CPDase display similarly elongated conformations, with BfPNK/CPDase being more extended (Fig. 8a, e, f and Table 3). The estimated Porod volume is in agreement with a monomeric form of both proteins. The p(r) profile displays a single peak with shoulders for both yeast and lancelet PNK/CPDase proteins (Fig. 8c), indicating an asymmetric scattering particle, possibly containing two domains bridged by a flexible linker. Dimensionless Kratky plots (Fig. 8d) further show that BfPNK/CPDase is the most flexible of the studied samples, in line with its relatively large volume-to-mass ratio (Table 3).
The 3D models for both PNK/CPDases are consistent with an arrangement, where the two domains form an elongated assembly. The model suggests that the active sites of the PNK domain and CPDase domain are accessible and not being blocked by the neighbouring domain. An open active site facilitates substrate binding without the need for large conformational changes. If the domains would cover each other more extensively, the shape of the molecule would be more globular.
The solution structure of full-length mouse CNPase is similar to that of the monomeric ScTrl1 PNK/CPDase and BfPNK/CPDase , although somewhat more compact, reflecting the presence of large flexible insertions in the yeast and lancelet enzymes. Mouse CNPase is a monomer with an elongated conformation in solution . Although the molecular masses of monomeric ScTrl1 PNK/CPDase [50.52 kDa] and monomeric mouse CNPase [44.8 kDa] are close to each other, the D max and R g values are comparatively higher for the ScTrl1 PNK/CPDase (Table 3). This could be due to the presence of a flexible linker region between the N-terminal PNK and C-terminal CPDase domains in ScTrl1 PNK/CPDase: this segment is absent in mouse CNPase (Fig. 9a). Based on SAXS data and molecular modelling, the C-terminal membrane-anchoring tail of the mouse CNPase lies in the middle region of the molecule, enabling the association of the active site to the vicinity of the lipid bilayer .
CPDase is compact in solution
The molecular size determined from the SAXS profile of ScTrl1 CPDase indicates that the molecule is more compact in solution than ScTrl1 PNK/CPDase and exists as a monomer (Fig. 8a, c, d, f & g, Table 3). Molecular mass determination shows the presence of a monomer. The Porod volume is also in line with the presence of monomeric ScTrl1 CPDase. The SAXS model shows the presence of an extension presumably in the N terminus of the protein; this tail possibly arises from the interdomain linker at the N terminus of the construct, suggesting the requirement for a shorter construct of the protein without the N terminal insertion to facilitate crystallization.
The solution structure of ScTrl1 CPDase resembles that of the catalytic domain of mouse CNPase . The mouse CNPase catalytic domain is also monomeric in solution. Both ScTrl1 CPDase and the catalytic domain of mouse CNPase adopt an elongated conformation possibly due to the opening of the active-site cleft, as found in the NMR structure of rat CNPase catalytic domain [25, 27]. The mouse CNPase catalytic domain displays a more compact structure in the presence of citrate and the ligand also reduces aggregation of the protein . It has been suggested that the overall solution conformation of fully active CNPase might be more open than that seen in the crystal structures, in which the enzyme is bound to citrate and sulphate . Different lengths of the mouse CNPase catalytic domain have been studied by SAXS . The presence of a C-terminal tail does not alter the conformation of the catalytic domain in solution, whereas N-terminally extended CNPase catalytic domain has a higher radius of gyration. Both N- and C-terminally extended protein is remarkably more elongated . The symmetric peak of the p(r) for ScTrl1 CPDase shows that the protein, similarly to the mouse CNPase catalytic domain, is more compact than the protein variants that also contain the N-terminal domain (Fig. 8c).
The solution shapes and behaviour of the PNK/CPDases and the CPDase are in many aspects comparable to those of the full-length mouse CNPase and its catalytic domain, respectively. The differences are likely caused by the relatively large insertions compared to CNPase: the yeast enzyme has a long insertion between the PNK and CPDase domains, while the lancelet enzyme has insertions in two loop regions of the CPDase domain (Fig. 9a, b). The SAXS model of full-length mouse CNPase indicates that the N-terminal PNK-like and the C-terminal phosphodiesterase domains form an elongated assembly , and as shown here, the same is true for the yeast and lancelet enzymes. The biphasic model built for the yeast enzyme (Fig. 8h) compares well to the mouse CNPase solution structure. Thus, the domain arrangements of the yeast and lancelet PNK/CPDase in solution resemble that of CNPase. An open elongated conformation of the full-length protein could be crucial for substrate binding . Superposition of a homology model of the CPDase domain of BfPNK/CPDase and the catalytic domain of MmCNPase shows that the active site consisting of the 2H motifs is structurally conserved [27, 53] (Fig. 9b). The homology model of the lancelet CPDase domain contains a symmetrical bilobed tertiary structure, which is similar to the crystal structures of the catalytic domain of mouse and human CNPase [26,27,28,29]. The inserted loops within the CPDase domain of BfPNK/CPDase (Fig. 9b) were also the most flexible in the recently-determined structure of a bacterial 2H phosphoesterase, LigT from E. coli .
Comparison of the SAXS models of ScTrl1 PNK/CPDase, BfPNK/CPDase, and ScTrl1 CPDase with the SAXS structural parameters of the homologous mouse CNPase indicates that the structural and hydrodynamic properties of the different constructs can be explained by the presence of sequence insertions in the yeast and lancelet enzymes. A structural model of the PNK/CPDase proteins (Fig. 9c) can be proposed based on the determined 3-dimensional shapes. Yeast PNK/CPDase exhibits an elongated conformation with the presence of a possibly flexible region between the N-terminal PNK and C-terminal CPDase domains. Removal of the flexible region, which is not present in mouse CNPase, might aid in crystallization of the domains separately. The model also indicates that the two loop insertions in the CPDase domain in lancelet PNK/CPDase result in an even more elongated molecule, and these long loops point away from the PNK domain at the other end of the molecule. It is possible that BfPNK/CPDase also has a flexible region at its N terminus.
The open elongated conformation of the PNK/CPDases might play a role in RNA substrate binding; however, this needs to be confirmed by high-resolution structures of the enzymes with bound substrates. Based on currently available structural data, we can conclude that both enzymes display structural similarities to other members of the 2H-phosphoesterase family - not only at the level of conserved sequence motifs, but also in their respective domain arrangements. The yeast and lancelet healing enzymes possess both PNK and CPDase activities, whereas the mouse CNPase, which also binds RNA, has no PNK activity. This suggests a loss of function at the CNPase N-terminal PNK-like domain over the course of evolution. On the other hand, BfPNK/CPDase may have a unique substrate specificity within the family with respect to its PNK activity. Thus, our results on the structure and function of tRNA healing enzymes from yeast and lancelet provide evidence supporting the evolution of vertebrate CNPase from functional tRNA healing enzymes.
- Bf :
Bovine serum albumin
cyclic cytosine monophosphate
2’,3’-cyclic nucleotide 3’-phosphodiesterase
Central nervous system
- D max :
maximum dimension of a particle
- E. coli :
European Molecular Biology Laboratory
Hydroxyethyl piperazineethanesulfonic acid
Human immunodeficiency virus
- Hs :
- Mm :
- NADP+ :
Oxidized form of nicotinamide adenine dinucleotide phosphate
Reduced form of nicotinamide adenine dinucleotide phosphate
Polymerase chain reaction
Protein Data Bank
Peripheral nervous system
- R g :
Radius of gyration
Small-angle X-ray scattering
- Sc or S. cerevisiae :
Sodium dodecyl sulphate-polyacrylamide gel electrophoresis
Size exclusion chromatography
Sequence and ligation independent cloning
Synchrotron radiation circular dichroism
tRNA ligase 1
transfer ribonucleic acid
Englert M, Beier H. Plant tRNA ligases are multifunctional enzymes that have diverged in sequence and substrate specificity from RNA ligases of other phylogenetic origins. Nucleic Acids Res. 2005;33:388–99.
Greer CL, Peebles CL, Gegenheimer P, Abelson J. Mechanism of action of a yeast RNA ligase in tRNA splicing. Cell. 1983;32:537–46.
Abelson J, Trotta CR, Li H. tRNA Splicing. J Biol Chem. 1998;273:12685–8.
Wang LK, Shuman S. Structure-function analysis of yeast tRNA ligase. RNA. 2005;11:966–75.
Wang LK, Schwer B, Englert M, Beier H, Shuman S. Structure-function analysis of the kinase-CPD domain of yeast tRNA ligase (Trl1) and requirements for complementation of tRNA splicing by a plant Trl1 homolog. Nucleic Acids Res. 2006;34:517–27.
Filipowicz W, Shatkin AJ. Origin of splice junction phosphate in tRNAs processed by HeLa cell extract. Cell. 1983;32:547–57.
Zillmann M, Gorovsky MA, Phizicky EM. Conserved mechanism of tRNA splicing in eukaryotes. Mol Cell Biol. 1991;11:5410–6.
Englert M, Sheppard K, Gundllapalli S, Beier H, Söll D. Branchiostoma floridae has separate healing and sealing enzymes for 5’-phosphate RNA ligation. Proc Natl Acad Sci U S A. 2010;107:16834–9.
Nasr F, Filipowicz W. Characterization of the Saccharomyces cerevisiae cyclic nucleotide phosphodiesterase involved in the metabolism of ADP-ribose 1“,2”-cyclic phosphate. Nucleic Acids Res. 2000;28:1676–83.
Mazumder R, Iyer LM, Vasudevan S, Aravind L. Detection of novel members, structure-function analysis and evolutionary classification of the 2H phosphoesterase superfamily. Nucleic Acids Res. 2002;30:5229–43.
Walker JE, Saraste M, Runswick M, Gay NJ. Distantly related sequences in the alpha- and beta-subunits of ATP synthase, myosin, kinases and other ATP-requiring enzymes and a common nucleotide binding fold. EMBO J. 1982;1:945–51.
Wang LK, Shuman S. Domain Structure and Mutational Analysis of T4 Polynucleotide Kinase. J Biol Chem. 2001;276:26868–74.
Galburt EA, Pelletier J, Wilson G, Stoddard BL. Structure of a tRNA repair enzyme and molecular biology workhorse: T4 polynucleotide kinase. Structure. 2002;10:1249–60.
Wang LK, Lima CD, Shuman S. Structure and mechanism of T4 polynucleotide kinase: An RNA repair enzyme. EMBO J. 2002;21:3873–80.
Wang LK, Shuman S. Mutational analysis defines the 5’-kinase and 3’-phosphatase active sites of T4 polynucleotide kinase. Nucleic Acids Res. 2002;30:1073–80.
Sawaya R, Schwer B, Shuman S. Genetic and biochemical analysis of the functional domains of yeast tRNA ligase. J Biol Chem. 2003;278:43928–38.
Westaway SK, Belford HG, Apostol BL, Abelson J, Greer CL. Novel activity of a yeast ligase deletion polypeptide: Evidence for GTP-dependent tRNA splicing. J Biol Chem. 1993;268:2435–43.
Myllykoski M, Seidel L, Muruganandam G, Raasakka A, Torda AE, Kursula P. Structural and functional evolution of 2′,3′-cyclic nucleotide 3′-phosphodiesterase. Brain Res. 2016;1641(A):64–78.
Bifulco M, Laezza C, Stingo S, Wolff J. 2’,3’-Cyclic nucleotide 3’-phosphodiesterase: a membrane-bound, microtubule-associated protein and membrane anchor for tubulin. Proc Natl Acad Sci U S A. 2002;99:1807–12.
Lappe-Siefke C, Goebbels S, Gravel M, Nicksch E, Lee J, Braun PE, et al. Disruption of Cnp1 uncouples oligodendroglial functions in axonal support and myelination. Nat Genet. 2003;33:366–74.
Muraro PA, Kalbus M, Afshar G, McFarland HF, Martin R. T cell response to 2′,3 ′ -cyclic nucleotide 3′ -phosphodiesterase (CNPase) in multiple sclerosis patients. J Neuroimmunol. 2002;130:233–42.
Peirce TR, Bray NJ, Williams NM, Norton N, Moskvina V, Preece A, et al. Convergent evidence for 2’,3’-cyclic nucleotide 3’-phosphodiesterase as a possible susceptibility gene for schizophrenia. Arch Gen Psychiatry. 2006;63:18–24.
Wilson SJ, Schoggins JW, Zang T, Kutluay SB, Jouvenet N, Alim MA, et al. Inhibition of HIV-1 particle assembly by 2’,3’-cyclic-nucleotide 3’-phosphodiesterase. Cell Host Microbe. 2012;12:585–97.
Drummond GI, Iyer NT, Keith J. Hydrolysis of Ribonucleoside 2’,3’-Cyclic Phosphates by a Diesterase from Brain. J Biol Chem. 1962;237:3535–9.
Kozlov G, Lee J, Elias D, Gravel M, Gutierrez P, Ekiel I, et al. Structural Evidence That Brain Cyclic Nucleotide Phosphodiesterase Is a Member of the 2H Phosphodiesterase Superfamily. J Biol Chem. 2003;278:46021–8.
Sakamoto Y, Tanaka N, Ichimiya T, Kurihara T, Nakamura KT. Crystal structure of the catalytic fragment of human brain 2′,3′ -cyclic-nucleotide 3′ -phosphodiesterase. J Mol Biol. 2005;346:789–800.
Myllykoski M, Raasakka A, Han H, Kursula P. Myelin 2′,3′ -cyclic nucleotide 3′ -phosphodiesterase: Active-site ligand binding and molecular conformation. PLoS One. 2012;7(2):e32336.
Myllykoski M, Raasakka A, Lehtimäki M, Han H, Kursula I, Kursula P. Crystallographic analysis of the reaction cycle of 2′,3′-cyclic nucleotide 3′-phosphodiesterase, a unique member of the 2H phosphoesterase family. J Mol Biol. 2013;425:4307–22.
Raasakka A, Myllykoski M, Laulumaa S, Lehtimäki M, Härtlein M, Moulin M, et al. Determinants of ligand binding and catalytic activity in the myelin enzyme 2’,3’-cyclic nucleotide 3’-phosphodiesterase. Sci Rep. 2015;5:16520.
Myllykoski M, Kursula P. Expression, purification, and initial characterization of different domains of recombinant mouse 2’,3’-cyclic nucleotide 3’-phosphodiesterase, an enigmatic enzyme from the myelin sheath. BMC Res Notes. 2010;3:12.
Stingo S, Masullo M, Polverini E, Laezza C, Ruggiero I, Arcone R, et al. The N-terminal domain of 2’,3’-cyclic nucleotide 3’-phosphodiesterase harbors a GTP/ATP binding site. Chem Biol Drug Des. 2007;70:502–10.
Myllykoski M, Itoh K, Kangas SM, Heape AM, Kang SU, Lubec G, et al. The N-terminal domain of the myelin enzyme 2’,3’-cyclic nucleotide 3’-phosphodiesterase: Direct molecular interaction with the calcium sensor calmodulin. J Neurochem. 2012;123:515–24.
Snaidero N, Velte C, Myllykoski M, Raasakka A, Ignatev A, Werner HB, et al. Antagonistic Functions of MBP and CNP Establish Cytosolic Channels in CNS Myelin. Cell Rep. 2017;18:314–23.
Mooij WTM, Mitsiki E, Perrakis A. ProteinCCD: Enabling the design of protein truncation constructs for expression and crystallization experiments. Nucleic Acids Res. 2009;37:W402–5.
Lees JG, Smith BR, Wien F, Miles AJ, Wallace BA. CDtool—an integrated software package for circular dichroism spectroscopic data processing, analysis, and archiving. Anal Biochem. 2004;332:285–9.
Lobley A, Whitmore L, Wallace BA. DICHROWEB: an interactive website for the analysis of protein secondary structure from circular dichroism spectra. Bioinformatics. 2002;18:211–2.
Lee J, Gravel M, Gao E, O’Neill RC, Braun PE. Identification of essential residues in 2’,3’-cyclic nucleotide 3’-phosphodiesterase. Chemical modification and site-directed mutagenesis to investigate the role of cysteine and histidine residues in enzymatic activity. J Biol Chem. 2001;276:14804–13.
Sogin DC. 2’,3’-cyclic NADP as a substrate for 2’,3’-cyclic nucleotide 3’-phosphohydrolase. J Neurochem. 1976;27:1333–7.
Ericsson UB, Hallberg BM, DeTitta GT, Dekker N, Nordlund P. Thermofluor-based high-throughput stability optimization of proteins for structural studies. Anal Biochem. 2006;357:289–98.
Konarev PV, Petoukhov MV, Volkov VV, Svergun DI. ATSAS 2.1, a program package for small-angle scattering data analysis. J Appl Crystallogr. 2006;39:277–86.
Majava V, Petoukhov MV, Hayashi N, Pirilä P, Svergun DI, Kursula P. Interaction between the C-terminal region of human myelin basic protein and calmodulin: analysis of complex formation and solution structure. BMC Struct Biol. 2008;8:10.
Majava V, Wang C, Myllykoski M, Kangas SM, Kang SU, Hayashi N, et al. Structural analysis of the complex between calmodulin and full-length myelin basic protein, an intrinsically disordered molecule. Amino Acids. 2010;39:59–71.
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.
Svergun DI. Determination of the regularization parameter in indirect-transform methods using perceptual criteria. J Appl Crystallogr. 1992;25:495–503.
Franke D, Svergun DI. DAMMIF, a program for rapid ab-initio shape determination in small-angle scattering. J Appl Crystallogr. 2009;42:342–6.
Volkov VV, Svergun DI. Uniqueness of ab initio shape determination in small-angle scattering. J Appl Crystallogr. 2003;(36):860–4.
Svergun DI, Petoukhov MV, Koch MH. Determination of domain structure of proteins from X-ray solution scattering. Biophys J. 2001;80:2946–53.
Svergun DI. Restoring Low Resolution Structure of Biological Macromolecules from Solution Scattering Using Simulated Annealing. Biophys J. 1999;76:2879–86.
Bekesi A, Pukancsik M, Haasz P, Felfoldi L, Leveles I, Muha V, et al. Association of RNA with the uracil-DNA-degrading factor has major conformational effects and is potentially involved in protein folding. FEBS J. 2011;278:295–315.
Gravel M, Robert F, Kottis V, Gallouzi I-E, Pelletier J, Braun PE. 2′,3′-Cyclic nucleotide 3′-phosphodiesterase: A novel RNA-binding protein that inhibits protein synthesis. J Neurosci Res. 2009;87:1069–79.
Myllykoski M, Kursula P. Structural aspects of nucleotide ligand binding by a bacterial 2H phosphoesterase. PLoS One. 2017;12:e0170355.
Schwer B, Aronova A, Ramirez A, Braun P, Shuman S. Mammalian 2’,3’ cyclic nucleotide phosphodiesterase (CNP) can function as a tRNA splicing enzyme in vivo. RNA. 2008;14:204–10.
Kelley LA, Sternberg MJE. Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc. 2009;4:363–71.
Notredame C, Higgins DG, Heringa J, Notredame C, Higgins DG, Heringa J. T-coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302:205–17.
We extend our special thanks to the beamline staff at EMBL-Hamburg/DESY, MAX-Lab, and ANKA, as well as to staff at the Biocenter Oulu Proteomics and Protein Analysis Core Facility.
This study has been supported by the Academy of Finland, the German Federal Ministry for Education and Research, and the Research and Science Foundation of the City of Hamburg (Germany). The funders had no role in the design of the study, the collection, analysis, and interpretation of data, or in writing the manuscript.
Availability of data and materials
All data generated or analysed during this study are included in the published article and its supplementary information files. Raw data files are available from the corresponding author on reasonable request.
GM planned and carried out experiments, analysed data, and was a major contributor in writing the manuscript. AR planned and carried out experiments, analyzed data, and participated in writing the manuscript. MM planned and carried out experiments, analysed data, and participated in writing the manuscript. IK obtained funding for the study, planned experiments, analysed data, and participated in writing the manuscript. PK obtained funding for the study, planned experiments and analysed data, and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers used for PCR amplification of target DNA fragments. (PDF 54 kb)
The yeast PNK/CPDase and CPDase proteins purified by Ni-NTA affinity chromatography. A) SDS-PAGE from a Ni-NTA purification of ScTrl1 PNK/CPDase. M: marker, E1-E8: fractions eluted with 500 mM imidazole. The calculated size of the protein was 52.7 kDa. B) SDS-PAGE from Ni-NTA purifications of ScTrl1 CPDase and ScTrl1 (N+) CPDase. M: marker, SF1 and W1 are supernatant and wash flow-through samples from the purification of ScTrl1 CPDase, E1-E5: ScTrl1 CPDase fractions eluted with 500 mM imidazole. The calculated size of the protein was 32.2 kDa. SF2 and W2 are supernatant and wash flow-through samples from the purification of ScTrl1 (N+) CPDase, E6-E10: ScTrl1 (N+) CPDase fractions eluted with 500 mM imidazole. The calculated size of the protein was 32.4 kDa. C) SDS PAGE from a Ni-NTA purification of ScTrl1 (N-) CPDase. SF: supernatant flow-through, W: wash flow-through, E1-E10: fractions eluted with 500 mM imidazole, M: marker. The calculated size of the protein was 29.6 kDa. (PDF 220 kb)
Polynucleotide kinase activity assay. PNK reaction mixtures of A) T4 PNK, B) MmCNPase C) MmCNP_N and D) MmCNP_C. The reactions were carried out in the presence of Mg2+, with ATP as phosphate donor and with either A20, dA20, or both. T4 PNK was used as a positive control to validate the assay setup. The MmCNP domains were tested separately, with a negative result in the chosen assay conditions. The arrow represents the direction of electrophoresis. (PDF 387 kb)
About this article
Cite this article
Muruganandam, G., Raasakka, A., Myllykoski, M. et al. Structural similarities and functional differences clarify evolutionary relationships between tRNA healing enzymes and the myelin enzyme CNPase. BMC Biochem 18, 7 (2017). https://doi.org/10.1186/s12858-017-0084-2
- Cyclic phosphodiesterase
- Polynucleotide kinase
- Substrate specificity
- Protein structure
- 2H family
- tRNA splicing