Cleavage of Model Substrates by Arabidopsis thaliana PRORP1 Reveals New Insights into Its Substrate Requirements

Two broad classes of RNase P trim the 5' leader of precursor tRNAs (pre-tRNAs): ribonucleoprotein (RNP)- and proteinaceous (PRORP)-variants. These two RNase P types, which use different scaffolds for catalysis, reflect independent evolutionary paths. While the catalytic RNA-based RNP form is present in all three domains of life, the PRORP family is restricted to eukaryotes. To obtain insights on substrate recognition by PRORPs, we examined the 5' processing ability of recombinant Arabidopsis thaliana PRORP1 (AtPRORP1) using a panel of pre-tRNASer variants and model hairpin-loop derivatives (pATSer type) that consist of the acceptor-T-stem stack and the T-/D-loop. Our data indicate the importance of the identity of N-1 (the residue immediately 5' to the cleavage site) and the N-1:N+73 base pair for cleavage rate and site selection of pre-tRNASer and pATSer. The nucleobase preferences that we observed mirror the frequency of occurrence in the complete suite of organellar pre-tRNAs in eight algae/plants that we analyzed. The importance of the T-/D-loop in pre-tRNASer for tight binding to AtPRORP1 is indicated by the 200-fold weaker binding of pATSer compared to pre-tRNASer, while the essentiality of the T-loop for cleavage is reflected by the near-complete loss of activity when a GAAA-tetraloop replaced the T-loop in pATSer. Substituting the 2'-OH at N-1 with 2'-H also resulted in no detectable cleavage, hinting at the possible role of this 2'-OH in coordinating Mg2+ ions critical for catalysis. Collectively, our results indicate similarities but also key differences in substrate recognition by the bacterial RNase P RNP and AtPRORP1: while both forms exploit the acceptor-T-stem stack and the elbow region in the pre-tRNA, the RNP form appears to require more recognition determinants for cleavage-site selection.


Introduction
Most tRNA genes are transcribed as precursor RNAs (pre-tRNAs) with both the 5' and 3' ends having additional residues that need to be removed to generate functional, mature tRNAs.The ubiquitous ribonucleoprotein (RNP) ribonuclease P (RNase P) is responsible for removing the 5' leader from pre-tRNAs.In Bacteria, RNase P is composed of one RNA subunit and one protein subunit, while in Archaea and Eukarya four or more proteins associate with the sole RNA [1,2].Irrespective of origin, the catalytic activity resides in the RNase P RNA (RPR) as evident from its ability, even in the absence of associated protein cofactor(s), to mediate cleavage of pre-tRNA as well as various other natural (e.g., pre-4.5SRNA) and artificial (e.g., model hairpin loop) substrates [1][2][3][4][5][6][7].
In several eukaryotes, there also exists an RNA-free RNase P that is composed solely of proteins [8].PRORP (proteinaceous RNase P) cleaves pre-tRNAs at the same site as the RNP variants, and is also involved in tRNA 5'-maturation.In Arabidopsis thaliana, three distinct PRORPs (AtPRORP1, 2 and 3) are present, but an RPR has not been identified [9].AtPRORP1 is localized to the mitochondria and chloroplasts, while AtPRORP2 and AtPRORP3 are targeted to the nucleus [9].Single-polypeptide PRORPs from A. thaliana nucleus/organelles have been characterized and shown to be active as individual entities, while the human mitochondrial native variant was purified as a complex with two other proteins [8,9].RNAi-mediated knock-down of AtPRORP1 showed protein synthesis defects in chloroplasts and mitochondria, although only photosynthesis was defective and respiration was unaffected; interestingly, the effects on 5' processing of individual organellar tRNAs were not uniform [10].To better understand these phenotypic effects and, more broadly, appreciate the choice of RNP-and proteinbased RNase P for pre-tRNA/RNA processing, it is important to understand how the two variants recognize and process their substrates [10,11], the motivation for this study.
By examining cleavage of pre-tRNAs and model substrates, residues at and near the cleavage site have been demonstrated to influence both cleavage-site recognition and cleavage efficiency of bacterial ribonuclease P (for a review, see [12]).Specifically, the residue N -1 , the discriminator base and the two C residues at the pre-tRNA 3' end, and the T-loop have key roles [7,[13][14][15]; for reviews, see [1,16].In contrast, we have little information about either the impact of individual substrate residues and chemical groups on cleavage or if members of the PRORP family process small model substrates.
Given the ability to chemically synthesize short RNAs (~50 nts), especially with desired chemical modifications, we previously invested considerable effort into design and validation of short hairpin model substrates for the RNP version of RNase P. We have now used this approach to investigate for the first time the effect of certain site-specific replacements (e.g., guanosine with inosine or a 2'-OH with a 2'-H) on substrate recognition and cleavage by AtPRORP1.Our data show that recombinant AtPRORP1 cleaves model hairpin loop substrates with at least a 1000-fold lower single-turnover rate than that observed for cleavage of the parental pre-tRNA (pSu1, the Escherichia coli tRNA Ser Su1 precursor).We also found a dramatic decrease in the cleavage rate upon replacement of either the 2'-hydroxyl at -1 or the seven-bp T-loop equivalent with a GAAA-tetraloop in the model substrate.Moreover, like the bacterial RPR, the -1 identity is an important cleavage-site determinant in the context of both pre-tRNA and model substrates, irrespective of whether the -1 residue is paired or not with the residue at the discriminator position.These results led to some predictions in terms of disfavored sequences for processing by AtPRORP1.We gained support for these predictions by examining the sequences of all mitochondrial and chloroplast tRNA genes from eight different green algae and plants, an analysis not reported before.Together, these findings provide new insights into AtPRORP1-mediated catalysis and offer possibilities to dissect the role of individual residues and chemical groups important for cleavage.We have also integrated our findings with two very recent studies on PRORP-mediated substrate recognition [17,18] that appeared during preparation of this manuscript.

Preparation of substrates for binding studies
The DNA template for in vitro transcription of pSu1 with a 5-nt trailer was generated by PCR using primers FWD (5'-taatacgactcactatagatctgaatggagag-3'; the italicized g was added to facilitate transcription) and REV (5'-ggtgtcggagagagggggattt-3'; the trailer sequence added is italicized).The DNA template was the plasmid pUC19-pSu1 [20].
3'-Labeling of pSu1 and the pATSer derivatives was performed with some modifications of a previously described procedure [26][27][28].For each substrate, 130 μM of in vitro transcribed RNA in 100 μL 100 mM NaOAc (pH 4.5) was oxidized by addition of 10 mM NaIO 4 , and incubated at 22°C for 1.5 h in the dark.The RNAs were then ethanol precipitated and re-suspended in 500 μL of 100 mM NaOAc (pH 5.2) using a 20:1 molar ratio of fluorescein-5-thiosemicarbazide (FTSC):RNA; FTSC was a generous gift of Prof. Edward Behrman, Ohio State University (OSU).The labeling reactions were carried out at 4°C for 16 h in the dark.Excess, unincorporated FTSC was removed by sequential phenol-chloroform and charcoal extractions, followed by purification using a 8% (w/v) polyacrylamide/7 M urea gel.The excised RNA was eluted at 4°C for 16 h into 1 M NaOAc (pH 4.9), and then subjected to ethanol precipitation.The 3'-labeling efficiency was typically >90%, as assessed by Abs 260 (RNA) and Abs 492 (fluorescein) values.

Cleavage assays and determination of k app
The cleavage reactions with AtPRORP1 (purified as described in ref. 28) were performed in buffer containing 20 mM HEPES-KOH (pH 7.4), 100 mM NH 4 OAc, 4 mM DTT, 10 mM Mg (OAc) 2 and 0.8 mM spermidine.To determine the optimal Mg 2+ concentration for cleavage, Mg(OAc) 2 was added separately to give the final concentration as indicated.All assays were performed at 37°C.The reactions were terminated by adding twice the assay volume of stop solution (10 M urea, 100 mM EDTA), and the products were separated on 25% (w/v) polyacrylamide/7 M urea gels.
The rate constant k app was determined under single-turnover conditions at pH 7.4 in the presence of 10 mM Mg 2+ , which was determined to be optimal for AtPRORP1-mediated cleavage of pSu1 and pATSerUG.The concentration of AtPRORP1 used was 0.37 μM for assays with pSu1 [except 1.1 μM for pSu1(-1C)] and 5.6 μM for assays with pATSer derivatives (except 4 μM for pATSer 3' truncated variants).The concentrations of AtPRORP1 used to generate the data are specified in the respective figure legends.The concentration of pSu1 and model substrates was 0.02 μM.For rate calculations, we used the 5' cleavage fragment as a measure of product formed.In each assay, the time of incubation was adjusted to ensure that the Secondary structures of substrates used in this study.Secondary structures of pSu1 and pATSer.The highlighted regions/residues were substituted to generate the different variants as indicated, A, adenosine, G, guanosine, U, uridine, I (Ino), inosine and D (DAP), 2,6-diaminopurine; dC, deoxycytosine; and dU, deoxyuridine.The canonical RNase P cleavage sites between residues N -1 and N +1 (correct cleavage denoted C 0 ), and the alternative cleavage sites between residues N -2 and N -1 (miscleavage denoted as M -1 ) are marked with black and grey arrows, respectively.The N +73 position, which immediately precedes the 3'-terminal CCA-motif, corresponds to the discriminator base.doi:10.1371/journal.pone.0160246.g001velocity measurements were in the linear range (typically 10% but never exceeding 40%).Each k app value is reported as a mean ± standard deviation of this value, which were calculated using data (six time points) from at least three independent experiments.

Fluorescence polarization binding assays and determination of K D values
Defined amounts of AtPRORP1, as indicated, were incubated individually with either 2 nM pSu1 or 20 nM pATSer derivatives that had been 3'-labeled with fluorescein [28].The binding reactions were performed in 20 mM HEPES (pH 7.2), 10 mM Ca(OAc) 2 , 100 mM NH 4 OAc, 4 mM DTT, and 5% (v/v) glycerol.The reactions were carried out for at least 10 min at 22-25°C in a 384-well plate (Corning Costar black round bottom).The fluorescence polarization values were then obtained using infinite M1000 PRO (Tecan), with the G factor set to 1.2.Polarization (P) observed in the presence of different AtPRORP1 concentrations were subtracted from that observed with the respective substrate alone to obtain ΔP at each protein concentration tested.The dissociation constants were then calculated by fitting to DP ¼ DP max x ½AtPRORP1 K D þ ½AtPRORP1 using Kalei-daGraph (Synergy).The curve-fit errors for each measurement did not exceed 26%, with R 2 values 0.96.Each K D value is reported as a mean ± standard deviation, which were calculated using data from at least three independent experiments.

Results
The identity of N -1 in pre-tRNA Ser (pSu1) influences cleavage by AtPRORP1 Studying the recognition and cleavage of a suite of model substrates (pATSer series) derived from Eco pre-tRNA Ser Su1 (pSu1) by the bacterial RNase P RNP has been gainful [6,7,14,15,[20][21][22][29][30][31][32][33][34][35].To facilitate a direct comparison of substrate recognition by the RNP and proteinaceous forms of RNase P, we therefore chose to exploit the same pAT series of model substrates.Moreover, compared to other pre-tRNAs used to study PRORP-mediated cleavage [17,18,36], Eco tRNA Ser is equipped with a longer variable loop thus enabling a comparison of structurally distinct pre-tRNAs.Towards this overall objective, we first investigated if a recombinant AtPRORP1 could cleave pSu1 [Fig 1 ; wild type pSu1 referred hereafter as pSu1(-1C)].
Eco RPR cleaves pSu1(-1C) predominantly at the canonical correct position between N -1 and N +1 (termed C 0 ), but also miscleaves between N -2 and N -1 (termed M -1 ; Fig 1) [13].In contrast, AtPRORP1 cleaved pSu1(-1C) mainly at M -1 but also at Together, these findings suggest that the identity of N -1 and/or pairing between N -1 and the discriminator base (as in C -1 :G 73 ) play an important role in cleavage-site selection.
Because an examination of the single-and multiple-turnover rates indicated that cleavage (or a preceding step) is likely to be rate limiting for AtPRORP1 [37], we determined the apparent rates (k app ) of cleavage for the pSu1 "-1 variants" under single-turnover conditions.We first determined that the optimal Mg 2+ concentration for cleavage of pSu1(-1U) by AtPRORP1 was 10 mM Mg 2+ (Fig 4A ); we found that the choice of cleavage site did not change with increasing Mg 2+ .Hence, we chose 10 mM Mg 2+ for the kinetic studies.
When we examined the different model substrates for cleavage at C 0 and M -1 , k app showed a three-fold variation with pSu1(-1C) being the weakest substrate.In contrast, k app for cleavage at M -1 (the incorrect site) was roughly 20-fold higher for pSu1(-1C) compared to the other three N -1 variants (Table 1) consistent with its miscleavage propensity.Irrespective of the substrate tested, the frequency of cleavage at M -1 and C 0 did not change as a function of time (not shown).
The 3'-CCA in pre-tRNA Ser (pSu1) is not a major determinant for cleavage by AtPRORP1 Eco pSu1 has a 3' terminal CCA-motif (Fig 1).However, eukaryotic and organellar tRNA genes in general do not encode CCA (see e.g.http://trna.ie.niigata-u.ac.jp/cgi-bin/trnadb/index.cgi.).When we analyzed the organellar tRNA sequences for 8 algal and plant species (available at http://plantrna.ibmp.cnrs.fr.),only 0.5% (2 out of 423) tRNA-encoding genes have a 3'-CCA: a choloroplast tRNA Ala in Cyanophora paradoxa and a mitochondrial tRNA Ile in Solanum tuberosum (potato).Thus, AtPRORP1-localized to the mitochondria and chloroplasts-may not encounter pre-tRNAs with 3'-CCA.Although the inference was drawn from a single end-point measurement, it was previously reported that the presence of the 3'-CCA in pre-tRNA decreases AtPRORP1 cleavage and might therefore serve as an anti-determinant [9].We therefore generated truncated pSu1 "-1 variants", which lack this CCA-motif (  1).With respect to cleavage site-selection, we did not observe any major difference with and without the 3'-CCA, if anything a small increase in cleavage at M -1 for pSu1(-1C) and pSu1(-1U) in the absence 3'-CCA (Fig 3A).Upon deletion of the 3'-CCA motif, we noted a modest increase in k app (at C 0 ) for substrates having A -1 or G -1 , while a decrease was detected for those with C -1 or U -1 .The most striking effect was a 3.5-fold decrease in k app for cleavage of pSu1(-1U) at C 0 (Table 1).A simple classification that the 3'-CCA motif acts as a positive or negative determinant is not possible given the substrate-context effects.

Cleavage of model hairpin loop substrates by AtPRORP1
We next investigated whether AtPRORP1 cleaves the model hairpin loop substrate pATSerCG, which is composed of the 5' leader, the amino acid acceptor-stem (with the 3'CCA-motif and a dangling 3'C), and the T-stem and loop of pSu1(-1C) (  The rates of cleavage (k app ) of pATSerCG and pATSerUG at 10 mM Mg 2+ were dramatically lower than their pSu1 counterparts (Table 1).For pATSerCG, the C 0 and M -1 rates are 2500and almost 900-fold lower, respectively, while for pATSerUG cleavage at C 0 was three orders of magnitude lower.Deleting the 3'-CCA-motif resulted in a modest increase in k app for both pATSerCG and pATSerUG.In this context, note that deletion of both C's and the 3' terminal A is needed to elicit a modest increase in k app (Table 1).
Despite weak cleavage of the model substrates, compared to the parental pre-tRNA, the qualitative trends with respect to cleavage-site selection are similar for pSu1 and the pATSer N -1 variants (Fig 3).For example, comparison of pSu1(-1C) and pATSerCG (both without 3'-CCA) reveals that the k app for cleavage at M -1 relative to C 0 is two-fold greater in each case (Table 1).For the same cohort with 3'-CCA, the k app for cleavage at M -1 relative to C 0 is 1.4-fold higher for pSu1(-1C) and four-fold for pATSerCG (Table 1).Hairpin loop substrate binds with lower affinity than pre-tRNAs to AtPRORP1 We next used a previously described fluorescence polarization assay [38] to determine the dissociation constants (K D values) for the binding of 3'-CCA-containing pATSerUG and pSu1 (-1U) to AtPRORP1.These binding reactions were performed in the presence of Ca 2+ , because  AtPRORP1 shows tight pre-tRNA binding but no detectable cleavage when Mg 2+ is substituted with Ca 2+ [28,38].The K D value for pATSerUG increased by almost 200-fold relative to pSu1 (-1U) (Table 2; see also C -1 and N -1 :N +73 pairing influence cleavage by AtPRORP1.Our results show that N -1 identity influences cleavage by AtPRORP1, as is particularly evident from results obtained with C -1 substrates that were cleaved preferentially at the alternative site M -1 (miscleavage).Both pSu1 and pATSer have G +73 as the discriminator base, and therefore have the possibility of C -1 : G +73 pairing.Thus, the bp just upstream of the correct cleavage site could affect fidelity and rate.To investigate this possibility, we next generated pATSer variants with different N -1 :N +73 options (Fig 1).These substrates are referred to as pATSerAG, pATSerGG, pATSerCIno (inosine at +73 can potentially form two H-bonds with C -1 ), pATSerDAPG (2,6-diamino purine at -1), pATSerInoG (inosine at -1), pATSerUDAP (2,6-diamino purine at +73 can potentially form three H-bonds with U -1 ) and pATSerCU.
The cleavage data (Figs 3 and 5) showed that pATSer acts as an AtPRORP1 substrate irrespective of the identity of residue -1.While substrates having A -1 , G -1 , DAP -1 , Ino -1 and U -1 (except for pATSerUDAP) were cleaved preferentially at the correct site, C -1 resulted in cleavage at M -1 even when it is not engaged in pairing with N +73 as evident from cleavage of pAT-SerCU (Fig 3B).Moreover, formation of a N -1 :N +73 pair with three H-bonds resulted in cleavage mainly at the alternative site M -1 (see pATSerCG and pATSerUDAP).Together, these data suggest that C -1 as well as the presence of a N -1 :N +73 pair with three H-bonds in a pATSer context influence the choice of cleavage site by AtPRORP1.Consistent with these findings, the k app for pATSerCU (absence of the -1/+73 pair) was ten-fold higher for cleavage at the correct site compared to pATSerCG while it was two-fold higher for pATSerCIno (Table 1).
Since our findings indicated that N -1 identity and the strength of the N -1 :N +73 pair play important roles in determining the rate and fidelity of cleavage, and that some combinations result in adverse effects with respect to AtPRORP1 catalysis, we postulated that a bias might become apparent from an analysis of the N -1 :N +73 sequence information among the organellar tRNAs in eight different green algae and plants (Fig 6 ; Table A in S1 File).From examining these 423 tRNAs, we observed the following features.First, C -1 :G +73 was present in only 0.7% of the tRNAs (3 instances) even though C -1 displayed a 10-fold higher incidence (7.6%, 32 out of 423).Second, there were eight examples of G -1 :C +73 (~2%), but seven of these were tRNA His ; a universal identity determinant of tRNA His for histidyl-tRNA synthetase is the presence of G -1 and an 8-bp acceptor stem.Last, there is a variable distribution of other pairing possibilities: A -1 :U +73 (6.4%), U -1 :A +73 (18.2%) and G -1 :U 73 (2.4%).Although an in-depth analysis is needed to draw firm conclusions, it appears that two hydrogen bonds in the N -1 :N +73 pair alone might not engender miscleavage, especially when N -1 is not a C. A better understanding of cleavage-site selection as well as the biological specificity of PRORP in vivo requires determining the k cat /K m for cleavage (at C 0 and M -1 ) of different pre-tRNAs that exemplify the natural variations.

Role of the 2'-OH at N -1 in cleavage by AtPRORP1
The 2'-OH at N -1 in pATSer plays an important role for Eco RPR-mediated cleavage at the correct site [32,33].For comparison, we therefore decided to study AtPRORP1 cleavage of pAT-Ser variants in which the 2'-OH at N -1 was replaced with 2'-H (deoxy).These variants,

Role of the pATSer loop in cleavage by AtPRORP1
Since the structure of the pATSerCG loop (corresponding to the T-loop in pSu1; Fig 1) influences cleavage efficiency and cleavage site-selection in Eco RPR-mediated cleavage [21,22], we tested the significance of the loop for AtPRORP1 catalysis.Indeed, replacing the loop in both pATSerCG and pATSerUG with a GAAA-tetraloop reduced cleavage efficiency dramatically (Fig 5 lanes 12 and 13).We detected cleavage at M -2 (between residues -1 and -2) and no cleavage at either C 0 or M -1 for pATSerCG GAAA , and very little cleavage at M -2 (if any) or any other position for pATSerUG GAAA .For both these substrates we were unable to determine k app .Interestingly, comparing the K D values for pATSerUG and pATSerUG GAAA revealed that AtPRORP1 binds both these substrates with roughly equal affinity (Table 2).Collectively, these data suggested that the structure of the loop in pATSer influences cleavage efficiency and site recognition but not binding.

Requirements for efficient and accurate cleavage by AtPRORP1
In addition to pre-tRNAs, PRORPs from various sources are capable of processing mRNAs, tRNA-like molecules called t-elements, and snoRNAs [9,41,42].Here, we investigated processing of mutant derivatives of pre-tRNA Ser (Su1)-based substrates by AtPRORP1 and interpret here the experimentally observed nucleobase preferences with identity biases in the sequences of organellar tRNAs.We have also drawn collectively from two recent complementary reports: Brillante et al. used AtPRORP3 and Thermus thermophilus (Tth) pre-tRNA Gly , and Howard et al. compared AtPRORP1, AtPRORP2 and AtPRORP3 for their ability to process both nuclear and organellar pre-tRNAs (including Arabidopsis mitochondrial pre-tRNA Cys ) [17,18].Overall, these results are expected to contribute to an understanding of the First, we showed that AtPRORP1 cleaves Eco pre-tRNA Ser (pSu1) that has a large variable loop, a structural feature that is known to affect the structural topography in the vicinity of the T-/D-loop region [43].We investigated pSu1 since previous studies with PRORP variants used pre-tRNAs with smaller variable loops (e.g., pre-tRNA Tyr , pre-tRNA Phe , pre-tRNA Gly , pre-tRNA Cys [8,9,17,18,28,29]. Second, we demonstrate here that AtPRORP1 cleaves model hairpin-loop substrates (45-nt pATSer variants) at least 1000-fold slower than the parental pre-tRNA Ser (Table 1).Our finding is consistent with the >1000-fold decrease reported for AtPRORP1-mediated processing of an Arabidopsis mitochondrial pre-tRNA Cys -derived stem-loop substrate compared to pre-tRNA Cys [18]; this decrease was less pronounced with AtPRORP2 and AtPRORP3 (30and 67-fold, respectively).Brillante et al. independently reported a 26-fold lower rate for AtPRORP3-mediated cleavage of a Tth pre-tRNA Gly -derived stem-loop substrate relative to pre-tRNA Gly [17].Given the striking similarity of the tertiary structures of AtPRORP1 and AtPRORP2 [38,44], their differences in processing stem-loop substrates is surprising.Our observations on AtPRORP1 contrast with bacterial RNase P, where the RPR with or without its protein cofactor exhibits only a two-to ten-fold lower activity with a "pAT-type" model substrate compared to its corresponding parental pre-tRNA or even with pre-4.5SRNA [4,21,45].Although these findings suggest that AtPRORP1 might not be capable of efficiently processing substrates such as Eco pre-4.5SRNA [45,46], which resembles pATSer, expression of AtPRORP1 in an E. coli strain that is temperature sensitive (ts) for RNase P activity resulted in growth at the non-permissive temperature [9].
Third, our data suggest that N -1 in the substrate contributes to cleavage efficiency of and site selection by AtPRORP1 (Table 1).Specifically, C -1 decreased the cleavage frequency at C 0 both in the context of pSu1 and pATSer substrates (Fig 3).We consider two possibilities why C -1 might interfere with correct cleavage: (i) the exocyclic amine in C -1 base results in unfavorable positioning in the AtPRORP1 active site; and (ii) formation of a C -1 :G +73 bp imposes a barrier for exposing the C 0 cleavage site, as has been suggested for bacterial RPR [14,15,29].To evaluate these postulates, it is instructive to compare the frequency and rates of miscleavage of pAT-SerCG (C -1 :G +73 ), pATSerCIno (C -1 :I +73 ) and pATSerCU (C -1 :U +73 ).With these three C -1 substrates, we notice a trend towards increasing correct cleavage and a higher overall rate as we transition from three to two to zero hydrogen bonds between N -1 and N +73 ; the four-fold higher preference for miscleavage with pATSerCG shifts to a 1.4-fold preference for correct cleavage with pATSerCU (Fig 3B; Table 1).Thus, both the identity and the strength of the bp at N -1 :N +73 are important in cleavage-site selection.A few additional comments in this regard: AtPRORP1 cleaves chloroplast pre-tRNA Phe with a C -1 :A +73 at C 0 with >95% and Tth pre-tRNA Gly with a C -1 :U +73 only at C 0 [18,29].In contrast, we find miscleavage (40% of total; Fig 3B ) of pATSerCU; akin to the other reports, we find a bias towards correct cleavage.However, it is clear that the impact of C -1 appears to be dependent on context and other structural elements (for instance, shorter D and variable loops in chloroplast pre-tRNA Phe and Tth pre-tRNA Gly , and a larger variable loop in Eco pSu1).Further support for this postulate stems from our sequence analyses (Fig 6; Table B in S1 File).While we noticed a negative bias for C -1 :G +73 in that there were only 0.7% organellar tRNAs from eight different algae/plants, nearly 8% of the total suite have C -1 .We recognize that tRNA nucleobase identities coevolve with a suite of tRNA processing and modification enzymes, including RNase P. As far as AtPRORP1 is concerned, while C -1 :G +73 is clearly not preferred, C -1 alone might be tolerated depending on the N +73 identity and other structural elements (Fig 6 ; Table A in S1 File; see below).
Fifth, we discovered that there is little interplay between N -1 and N +1 in cleavage-site selection by AtPRORP1, a notable difference compared to bacterial RNase P. G -1 -containing pSu1 and pATSer variants are cleaved with a high frequency at M -1 by bacterial RPR.This is particularly true for substrates having G -1 :C 73 (e.g.pre-tRNA His ) [14, 15, 31, 53-57; Mao and Kirsebom, unpublished data].G +1 , which has been suggested to help position the nucleophile during RNase P-mediated cleavage, is indeed present in a majority of bacterial tRNAs (see e.g., http://trna.ie.niigata-u.ac.jp/cgi-bin/trnadb/index.cgi.)[34,49,58].Thus, the presence of G -1 leads to increased cleavage at M -1 [15], likely due to metal-ion or other anchoring determinants now being present at both G +1 and G -1 .Although organellar tRNAs from the eight green algae and plants that we analyzed also favor G +1 (nearly 75% bias; Table B in S1 File), it appears that AtPRORP1 might not rely on G +1 as a guide for cleavage-site selection.Unlike bacterial RPR, which cleaves pSu1 and pATSer (having G -1 and G +1 ) at the incorrect site with either higher or similar frequencies as substrate counterparts with C -1 [14, 15; Mao and Kirsebom, unpublished data], AtPRORP1 cleaved a G -1 substrate mainly at the correct site C 0 [for example, see pSu1 (-1G), Fig 2].However, there is a hierarchy in cleavage-site selection by PRORPs, with contributions from multiple factors such as N -1 identity and the N -1 :N +73 bp (especially, a G -1 :C +73 bp) as evident from the following observations.In spinach chloroplast pre-tRNA His , G -1 is encoded in the gene; 5' processing of this precursor using a spinach S100 extract results in a 5'matured tRNA His with G -1 [59]; similarly, recombinant AtPRORP1 cleaved potato tRNA His predominantly between G -2 and G -1 [60].This scenario with plants contrasts with yeast tRNA His , where G -1 is added after RNase P processing [61].Also, AtPRORP1, AtPRORP2 and AtPRORP3 mis-cleave (at a frequency ranging from 28% to 72%) A. thaliana nuclear pre-tRNA Phe with U -1 :A 73 [18].Swapping the native C -1 :U 73 in pre-tRNA Gly to G -1 :C 73 led to 100% mis-cleavage at M -1 by AtPRORP3 [17].
Sixth, comparing the K D and k app values, respectively, for binding and cleavage of pre-tRNA Ser Su1 and pATSerUG by AtPRORP1 revealed the importance of the D-loop, the variable loop and the anticodon stem and loop (  2).Replacement of the native T-loop (seven nt) with a GAAA tetraloop in pATSer did not affect binding but eliminated cleavage at the correct position C 0 for both the C -1 and U -1 variants (Table 2).Hence, at least for cleavage of model hairpin-loop substrates, the T-loop equivalent contributes to the rate and fidelity but not binding.Our observations, which emphasize the importance of the T-/ D-loop region for binding and processing by AtPRORP1, are consistent with findings from earlier studies.Substitution of residues at positions 18 or 19 in the D-loop, or 56, 57 or 58 in the T-loop influenced the cleavage efficiency of AtPRORP1 [9,62].A substrate in which the anticodon stem and loop is deleted was cleaved with high efficiency, whereas removal of the Dloop resulted in an RNA for which no detectable cleavage was observed [9].Footprinting analysis of pre-tRNA Cys further indicated that U16, G18, G19 (D-loop) and C56 (T-loop) are protected when bound to AtPRORP1 [11].Moreover, the K M(STO) and k react (kinetic constants determined under single turn over) of AtPRORP3-mediated processing of pre-tRNA Gly decreased by 1200-and 26-fold, respectively, upon deletion of the D-stem-loop and anticodon stem-loop [17].Taken together, it is clear that efficient and correct cleavage depends on a productive interaction between the T-/D-loop region and AtPRORP1, as has been shown for bacterial RPR [7,21,22; see also 51,[63][64][65].
Last, we find that the absence or presence of 3'-CCA in either pSu1 or pATSer does not affect cleavage-site selection by AtPRORP1 (Fig 3).The rate of cleavage, however, does change by two-to three-fold but not in any predictable fashion with the various substrates that we studied (Table 1).AtPRORP1-mediated cleavage of a plant mitochondrial pre-tRNA Cys was shown to be inhibited by the presence of a 3'-RCCA motif [11].Together, these results with AtPRORP1 emphasize an important difference compared to bacterial RNase P, where the rate and fidelity of cleavage are dramatically affected when the 3'-CCA is deleted from either a pre-tRNA or "pAT-type" substrate; these results are expected due to the base pairing between the 3'-RCC of the pre-tRNA and a conserved GGU-motif in the RPR [13,66].Unlike bacterial pre-tRNAs, a 3'-CCA was predicted to be present in the initial pre-tRNA transcript for only 0.5% of the total suite of 423 organellar genome-encoded pre-tRNAs in eight plant/algal species.Thus, the 3'-CCA is unlikely to be a major contributor to AtPRORP1 catalysis (see also [18]).

Substrate recognition by the bacterial ribozyme variant and AtPRORP1
For bacterial RNase P, biochemical and genetic studies have provided insight into substrate recognition features, and these were confirmed and extended by the crystal structure of the bacterial RNase P-tRNA complex [51]: (i) N -1 in the pre-tRNA has a key role and might interact with a specific base in the RPR; (ii) the 2'-hydroxyl N -1 is used to coordinate metal ions essential for catalysis; (iii) G +1 in the pre-tRNA acts as a guidepost in the RNase P-substrate complex; (iv) the T-loop in the pre-tRNA is specifically recognized by an architectural motif of two inter-digitated T-loops in the RPR; and (v) 3'-RCC sequence of the pre-tRNA pairs with a conserved GGU sequence in the RPR [51]; for reviews see e.g.[1,11].The anticodon stem-loop was shown to be dispensable, which is expected given that all tRNAs are processed by RNase P. Two previous models [11,42] show how AtPRORP1 might use the "acceptor-T-stem" stack as the main recognition determinant, an idea that is supported by our finding here that the pAT-Ser-type variants, which possess only the acceptor-T-stem stack element, are cleaved by AtPRORP1 with the same fidelity as the parental tRNA counterparts (Fig 3 and Table 3).
With respect to pre-tRNAs, however, AtPRORP1 is likely to interact with the amino acceptor-stem and the T-/D-loop region (Table 3).It is possible that the distance between T-/Dloop region and cleavage site determines metal-ion binding and cleavage-site selection by AtPRORP1 [11], but this notion might need if AtPRORP1 (like the bacterial RPR) accepts substrates with shorter acceptor-/T-stem stacks.Additional insights are needed to ascertain whether this amounts to a measuring mechanism that has been suggested for bacterial and eukaryotic RNase P [58,63,67,68].Since we observe binding of both pATSerUG and pATSerUG GAAA but cleavage of only the former (Table 2; Fig 5), an induced-fit mechanism based on T-loop recognition is likely with AtPRORP1, again mirroring a proposal for bacterial RNase P [7,21,22].Bacterial RPR/RNase P uses multiple determinants to define its cleavage site whereas AtPRORP1 appears to employ fewer elements and differs notably in not using either the 5'leader or the 3'-trailer (Table 3) [17,18].While this difference might signify how binding energy and cleavage-site selection are accomplished by nucleic acid-versus protein-based RNase P, it likely reflects the culmination of a catalytic strategy based on the co-evolution of each catalyst with its entire suite of substrates not just pre-tRNAs.Both forms of RNase P have honed in on the common denominators in all pre-tRNAs: the acceptor-T-stem stack and the T-/D-loop interaction [69], which incidentally is used as a recognition determinant by other RNAs and proteins that act on tRNA [70].

Fig 1 .
Fig 1.Secondary structures of substrates used in this study.Secondary structures of pSu1 and pATSer.The highlighted regions/residues were substituted to generate the different variants as indicated, A, adenosine, G, guanosine, U, uridine, I (Ino), inosine and D (DAP), 2,6-diaminopurine; dC, deoxycytosine; and dU, deoxyuridine.The canonical RNase P cleavage sites between residues N -1 and N +1 (correct cleavage denoted C 0 ), and the alternative cleavage sites between residues N -2 and N -1 (miscleavage denoted as M -1 ) are marked with black and grey arrows, respectively.The N +73 position, which immediately precedes the 3'-terminal CCA-motif, corresponds to the discriminator base.

Fig 3 .
Fig 3. Frequencies of cleavage-site selection by AtPRORP1.Histograms summarizing cleavage-site selection frequencies (in %) during AtPRORP1-mediated cleavage of pSu1 "-1" (A) and pATSer (B) variants.Mean and standard deviation values were calculated using data from at least three independent experiments.doi:10.1371/journal.pone.0160246.g003 Fig 1), and assessed their fidelity and rate of cleavage by AtPRORP1 (Fig 2, lanes 14 to 17; Fig 3A and Table Fig 1).Indeed, pATSerCG acts as a substrate for AtPRORP1 (Fig 5, lane 9), and as expected based on the fidelity of processing of pSu1(-1C), pATSerCG was also cleaved mainly at M -1 (Fig 3B).Substitution of C at -1 with

Fig 4 .
Fig 4. Effect of varying Mg 2+ concentration on AtPRORP1-mediated cleavage.AtPRORP1-mediated cleavage of the pSu1(-1U) (A) and pATSerUG (B) as a function of Mg 2+ concentration at 37°C.Mean and standard deviation values were calculated using data from at least three independent experiments.doi:10.1371/journal.pone.0160246.g004

#C
and U correspond to residue identity at the -1 position while G, Ino (inosine) and U refer to residue identity at the discriminator position "+73" (numbering same as in tRNA; Fig1).## k app values determined at 25 mM Mg 2+ for these substrates.While these experiments were performed prior to our establishing 10 mM Mg 2+ as being optimal, the rate and fidelity of cleavage is largely unchanged between 10 to 25 mM Mg 2+ .Δ3'AC, Δ3'CAC and Δ3'CCAC indicates residues in the 3'CCAC motif that were deleted.ND, not determined; NA, not applicable.doi:10.1371/journal.pone.0160246.t001

Fig 5 .
Fig 5. AtPRORP1-mediated cleavage of pATSer variants.Representative gel showing AtPRORP1-mediated cleavage of 3' CCA-motif-containing pATSer variants.Lanes 1 to 6 and 14 to 17 are negative controls (loaded in the same order as the reactions with AtPRORP1 in lanes 8 to 13 and 19 to 22, respectively); and lanes 7 and 18 (size marker) refers to cleavage of pATSerUG by Eco RPR.The final concentration of AtPRORP1 was 6.6 μM and the reactions were performed at 37°C for 60 min in the presence of 10 mM Mg 2+ .The position of each 5' cleavage fragment (5' CL Frags) generated after cleavage is indicated.The two lower panels represent overexposure to better highlight the 5'-cleavage products in the upper panels.(Note: Fig A in S1 File shows cleavage of pATSer derivatives without the 3' CCA-motif.)doi:10.1371/journal.pone.0160246.g005 Fig B in S1 File).This change, which corresponds to a loss of 3.2 kcal/ mol in binding, reflects the importance of the D stem-loop, and perhaps the T-/D-loop tertiary contacts, for tight substrate binding by AtPRORP1.The model substrate also lacks the anticodon stem-loop, but this structural element has been shown to be dispensable for substrate recognition and cleavage by the RNP and AtPRORP forms of RNase P[4,6,9,13,20,39].
pATSerU deoxy G and pATSerC deoxy G (Fig 1), were subjected to cleavage by AtPRORP1.For the variants carrying 2'-H at N -1 , we did not detect any cleavage at any position (Fig C in S1 File).This result is in contrast to what has been reported for Eco RPR-mediated cleavage of the same substrates (see Discussion).

Fig 6 .
Fig 6. Analysis of N -1 :N +73 identities in mitochondrial and chloroplast tRNAs.Analysis of N -1 :N +73 identities in 423 mitochondrial and chloroplast tRNAs from eight different green algae and plants (sequences obtained from http://plantrna.ibmp.cnrs.fr/).Table A in S1 file lists the individual distributions in each species.doi:10.1371/journal.pone.0160246.g006 Fig 5 and Fig A in S1 File; Table

Table 1 .
Rate of cleavage (k app ) of pSu1 and pATSer variants at 10 mM Mg 2+ .
Each value listed is a mean ± standard deviation determined from three or more independent experiments.