|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 283, Issue 6, 3655-3664, February 8, 2008
Structural and Functional Analyses of the Severe Acute Respiratory Syndrome Coronavirus Endoribonuclease Nsp15*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
10% (as shown by the World Health Organization, www.who.int/csr/sars/country/en/country2003_08_15.pdf). In addition to their medical importance, coronaviruses are of interest for their large
30-kb positive-strand genome and novel mechanisms that have evolved to replicate and transcribe this large RNA (5, 6). In keeping with the novel strategies used, coronaviruses encode several unusual RNA-processing enzymes, including an RNA endoribonuclease, an RNA methyltransferase, a second RNA-dependent RNA polymerase that generates primers for coronavirus replication (7), and a mechanism to decrease replication errors (8).
Coronavirus subgenomic RNAs are particularly interesting in that they all have the same 5' leader sequence derived from the 5' end of the genomic RNA. This organization requires recombination as part of transcription. Various mechanisms have been proposed for subgenomic RNA production, but a discontinuous transcription mechanism is increasingly favored (5). This model proposes that transcription regulatory sequences in the minus-strand RNA direct translocation of the ternary complex to the 5' leader sequence, where transcription resumes. The minus-strand RNAs serves as the template for subgenomic RNA transcription. Ribonucleases that process the RNA intermediates for transcription have been proposed (9-11), but the mechanism for the process is still not completely understood.
Nsp15 (nonstructural protein 15) was predicted to be an RNA endoribonuclease as part of a bioinformatics analysis of the SARS-CoV genome (7). Substitutions of critical residues built into infectious clones of the human Coronavirus resulted in decreased infectivity and reduced genome RNA levels (12). A more thorough analysis of mutations in the Nsp15 ortholog of the related Arterivirus revealed that several residues, including those in the putative active site, reduced viral plaque size and decreased viral titers by up to five logs (13). Similar results were observed with Nsp15 of mouse hepatitis virus A59 (14). In both arterivirus and mouse hepatitis virus, reduced subgenomic RNA levels were associated with mutations in Nsp15. Other properties of virus infection were also affected. Therefore, it is likely that Nsp15 and its orthologs play multiple roles in viral infection (14).
To better understand how the activity of Nsp15 affects processing of viral RNAs, Bhardwaj et al. (15) and Ivanov et al. (12) independently produced recombinant Nsp15 of the SARS-CoV in Escherichia coli and demonstrated endoribonuclease activity that preferentially cleaves RNAs at uridylates. This activity was stimulated by Mn2+ but not by other divalent metals such as Mg2+ (15). More recently, Bhardwaj et al. (16) demonstrated that Nsp15 cleaves 3' of uridylates through the formation of a 2'-3' cyclic phosphodiester product.
Ribonucleases are generally grouped on the basis of whether they cleave 5' or 3' of the cognate phosphodiester (17). The former includes RNase H and RNase III, whereas the latter includes members of the RNase A family and RNase T1 (18). The two types of enzymes also differ in their metal ion dependence. The RNase H and III enzymes are metal-dependent, whereas RNase A and T1 are metal-independent (19). Thus, Nsp15 does not obviously belong to either group. Its cleavage is similar to RNase A, but its activity is enhanced by metal ion. This combination of properties is shared by XendoU, an endoribonuclease involved in small nucleolar RNA biosynthesis (19).
An unusual feature of Nsp15 is that it assembles into a hexamer in solution and that disruption of hexamer formation reduces endoribonuclease activity (20). Transmission electron microscopy analysis showed that Nsp15 consists of a dimer of trimers that interact in an end-to-end fashion (16). A low resolution co-crystal of Nsp15 and RNA suggested that the RNA binds to the outside of the hexamer and interacts with both trimers (16).
Two Nsp15 structures have been determined. Ricagno et al. (21) determined the x-ray structure of the SARS-CoV Nsp15 protein from crystals with one protomer in an asymmetric subunit. The crystal structure of mouse hepatitis virus Nsp15 was solved by Xu et al. (22), also with one subunit in the asymmetric unit. When arranged into a hexamer, the six active sites are exposed at the surface of the hexameric structure. Furthermore, the active site of Nsp15 can be modeled according to the residues of RNase A that share a common mechanism of catalysis (21, 22). In this structure, the C-terminal tail of Nsp15 folds back on the active site, suggesting that it may play a role in specific recognition of the cognate uridylate. These crystal structures, however, lack RNA, and the details of RNA interaction and many requirements of the subunit interactions remain unclear.
Using the coordinates of Ricagno et al. (21), we solved the structure of an active site mutant of the SARS-CoV Nsp15 by molecular replacement. Notably, this structure has six subunits within the asymmetric unit. Predictions of the structure were examined, leading to elucidation of the requirements for specificity in RNA cleavage, subunit interaction, and RNA binding.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
|
Molecular Docking—Dock 6.0 and its accessory programs (31) were used to carry out the flexible docking of 3'-UMP or CMP (from Protein Data Bank code 4RSK [PDB] -RNASE A:3'UMP complex structure) into the Nsp15 crystal structure (Protein Data Bank code 2H85), which was kept rigid. DMS (32), a software that computes molecular surface, was used to create the molecular surface of the receptor. The negative image of the binding site was defined using Sphgen (31) within the 10 Å radius of the catalytic site residues His234, His249, Ser293, Lys289, and Tyr342 to adopt the sphere-matching algorithm of Dock 6.0. An incremental construction (Anchor-and-Grow method) was used to allow for flexibility of the ligand. Automatic matching mode was used with the 20 configurations/ligand building cycle. Interaction between the ligand and the catalytic site was evaluated by the grid score (33, 34) followed by visual inspection.
Constructing the Nsp15-S293A Model—The atomic coordinates of Nsp15 crystal structure (Protein Data Bank code 2H85) were used to construct a model of Nsp15-S293A using Modeler (Insight II, www.accelrys.com). Residues within a 10 Å radius of A293 were subjected to a medium level simulated annealing optimization with respect to the variable target function of Modeler.
Site-directed Mutagenesis and Protein Purification—Amino acid substitutions were introduced in the Nsp15 expression plasmid using the QuikChange site-directed mutagenesis protocol as recommended by Stratagene. The entire open reading frame was sequenced to confirm the presence of directed mutation and absence of unintended mutations. WT and mutant proteins with a His tag at their N termini were expressed in E. coli and purified by metal ion affinity chromatography, followed by Mono Q ion exchange chromatography as described previously (15). The purified proteins were stored in 50 mM Tris (pH 7.9), 300 mM NaCl, 1 mM dithiothreitol, 50% (v/v) glycerol at -20 °C. The protein concentrations were determined by absorbance at 280 nm.
Size Exclusion Chromatography—WT or mutant protein samples (
100 µg) were filtered through a Superdex 200 column using a Pharmacia fast protein liquid chromatography system calibrated with the following molecular mass standards (kDa in parentheses): blue dextran (2000); thyroglobulin (660); ferritin (440); catalase (232); aldolase (140); bovine serum albumin (67); and chymotrypsinogen (25) according to the protocol described in Guarino et al. (20).
Reversible Protein-RNA Cross-linking and Peptide Mapping—H249A (10 µM) was incubated with or without 20 µM biotinylated RNA, Bio-rUC10 in 20 mM Hepes (pH 7.5), 5 mM MnCl2, and 1 mM dithiothreitol at room temperature for 5 min followed by the addition of formaldehyde to a final concentration of 0.1% in a 25-µl reaction. The cross-linking reaction was stopped by the addition of glycine at 0.2 M final concentration. After 5 min, sequencing grade trypsin (Trypsin Gold; Promega, Madison, WI) was added at a protease:substrate ratio 1:50 (w/w) in 100 mM NH4HCO3 (pH 7.8) at 37 °C and allowed to digest overnight. Streptavidin magnetic beads (New England Biolab, Beverly, MA) were used to capture biotinylated RNA and RNA-peptide conjugates as described previously by Kim et al. (4). RNA-peptide conjugates were reversed by incubating the samples at 70 °C for 1 h. The samples were centrifuged at 3000 x g for 5 min, and the supernatants containing the peptides were desalted using a Ziptip (Millipore, Bedford, MA). The bound peptides were eluted in 2.5 µl of 70% acetonitrile and 0.1% trifluoroacetic acid. The eluted samples were analyzed by MALDI-ToF.
Affinity Chromatography—10 µg of WT or mutant proteins were incubated on ice with
2 mg polyU-agarose in buffer A (50 mM Tris, pH 7.5, 5% glycerol, 1 mM β-mercaptoethanol, 10 mM MnCl2) containing 30 mM NaCl. After 90 min, the agarose beads were washed by centrifugation at 6,000 rpm for 10 min to remove unbound protein. Bound protein was eluted with buffer A containing increasing concentrations of 50, 100, 200, 300, 400, and 500 mM NaCl. The fractions were analyzed by SDS-polyacrylamide gel electrophoresis followed by staining with Coomassie Brilliant Blue. Protein was quantitated using the Coomassie Plus Protein Assay kit from Pierce.
Endoribonuclease Assays—The real time endoribonuclease assay used a substrate purchased from Integrated DNA Technologies, Inc. (Coralville, Iowa). The substrate has a 5' carboxyl fluorescein fluorophore at the 5' end and tetramethylrhodamine at the 3' end, which quenches 5' carboxyl fluorescein fluorescence. Of the four nucleotides, only the cognate cleavage nucleotide is a ribonucleotide, whereas the others are deoxyribonucleotides. Protein and RNA substrate were incubated in 20 mM Hepes (pH 7.5), 5 mM MnCl2, 1 mM dithiothreitol, as indicated. Substrate cleavage was monitored with an excitation wavelength of 492 nm and an emission wavelength of 518 nm (16). Changes in fluorescence were measured over time using an LS55 spectrometer (PerkinElmer Life Sciences).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
A single crystal of H234A protein, hereafter named Nsp15H, diffracted up to 2.64 Å with unit cell parameters a = b = 305.757 and c = 88.740 in the R3 space group. The asymmetric unit contained six molecules, and the positions of each subunit of Nsp15H were solved by molecular replacement using the structure of the catalytically active Nsp15 V291G-D300G double mutant (Protein Data Bank code 2H85) (21) as the search model. The Ricagno et al. (21) model, hereafter named Nsp15VD, has only one molecule/asymmetric unit. The asymmetric unit of Nsp15H provides additional insight into the heterogeneity among individual subunits of the hexamer and allows comparison with the WT active site in Nsp15VD.
The overall structure of Nsp15H is very similar to that of Nsp15VD, with a root mean square deviation difference of 0.53-0.73 Å (344C
atoms of each subunit of Nsp15 H were super-imposed on Nsp15VD structure using the LSQKAB routine of CCP4 (26)). Both Nsp15H and Nsp15VD consist of a dimer of trimers, which is consistent with the solution state and with single Nsp15 particles analyzed by transmission electron microscopy and cryoelectron microscopy (16).
Each subunit contains nine
-helices and twenty-one β-strands. A subunit is further organized into three domains: an N-terminal domain (residues 1-61), a middle domain (residues 62-181), and a C-terminal domain (residues 182-345) (Fig. 1A). Alanine substitutions of highly conserved residues have demonstrated that the C-terminal domain contains the active site (20), which faces away from the center of the hexamer and contains the extreme C-terminal residues, consistent with the observations of Ricagno et al. (21). When viewed from the top of the hexamer, a pore through the trimer is evident, with the N-terminal domains of the trimer lining the bottom of the pore (Fig. 1B). This pore has an inner diameter of
12 Å and does not interact with the substrate RNA according to the cryoelectron microscopy structure (16). The middle and the N-terminal domains have extensive contacts with the other subunits of the hexamer. These interactions are evident in the side view of the Nsp15 hexamer (Fig. 1C).
|
Heterogeneity is most obvious in the subunit catalytic pockets formed by residues His234, His249, Trp332, Tyr342, and Lys289. Superposition of the active site residues showed the degree of conformational flexibility between the subunits (Fig. 1E). Although the positions of Lys289 and His234 are highly similar in each of the subunits, Tyr342 is especially flexible and adopted different rotomer conformations in the different subunits (Fig. 1E). In addition, Trp332 is disordered in subunit T2c and not visible in the electron density map. The distance between the side chain oxygen atom of Tyr342 to
1-nitrogen atom of Trp332 varied between 3.0 and 4.1 Å in Nsp15H.
Subunit Interactions—Based on the examination of a panel of mutants, Guarino et al. (16, 20) proposed that hexamer formation is a prerequisite for enzymatic activity and RNA binding. A key mutant in this analysis was E3A, located in the N-terminal domain of Nsp15. It was inactive at lower concentrations but recovered some catalytic function at higher concentrations that should favor oligomerization (20). Structural data suggests that E3 and the adjacent L2 residue of one subunit are involved in trimer-trimer interactions with L2 and E3 of a subunit in a different trimer (Fig. 2A, inset). This interlocking structure thus provides a reasonable explanation for the effects of the E3 mutation on Nsp15 oligomerization (20). Furthermore, because oligomerization is required for activity, it is important to examine oligomerization in greater detail.
|
In interaction surface B, the backbone oxygen atom of residue Asp39 of subunit T1a forms an ionic bond with the H1-nitrogen atom of Arg90 in T1c (Fig. 2B). In addition, the
2-oxygen atom of Asp39 forms a H-bond with the
1-oxygen atom of Thr48 in T1c. To test the functional relevance of these interactions, we analyzed hexamer formation and enzymatic activity of the mutant versions D39A and R90A. Oligomerization state was analyzed using a rapid and reproducible gel filtration assay. The validity of this assay has been previously confirmed by analytical centrifugation assays and also a native gel analysis (20). Both D39A and R90A eluted in a gel filtration column as monomers (Fig. 2B), confirming the importance of these residues in hexamer formation.
A real time fluorescent assay previously described in Bhardwaj et al. (16) was used to monitor endoribonuclease activity. The substrate contains three deoxyribonucleotides, one cleavable ribonucleotide (rU), a 5' fluorophore, and a 3' quencher. Cleavage of the substrate results in an increase in fluorescence, which can be measured in real time. Functional analysis of D39A and R90A, using the real time endoribonuclease assay, revealed that the activity of these mutants was less than 1% of the WT protein.
In surface C, the backbone oxygen atom of N163 in T1b interacts with the
-nitrogen atom of K264 in T1a. In addition, the
2-nitrogen atom of N163 in T1b bonds the
2-oxygen atom of T281 in T1a (Fig. 2C). Alanine substitution of Asn163 should disrupt the latter interaction and resulted in a nearly equal distribution between the monomer and the hexamer, suggesting that hexamer formation was affected but not abolished. In the real time endoribonuclease activity assay, N163A retained 29% of the activity of the WT protein.
In surface D, the H1-nitrogen atom of Arg61 in subunit T1a interacts with the
1-oxygen atom of Glu266 in T1c. Gel filtration analysis of the mutant protein E266A showed that it eluted primarily as a monomer. Furthermore its endoribonuclease activity was less than 0.4% of WT (Fig. 2D). With the three sets of mutants that affected intratrimer interactions, it is interesting to note that only mutant N163A retained some ability to form hexamers and also retained some endoribonuclease activity.
|
1-oxygen and
1-nitrogen atoms of Asn52 in T1a with backbone nitrogen atom of Ile26 and oxygen atom of Val24, respectively, in T2b. When mutant N52A or I26A,I27A were tested, both proteins existed primarily as monomers. Again, both mutants were severely debilitated for endoribonuclease activity, with cleavage rate of rU being less than 1% of WT (Fig. 2E). Interestingly, mutants E3A, N52A, and I26,27A, which were predicted to form trimers because these residues are involved in trimer-trimer interactions, all eluted as monomers and not trimers as we expected (Fig. 2, A, inset, and E). These results suggest that a trimer is not a stable intermediate in the formation of hexamers and that the equilibrium favors monomers and hexamers. Some trimer (<10%) was seen previously in gel filtration columns (20). It is possible that these reflect a dissociation event.
Functional Analysis of Active Site Residues—An interesting aspect of the Nsp15 active site is that the C-terminal tail contributes to its formation. To obtain an initial model for substrate recognition at the Nsp15 active site, we performed a molecular docking analysis with UMP. The results suggested that some of the active site residues might contribute to specific interactions with U and help discriminate against a C (Fig. 3A). First, Ser293 could potentially H-bond with the N3 imino hydrogen and the C2 carbonyl oxygen of the uracil. This information is consistent with the proposal of Ricagno et al. (21), who postulated that Ser293 was comparable in spatial location to Thr45, the residue in RNase A that hydrogen bonds to the base and specifies cleavage. Second, the phenyl side chain of Tyr342 could stack against the uracil, resulting in Pro343 and Leu345 occupying the space above the C4 carbonyl of uracil. In addition, three basic residues, Lys289, His234, and His249 that have been previously identified by bioinformatics analysis to participate in catalysis (7) could form an ionic surface that interacts with the ribose portion of the uridylate and the adjacent phosphate. Consistent with their proposed roles in catalysis, these three residues are positioned to act on the ribose 2' OH and the adjacent phosphodiester group (Fig. 3A). It is known that mutations in any of these three residues reduced RNA cleavage to background (20).
To examine whether and how residues in the catalytic site, including the C-terminal tail, affect endoribonuclease activity, we made mutant versions of Nsp15 that changed Ser293 and Tyr342 to alanines to disrupt H-bond or stacking interactions, which these residues might have with substrate. Pro343 and Leu345 were mutated to glycines to widen the catalytic pocket. Endoribonuclease analyses with the rU substrate revealed that the WT protein had a cleavage rate of
1.5 µmol-1. The rate for mutant S293A was
0.25 µmol-1, P343G and L345G cleaved the substrate at 0.02 and 0.055 µmol-1 respectively, whereas the activity of Tyr342 was undetectable (Fig. 3B). These results confirm that the residues identified in the catalytic pocket are involved in catalysis.
Based on the structure, we did not expect oligomerization to be affected by the changes in proteins S293A, Y342A, P343G, and L345G. To test this directly, however, we examined the oligomerization state of these proteins by gel filtration. WT Nsp15 and all four of these mutants eluted in a position corresponding to hexamers (
240 kDa) (Fig. 3C), confirming that the C-terminal tail does not affect oligomerization but is required for endoribonuclease activity.
|
Specificity for Uridylate—Molecular docking of UMP into the Nsp15 active site (Fig. 4A) suggests that N3 and O2 of the uridylate could form H-bonds with
-oxygen atom and backbone nitrogen atoms, respectively, of Ser293, whereas O4 could interact with either Leu345 or Pro343 via H-bonds (Fig. 4B). This possibility is supported by our earlier results showing that modification of either N3 or O4 of uridylate affected cleavage by Nsp15 (16). Cytosine should not form equally strong interactions with the active site residues. Although the O2 of cytosine could interact with backbone nitrogen atom of Ser293, the
-oxygen atom of Ser293 would need to reorient to form a second H-bond with N3 of cytosine. Additionally, the bulkier NH2 group at C4 could sterically hinder the fit of cytosine in the pocket (Fig. 4B, inset).
Molecular modeling and energy minimization of mutant S293A with CMP or UMP indicate that in absence of Ser293, Pro343 could interact with either O4 of uridylate or the C4 amino group of cytosine to stabilize binding (supplemental Fig. S1A). In addition, the side chains of Lys289, Leu248, and Leu345 would adjust to better complement the electronics of a cytosine in the active site (supplemental Fig. S1B). The S293A mutation should also result in a decreased anionic nature of the active site cleft and a narrowing of the active site cavity, which could accommodate a uridylate or a cytidylate equally well, although with a lower affinity as compared with the affinity of WT protein to a uridylate. Taken together, our observations from the models is that uridylate should be a preferred substrate over purines or a cytidylate because it can form stronger specific interactions and has a better steric fit in the catalytic pocket.
To test these predictions, first we analyzed the effect of mutations that would potentially alter the interaction between uridylate and the catalytic residue Ser293. In addition to the already characterized S293A, we changed Ser293 to asparagine (S293N) or threonine (S293T). All three mutant proteins eluted in a gel filtration column as hexamers (data not shown). To examine the specificity for RNA cleavage with greater sensitivity, we used the fluorescent real time assay and with 0.25 µM of S293A or WT Nsp15, a concentration that is 5-fold higher than our typical assay. We also compared the normal rU substrate and a similar substrate having cytidylate as the only ribonucleotide (named rC), as well as substrates with guanosine (rG) and adenosine (rA). WT reproducibly cleaved rU with a rate of 4.8 µmol-1 and rC with a rate of 0.285 µmol-1 (Fig. 4C). The ratio of rU to rC cleavage, or the specificity factor, is 16.8 for WT. There was no cleavage by the WT protein with rG or rA (Fig. 4D and data not shown).
Mutant S293N was severely decreased in activity, with rates of cleavage for rU and rC at less than 0.017 µmol-1. S293T cleaved rU and rC with rates of 4.6 and 0.28 µmol-1, respectively, thus generating a specificity factor of 16.4, which is comparable with that of WT (supplemental Fig. S2). The fact that S293T is not affected in endoribonuclease activity is not surprising given that the comparable residue in the mouse hepatitis virus Nsp15 is a threonine, and the infectivity of EAV was not significantly affected either residue at the comparable position in its endoribonuclease ortholog (13).
S293A had a more interesting activity profile. It had reduced cleavage rate, with rU being cleaved at 2.35 µmol-1 (Fig. 4C). Unlike WT enzyme, however, S293A cleaved rU and rC with comparable rates (2.55 and 2.35 µmol-1), resulting in a specificity factor of 1.1. Furthermore, S293A has comparable Km values for rU (26.7 + 8.5 µM) and rC (28 + 2.6 µM), values similar to that of WT Nsp15 for rU (30.3 ± 9.5 µM). Like WT, S293A did not cleave rA or rG, thus discrimination against purine substrates was not altered. These results support our molecular docking predictions and the idea proposed by Ricagno et al., (21). These modeling results suggest that Ser293 participates in specific recognition of the substrate by Nsp15 and further extends the similarities in the cleavage mechanism for Nsp15 and RNase A.
Next, we analyzed the mutants P343G and L345G. We predicted that these substitutions could potentially widen the pocket and/or increase flexibility in the active site. P343G was reduced in catalytic activity to 1.3% of WT and had a specificity factor of 10.4 ± 2 for rU/rC. The more interesting result was with L345G. Leucine at this position should sterically block the binding of a purine in the WT protein (Fig. 4A). Substitution of glycine at this position decreased overall cleavage activity and resulted in a lower rU/rC specificity factor of 2.3 ± 0.2 relative to WT. Interestingly, L345G cleaved purine substrates at a significantly higher rate than WT or other catalytic mutants (Fig. 4D). Because leucine 345 is the penultimate residue, these results suggest that the C terminus of Nsp15 contributes to the specific recognition of uridylates in part by preventing the access of other nucleosides into the active site.
Nsp15-RNA Interactions—Images of Nsp15 in 2D crystals in complex with RNA indicates that an RNA of 16 nucleotides binds across the interface between two trimers (16). Furthermore, only hexamers bind RNA in solution, indicating that RNA binding follows oligomerization (16). The catalytic pockets are surface-exposed in both Nsp15VD and Nsp15H crystal structures. These pockets, however, are too narrow in both structures to accommodate an RNA molecule, and a conformational change is necessary for RNA binding. The flexibility observed in this region may be a part of the mechanism to allow RNA docking into the catalytic pocket.
Three different RNA access paths for the catalytic site of Nsp15 hexamer were identified based on surface complementarity and were finalized by manual docking of a single-stranded RNA of sixteen uridylates (built using the biopolymer module of Insight II (Accelrys Inc.)). These dockings were done with some steric clashes that needed to be resolved by slight conformational changes in the protein. The three potential paths for RNA binding are shown in Fig. 5A. To examine whether any or all of these potential paths contact RNA, we used a method that does not contain built-in bias in selecting RNA-contacting residues, a reversible cross-linking peptide fingerprinting method that was previously used to identify the template-binding portions of the hepatitis C virus RNA-dependent RNA polymerase (4). Briefly, a biotinylated RNA is cross-linked to a protein by use of formaldehyde. The reaction is then treated with trypsin to digest the protein after lysines and arginines, and the RNA and RNA-peptide complexes are enriched by streptavidin affinity chromatography. After reversing the cross-links, the peptides are purified from the RNA and subjected to MALDI-ToF spectroscopy. The peptides are identified from peptides generated by a virtual cleavage.
The biotinylated RNA used in the reversible cross-linking peptide fingerprinting contained multiple uridylates to facilitate binding by Nsp15 (16) and several cytidylates that contain amines that can be cross-linked with formaldehyde. The catalytically inactive mutant H249A was used in these experiments to maximize the production of RNA-protein complexes by avoiding catalytic cleavage of the substrate. A typical mass spectrum for the H249A-RNA cross-linked sample is shown in the lower panel of supplemental Fig. S3, whereas the result from a control reaction performed in the absence of RNA is in the upper panel. After subtracting the peaks corresponding to streptavidin and/or trypsin, eight ions that could be attributed to H249A were identified (Fig. 5B). Peptides containing over-lapping sequences to these eight were also observed when the samples were treated with chymotrypsin (data not shown), providing a confirmation for their assignment.
When the peptides were mapped back to the Nsp15 structure, the pattern shown in Fig. 5C was observed. Notably, the affinity-purified peptides border the predicted paths identified by bioinformatics analysis. In addition, peptides corresponding to the active site were found. The location of the peptides, however, did not allow us to identify a single path, suggesting either that all three are used at some level or that this analysis is insufficient to discriminate among the possibilities.
To obtain a better resolution of the RNA paths, we used the locations of the RNA-interacting peptides in the crystal structure and sequence conservation to guide the design of additional amino acid substitutions. This analysis required assumptions of the more likely residues based on the probability of particular side chain atoms to participate in formaldehyde cross-linking as well as orientation of the side chains within the Nsp15 hexamer structure. Formaldehyde preferentially cross-links residues that have basic, aromatic, or sulfur-containing side chains (29). Based on the orientation of the side chains, alanine substitutions of residues His18, Asn73, Asn82, Asn83, Arg181, Phe213, and Lys256 were constructed. In addition, based on their conservation and presence in the affinity purified tryptic peptides, alanine substitutions at Asn28, Asn29, and Thr47 were also included in the mutational analysis. S293A, P343G, Y342A, and L345G were previously made in our analysis of the active site, and they are known to have effects on RNA binding, respectively (Fig. 3D). Because of their proximity, Asn28 and Asn29 were mutated together, as were Asn82 and Asn83. Because hexamer formation is a prerequisite for RNA binding, all of the mutants were examined for RNA binding and were found to be hexameric (data not shown).
Next, we used polyU affinity resin to compare the stability of binding of mutant Nsp15 proteins with that of WT. Mutants K181A, F213A, and N82A,N83A bound to polyU resin at levels comparable with WT, with the peak of the proteins eluting at 400 mM NaCl (Fig. 5D). Mutants H18A, T47A, N28/29A, N73A, and K256A were significantly decreased for RNA binding, with the peak of the proteins eluting at 200 mM NaCl or lower, suggesting that these residues directly or indirectly contribute to RNA binding. When located onto the Nsp15 structure, the effects of the mutants indicate that Path III is most likely involved in stable binding to RNA (Fig. 5, A and E).
|
| FOOTNOTES |
|---|
* This work was supported by National Institutes of Health Grant 1R56AI06164 and a Small Grant for Exploratory Research from the National Science Foundation. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
The on-line version of this article (available at http://www.jbc.org) contains supplemental Table S1 and supplemental Figs. S1-S3. ![]()
1 Both authors contributed equally to this work. ![]()
2 To whom correspondence should be addressed. Tel.: 979-458-2235; Fax: 979-845-9274; E-mail: ckao{at}tamu.edu.
3 The abbreviations used are: SARS, severe acute respiratory syndrome; CoV, coronavirus; WT, wild type; MALDI-ToF, matrix-assisted laser desorption ionization time-of-flight. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. M. Gustin, B.-J. Guan, A. Dziduszko, and D. A. Brian Bovine Coronavirus Nonstructural Protein 1 (p28) Is an RNA Binding Protein That Binds Terminal Genomic cis-Replication Elements J. Virol., June 15, 2009; 83(12): 6087 - 6097. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. D. Nedialkova, R. Ulferts, E. van den Born, C. Lauber, A. E. Gorbalenya, J. Ziebuhr, and E. J. Snijder Biochemical Characterization of Arterivirus Nonstructural Protein 11 Reveals the Nidovirus-Wide Conservation of a Replicative Endoribonuclease J. Virol., June 1, 2009; 83(11): 5671 - 5682. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Yi, E. Letteney, C.-H. Kim, and C. C. Kao Brome mosaic virus capsid protein regulates accumulation of viral replication proteins by binding to the replicase assembly RNA element RNA, April 1, 2009; 15(4): 615 - 626. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. T. Ranjith-Kumar, A. Murali, W. Dong, D. Srisathiyanarayanan, R. Vaughan, J. Ortiz-Alacantara, K. Bhardwaj, X. Li, P. Li, and C. C. Kao Agonist and Antagonist Recognition by RIG-I, a Cytoplasmic Innate Immunity Receptor J. Biol. Chem., January 9, 2009; 284(2): 1155 - 1165. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. T. Ranjith-Kumar, K. E. Duffy, J. L. Jordan, A. Eaton-Bassiri, R. Vaughan, S. A. Hoose, R. J. Lamb, R. T. Sarisky, and C. C. Kao Single-Stranded Oligonucleotides Can Inhibit Cytokine Production Induced by Human Toll-Like Receptor 3 Mol. Cell. Biol., July 15, 2008; 28(14): 4507 - 4519. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |