|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 15, 11410-11426, April 13, 2007
In Vitro Comparative Kinetic Analysis of the Chloroplast Toc GTPases*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
-proteobacteria or cyanobacteria. During the course of evolution, in both cases, the vast majority of the endosymbiont genes was relocated to the host nucleus (1). Moreover, in higher plants, most of the proteome from both organelles is derived from nuclearly encoded proteins that are translated on free cytosolic ribosomes and then post-translationally translocated into the organelle via distinct protein complexes located in each of the two membranes that enclose these organelles (2, 3). These mitochondrial and chloroplast translocators are denoted as Tim/Tom and Tic/Toc, respectively (Translocator of the inner/outer membrane of mitochondria/chloroplast). The Toc5 complex consists of three key proteins denoted by their apparent molecular masses, Toc34, Toc75, and Toc159, and most likely exist with a stoichiometry of 4:4:1, respectively (4). Proteins targeted from the cytosol to these mitochondrial and chloroplast translocators require additional "information" in the form of an N-terminal targeting sequence known as a presequence and transit peptide, respectively (1). This additional targeting sequence has been added during evolution, yielding a larger precursor protein, and is cleaved once the precursor is imported into the organelle (5).
In plant cells, despite colocalization of these two post-translational targeting processes, very little mis-sorting occurs in vivo or in vitro (6, 7). The high fidelity of these two processes is maintained, in part, by the specific recognition of the targeting sequences by the import receptors associated with the translocator. In mitochondria, it has been shown that the outer membrane receptor, Tom20, recognizes and binds the presequence (8, 9). In chloroplasts, both Toc34 and Toc159 contribute to import specificity through the coordinate or sequential recognition of the transit peptide (10-12). Although both the Tom and Toc receptors function to promote preprotein binding specificity, only the Toc receptors perform as a molecular switch and/or molecular motor (13). These additional activities are derived from the intrinsic GTPase activity inherent in Toc34 and Toc159 (14). Further regulation of protein translocation has been suggested to occur via specific phosphoryation/dephosphorylation of the Toc GTPases (15, 16); however, the in vivo significance of this phosphorylation has not been established (17).
From sequence analysis of Arabidopsis, the Toc GTPases fall into two gene families as follows: the small Toc GTPases (Toc33 and Toc34) (18) and the large acidic Toc GTPases (Toc90, Toc120, Toc132, and Toc159) (19, 20). However, in pea, only Toc34 and Toc159 isoforms have been detected. The small GTPases are inserted into the chloroplast outer membrane via their C-terminal transmembrane domain and contain a globular GTP-binding domain that contains the canonical G1 and G3-G5 domains involved in GTP/GDP binding similar to the Ras family of small GTPases (21). Arabidopsis Toc159 contains a central GTP-binding domain (residues 728-1092) with high sequence similarity (
63%) to the small Toc GTPases; however, Toc159 also contains two additional domains as follows: a large proteolytically sensitive N-terminal acidic domain (residues 1-727) and a C-terminal membrane-anchoring domain (residues 1093-1503) (22). The other acidic Toc GTPases share this domain structure, yet the acidic regions are reduced in length (23). Very little is known about the structure and function of the acidic and membrane domains.
Although it has been shown in pea chloroplasts that an early preprotein binding intermediate is mediated in a GTP-dependent manner by the Toc complex (24), the details of how the two different GTPases, Toc34 and Toc159, participate in binding/translocation is debated. Several reports on the interaction of the Toc complex with preproteins have been published recently and result in the following two models: one in which a soluble form of Toc159 interacts with the preprotein in the cytoplasm and is subsequently targeted to Toc34, and Toc75 at the outer envelope to begin protein translocation (25); however, these results have been recently challenged (26). The second model proposes that the preprotein first interacts with Toc34 at the chloroplast outer membrane surface and is then "handed over" to translocon-associated Toc159 (15). After the recognition event involving transit peptide and/or preprotein interaction with the cytosol-exposed portions of Toc34 and Toc159 and concurrent rounds of GTP hydrolysis, the preprotein is inserted into the
-barrel protein Toc75 and translocated across the outer membrane of the chloroplast (27, 28).
Despite multiple reports investigating the interaction of the Toc34 and Toc159 homologues with preproteins and transit peptides, the role of GTP binding and hydrolysis is far from understood (1, 4, 5, 10-12, 15, 16, 18, 21, 26, 29-33). In fact, the few reports on GTPase activity have been primarily focused on the small GTPases and report enzymatic properties such as Vmax and KM values that disagree by nearly 2 orders of magnitude (11, 21). In Arabidopsis, the two related homologues, Toc33 and Toc34, have had their nucleotide binding properties investigated, yet there have not been systematic investigations on their potentially different enzymatic properties (33). Finally, atToc159G has had only limited analysis performed on its nucleotide binding and/or hydrolysis properties (31, 34). Mechanistically, from both the crystal structure (21) and the low resolution EM structure of the Toc translocon (4), it is clear that the Toc GTPases may be able to form either a functional homo- or heterodimer. The significance of this dimerization is still debated (35), yet it has been proposed that one possible explanation is that the two subunits can function as GTPase-activating proteins that may reciprocally stimulate and or coordinate the hydrolysis of GTP (21). This dimeric, self-stimulatory interaction was first reported in the crystal structure where the presence of Arg-133 from one monomer interacting deep within the GTP-binding pocket of the other monomer can be observed. Reciprocal GAP activity is not a novel idea; for instance, reciprocal GAP activity has become an accepted model for the function of SRP54/Ffh and SR
/FtsY in protein secretion (36), the only other family of small GTPases associated with protein translocation (37).
Therefore, to enable a careful analysis of how preprotein binding, membrane association, and homo- and heterodimerization influence the GTPase activity of the various Toc proteins, we first present a systematic comparative kinetic analysis of the Toc34 subfamily and atToc159G. We report the Vmax, KM, Ea, and kcat values for GTP hydrolysis, as well as the Kd and Ki values of nucleotide binding for each protein. Using psToc34 as a model, we also evaluate how concentration-dependent homodimerization, an "arginine finger" mutation, and transit peptide interaction affect subunit dimerization and GTP hydrolysis for psToc34. Finally, using the crystal structure of psToc34, we have developed molecular models for atToc33, atToc34, and atToc159G. Comparing these models to each other and to p21 Ras and other small GTPases may provide new insights on possible mechanisms of how the Toc GTPase cycle may be modulated during chloroplast protein import.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
-D-galactopyranoside for 4 h at 37°C. The E. coli cells were lysed by three passages through a French press (Thermo Electron Corp., Waltham, MA) at 20,000 p.s.i. in the presence of lysis buffer (20 mM NaPO4 (pH 8.0), 0.5 M NaCl, 1 mM MgCl2, 0.1% Triton X-100, 1 mM
-mercaptoethanol). The soluble protein was separated from inclusion bodies and cell debris by centrifugation at
45,000 x g for 30 min. The His-tagged proteins were purified under nondenaturing conditions using Talon Co2+ IMAC (BD Biosciences). The bound proteins were washed with 15 mM imidazole and eluted with 250 mM imidazole. Eluent was immediately diluted with 2.5x GBS and 50% glycerol yielding a storage condition of 1x GBS (20 mM Tricine-KOH, pH 7.65, 1 mM MgCl2, 50 mM NaCl, 1 mM
-mercaptoethanol) (16) with 20% glycerol and frozen at -80 °C until used. The concentrations of the purified recombinant proteins were determined using the Bradford assay (Bio-Rad). The purity of the proteins was estimated by scanning densitometry of a Coomassie Brilliant Blue-stained SDS-polyacrylamide gel to be greater than 95%.
Purification of the Full-length Transit PeptidesThe subcloning of the sequence encoding the transit peptide of the small subunit of ribulose-bisphosphate carboxylase/oxygenase from Nicotiana tobacum (SStpNt) into PTYB2 vector (New England Biolabs, Beverly, MA) was described previously (10). The growth and purification of SStpPs were performed as described previously for SStpNt (38). The eluted peptide was lyophilized to remove the
-mercaptoethanol and concentrate the protein. Upon dissolving the lyophilized peptide in water, the contaminating DnaK remained insoluble and could be removed by centrifugation (Fig. 1C). The purity of the peptide was initially determined by SDS-PAGE on a 19.2% Tris/Tricine gel and staining with Coomassie Brilliant Blue. The peptide was quantified using the BCA protein assay (Pierce) and then realiquoted and lyophilized in small, 0.5-mg aliquots. The transit peptide was then stored as a lyophilized powder at -20 °C until used. The peptide was resuspended in double distilled H2Oto the appropriate concentration for enzyme assays.
MALDI-TOF Mass SpectrometryMALDI-TOF MS was performed on a Bruker Daltonics MicroflexTM mass spectrometer, with apomyoglobin used as an internal standard for calibration. Lyophilized standard protein and transit peptide were resuspended to a final concentration of 50 pmol/µlin10 µl of 50% acetonitrile (v/v) with 0.1% trifluoroacetic acid (v/v). The samples were mixed with 10 µl of
-cyano-4-hydroxycinnamic acid matrix (10 mg/ml in 50% ACN (v/v), 0.1% trifluoroacetic acid (v/v)), spotted, and dried thoroughly. Mass spectra were acquired in positive ion mode.
Phosphate Release Assay for GTP HydrolysisGTP hydrolysis of psToc34, atToc33, atToc34, and atToc159G were all performed in the same manner using a modification of an activated charcoal, phosphate-release assay originally developed by our laboratory (39). In a flat-bottom 96-well microplate, reactions were carried out in GBS, 10 nM [
-32P]GTP, and varying amounts of cold GTP (and any additional substrate) in a final volume of 100 µl. The addition of Toc protein started the reaction, and at various times, as indicated in the figures, a 12.5-µl aliquot was removed and added to 200 µl of 10% w/v activated charcoal in 50 mM HCl, 5 mM H3PO4. Vigorous pipetting with a multichannel pipette mixed the charcoal/aliquot mixture. At the completion of the reaction, the charcoal suspension was transferred to a 96-well 0.2-µm filter plate (Pall Corp., Ann Arbor, MI) and vacuum-filtered into a receiving 96-well plate. An 8-µl aliquot was taken of the filtrate and added to 150 µl of MicroScint 40 scintillation fluid (PerkinElmer Life Sciences). Again, vigorous pipetting was performed to ensure even distribution in the scintillation fluid. Counting was performed on a 96-well microplate scintillation counter (TopCount NXT Microplate Scintillation Counter, PerkinElmer Life Sciences). In all cases, minimums of three independent experiments were performed, counted in duplicate for each data point, and crosstalk filters applied. Counts/min data were graphed as counts/min versus time using the Prism software (Prism 4c, GraphPad Software, San Diego).
Determination of Enzyme Velocity and Substrate AffinityThe conditions were designed to ensure that the reaction rate was linear over the 30-min time course of the assay. By using multiple time points,
cpm/min values were calculated from the slope of the line generated by counts/min versus time. With the known specific activity of the isotope, the isotope purity, the enzyme concentration, and the counting efficiency, the slope was then transformed into a velocity expressed as nanomoles of GTP hydrolyzed/min/µmol of Toc protein. To determine the Vmax and KM values for GTP hydrolysis, the GTP concentration was varied from 1 to 100 µM. Once the rate values had been calculated for each GTP concentration, these rates (y axis) are graphed against the substrate concentration (x axis) and fit using the Michaelis-Menten equation (Equation 1) to determine the Vmax and KM values using the GraphPad Prism program,
![]() |
where v is the velocity; Vmax is the maximum velocity; KM is the Michaelis-Menten value, and [S] is the substrate concentration.
Optimization of Buffer Composition, pH, and Storage Conditions for Toc GTPasesBuffers used to determine optimal storage conditions for the Toc proteins were as follows: GBS, HMK (50 mM HEPES-KOH (pH 7.5), 2 mM MgCl2, 40 mM KOAc, 12.5 mM imidazole, 10% glycerol, 1 mM
-mercaptoethanol) (34); KSS (50 mM Tris-HCl (pH 8), 50 mM KCl, 5 mM MgCl2, 5 µM GTP, 1 mM dithiothreitol) (25); and TTP (20 mM Tris-HCl, pH 8.0, 100 mM NaCl, 1 mM MgCl2, 10 mM imidazole, 10% glycerol, 1mM
-mercaptoethanol), a variation of the buffer used during the crystallization of psToc34 (21). In addition, psToc34 was subjected to the following freezing conditions: -80 °C slow, -80 °C fast, -20 °C, and storage in liquid nitrogen; "fast" and "slow" refer to the relative speed at which the samples reached -80 °C. Fast-treated samples were snap-frozen in liquid nitrogen for 30 s prior to being placed in a -80 °C freezer, and slow-treated samples were placed directly into a -80 °C freezer in a foam insulating container. Furthermore, psToc34 was incubated in GBS titrated through a pH range that extended from 6.1 to 11.9. GBS was made with the following buffering agents for a given pH range: MES (pH 6.1-6.7), HEPES (pH 6.8-7.5), Tricine (pH 7.6-8.4), Tris (pH 8.5-9.1), and CAPS (pH 9.2-11.9). At each pH titration step, psToc34 was assayed for maximal velocity as described earlier.
Effect of Cation Replacement on psToc34 ActivitypsToc34 was incubated in 1x GBS buffer in the presence of 5 mM (saturating conditions) selected metal salts as follows: MgCl2, CaCl2, CdCl2, CoCl2, CuSO4, MnCl2, NiSO4, and ZnCl2 for 10 min at 25 °C to allow for cation replacement. A mixture of both [
-32P]GTP (10 nM) and variable amounts of cold GTP was added to start the reactions. The hydrolysis activity of psToc34 was measured as described above and graphed as a percentage of psToc34 activity in GBS buffer.
Effect of Protein Concentration on the Hydrolysis Rate of the Toc GTPasesThe four Toc proteins were assayed as described above in GBS at 1-3, 5, 10, 15, and 20 µM protein concentration for hydrolysis activity. The resulting data are plotted as rate of hydrolysis as a function of protein concentration in micromolars.
Construction of psToc34(R133A) MutantMutation of Arg-133 was performed with the following mutagenesis primers from 5' to 3': GATGCCTACGCGGTGGACAACCTGGAC forward and GGTTGTCCACCGCGTAGGCATCCAGC. All mutagenesis was carried out with QuickChange PCR mutagenesis (Stratagene).
Circular Dichroism Spectroscopic Analysis of psToc34 (R133A)Circular dichroism spectroscopy was performed on an Aviv Series 202 circular dichroism spectrometer (Aviv Biomedical, Lakewood, NJ). After dialyzing extensively into 50 mM NaH2PO4, pH 7.8, 10 mM NaF, the final concentration of the psToc34(R133A) mutant protein was determined to be 7.58 µM using the Bradford assay. Circular dichroism spectra were collected at 25 °C as 2-s averages at 1 nm intervals from 185 to 285 nm. Five spectra were averaged, corrected for buffer contributions, smoothed, and converted to molar ellipticity using the Aviv software, version 2.71. Deconvolution was performed as described previously (40) using the CDPro software with IBase3, a reference set of 37 soluble proteins, available on line. From these data, the secondary structural features of psToc34(R133A) were determined. The Stride program (41) was employed to calculate the percent secondary structural elements of crystallized wild-type psToc34, chain B, PDB code 1H65 (21).
Sedimentation Velocity Analysis of psToc34Sedimentation velocity analyses were performed using a Beckman Optima XL-I analytical ultracentrifuge equipped with both absorbance and interference optical detection systems. After dialyzing samples extensively into GBS as described above, 400 µl was loaded into Epon double sector centerpieces. For experiments at different concentrations, all dilutions were made using the dialysate, and all experiments used dialysate as the reference buffer. Sample cells were loaded into the rotor and allowed to equilibrate to 25 °C under vacuum for at least 1 h. After being fully equilibrated, samples were centrifuged at 50,000 rpm at 25 °C using the An50Ti rotor, and data were collected using interference optics. Analysis of the sedimentation velocity data was performed using the program Sedfit, which uses Lamm equation modeling to carry out size distribution analyses of velocity profiles (42). Velocity profiles were analyzed by direct boundary modeling using distributions of Lamm equation solutions, c(s) (43, 44). The solvent viscosity,
, and the density,
, were experimentally determined to be 0.00896 g cm-1 s-1 and 1.0003 g ml-1, respectively, at 25 °C using an Anton Parr densitometer by Prof. Jochen Weiss, University of Massachusetts, Amherst. The partial specific volume of the protein,
, was calculated to be 0.7441 ml g-1 by inputting the protein sequences into the program Sednterp (45). The nonlinear regression using distributions of Lamm equation solutions, c(s), was combined with the determinations of the best fit weight average frictional ratio, the best fit meniscus position, and the algebraic calculation of systematic time invariant and radial invariant noise components and maximum entropy regularization using a confidence level of >0.8 (42). The resulting sedimentation coefficient distributions, expressed in Svedberg units (S), were exported, normalized to respective protein concentrations, and graphed using GraphPad Prism 4c (GraphPad Software, Inc.).
Activation Energies of the Toc GTPasesReactions were performed ranging in temperature from 4 (277 K) to 45 °C (318 K) in 0.6-ml microcentrifuge tubes. Prior to initiation of reaction by addition of Toc protein, the temperature of the reaction buffer and nucleotide was allowed to equilibrate for 10 min. Rates were determined as above in hydrolysis assays and graphed as the natural log of the rate (lnk) versus the inverse of the absolute temperature, expressed in Kelvin (103/T, K-1). Ea value was determined using the Arrhenius equation (Equation 2),
![]() |
where k is rate in nanomoles of GTP hydrolyzed per min per µmol of Toc protein; R is the universal gas constant; T is temperature in Kelvin; and Ea is the energy of activation expressed in kcal/mol.
Nucleotide Binding Studies Using NTA-Magnetic AbsorptionUsing iron-incorporated agarose beads with Ni-NTA chemistry (Qiagen Inc., Valencia, CA), we were able to recapture the various Toc proteins via their C-terminal hexahistidinyl tags in a 200-µl reaction volume. This approach affords us the ability to selectively immobilize and extensively wash each protein in the presence of a 96-well hybrid plate magnet. Competition assays for nucleotide binding were performed in GBS buffer with each purified Toc GTPase. The assay contained 20 nM of the GTPase, [
-32P]GTP (MP Biomedicals, Irvine, CA) with a specific activity of 3000 Ci/mmol, an excess of Ni-NTA-agarose magnetic particles, and increasing concentrations of competing, unlabeled nucleotide as indicated in the Fig. 6 legend. The conditions were optimized to cover >3 orders of magnitude of the competitor with sufficient data points to permit an accurate nonlinear fit. The reaction was carried out in a 96-well microplate in a total volume of 100 µl. First, a serial dilution of unlabeled competitor was carried out in a separate plate leading to a concentration 4-fold greater than indicated on the x axis. Then 25 µl of the resulting solution was transferred with a 12-well multichannel pipette into each row of wells in the reaction plate. The [
-32P]GTP was diluted and added to each well of the reaction plate, also as a 25-µl addition. Finally, the Ni-NTA-agarose magnetic particles and enzyme were added to start the reaction in a volume of 50 µl. The reaction proceeded for 30 min on ice before the plate was transferred to a novel high performance hybrid plate magnet. We used a new class of magnet plates recently developed at the Department of Energy Joint Genome Institute at Lawrence Berkeley National Laboratory for biological and industrial applications (46). They utilize hybrid technology that combines linear permanent magnetic and ferromagnetic material to produce significantly higher fields and gradients than currently available on the commercial magnetic plates. This hybrid structure exhibited maximum fields in excess of 9000 G, which allowed for greater holding forces on magnetized targets that are being processed as well as permitting faster draw down. Following immobilization in high field strength, the reaction mixture could be easily removed by aspiration and replaced with 300 µl of GBS. Each well could then be washed continuously with GBS for 2-3 min with a microtiter plate washer (Skatron Microwash II, Norway) to remove all nonspecific background. The wash buffer was then aspirated and replaced with 100 µl of elution buffer (GBS with 500 mM imidazole) and removed from the magnet for 5 min with periodic agitation. After replacement on the magnet for 5 min to ensure immobilization of the beads, 20 µl of the elution buffer was removed for scintillation counting. These 20-µl samples were removed from each well in triplicate, and each competition reaction was performed in duplicate. In addition, to avoid the ambiguous results associated with homologous binding competitions performed at a single concentration of labeled nucleotide, at least two concentrations of [
-32P]GTP were used in each case. Thus, the entire competition experiment was repeated using a 3-fold higher concentration of the labeled nucleotide. The samples were counted in a 96-well format using the Top Count scintillation counter (PerkinElmer Life Sciences) with Microscint fluid, and cross-talk filters were applied before the data were analyzed.
Global Fit Analysis of Nucleotide Binding DataUsing the experimental procedure described above, there were a total of 12 data sets per global fit analysis (i.e. two replicates, each counted in triplicate at two different nucleotide concentrations). These values were fit to Equation 3 using Graphpad Prism 4.0c with the log Kd value shared for global fitting and the amount of hot and cold nucleotide constrained according to the experimental design.
![]() |
The remaining variables were given best fit values by nonlinear regression. For heterologous binding competitions, the experiments were performed in the same manner, but the analysis was fit by Equation 4,
![]() |
In this case the Kd value measured in the previous experiment was substituted into the equation, and the log Ki value was shared for all data sets for global fitting. The amount of labeled and nonlabeled GTP was held constant, and the remaining variables were given best fit values by nonlinear regression. For data analysis, a confidence of 95% was used.
Peptide Stimulation of Hydrolysis ActivityHydrolysis experiments were performed with psToc34 and the addition of increasing concentrations of GTP, as described above, in the presence of a fixed concentration of SStpPs, SStpNt, or the negative control RCMLA to generate linear counts/min versus time data. Linear regressions were applied generating a series of slopes for each concentration of GTP used. This was repeated for multiple concentrations of all peptides as follows: 2.5, 4, 6.3, 10, 15.8, 25, 40, 63, and 100-fold molar concentrations to psToc34. The slopes from the counts/min data were mathematically transformed to rate data and fit to Equation 1 to determine Vmax and KM values.
Effect of Transit Peptide on Nucleotide ExchangepsToc34 (0.5 mg) was incubated with [
-32P]GTP (0.5 mCi of 800 Ci/mmol, MP Biomedicals) and EDTA in GBS for 1 h at room temperature in a 500-µl reaction volume. The unbound nucleotide, EDTA, and imidazole were removed by loading the mixture onto a PD-10 desalting column (Amersham Biosciences) and eluted in GBS;
2% of the psToc34 was routinely labeled. The eluent was collected in 60-µl fractions; psToc34-containing fractions began after
3.2 ml of elution, and the concentrations of these fractions were determined using the Bradford assay. The requisite amount of psToc34 was added to create a 200 nM solution, and 50 µl was aliquoted into each well of a microtiter plate. A time course of the dissociation of GTP label was started with a 50-µl addition of 100-fold GTP or 100-fold GTP and transit peptide at 0.3-, 1.5-, 7.5-, or 37.5-fold molar excess over psToc34; GBS was used as a negative control. After each of the wells had incubated for the indicated time, the plate was placed on the hybrid magnet and then washed 12 cycles in a microtiter plate washer. The remaining GBS was aspirated away and replaced by GBS with 50 µl of 500 mM imidazole. Half of the elution was transferred into 100 µl of MicroScint scintillation fluid and counted in the TopCount (Packard Instrument Co.). Data were fit using a single-phase exponential decay in Graphpad Prism 4.0c.
Homology Modeling of Toc ProteinsatToc33, atToc34, and atToc159G were modeled using the Modeler modeling package (47) with the crystal structure of psToc34 as a template. The coordinates were obtained from the Brookhaven Protein Data Bank (PDB code 1H65). The sequence of each homologous Toc protein was aligned using the Modeler algorithm Salign, which uses a variable gap penalty to avoid artificially introducing insertions into conserved areas of secondary structure. Because the percent identity was 62.9, 66.0, and 32.0 for atToc33, atToc34, and atToc159G, respectively, the automodel scripts of Modeler were used to generate 20 candidate models for each Toc protein ignoring out-gaps. These models were evaluated using the discrete optimized potential energy function from Modeler, and all amino acids aligned with psToc34 were compared summing the square of the differences relative to the discrete optimized potential energy function of the psToc34 crystal structure. All candidate models agreed closely with the discrete optimized potential energy function of psToc34, but for each set of candidate models, the lowest global score was chosen for further analysis. All structural models were rendered using Virtual Molecular Dynamics (48), and structural alignments were performed with Stamp (49). r.m.s.d. values of the
-carbons were calculated relative to psToc34 using the r.m.s.d. plug-in of VMD, and per residue r.m.s.d. values were calculated by the r.m.s.d. calculator of Multiseq (50). All figures were ray traced in POV-Ray 3.6.
| RESULTS |
|---|
|
|
|---|
TM and atToc33
TM (E. coli TG1) and psToc34
TM and atToc159G (E. coli BL21[DE3]) as soluble protein and purified via IMAC using the C-terminal hexahistidinyl tag (Fig. 1A). To verify purity, 100 ng of each protein was run on 15% acrylamide SDS-PAGE and concluded to be greater than 95% pure (Fig. 1B).
|
-mercaptoethanol (Fig. 1C). Success of cleavage was assessed by removing chitin matrix, boiling with sample solubilizing buffer, and running on SDS-PAGE (Fig. 1C, lane 7). The upper band in Fig. 1C, lane 7, is due to uncleaved transit peptide still bound to the chitin matrix via the chitin-binding domain. Cleaved protein was purified to homogeneity, as confirmed by gel electrophoresis (Fig. 1C, lane 9) and MALDI-TOF mass spectrometry (Fig. 1D). MALDI-TOF MS of transit peptides was performed with apomyoglobin used as an internal standard on a Bruker-Daltonics MicroflexTM mass spectrometer shown in Fig. 1D. The +2 and +3 charged apomyoglobin species can be seen at 8477 m/z and 5651.7 m/z, respectively; the average molecular weight of apomyoglobin is 16,952. The +1 and +2 charged SStpNt species are found at 6125.9 and 3062.6 m/z. The FindPept tool was used to identify the peptides that correspond to the main peaks. The largest peak at 6125.9 m/z corresponds to the full-length transit peptide, whereas the minor peak at 4152.6 m/z represents either an alternative translation initiation at Met-22 or alternatively a degradation product where short N- and C-terminal regions were proteolytically removed during expression and purification (Fig. 1C). Similarly, MALDI-TOF analysis of the purified Pisum sativum SS-tp indicated that the full-length peptide had the expected mass of 6182.2 (data not shown). Effects of Buffers, Storage Conditions, Enzyme Concentrations, and Metal IonsOne possible explanation for the variation in reported activities observed for the Toc GTPases is their sensitivity to various buffers, storage conditions, and metal ion contaminants. As such, it was important to test how different buffer systems affected the hydrolysis activity measurements. We tested four buffers for their effect on the GTPase activities of psToc34. A comparison of the enzymatic activity of psToc34 in these four buffers is shown in supplemental Fig. 1A. Despite the differences in their composition, all of these buffers yielded relatively similar activities of psToc34; however, GBS appeared to have slightly higher levels of GTPase activity after background subtraction.
Upon determining the best buffer for enzymatic activity, an additional variable is how the protein is stored prior to use; therefore, we also tested several freezing and storage conditions as follows: -80 °C slow, -80 °C fast, liquid nitrogen storage, and -20 °C as shown in supplemental Fig. 1A. We found that -80 °C slow freeze afforded the greatest maintenance of activity, and under optimal storage conditions, we were able to freeze and thaw the protein >3 times before activity became impaired (data not shown). Furthermore, GBS buffer was able to maintain the GTPase activity of the Toc proteins at -80 °C in 20% glycerol for up to 4 weeks, after which the activity began to diminish (data not shown).
The activity of psToc34 was assayed as a function of the pH in supplemental Fig. 1B. We varied several buffers (each buffering agent at 20 mM) from pH 6.0 to 12.0 and tested the maximal velocity as described above. Consistent with previous publications reporting the activity of psToc34 in GBS buffer, pH 7.6 was chosen (16). According to supplemental Fig. 1B, the pH optimum for psToc34 is between 7.5 and 8.2, and we therefore maintained the pH at 7.6 in accordance with the published GBS buffer.
Because it has been reported previously that Ras and other small GTPases are affected by the metal cofactor in the nucleotide-binding site, we attempted to remove Mg2+ from the active site of the Toc proteins by dialysis with EDTA in Mg2+-free GBS. This treatment consistently resulted in protein precipitation (data not shown). Alternatively, we were able to add various metal salts at a large molar excess to determine whether these metal ions could possibly exchange with the Mg2+ ion in the active site. The results of how MgCl2, CaCl2, CdCl2, CoCl2, CuSO4, MnCl2, NiSO4, and ZnCl2 affect activity is shown in supplemental Fig. 1C. Suppression of activity was seen with all metals but two, CaCl2 and ZnCl2, of which the ZnCl2 stimulated the hydrolysis activity nearly 2.7-fold. As a control, each metal salt was incubated with nucleotide (in the absence of enzyme) showing no effect on the spontaneous rates of GTP hydrolysis (data not shown).
|
-32P]GTP. This assay employs liquid scintillation counting to quantitate the amount of Pi released by GTP hydrolysis. In this assay, we incubated the
TM Toc proteins with [
-32P]GTP at 25 °C. Aliquots of the reaction were added to 10% (w/v) activated charcoal where only the hydrolysis product 32PO3 along with H2O, HCl, and H3PO4 remain in the filtrate after a 0.2-µm vacuum filtration step. Aliquots of this filtrate were counted in a scintillation counter. A bovine serum albumin control yielded counts analogous to spontaneous hydrolysis (data not shown). Counts/min data from at least three independent experiments were obtained and graphed as nanomoles of GTP hydrolyzed per min per µmol of Toc protein (Fig. 2, A-D). This analysis indicates that atToc33, atToc34, and psToc34 share very similar KM and Vmax values, whereas atToc159G does not. Values from 3.5 to 6.3 µM for the Michaelis-Menten constant, KM, and 35-37 nmol of GTP hydrolyzed/min/µmol Toc protein for the maximal velocity, or Vmax, were obtained for the smaller GTPases (Table 1). The G domain of atToc159 had a KM and a Vmax of 21.16 µM and 10.13 nmol GTP hydrolyzed/min/µmol Toc protein, respectively (Table 1). In the case of the Toc34 subfamily of proteins, all three displayed a sharp substrate-velocity curve, whereas atToc159G did not, representative of its relatively low maximal velocity and decreased affinity for substrate.
|
12.5 µM, suggesting that this is the equilibrium midpoint of psToc34. At concentrations below this threshold value, psToc34 is more likely to exist at a monomer, whereas higher concentrations increase the population of psToc34 dimers. Interestingly, the GTP hydrolysis rate of the four Toc GTPases, when corrected for protein concentration, do not increase in rate as a function of protein concentration (Fig. 3C). The dashed and dotted lines in Fig. 3C represent 1.7 and 16.9 µM, respectively, corresponding to concentrations where AUC analysis revealed psToc34 to be monomeric and mostly dimeric, respectively (Fig. 3D). At 33.8 µM psToc34, the protein is predominantly dimeric (Fig. 3D). Furthermore, the monomeric species decreases concomitantly with the increase of the dimeric species as protein concentration is increased as is shown graphically in Fig. 3E. Therefore, the maximal rate of the psToc34 is not affected by protein dimerization in vitro, but this does not preclude the possibility that dimerization in vivo may affect the hydrolytic rate and thus protein import.
|
Wild-type psToc34, as shown in Fig. 3, has a defined profile in the analytical ultracentrifuge indicative of a concentration-dependent monomer/dimer equilibrium. However, AUC analysis of the mutant R133A reveals a single sedimenting species with an S value of 2.8. Moreover, this single species is similar in sedimentation behavior to the wild-type monomer (Fig. 3D). However, unlike wild-type psToc34, the R133A mutant remains monomeric at all three protein concentrations tested: 3.4, 33.8, and 170 µM (Fig. 4A). Comparatively, dimerization of wild type was observed as low as 16.9 µM (Fig. 3D). This AUC analysis indicates that the R133A mutation severely disrupts the psToc34 homodimeric interface.
In addition to eliminating concentration-dependent dimerization, the mutation of R133A also abolishes the in vitro GTP hydrolytic activity (Fig. 4B). The low level of hydrolysis is essentially indistinguishable from the rates of spontaneous GTP hydrolysis, whereas wild-type psToc34 exhibits a sharp substrate-velocity curve, also shown for comparison in Fig. 4B. This very low activity precludes determination of the Vmax or KM values with any statistical certainty. These two observations taken together indicate that the R133A mutation disrupts both GTP hydrolysis in the monomer and prevents homodimerization.
To rule out the possibility that the R133A mutant resulted in misfolding or some global structural change, we performed circular dichroism spectroscopy of purified psToc34(R133A) shown in supplemental Fig. 2. Comparison of the data to the recently published CD spectra of wild-type psToc34 (10) reveals high spectral similarity. To provide more detailed insight into the structural properties of the mutant protein, we deconvoluted the CD spectra using the CDPro program. The secondary structural content predicted by this program is shown in supplemental Table 1. Furthermore, comparison of the data with the secondary structure content derived (using Stride, supplemental Table 1) from the crystal structure of wildtype psToc34 (21) is in close agreement with our measurements for R133A, indicative of a very similarly folded protein. Thus, the loss of activity of psToc34(R133A) is not because of protein misfolding.
|
|
Nucleotide Binding Activity of the Toc GTPasesWe used a novel approach to determine the equilibrium dissociation constants of the various Toc proteins for GTP, GDP, and other nucleotides. The particular Toc protein is incubated with a known concentration of radiolabeled [
-32P]GTP as well as increasing concentrations of unlabeled GTP in the case of homologous competition. Homologous and heterologous competition assays were carried out on the Toc GTPases using a variety of substrates to determine affinities and to gain insight into the significance of the interacting moieties, as shown in Fig. 6, A-F. All of the Toc GTPases in this investigation bound GTP but with varying degrees of affinity, ranging from 29 nM to 1.2 µM (Table 2). Although these enzymes have already been implicated as GTPases, no systematic enzymatic comparison has been made. Unlike the hydrolysis data, psToc34 and atToc33 both bind GTP quite strongly with Kd values of 50 and 29 nM, respectively, whereas atToc34 and atToc159G have considerably lower affinities. The difference in affinity of the psToc34 family for GTP and GDP was not great, showing only a 2-3-fold difference. In contrast, whereas atToc159G does not bind GTP as tightly as the small Tocs, it appears to be the most effective discriminator of GTP and GDP with nearly a 10-fold difference (Table 2).
|
-S and GMP-PNP (Fig. 6E). These results show that, although they fail to be hydrolyzable, the GTP-like properties of GTP-
-S and GMP-PNP both compete for binding; however, GTP-
-S with its lower Ki value suggests it is the better nonhydrolyzable GTP analogue for psToc34. In this assay, psToc34 did not have any significant affinity for the cGMP analogues, the purine ATP, nor the pyrimidine TTP (Fig. 6F). However, most interestingly, it was able to slightly bind the pyrimidine UTP, which shares an analogous placement of a carbonyl oxygen with the purine GTP, suggesting that this oxygen is critical for recognition by psToc34. This may also explain the ability of psToc34 to bind XTP. Homology Modeling of the Toc GTPasesatToc33, atToc34, and atToc159G were aligned with the protein sequence from the crystal structure of psToc34 (PDB code 1H65, chain B) using the Salign algorithm. Using this alignment, molecular models were made with the homology modeling software Modeler. The high degree of identity/similarity and the small number of inserts between the proteins were well suited to this homology modeling technique and did not require loop-modeling refinement. The overall fold of each of the proteins was very similar according to the models, and in fact, when superimposed they are indistinguishable from the majority of each structure as shown in Fig. 7A. Not including the inserts, the total average r.m.s.d. value for each of the modeled proteins was 0.128, 0.142, and 0.193 for atToc33, atToc34, and atToc159G, respectively. Fig. 7A is the structural alignment of atToc33, atToc34, and atToc159G models with the crystal structure of psToc34, and the color scale denotes per residue r.m.s.d. values. There are a few small insertions, mostly in atToc159G. Fig. 7B displays the per residue r.m.s.d. as a function of amino acid position of psToc34; at each gap the inserted residues are denoted below. The residues with large r.m.s.d. values are for the most part near insertions because Modeler satisfies spatial restraints based not only on the template but also creates bonds with reasonable chemical properties, and therefore residues near gaps are moved to allow a feasible model of the insert.
|
-phosphate of GTP is nearly identical with the placement of the conserved Arg residues that have been implicated as arginine fingers found in several GAPs that have been co-crystallized in complex with their cognate GTPases (Fig. 7E and references within the Fig. 7 legend). Although the primary sequence alignment of atToc159G would suggest that that Arg-953 may perform the equivalent function as an arginine finger upon hetero- or homodimerization, our model suggests that this is probably not the case. As shown in Fig. 7D, Arg-953 from atToc159G is displayed in cyan, and although still modeled to be within an
-helical segment with a similar backbone position, its side chain is clearly projecting away from the bound GDP, thus preventing the amino group from stabilizing the
-phosphate as has been postulated for other GAP arginine fingers.
|
35 nmol of GTP per min per µmol of Toc (Fig. 2B). Furthermore, when psToc34 is incubated with a 25:1 molar excess of reduced carboxymethylated lactalbumin (RCMLA), an inherently unstructured peptide, the rate of 32PO3 release remains constant and is very similar to psToc34 alone (Fig. 8A), and the velocity reflects that of "unstimulated" psToc34 (Fig. 8B). Furthermore, when psToc34 is incubated in the presence of a similar molar excess of either SStpPs or SStpNt (the transit peptide from the small subunit of ribulose-bisphosphate carboxylase/oxygenase from either pea or tobacco, respectively), the amount of 32PO3 release and the rate of GTP hydrolysis increases (Fig. 8, A and B). Fig. 8A represents a single concentration of GTP (25 µM) with SStpPs, SStpNt, and RCMLA all present at 25-fold molar excess to Toc protein. This was repeated with multiple concentrations of GTP, aliquots removed at 0, 10, 15, and 20 min, and scintillation counted yielding counts/min data that were fit to linear regression (example shown in Fig. 8A). Each linear regression yielded a slope that was mathematically transformed into a rate expressed in nanomoles of GTP hydrolyzed per min per µmol of psToc34, as seen in Fig. 8B (data shown represents psToc34 with one concentration of peptide, 25 µM, and various GTP concentrations). A substrate-velocity plot was created for each of 11 peptide concentrations ranging from 0.58- to 100-fold molar excess (data not shown). No concentration of RCMLA stimulated the velocity of psToc34 at any of the GTP concentrations assayed (Fig. 8B and data not shown). SStpPs and SStpNt both stimulated the maximal velocity as well as reduced the KM (Fig. 8, B and D). Using the nonstimulated velocity curve in Fig. 8B (solid squares, solid line), we determined the fold increase in rate at each of the 11 peptide concentrations of SStpPs, SStpNt, and the negative control RCMLA (Fig. 8C). Each data point for each treatment in Fig. 8C represents the Vmax value obtained from nonlinear regression for each treatment (culmination of all data points for a given treatment) divided by the Vmax value of psToc34 with no treatment (Fig. 8B). Both SStpPs and SStpNt increase the rate of GTP hydrolysis by psToc34 by
1.6-fold, whereas RCMLA does not affect the rate (Fig. 8C). Furthermore, the molar excess of SStpPs and SStpNt required to give half the maximal rate increase is 3.5 and 3.8 molar fold, respectively. SStpPs and SStpNt also alter the KM value or the affinity of psToc34 for the substrate GTP (Fig. 8D). Each data point in Fig. 8D represents the KM value calculated from the nonlinear fit of each molar excess of SStp to psToc34. The molar excess of SStpPs and SStpNt required to reduce the KM by half is 4.6 and 4.2, respectively.
|
-32P]GTP, and the dissociation of labeled nucleotide was monitored by scintillation counting as a loss of counts as a function of time under several conditions. Fig. 9, dashed line with closed squares, represents the basal rate of GDP loss with psToc34 alone in GBS. The remaining data indicate that there is no significant difference between the rate of exchange with a 100-fold excess of GTP and no SStpNt (open triangles) relative to a 100-fold excess of GTP at any amount of SStpNt, as indicated in Fig. 8 legend. | DISCUSSION |
|---|
|
|
|---|
TM forms of psToc34, atToc33, atToc34, and atToc159G. Although it would not be surprising for these rates to vary between the different homologues, the substantial variation in the reported kinetic properties between different groups, and on occasion even within a given group for the same protein, is troubling. For example, with psToc34, careful analysis of the published data reveals hydrolysis rates within the same article that vary from 0.588 to 26.5, when expressed as nanomoles of GTP/min/µmol Toc protein (21). Moreover, substantial variability is seen when comparing the hydrolytic rates (in nmol/min/µmol) of homologous proteins between pea and Arabidopsis, where values range from 0.588 for psToc34 (21) to
300 for atToc33 and atToc34 (16). Upon comparing the reported values for rates of hydrolysis in the literature, a lack of standardization becomes apparent. Also, there were experimental design problems, including reporting rates of GTP hydrolysis that are nonlinear with time (14, 26, 33), kinetic analyses performed with subsaturating levels of GTP (33, 34), and application of nonlinear fits to poor or limited data with no error analysis (35). To help reconcile these conflicting reports, we have concurrently determined the enzymatic properties of these four Toc GTPases in parallel, thus enabling a direct comparison of their GTPase properties.
|
30% change in activity; those reported for the Toc GTPases vary in some cases by 100-fold (16).
It has been shown for other small GTPases that require Mg2+ for their activity (30-34) that very few other cations can substitute for Mg2+ in the catalytic site. We attempted to replace Mg2+ in the active site of psToc34 with several different cations and found that only Zn2+ could support hydrolysis activity of psToc34 at or above the Mg2+ level. Although zinc has a smaller radius than magnesium (1.60 versus 1.34 Å) it may still be able to fit into the metal-binding site of psToc34, thereby stabilizing GTP. Indeed, other divalent metal cations have been shown to bind to other GTPases, for example the replacement of the Mg2+ by Mn2+ in p21 Ras results in a higher GTPase activity (55). This stimulation has been suggested to involve an increase of the pKa of the
-phosphate of GTP resulting in an increase of reaction rate at neutral pH. Schweins et al. (55) determined that this relationship between pKa of nucleotide and 1H abstraction potential indicated that the GTP itself acts as a general base in the reaction. This divalent metal cation-mediated rate increase was also shown to be the case with two other GTPases Rap1A and Ran. The lack of elevated GTP hydrolysis by psToc34 with Mn2+ is puzzling; however, the non-Ras-like G2 motif and its interactions with the
-phosphate of bound nucleotide may account for these observations.
As a result of our systematic comparison, we observed a number of differences between the enzymatic properties of atToc159G and the small Toc GTPases; atToc159G demonstrated a lower Vmax, a higher KM, a lower Ea, a higher kcat, and a relatively lower affinity for GDP and GTP. One possibility that may account for these enzymatic differences are the residues coordinating the Mg2+ ion. We took advantage of the molecular model that we built to explore this possibility. From the crystal structure of psToc34, the Mg2+ is coordinated between residues Glu-73 and Ser-53; however, from our homology-based model of atToc159G, this glutamate residue is replaced by a serine (Fig. 7D). The GDP and Mg2+ are shown with residues from psToc34 in pink and equivalently aligned residues (Ser-889 and Ser-869) of atToc159G in cyan (Fig. 7D). Whether these two serines are capable of coordinating the Mg2+ is not known; however, analogous small GTPases (e.g. Rho) have been reported to contain a threonine-threonine Mg2+ coordination site (56). This Ser-Mg2+-Ser scenario is not likely because the participation of dual serines has not been observed in any GTPases to date nor have two serines been reported in the architecture of >200 other Mg2+ coordination sites (57). Alternatively, our model suggests that another acidic residue, Asp-884, may be positioned such that it could provide an alternative Mg2+ coordination site for GTP/GDP binding to atToc159G (Fig. 7D). Altering the coordination site of the Mg2+ ion may thus change the conformation of the bound nucleotide, leading to the enzymatic differences observed between atToc159G and the other small Toc GTPases. In this study we have primarily utilized the homology models to provide some probable structural insight for our in vitro observations. As a future direction we intend to use these models to introduce mutations that modulate individual functions of Toc GTPases, e.g. dimerization, hydrolysis, nucleotide specificity, and/or peptide affinity. These constructs will allow further examination of the interdependence of these functions as well as indicate the legitimacy of these models.
Although the enzymatic characterization of psToc34 clearly shows that the rate of GTP hydrolysis is insensitive to enzyme concentration, analysis by AUC reveals that the oligomeric properties of psToc34 over a similar concentration range leads to an increased population of the dimer. The behavior of psToc34 during centrifugation is suggestive of a fast equilibrium between the monomer and dimer, where the equilibrium is shifted toward the dimeric form as the protein concentration is increased. Whereas dimerization as a function of protein concentration is not a novel idea (58, 59), the dimeric form of psToc34 is no more active than its monomeric form. Interestingly, when we characterized the oligomeric behavior of a mutation, R133A, we observe a highly soluble, exclusively monomeric species at all protein concentrations tested, including those that were shown to be primarily dimeric for the wild type protein. Furthermore, the R133A mutant was unable to hydrolyze GTP at any appreciable rate. In light of evidence that psToc34(R133A) is correctly folded, the structural basis for why this mutant is unable to undergo dimerization or perform GTP hydrolysis is unknown. Nonetheless, the causal relationship between loss of dimerization and GTP hydrolytic activity was not resolved by our current experimentation.
Although the precise relationship between oligomerization and GTPase activity is complex, one possible explanation of these results is that when psToc34(R133A) is expressed in E. coli, it is loaded with Mg2+ GTP during expression and folding. However, because this mutation is catalytically inactive, there is no opportunity to hydrolyze this GTP to GDP + Pi, thus trapping this protein in the GTP-bound form. This is unlike the wild type protein which, when isolated and crystallized, was found to be exclusively in the GDP-bound form (21); there is no GTP-bound crystal structure to date. Our data support the proposal that Arg-133 functions as the arginine finger that contributes to the stabilization of the dimeric form of psToc34 (21). This is supported by our molecular models, which consistently place a conserved Arg at this corresponding residue, as spanning across the dimeric interface functioning to stabilize the
-phosphate associated with the GDP in the trans subunit (Fig. 7).
Because this work has been limited to the soluble forms of the Toc GTPases, how their transmembrane domains may influence either their oligomeric behavior or their catalytic activity remains unclear. However, using the soluble form of psToc34, we have clearly demonstrated the ability of this domain to form homodimers, yet we failed to see any evidence of its ability form homotetramers. This observation is consistent with the recent crystal structure, which also reported a dimeric structure with a 2-fold axis of symmetry. The ability of psToc34 to oligomerize into only dimers is especially intriguing in light of the recent low resolution structure of the pea Toc translocon that contains four Toc34 subunits with a single Toc159 subunit (4). This may suggest that the native Toc complex contains a pair of co-existing Toc34 dimers. Although speculative, the homodimerization that we observe in vitro may also exist in vivo and may become disrupted by either precursor binding or, alternatively, through interaction with Toc159, to yield a transient heterodimer. Our future work will be to dissect the interplay by which both nucleotides and transit peptides modulate the oligomeric properties of the Toc GTPases.
The comparison of activation energies for GTP hydrolysis between the Toc proteins shows values of
10-15 kcal/mol (Fig. 5). This is similar to the activation energy of another small GTPase, EF2, of
20 kcal mol-1 (60). It was, however, interesting to observe that although atToc159G has a lower Vmax and a higher KM value, it appeared to have a significantly lower activation energy (10.1 versus
15). This would suggest that although atToc159G has a lower affinity for productive binding of GTP, upon binding, the orientation of the nucleotide was in some ways more favorable for hydrolysis, thus a lower energy transition state. The lower nucleotide affinity reflected in the KM could be the result of some steric constraints imparted by the additional N and C termini residues found in atToc159G. It is also possible that the addition of these flanking sequences may somehow contribute to its distinct kinetic and thermodynamic parameters. Unfortunately, we have no structural insight into these regions, because when building our molecular models, we necessarily truncated these sequences because they were not present in psToc34 and would have disrupted the structural alignment. How the acidic N-terminal and membrane-associated C-terminal residues influence the enzymatic properties as well as homo/heterodimerization of atToc159 will be difficult to ascertain because the full size protein is insoluble when expressed in a heterologous system.
Heterodimerization and reciprocal activation of GTP hydrolysis is not a novel idea; in fact, Sun et al. (21) has proposed homodimerization between psToc34 monomers and possible heterodimerization between atToc33 and atToc34. This has been well studied in another protein transport system that involves GTPases, the SRP54/Ffh from E. coli. In this system, it has been shown that heterodimerization elicits a reciprocal GAP activity from each monomer, thus coordinating GTP hydrolysis, completing the GTPase cycle, and thereby facilitating protein transport. Although the interactions between the Toc GTPases may mechanistically be different from the SRP54/Ffh GTPases, the possibility of homodimerization and/or heterodimerization (21, 35) may still modulate the activities of the Toc GTPases and influence the functionalities of these proteins in vivo.
Furthermore, in the crystal structure, psToc34 appears to be a homodimer with an arginine finger extending from one monomer into the other (21). Fig. 7E is a stereo view of several small GTPase·GAP complexes showing the proximity of the arginine finger to the
-phosphate of GDP; Arg-133 of psToc34 is colored according to element name, and other proteins (Ras, Rho, Rab, CDC42 and psToc34) are in solid colors; see Fig. 7 legend for details. These complexes were aligned about the GDP and Mg2+ molecules; psToc34 as well as all of these canonical small GTPases appear to bind GDP and Mg2+ in a very similar conformation. Both atToc33 and atToc34 have equivalent arginine residues that appear to coincide with the position of the Arg-133 of psToc34. The arginine finger of each GAP can be seen in the area near the
/
-phosphate, and psToc34 fits well with the terminal portion of the arginine side chain. Although many GAPs do in fact have a tyrosine adjacent to the arginine finger, the tyrosine of psToc34 is likely not serving the same function as implicated previously (21).
Unfortunately, investigations into the ability of preproteins and/or transit peptide-derived synthetic peptides to stimulate the GTP hydrolysis rates of the Toc homologues is another area of some conflict and confusion in the literature. For instance, one group reported that there is an
30-fold stimulation of hydrolysis by Toc34 when in the presence of preprotein (11); however, this same group subsequently reported no effect upon the addition of synthetic peptides comprising the N and C termini of the prSSU transit peptide (26). Furthermore, an N-terminal peptide from prSSU had the equivalent stimulatory activity as the full-length preprotein with atToc159G alone and with the reconstituted complex. However, in this more recent paper, they reported only
4-fold stimulation and not the previously reported
30-fold effect. Although the enzymatic properties of the native or reconstituted Toc complex may certainly be different from the individual components, our results (Fig. 3) clearly show that the oligomeric state of psToc34 and the other Toc homologues does not stimulate the hydrolysis rates in vitro. A further concern when comparing synthetic peptides to the full-length preproteins is that E. coli expressed preproteins are insoluble and require chemical denaturation for solubilization. This may interfere with the GTPase assay because of indirect effects of the denaturant. Additionally, the use of short synthetic peptides could have unintentionally split a continuous recognition element present in the full-length transit peptide. Finally, short synthetic peptides would abrogate any structural contributions from the mature domain as has been clearly demonstrated for prSSU (61).
Nevertheless, it is clear from our results using a soluble, fulllength transit peptide with psToc34 that the transit peptide alone does in fact stimulate the rate of GTP hydrolysis. We further show that this stimulation increases not only the Vmax of the enzyme but also lowers the KM value for GTP. The
2-fold stimulation that we observe is very reproducible and is not the result of nonspecific interactions. Although this is much more activity than the synthetic peptides were capable of stimulating psToc34 as reported by Becker et al. (26), it is still significantly less than the initial activity reported value by prSSU (11). Moreover, we have also tested the full-length prSSU and found that it does demonstrate a higher level of stimulation yet plateaus at
3-fold (data not shown). Causes for the variability between laboratories in the level of peptide and preprotein stimulation of the hydrolysis of GTP by the Toc GTPases remain unclear. All the same, the ability of the full-length transit peptide to stimulate the rate of GTP hydrolysis suggests that this sequence may somehow regulate the activity/structure of the GTPase during Toc-mediated precursor recognition. The fact that this stimulation was not observed upon addition of the unstructured protein, RCMLA, suggests that the GTPase stimulation is a result of specific interaction(s) being mediated by one or more regions of the transit peptide and is not simply the interaction of the exposed peptide backbone as has been suggested for RCMLA stimulation of ATPase activity of the chaperone hsp70 (39). These observations suggest that some region or domain of the transit peptide somehow functions to regulate the activity/structure of the GTPase during Toc-mediated precursor recognition.
Transit peptide-mediated stimulation of the GTPase activity of Toc34 and Toc159 has been reported (10, 15, 26, 31). Although others have even used the term GTPase-activating protein, we feel that this was premature. The regulatory mechanism of switch GTPases involves both activating proteins and exchange factors. If the only observable quantity is a hydrolysis product, then both exchange factors and activating proteins will lead to an apparent increase in rate. We have carefully examined the stimulatory effect on hydrolysis by titrating the transit peptide and saturating that effect as seen in Fig. 8. In a separate assay designed to follow the rate of exchange, no significant differences are observed with GTP alone and the addition of nearly a 40-fold excess of transit peptide (Fig. 9), which was well into the saturating concentration according to Fig. 8C. Therefore, taken together, these data suggest that the transit peptide is in fact functioning as a GAP and not a guanine nucleotide exchange factor or GEF. The ability of two different preproteins (prSSU and prOE33) has already been shown to stimulate the rate of GTPase hydrolysis of the purified translocation complex (11), suggesting that this may be a universal feature of chloroplast precursors and transit peptides.
The evolutionary process of how transit peptides acquired the ability to function as a GTPase-activating protein is unknown. Moreover, there are no structural details known on how transit peptides stimulate the GTPase activity. It has become clear that for many of the well characterized, small GTPases (Rho, Rac1, and Cdc42), there are several divergent classes of proteins that are capable of stimulating the hydrolysis of GTP (62). The Rho family of small GTPases has been found to be the target of several bacterial cytotoxins that transiently down-regulate these GTPases. This can be done by both covalent modifications such as ADP-ribosylation and glycosylation (63), as well as through cytotoxins or virulence factors, such as YopE, ExoS, and SptP, that mimic the activity of eukaryotic cytosolic GAPs (64). Interestingly, although these proteins have very little sequence similarity with each other, they have been shown to have very similar three-dimensional structures (65). Crystal structures of eukaryotic GAPs with their cognate GTPases show striking similarity in observed function with the bacterial cytotoxins, especially ExoS from Pseudomonas aeruginosa with Rac (66). Of significant note is the critical Arg residue (Arg-144) that projects outward interacting with the bound GTP molecule of the GTPase. Comparison of this model with the structure of other eukaryotic GAPs is remarkable because outside of the conserved Arg residue, there are no other structural similarities of these bacterial cytotoxins to the cognate eukaryotic GAPs. This would suggest that the bacterial GAPs arose via convergent evolution with the functional analogues found in eukaryotes. This observation also suggests that multiple GAPs with highly divergent structures can perform the equivalent functional role. This common activity in the absence of a conserved structure may also be observed in how transit peptides function as GAPs.
Interestingly, the GAP activity associated with the bacterial cytotoxins is the disruption of the actin cytoskeleton in the host cell thereby avoiding engulfment and destruction by phagocytosis (67). To accomplish this activity they rely on a type III secretory pathway to inject these virulence factors into the host cell. It has been proposed earlier that transit peptides could have evolved from cyanobacterially expressed virulence factors (1, 68). Possibly these ancestral transit peptide progenitors were secreted into the host cytoplasm via the SynToc75 complex. Upon secretion, they function as virulence factors able to alter the actin cytoskeleton of their endocytotic host via the GAP mechanism, similar to YopE/ExoS/SptP. Although speculative, it is possibly through this GAP activity of the secreted peptides that these cyanobacteria were able to thwart the host defense response that would normally lead to their degradation. This GAP activity thus enables the transition from endocytosed bacteria into a semi-autonomous organelle. Finally, through unknown evolutionary changes, these ancestral host-encoded soluble GTPases were ultimately recruited and redeployed to the chloroplast envelope where they became receptor/gatekeeper GTPases. This ancestral host origin of the Toc GTPases is consistent with the fact that to date there are no clear cyanobacterially encoded Toc34 or Toc159 homologues.
Clearly, the mechanism of how transit peptides function as GAPs is an important aspect to their role in chloroplast protein import. However, in the absence of any structural information on how transit peptides are recognized by the Toc receptors, the mechanism for how they stimulate GTP hydrolysis will be speculative. For example, although SStp has two relatively conserved Arg residues, whether they function as Arg-finger(s) is not known and will require more investigation. The significant enzymatic differences that we observe between the small and acidic Toc GTPases suggest that these two subfamilies may inherently function differently in chloroplast protein import. Whether the other acidic Toc GTPases (atToc90, atToc120, and atToc132) function similarly to atToc159 and/or display some form of cooperative or synergistic behavior in complexes will be interesting. Finally, our future work will help clarify how the dual substrates (GTP and transit peptides) coordinately modulate the organization and enzymatic properties of the Toc GTPases in the translocon.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Figs. 1 and 2 and supplemental Table I. ![]()
1 Both authors contributed equally to this manuscript. ![]()
2 Supported in part by a Science Alliance Award from University of Tennessee, Knoxville. ![]()
3 University of Tennessee, Knoxville undergraduate supported in part by Hope, Volunteer, and Merit Scholarships. ![]()
4 To whom correspondence should be addressed: Biochemistry, Cellular and Molecular Biology Dept., University of Tennessee, Knoxville, TN. Tel.: 865-974-4082; Fax: 865-974-6306; E-mail: bbruce{at}utk.edu.
5 The abbreviations used are: Toc, translocator of the outer membrane of the chloroplast; TIC, translocator of the inner membrane of the chloroplast; TOM, translocator of the outer membrane of the mitochondria; AUC, analytical ultracentrifugation; XTP, xanthosine 5'-triphosphate; GTP-
-S, guanosine 5-O-(3-thiotriphosphate); GMP-PNP, guanylyl-5-[
,
-imido]-triphosphate; MES, 4-morpholineethanesulfonic acid; Tricine, N-[2-hydroxy-1,1-bis(hydroxymethyl)ethyl]glycine; CAPS, 3-(cyclohexylamino)-propanesulfonic acid; NTA, nitrilotriacetic acid; r.m.s.d., root mean square deviation; MALDI-TOF MS, matrix-assisted laser desorption ionization time-of-flight mass spectrometry; PDB, Protein Data Bank; IMAC, immobilized metal affinity chromatography; RCMLA, reduced carboxymethylated lactalbumin; GAP, GTPase-activating protein. ![]()
| ACKNOWLEDGMENTS |
|---|
-32P]GTP and [
-32P]GTP. We also thank Dr. Steve Wright for helpful discussions. | REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Rahim, S. Bischof, F. Kessler, and B. Agne In vivo interaction between atToc33 and atToc159 GTP-binding domains demonstrated in a plant split-ubiquitin system J. Exp. Bot., January 1, 2009; 60(1): 257 - 267. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Bionda, P. Koenig, M. Oreb, I. Tews, and E. Schleiff pH Sensitivity of the GTPase Toc33 as a Regulatory Circuit for Protein Translocation into Chloroplasts Plant Cell Physiol., December 1, 2008; 49(12): 1917 - 1921. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Koenig, M. Oreb, K. Rippe, C. Muhle-Goll, I. Sinning, E. Schleiff, and I. Tews On the Significance of Toc-GTPase Homodimers J. Biol. Chem., August 22, 2008; 283(34): 23104 - 23112. [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 |