Importin-β (Impβ) is a major transport receptor for Ran-dependent import of nuclear cargo. Impβ can bind cargo directly or through an adaptor such as Importin-α (Impα). Factors involved in nuclear transport have been well studied, but systems analysis can offer further insight into regulatory mechanisms. We used computer simulation and real-time assays in intact cells to examine Impα–β-mediated import. The model reflects experimentally determined rates for cargo import and correctly predicts that import is limited principally by Impα and Ran, but is also sensitive to NTF2. The model predicts that CAS is not limiting for the initial rate of cargo import and, surprisingly, that increased concentrations of Impβ and the exchange factor, RCC1, actually inhibit rather than stimulate import. These unexpected predictions were all validated experimentally. The model revealed that inhibition by RCC1 is caused by sequestration of nuclear Ran. Inhibition by Impβ results from depletion nuclear RanGTP, and, in support of this mechanism, expression of mRFP-Ran reversed the inhibition.
The transport of proteins in and out of the nucleus is a major cellular process (Görlich and Kutay, 1999; Macara, 2001; Weis, 2002). Transport not only localizes proteins destined for the nucleus or cytoplasm, but also functions as a key step in signal transduction pathways and in the regulation of cell cycle progression. The nuclear pore complex is a large multi-protein assembly that penetrates the nuclear envelope and connects the cytoplasm to the nucleus (Fahrenkrog and Aebi, 2003; Suntharalingam and Wente, 2003). Although small proteins <20–40 kD can move through the NPC by passive diffusion, larger proteins require soluble transport receptors called karyopherins (Chook and Blobel, 2001). Importin-β (Impβ) is one of the predominant karyopherins that drive import. Although a small number of cargo proteins may bind Impβ directly, most cargoes require the adaptor protein importin-α (Impα).
Impα contains a COOH-terminal region that binds directly to proteins containing an NLS sequence (Conti et al., 1998; Herold et al., 1998). However, Impα also contains an NH2-terminal auto-inhibitory domain that blocks the NLS-binding site (Kobe, 1999). Binding of Impβ to Impα relieves this auto-inhibitory blockade and allows Impα to bind cargo proteins with high affinity (Fanara et al., 2000).
Translocation of the Impα–Impβ–cargo complex through the NPC is thought to be mediated by weak hydrophobic interactions between Impβ and nucleoporins (Bayliss et al., 2000). On the nuclear side of the NPC, RanGTP dislodges Impβ from the complex (Görlich et al., 1996). The Impα–cargo heterodimer then dissociates, assisted by the high affinity binding of CAS to NLS-free Impα. CAS is a karyopherin that serves to export Impα in association with RanGTP (Kutay et al., 1997). Both the RanGTP–Impβ and RanGTP–CAS–Impα complexes translocate back to the cytoplasm where RanGAP (assisted by RanBP1 and Impα) hydrolyzes the RanGTP to RanGDP (Bischoff and Görlich, 1997; Floer et al., 1997; Petersen et al., 2000). The export complexes disassociate and the transport receptors are recycled for another round of import (Fig. 1).
The gradient of RanGTP between the nucleus and cytoplasm provides the motive force to drive transport. This gradient is created by an asymmetric distribution of the nucleotide exchange factor RCC1 (also called RanGEF) and of RanGAP. RCC1 is a resident nuclear protein that promotes the exchange of RanGDP to RanGTP (Bischoff and Ponstingl, 1991a). Excluded from the nucleus, RanGAP acts to maintain Ran in the GDP-bound state in the cytoplasm (Fig. 1). The Ran gradient is highly dynamic, with RanGDP rapidly imported by transport factor NTF2 (Ribbeck et al., 1998; Smith et al., 1998).
Smith et al. (2002) conducted the first in silico analysis of Ran transport. A system of ordinary differential equations (ODEs) was created using the Virtual Cell software to describe reactions and fluxes of proteins involved in Ran Transport (Loew and Schaff, 2001; Smith et al., 2002). These proteins were set at known initial concentrations and allowed to come to steady state. A jump in cytoplasmic Ran concentration was used to simulate the micro-injection of recombinant Ran. Sensitivity analysis showed that RCC1 dominated changes in steady-state Ran flux, whereas NTF2 and the permeability of NTF2–RanGDP had smaller effects. Increases in RanBP1, RanGAP, karyopherins, and the permeability of the RanGTP–karyopherin complex had no effect on Ran flux at steady state. A similar analysis to determine rate-limiting components of the initial rate of Ran import after cytoplasmic injection showed that NTF2 was limiting, but that the initial rate of import was insensitive to RCC1.
To compare simulation results with Ran flux in live cells, a fluorescently tagged Ran (RanFl) was injected into BHK21 cells and imaged using confocal microscopy. Time-lapse images were used to quantify nuclear accumulation and import rate of the Ran. Co-injected NTF2 produced positive changes in the initial rate of import and steady-state N/C ratio of Ran, as predicted by the model. RanGAP produced no changes, which also is consistent with the model. RCC1 depletion was tested using tsBN2 cells. These cells possess a temperature-sensitive RCC1 mutation, so that at 39.5°C the mutant RCC1 is rapidly degraded (Uchida et al., 1990). Temperature-shifted tsBN2 cells showed reduced nuclear accumulation of Ran but only small affects on initial rate of Ran import, which is consistent with the model.
Görlich and colleagues (Görlich et al., 2003) developed a model to study how the nuclear–cytoplasmic Ran gradient couples to nuclear import through Impβ. The authors modified several kinetic parameters in their simulation of the Ran GTPase system. Most significantly, the Michaelis-Menten model of RCC1 used by Smith et al. (2002) was replaced by a four-step reversible series of reactions for GTP–GDP exchange, based on a detailed in vitro kinetic analysis by Klebe et al. (1995b). The kcat/Km for RanGAP was also adjusted upward to reflect enzyme activity at 37°C. This simulation predicted that decreases in nuclear RCC1 would reduce RanGTP accumulation, but that RCC1 is not limiting for the system. Conversely, the simulation suggested that RanGAP is limiting. Lastly, nuclear accumulation of RanGTP was shown to be heavily dependent on the GDP/GTP ratio used in the model. The model also predicted a close relationship between the RanGTP gradient and the nuclear/cytoplasmic ratio of cargo concentration at steady state.
In summary, the Smith and Görlich kinetic models show how the RanGTP gradient is established and used to drive a nuclear cargo gradient, and that even a simple system of interacting proteins can generate behavior that would difficult to predict using a static model.
In the current study, we extend the RanGTPase model to show how the Ran gradient couples to nuclear cargo import through the Impα–Impβ system. Using a compartmental model, we show how the system of known biochemical reactions involved in Impα–Impβ import can successfully recapitulate nuclear import seen in vivo. Predictions from the model were tested in vivo by monitoring nuclear import in whole cells. Unexpectedly, the model predicted that increasing the concentrations of RCC1 or Impβ would inhibit rather than stimulate the initial rate of cargo import, and these effects were validated experimentally.
We constructed a compartmental model of Impα–Impβ-mediated nuclear transport using the program Jarnac, a biochemical simulation package for Windows (Sauro et al., 2003). An overview of the model is shown in Fig. 1 and is described in Materials and methods. A more detailed schematic is provided in the Fig. S1 , together with a script file of the model than can be run in Jarnac (http://220.127.116.11/downloads/ISetup.exe).
Quantitative models of nuclear protein import are necessarily complex, because of the large number of factors involved in a single transport cycle. Moreover, there are multiple pathways for nuclear import and export, all of which use some factors such as Ran and NTF2 in common, but which use different karyopherins (Görlich and Kutay, 1999; Chook and Blobel, 2001; Macara, 2001). The concentrations and even the identities of all the cargoes for each of these karyopherins remain mostly unknown. However, the kinetics of the Ran cycle and of individual components of the Impα–β pathway have been well studied. Therefore, in our computer model we used experimentally determined parameters for Ran and Impα–β, and included terms for a generic karyopherin to account for other nuclear transport pathways. The same kinetic parameters were used for the generic karyopherin as for Impβ. This approach seems reasonable because karyopherins such as transportin bind Ran with an affinity similar to that of Impβ, while karyopherins that have a much lower affinity for Ran will not significantly perturb the system. Another complication is that there are probably hundreds of endogenous cargoes for Impα–β, most of which have not yet been characterized, and their cytoplasmic concentrations are unknown. Necessarily, therefore, we initially treated the concentration of endogenous cargo as a floating parameter. Endogenous cargo and labeled cargo reactions were modeled with identical but separate reaction pathways (Fig. S1) using the parameters listed in Tables I–III.
The computer model was first run in the absence of any perturbation for a period of 1,000 s, to bring the system to steady state. A RanGTP gradient becomes established, and the transport factors distribute between the cytoplasmic and nuclear compartments (Table IV).
To simulate the micro-injection of labeled cargo, the cytoplasmic concentration of labeled cargo was stepped up from zero to some chosen initial value, and the simulation was run for an additional 500 s to track the import of this cargo into the nuclear compartment. Initial rates of cargo import were calculated from the rates of nuclear accumulation during the first 60 s after the addition of labeled cargo. To compare the results of the model with import in vivo, we injected GST-GFP-NLS (GGNLS) into the cytoplasm of whole cells and imaged nuclear import by confocal microscopy (Fig. 2 A). Both the initial rate and steady-state concentration of the virtual cargo compare well with the nuclear transport of GGNLS observed in live cells (Fig. 2 B).
A key feature in any model of nucleo-cytoplasmic transport is the permeability of the karyopherins, cargo, and other components across the nuclear envelope. Rather than model the detailed structure of individual pores, one can take advantage of the fact that translocation between the nuclear and cytoplasmic compartments behaves to a first approximation like passive diffusion, and can therefore be modeled using a simple permeability term (Table III). The permeabilities for some components, such as NTF2, have been determined previously, using reconstituted transport in digitonin-treated cells (Ribbeck and Görlich, 2001). We have also measured permeabilities both using digitonin-treated cells and micro-injection into intact cells to measure permeabilities for Impα, Impβ, and GST-GFP through the NPC (Fig. S2). Impβ, labeled with Oregon green isothiocyanate (OGI), gave a permeability of 0.4 s−1 in whole cells, compared with 0.95 s−1 reported in a previous study in digitonin-permeabilized cells (Görlich et al., 2003). The leakage of free RanGTP from the nucleus to the cytoplasm was represented using a permeability of 0.03 s−1, taken from the previously measured permeability for RanGDP (Ribbeck and Görlich, 2001).
One important limitation of the model is that it does not include terms for competition between different transport pathways at the level of permeability across the nuclear envelope. Thus, permeabilities measured in vitro (using digitonin-treated cells) are likely to be rather higher than permeabilities measured for the same proteins in intact cells where the proteins are competing with other endogenous proteins for access to the pores. However, we know that the NPC has a high capacity and the import rate is, to a first approximation, linear with concentration over the concentration range of GGNLS used in this study. Competition is unlikely, therefore, to be a major source of error in the modeling. Moreover, the sensitivity analysis (see below) suggests that the cargo import rate is not limited by the permeabilities of many of the transport components in the model.
To identify possible rate-limiting steps in nuclear protein import, we performed a global sensitivity analysis on reactant concentrations in the model (Fig. 3). Reactant concentrations were varied individually over a five log scale from 0.0001 to 10 times the initial concentration. After allowing the system to come to steady state for 1,000 s, the cytoplasmic injection of GGNLS was simulated. Normalized initial rates of GGNLS import were then plotted as a function of initial concentration.
These simulations showed that GGNLS import responds most dynamically to the Impα and Ran concentrations, and is limited by the estimated cellular concentration of these factors. The slope of the curve for Impα is higher than that for Ran, suggesting that transport is more highly responsive to small deviations from the basal concentration of Impα than of Ran. Surprisingly, however, import is predicted to be relatively insensitive to the cellular concentration of CAS, which recycles Impα back to the cytoplasm.
The initial import rate is also highly sensitive to the concentration of NTF2, which mediates the import of Ran into the nucleus, although the dynamic capacity is not so large as for Impα or Ran. Perturbations in the concentration of NTF2 will alter the import rate because the RanGTP gradient across the nuclear envelope is coupled to the rate at which RanGDP is recycled back into the nucleus.
The sensitivity analysis showed that RanGAP is not limiting for cargo import, and that the concentration of its cofactor, RanBP1, can be reduced at least 10-fold before cargo import is affected. Cytoplasmic RanBP1 forms a RanBP1–Impβ–RanGTP complex but then is immediately released in the next reaction step when RanGAP triggers GTP hydrolysis, and the complex disassociates (Fig. 1). This rapid recycling of RanBP1 allows it to be effective at submicromolar levels. By contrast, the nucleotide exchange factor for Ran, RCC1, shows an unexpectedly complex relationship to cargo import, in that either a reduction or an increase in the cellular RCC1 concentration relative to the endogenous level strongly inhibits import. A reduction in RCC1 reduces the concentration of free nuclear RanGTP, which leads to a nuclear accumulation of Impβ–Impα–cargo complexes that cannot dissociate, and of Impα which cannot be exported back to the cytoplasm by CAS. But why would an increase in RCC1 also inhibit cargo? The explanation lies in the four-step kinetic model for nucleotide exchange, described by Wittinghofer and colleagues (Klebe et al., 1995b). In this model, RCC1 has a significant affinity for Ran, and high concentrations of the exchange factor sequester Ran in the nucleus as an RCC1–Ran complex. The free concentration of RanGTP available to dissociate import complexes and export Impα is therefore reduced, with a consequent inhibition of cargo import.
Another unexpected outcome of the sensitivity analysis was the inhibitory effect of high Impβ concentrations. Intuitively, one might have predicted that increasing the amount of the transport receptor for GGNLS cargo would stimulate import. However, the receptor and Ran levels are finely balanced in the model and, if in excess, the Impβ translocates into the nucleus in an unloaded state, then sequesters free RanGTP, and delivers it to the cytoplasm where it is hydrolyzed. This process is called futile cycling, and it depletes the nucleus of RanGTP, thereby inhibiting cargo import. A similar, but less dramatic effect is seen at high concentrations of CAS, which also reduce import (Fig. 3).
We also performed a sensitivity analysis of the NPC permeability constants for various transport components (Fig. S3), which predicts that the system is fairly robust to perturbations in many of these constants. Thus, the predictions of the computer simulations are not dependent on their exact magnitudes. This is important, because we do not have experimental data on the permeabilities of all the components of the model.
Experimental tests of this type of global sensitivity analysis are complicated by the need to change the endogenous concentration of each factor by a known amount, before the injection of the fluorescent import cargo. Ectopic expression from a plasmid varies from cell to cell, and cannot easily be measured until after the experiment has been completed. Micro-injection is also highly variable from cell to cell, because the injected volume is very sensitive to back pressure and needle diameter. Therefore, we chose in most cases to use co-injection of the GGNLS cargo with known concentrations of recombinant factors. Although the amount of injected cargo and factor will vary from cell to cell, the actual concentrations can be calculated from the total GFP fluorescence intensity within each cell. One advantage of computational modeling is the ease with which effects on transport rate can be calculated and compared with experimental data even when there are two variables for each data point (the cargo concentration and the factor concentration).
Impα is a limiting reactant for nuclear import
To test the predictions of the computer model for nuclear import, we performed co-injection experiments using 20 μM GGNLS as the labeled cargo. The GGNLS was injected together with recombinant transport factors into the cytoplasm of HeLa cells that were maintained at 37°C. The total amount of injected cargo and its initial rate of import were then quantified as described in Materials and methods. The variability in the injected volume from cell to cell generates a range of initial cytoplasmic concentrations that can be plotted against the initial rates. To provide a direct comparison with the model, we generated a series of simulations that directly mirrored the conditions used in the co-injection experiments, over the same range of initial cargo and factor concentrations.
The initial cytoplasmic concentration of cargo varied from ∼0.5 to 4 μM, and the initial import rate was linear with concentration over this range (Fig. 4 A), showing that, as predicted by our model, the Impα–β transport system is operating well below maximum capacity in HeLa cells. The concentration at half-maximal import rate (K0.5) for cargo in silico was calculated to be ∼12 μM.
We first tested the effect of co-injecting recombinant Impα. The model predicted that this factor is rate-limiting for cargo import (Fig. 4 B). To ensure that the injected factor was functional, we used a permeabilized cell assay. A portion of the protein was fluorescently labeled and incubated with Impβ, Ran, NTF2, RanBP1, and RanGAP. The label rapidly accumulated in the nuclei of permeabilized cells (Fig. 4 C). The Impα was also able to bind to GST-Impβ attached to glutathione–Sepharose beads (Fig. 4 D). An unlabeled portion of the same Impα preparation was then co-injected into intact HeLa cells with the GGNLS cargo. As shown in Fig. 4 A, the Impα induced a significant increase in the initial rate of cargo import, similar in magnitude to that predicted by the computer model.
Import is insensitive to elevated CAS concentration
During each transport cycle, Impα needs to be recycled back to the cytoplasm. It is exported in a complex with RanGTP and CAS. Therefore, we expected that increasing the concentration of CAS in the cell would increase the availability of cytoplasmic Impα and therefore increase the cargo import rate. However, the model predicts that CAS is not limiting for import (Fig. 3 and Fig. 5 B). We addressed this possibility using recombinant CAS. To test the functionality of the CAS, we used a permeabilized cell assay, and observed that the recombinant CAS could accumulate in the nuclei of HeLa cells (Fig. 5 C). Second, we conducted binding assays between CAS, Impα, and RanGTP, using a constitutively active mutant of Ran (RanQ69L) (Fig. 5 D). The recombinant CAS bound Impα and showed increased binding in the presence of RanGTP. However, when the CAS was co-injected with GGNLS into intact HeLa cells, it had no detectable effect on import. These data show that CAS is not limiting for the initial rate of α–β cargo import into the nucleus.
Factors in the RanGTPase system are limiting for nuclear import
Cargo import via Impα–β depends on the nucleo-cytoplasmic gradient not only for the recycling of Impα, but also for the release of cargo from the Impα–β complex. We have shown previously that NTF2 is limiting for Ran import, and we therefore expected that both NTF2 and Ran might be limiting for NLS cargo import. Indeed, the computer model supports this idea (Figs. 3 and 6). To test these predictions we co-injected GGNLS with either wild type His6-tagged Ran or with NTF2. In both cases, a significant increase in the initial rate of cargo import was observed (Fig. 6, A and C). The effect of Ran co-injection was slightly larger, experimentally, than predicted from the model. This difference arises in part because we make the assumption (because of scatter in the data) that the initial rate varies linearly with GGNLS concentration, over the experimental range used, whereas in the model the rate curve is hyperbolic. Additionally, secondary effects, such as competition for access to the nuclear pores, were not included in the model.
We next tested other components of the Ran system that are required for maintenance of a RanGTP gradient (Fig. 7). The Ran gradient has been predicted to be sensitive to the concentration of extra-nuclear RanGAP (Ribbeck and Görlich, 2001). However, our cargo import model shows that RanGAP is not limiting for cargo import in intact cells at 37°C (Fig. 3 and Fig. 7 B). Experimental assays confirmed this (Fig. 7 A): although co-injection of active RanGAP appeared to slightly increase import, this effect falls below the level of statistical significance (P > 0.05).
Next, we tested the effect of the RanGAP coactivator, RanBP1, on the system. The computer model predicted that RanBP1 is not limiting, but the co-injection of recombinant RanBP1 reproducibly increased the initial rate of cargo import (Fig. 7, C and D). Why is this? The global sensitivity analysis performed on the model shows that RanBP1 acts catalytically to stimulate RanGAP, and is needed at only substoichiometric levels to maintain an optimal transport rate. Large increases in RanBP1 concentration in the model have little effect on import compared with the baseline level (Fig. 3). The unexpected disparity between the simulation and experimental results, together with the observation that elevated RanGAP does not significantly increase import, suggests that the role of RanBP1 may not be completely understood.
Increased nuclear RCC1 inhibits cargo import
One unexpected prediction of the computer simulations was that the endogenous RCC1 concentration is finely balanced to provide the maximum RanGTP gradient. A reduction in the amount of RCC1 in the nucleus will of course reduce the production of RanGTP, but an increase in the amount of RCC1 also is predicted to reduce the free RanGTP concentration, through the formation of an RCC1–RanGTP complex.
A small decrease in nuclear RanGTP in response to elevated RCC1 was proposed previously (Görlich et al., 2003), but this decrease was not predicted to alter the dynamic capacity of the system for Ran-driven transport, and no analysis of transport rates was performed. We found that cytoplasmic co-injection of recombinant RCC1 with GGNLS in equimolar ratio did inhibit import (unpublished data) but this result could arise because the RCC1 generates RanGTP in the cytoplasm, which would disassociate import complexes, and it might compete with the GGNLS for binding to Impα (Nemergut and Macara, 2000). Therefore, we tested the effect of RCC1 by transfection of HeLa cells with an mRFP-RCC1 construct, and then compared the import rates for injected GGNLS in cells that expressed the RFP fusion protein versus those that did not. As shown in Fig. 8 A, GGNLS import was dramatically inhibited in those cells that express mRFP-RCC1, as compared with adjacent cells that were injected 10–20 s earlier and that do not express the fusion protein. Cells transfected with mRFP alone showed no inhibition of import (Fig. 8 B), confirming that the effect is dependent on RCC1 expression and is not a nonspecific response to the mRFP.
RCC1 possesses an NH2-terminal NLS that can bind Impα. It is conceivable, therefore, that elevated RCC1 levels could inhibit cargo import by sequestering Impα in the nucleus. This explanation is perhaps unlikely, because free Impα has a very low affinity for cargo (Fanara et al., 2000), and the total nuclear concentration of NLSs is probably high. Nonetheless, we tested this possibility by staining RFP-RCC1 transfected cells for endogenous Ran. Sequestration of Impα in the nucleus would not be expected to alter nuclear Ran, whereas sequestration of Ran by the RCC1 would be predicted to cause an increase in the nuclear/cytoplasmic ratio of Ran. This outcome is seen in Fig. 8 C.
Excess Impβ suppresses the initial rate of nuclear import
Counter-intuitively, the model predicts that the co-injection of Impβ with cargo suppresses rather than stimulates the initial rate of GGNLS import (Fig. 3 and Fig. 9 B). Why does this surprising effect occur? The simulation shows that excess Impβ undergoes futile cycling—that is, it translocates into the nucleus without cargo, where it sequesters free RanGTP. Without an available pool of free RanGTP, the Impα–Impβ–cargo heteromer cannot dissociate, so the cargo is not released, and will not accumulate in the nucleus.
Experimental observations support the prediction of the model. As shown in Fig. 9 A, co-injection of recombinant Impβ suppressed the initial rate of GGNLS import. However, fragments of Impβ that lack the Ran-binding domain can inhibit transport by binding with high affinity to the nuclear pore complex. Thus, proteolysis of the Impβ might produce inhibitory contaminants that inadvertently produce the result predicted by the model. To test the functionality of the recombinant Impβ, we performed both in vitro and in vivo tests. We first checked that the majority of the purified Impβ could bind GST-Impα and GST-Ran(Q69L) (Fig. 9 C). Second, we injected fluorescently labeled Impβ into the cytoplasm of live cells (Fig. 9 D). The rapid nuclear accumulation of Impβ in whole cells supports the conclusion that the Impβ is functional and should not inhibit nuclear transport by binding irreversibly to the pores.
The underlying mechanism for the inhibition of cargo import by increased Impβ is the sequestration of nuclear RanGTP, and the model predicts that addition of exogenous Ran can reverse this effect (Fig. 10 A). The additional Ran accumulates in the nucleus in the GTP-bound form, replenishing the free RanGTP pool that has been reduced by the excess Impβ, and hence stimulating import (Fig. 10 B). To test this prediction, we first transfected HeLa cells with a fusion protein of Ran and RFP (mRFP-Ran). We then co-injected GGNLS and Impβ into those cells that showed nuclear accumulation of mRFP-Ran. Compared with nontransfected cells from the same population, the mRFP-Ran cells showed significant reduction of the Impβ effect on initial rate (Fig. 10 B). This result strongly argues that the observed effect of Impβ on import rates is not a consequence simply of contamination by inhibitory fragments of Impβ, and also validates the computer model.
Computer models are relevant to cell biology only to the extent that they can generate testable predictions, and provide new insights into cellular processes. They can suggest unexpected relationships between steps in a process, or expose flaws in our understanding by failing to correctly simulate an experimental outcome. We believe that our model of nuclear protein import satisfies both of these requirements. The model nicely describes the general features of Impα–Impβ-mediated import, in terms of rates and the requirement for a Ran-gradient but, more importantly, it made several unanticipated predictions about the response of the system to changes in the concentrations of certain factors. Moreover, it fails to account for certain experimental observations, which illuminates our ignorance of the functions of some of these factors. Specifically, it did not predict the observed increase in cargo import in response to an elevation in RanBP1 concentration. This failure suggests that RanBP1 performs roles not included in the model. We know that RanBP1 can form a ternary complex with RanGDP and Impβ, for which no function has yet been ascribed, and this interaction was not included in the model (Chi et al., 1997; Plafker and Macara, 2002). However, such a complex would likely reduce rather than increase the rate of cargo import, by inhibiting RanGDP import by NTF2, and therefore reducing the RanGTP gradient. Additionally, RanBP1 can shuttle in and out of the nucleus, but this process would again reduce rather than increase cargo import, by depleting RanGTP from the nucleus (Plafker and Macara, 2000a). We speculate that RanBP1 might either facilitate the dissociation of transporters from the NPC, or alter the permeability of the NPC to transporters. These processes were not included in the model, but might be expected to increase transport rates. However, we have not detected any accumulation of RanBP1 at the nuclear envelope in intact cells (unpublished data). Finally, RanBP1 is known to form a complex with Mog1 in the presence of Ran (Steggarda and Paschal, 2000). Mog1 is necessary for efficient nuclear import, although its function remains unknown, and elevated RanBP1 might in some way alter Mog1 activity (Baker et al., 2001).
The model was successful in making two unexpected and accurate predictions about the behavior of the nuclear transport system: that elevations in the cellular concentrations of either impβ or RCC1 would inhibit rather than stimulate cargo import. The inhibitory effect of Impβ demonstrates that futile cycling of transport receptors across the nuclear envelope can significantly deplete the Ran gradient. The levels of Impβ and of similar karyopherins must therefore be carefully regulated to reduce the steady-state concentration of free receptors in the cell. We have suggested previously that transport cofactors such as RanBP3 and NPAP60 might also function to reduce futile cycling (Lindsay et al., 2001; Nemergut et al., 2002). We have not included these cofactors in the current model, but expect that they might reduce the inhibitory effect of high Impβ levels.
The effect of RCC1 is of interest because it is a consequence of the kinetic model for exchange that was used in the transport model. In our earlier study, virtual cell simulations of Ran transport, we used simple Michaelis-Menten kinetics that did not produce the same outcome (Smith et al., 2002). These results therefore suggest that the four-step in vitro kinetic analysis of exchange by RCC1 provides an accurate reflection of the in vivo kinetics. Our model also predicts that RanGAP is not rate-limiting, and this was validated experimentally.
We believe that our model provides a foundation on which to test other concepts in nuclear transport. For example, it is not apparent why the cell uses the Impα adaptor to bind cargo, rather than allowing all cargo to interact directly with the Impβ karyopherin. One speculation is that the coupling of cargo import to the export of Impα, mediated by Cas–RanGTP, might drive the creation of a higher nuclear–cytoplasmic cargo gradient. The potential effects of competition between karyopherins for binding to the same cargo, or to the NPC, are also complex. In the future we will test these ideas both experimentally and in silico, using a modified transport model.
Materials And Methods
We first constructed a compartmental model of the Impα–Impβ import system using the program Jarnac, a biochemical simulation package for Windows (Sauro et al., 2003). Biochemical reactions were localized to the cytoplasm or nucleus in an idealized cell. Flux rates between compartments were modeled as the product of a permeability factor and the concentration difference of a reactant between the nucleus and cytoplasm. The reactions were converted internally by Jarnac to a series of coupled ODEs. Jarnac was then used to solve the ODEs based on a set of initial values. Protein concentrations and reactions rates were taken where possible from published reports (Tables I and II). Permeability constants for Impβ, Impα, and GST-GFP were determined experimentally as described in Tables I–III. The binding of Impα to Impβ as well as Impα–Impβ to the cargo were modeled as two-step reversible reactions (Catimel et al., 2001). We also adopted GTP and GDP concentrations and a four-step reversible reaction for nucleotide exchange by RCC1 (Klebe et al., 1995b). Concentrations of endogenous cargo proteins carried by Impα–Impβ or Impβ alone were treated as floating parameters. To account for the reduction of the Ran gradient by transport receptors other than Impβ, we included a reaction pathway for a “generic” karyopherin. This pathway used reactions rates identical to those used for the Impβ. Evidence exists in yeast (Gilchrist et al., 2002) that RanGTP and CAS act cooperatively to actively displace Impα from its cargo. Therefore, we modeled this reaction as a single-step displacement.
After construction of the model, parameter optimization was conducted using a genetic algorithm implemented in the Jarnac scripting language (Bäck and Schwefel, 1993).
Constructs, cell culture, and transfections
The ORFs for human Ran and RCC1 were subcloned into a mammalian vector (pKmRFP1) to generate mRFP-Ran and mRFP-RCC1 fusion proteins (T. Chen, University of Virginia). The mRFP1 was provided by R. Tsien (University of California, San Diego, La Jolla, CA; Campbell et al., 2002) and recombinant Cas by M. Yamada (University of Virginia). HeLa cells were transfected with the mRFP vectors using Effectene (QIAGEN) according to the manufacturer's instructions. The medium was changed to Ringer's solution 12–24 h after transfection, and the cells were imaged using confocal microscopy.
HeLa cells were fixed and immunostained with anti-Ran and anti-Impα mAbs (Signal Transduction Laboratories) as described previously (Smith et al., 1998).
Protein binding assays
Proteins were purified as described previously (Plafker and Macara, 2000b; Brownawell and Macara, 2002; Smith et al., 2002). Recombinant Impβ was tested for its ability to bind both RanGTP and Impα. γ-Bind–Sepharose beads were first blocked with 3% BSA. GST-Impα and Impβ were then added in binding buffer and incubated for 1 h at 4°C (20 mM Tris, pH 7.4, 2 mM MgCl2, 200 mM NaCl, 0.1% BSA, 0.1% Tween, 2 mM DTT). The beads were washed with binding buffer lacking BSA. Proteins were eluted with Laemmli sample buffer, separated by SDS-PAGE and visualized by Coomassie brilliant blue staining. Binding of Impβ to the constitutively active Ran mutant, RanGTP(Q69L), and of CAS to GST–Impα and Ran(Q69L), were tested in a similar manner.
HeLa cells were cultured in Bioptek Delta-T dishes. The medium was exchanged to a physiological saline (25 mM Hepes, pH 7.2, 110 mM NaCl, 5 mM KCl, 2 mM CaCl2, 1 mM MgSO4, 10 mM glucose, 1 mM KH2PO4, 1 mg/ml BSA) before microinjection. A fusion protein of GGNLS was used as a cargo to quantify nuclear import. The GGNLS was diluted with microinjection buffer to a final concentration of 20 μM. Co-injected proteins were added to give a final concentration of 5 μM in the injection mixture. Samples were centrifuged at 10,000 g at 4°C for 5 min before injection.
A confocal microscope (model Meta LSM510; Carl Zeiss MicroImaging, Inc.) was used to collect time-lapse images. Cells were imaged using a 20× objective lens (NA = 0.65) at a zoom of 1.0. Detector gain was kept beneath 550 to maintain linearity of pixel intensity to concentration. Cells were kept at 37°C and the medium was changed every 15 min to avoid changes in concentration caused by evaporation. Cells were injected using Femtotips on an Eppendorf Injectman NI 2. After cytoplasmic injection, time-lapse images were taken every 20 s to quantify nuclear import for GGNLS. Average nuclear pixel intensity was related to protein concentration by the use of a standard dilution series of known concentration. Initial rates were calculated as the linear regression of the first three data points of nuclear accumulation. The significance of differences between control and experimental rates over a range of GGNLS concentrations was computed by a two-way ANOVA, using the R statistical package (http://cran.r-project.org/). Differences were assumed to be significant for P < 0.05. To quantify permeability factors, Impβ and Impα were first labeled with OGI. After cytoplasmic injection, images were acquired every second. The slope of initial nuclear concentration versus initial rate provided the permeability factor for each protein.
Online supplemental material
Computer model. A complete schematic of the computer model is provided in Fig. S1. A script file containing the model is also included and can be run using the program Jarnac, which can be downloaded from http://18.104.22.168/downloads/ISetup.exe.
Permeability assays. Fig. S2 shows import rate data for recombinant OGI-labeled Impα and Impβ, and for GST-GFP.
Sensitivity analysis of permeabilities. Fig. S3 shows the effects of varying permeabilities for various components of the model on the initial rate of GGNLS import.
We thank Roger Tsien for the mRFP1 plasmid, Masami Yamada for recombinant CAS, and Ting Chen for the pK-mRFP-RCC1 and pK-mRFP-Ran plasmids. Thanks also to Leslie Loew, Boris Slepchenko, and Jim Schaff of the Virtual Cell Group at the University of Connecticut for useful early discussions concerning compartmental simulation.
This work was supported by grants GM50526 and GM064437 from the National Institutes of Health and Department of Health and Human Services.
Abbreviations used in this paper: GGNLS, GST-GFP-NLS; Impα, importin-α; Impβ, importin-β; NPC, nuclear pore complex; ODE, ordinary differential equation; OGI, Oregon green isothiocyanate.