Crowding-induced phase separation of nuclear transport receptors in FG nucleoporin assemblies
Abstract
The rapid (<1 ms) transport of biological material to and from the cell nucleus is regulated by the nuclear pore complex (NPC). At the core of the NPC is a permeability barrier consisting of intrinsically disordered phenylalanine-glycine nucleoporins (FG Nups). Various types of nuclear transport receptors (NTRs) facilitate transport by partitioning in the FG Nup assembly, overcoming the barrier by their affinity to the FG Nups, and comprise a significant fraction of proteins in the NPC barrier. In previous work (Zahn et al., 2016), we revealed a universal physical behaviour in the experimentally observed binding of two well-characterised NTRs, Nuclear Transport Factor 2 (NTF2) and the larger Importin-β (Imp-β), to different planar assemblies of FG Nups, with the binding behaviour defined by negative cooperativity. This was further validated by a minimal physical model that treated the FG Nups as flexible homopolymers and the NTRs as uniformly cohesive spheres. Here, we build upon our original study by first parametrising our model to experimental data, and next predicting the effects of crowding by different types of NTRs. We show how varying the amounts of one type of NTR modulates how the other NTR penetrates the FG Nup assembly. Notably, at similar and physiologically relevant NTR concentrations, our model predicts demixed phases of NTF2 and Imp-β within the FG Nup assembly. The functional implication of NTR phase separation is that NPCs may sustain separate transport pathways that are determined by inter-NTR competition.
Editor's evaluation
This theoretical study describes the interaction of a planar brush or film of the resident unstructured components of the nuclear pore complex (NPC) called nucleoporins (FG-nups) and different nuclear transport receptors (NTRs). The authors describe impacts of competitive binding that give rise to enrichment of the NTRs, NTF2 and importin-β, at different depths of the FG-nup film, which could relate to experimental observations in other studies, as well as evidence that crowding could promote the rate of nuclear transport by modulating FG-NTR binding/unbinding. The conclusions were found to be generally supported by the data, relevant to the field of nuclear transport, and able to make specific predictions that can be experimentally tested in the future.
https://doi.org/10.7554/eLife.72627.sa0Introduction
Embedded in the nuclear envelope are nuclear pore complexes (NPCs), large hour-glass shaped channels (inner diameter ∼40 nm) that regulate biomolecular transport between the cytoplasm and nucleoplasm (Alberts, 1994; Wente, 2000). The NPC presents an exclusion barrier to inert molecules, with the degree of exclusion increasing with molecular size. This physical barrier arises from a dense (mass density 100–300 mg/ml) assembly of moderately cohesive intrinsically disordered phenylalanine-glycine nucleoproteins (FG Nups; Hoogenboom et al., 2021). In addition, the barrier contains relatively high numbers (∼20–100) of nuclear transport receptors (NTRs), globular proteins that facilitate the translocation of cargo by transiently binding to the FG Nups (Lowe et al., 2015; Kim et al., 2018; Hayama et al., 2018). The known roles of NTRs in nucleocytoplasmic transport include ferrying specific cargo in and/or out of the nucleoplasm, returning RanGTP to the cytoplasm, and enhancing the exclusion of inert molecules (Jovanovic-Talisman et al., 2009; Aitchison and Rout, 2012; Lowe et al., 2015; Jovanovic-Talisman and Zilman, 2017; Kapinos et al., 2017). However, it remains to be fully elucidated how different NTRs organise themselves within the permeability barrier and how this organisation affects transport (Stanley et al., 2017; Jovanovic-Talisman and Zilman, 2017; Hoogenboom et al., 2021).
Qualitatively, NTRs are required to facilitate cargo transport, yet their presence in high numbers poses a significant risk of jamming the transport channel due to crowding effects (Hoogenboom et al., 2021). These seemingly contradictory phenomena have inspired various propositions about more subtle roles of NTRs in the NPC, such as their being essential to maintaining the barrier to non-specific transport; and the existence of spatially segregated, separate transport pathways through the NPC, where different NTRs and/or cargoes may take different trajectories through the barrier (Shah and Forbes, 1998; Yang and Musser, 2006; Naim et al., 2007; Fiserova et al., 2010; Yamada et al., 2010; Ma et al., 2012; Kapinos et al., 2014; Lowe et al., 2015; Lim et al., 2015; Ma et al., 2016; Kapinos et al., 2017; Perez Sirkin et al., 2021). Furthermore, NTRs of varying size and affinity to the FG Nups may play different roles in maintaining efficient transport, e.g., some smaller and more cohesive NTRs may play the role of a cross-linker, modulating the distribution of FG Nup mass in the pore, thereby influencing the trajectories of larger NTRs or cargoes (Perez Sirkin et al., 2021). An alternative mechanism, involving the switching between import and export transport states, has also been proposed (Kapon et al., 2008). Another observation to keep in mind is that the apparent binding affinities of NTRs to FG Nups in vitro appear too tight to enable the fast transport as seen in native NPCs; this can be strongly modulated, however, by the presence of other cellular proteins that compete – non-specifically – with the NTRs close to the FG Nup mass (Tetenbaum-Novatt et al., 2012; Lennon et al., 2021).
It is difficult to test different hypotheses regarding how NTR crowding affects the NPC barrier in an in vivo setting, due to the complexity of probing multitudes of diverse proteins in a dense nanoscale channel, as is the case in the NPC. To circumvent this complexity, various studies have reverted to much simpler in vitro FG Nup and NTR assemblies that resemble the physical environment of the NPC, e.g., considering the behaviour of FG Nups in solutions or assemblies with a similar mass density of FG Nups as found in the NPC (∼200 mg/ml) (Ghavami et al., 2016; Davis et al., 2020). Particularly well-controlled model systems are polymer film assemblies, where copies of an FG Nup are anchored to a planar surface and where NTRs (typically of one type) are introduced in the bulk volume above the surface (Eisele et al., 2010; Schoch et al., 2012; Eisele et al., 2012; Kapinos et al., 2014; Zahn et al., 2016; Vovk et al., 2016). It merits emphasising that FG Nup polymer film assemblies are but a simplified model for the NPC in that, typically, only one type of FG Nup and one type of NTR are probed, whereas in the NPC there are multiple different types of FG Nups and NTRs in the inner channel; and in that – in polymer film assemblies – the geometry of the nanoconfinement of the FG Nups and NTRs arises from a single hard planar wall, as compared with the more complex cylindrical confinement in the NPC. Nonetheless, they provide a most suitable platform to discover general principles of FG Nup and NTR behaviour that may be obscured in experiments on in vivo NPCs; and provide a foundation for a step-by-step increase in the level of complexity towards that of real NPCs. For the behaviour of planar films containing one type of FG Nup and one type of NTR, the main findings thus far have been: (1) that NTRs of one type (such as NTF2 and Importin-β [Imp-β]) bind to FG Nups in a rather generic way, suggesting universal physical principles – based on negative cooperative binding – might govern their behaviour (Vovk et al., 2016; Zahn et al., 2016); (2) NTRs readily penetrate the FG Nup films, with only limited effects on the collective morphology of the FG Nups (little swelling or compaction) (Eisele et al., 2010; Kapinos et al., 2014; Wagner et al., 2015; Vovk et al., 2016; Zahn et al., 2016); (3) that such systems can replicate the basic selective mechanism in the NPC, i.e., inert proteins tend not to penetrate the collective FG Nup phases but NTRs do, consistent with in vivo NPC functionality and with experiments on bulk solutions of FG Nups and NTRs (Schmidt and Görlich, 2015; Schmidt and Görlich, 2016); (4) the number of transport proteins in the FG Nup films can vary by orders of magnitude as a function of NTR numbers in solution above the film, in a highly non-Langmuir manner, where complex many-body interactions preclude the use of simple one-to-one binding models (Eisele et al., 2010; Schoch et al., 2012; Kapinos et al., 2014; Wagner et al., 2015; Vovk et al., 2016; Zahn et al., 2016). With the caveat that only a subset of NTRs have been probed, investigations of planar assemblies of FG Nups and NTRs highlight the fine-tuned balance of FG Nup attachment density, FG Nup cohesion, FG Nup-NTR affinities, the amount of non-specific proteins in the system, and NTR concentrations, where minor changes in this balance can lead to qualitatively different binding scenarios (Vovk et al., 2016; Zahn et al., 2016; Stanley et al., 2017; Lennon et al., 2021).
Taking a next step towards the complexity of the in vivo NPC and considering the large number (∼20–100) of NTRs of different sizes and affinities in the NPC inner channel (Lowe et al., 2015; Peters, 2009a; Peters, 2009b; Lim et al., 2015; Kim et al., 2018; Hoogenboom et al., 2021), one may next ask how the binding affinity of a specific NTR to a planar assembly of FG Nups depends on the binding behaviour of other NTRs. Physiologically, the answers to this question may explain how binding and thereby transport of specific NTRs, and their associated cargoes, can be modulated (if at all) by the presence of other types of NTRs, since this would directly impact on the transport function of the NPC. More generally, such answers will aid our understanding of how the NPC maintains fast and efficient transport whilst hosting a dense milieu of FG Nup motifs, NTRs, and cargoes. Finally, a systematically probing of crowding effects provides a means to assess various conceptual models of NPC transport such as the ‘Kap-centric’ barrier model (Wagner et al., 2015) and reduction-of-dimensionality models (Peters, 2005).
Here, we aim to understand how the crowding by different NTRs may affect the organisation and thereby transport functionality of FG Nup assemblies. Focusing on planar FG Nup assemblies as a well-controlled model system (Zahn et al., 2016), we use physical modelling to probe how the binding of one type of NTR could be affected by other types of NTRs, in a way that can be tested by currently available experimental setups. To explore the effects of mixed NTR crowding, we model a ternary mixture containing two different NTRs and one type of FG Nup in a polymer film assembly. When modelling FG Nups and NTRs, one can take various coarse-grained approaches, for instance one can take an all-atom approach (Miao and Schulten, 2009; Gamini et al., 2014; Raveh et al., 2016), or account only for the amino acids (Ghavami et al., 2013, Ghavami et al., 2014; Ghavami et al., 2018), or work only with the generic patterning of FG/hydrophobic/hydrophilic/charged ‘patches’ (Tagliazucchi and Szleifer, 2015; Matsuda and Mofrad, 2021; Davis et al., 2021), or completely neglect sequence details altogether in a ‘homopolymer’ approach (Moussavi-Baygi et al., 2011; Osmanovic et al., 2012; Osmanović et al., 2013b; Vovk et al., 2016; Zahn et al., 2016; Timney et al., 2016; Davis et al., 2020). Each approach has its strengths and weaknesses. For instance higher-resolution modelling can account for greater molecular complexity, but with the difficulty in probing large time and length scales. In contrast, homopolymer modelling provides access to larger time and length scales with more robust parameterisation and simplicity of execution, at the expense of (sub)molecular detail. In this work, we build on our previous minimal modelling framework based on treating FG Nups as sticky and flexible homopolymers and NTRs as uniformly cohesive spheres (Osmanović et al., 2013b; Zahn et al., 2016; Davis et al., 2020).
Results
Minimal physical modelling facilitates the understanding of many aspects of NPC functionality in terms of general principles, but it requires quantitatively accurate parameter settings to make meaningful predictions (Osmanović et al., 2013a; Jovanovic-Talisman and Zilman, 2017; Hoogenboom et al., 2021). In this work, the minimal modelling framework we employ is that of coarse-grained classical density functional theory (DFT), that has been previously used to model aspects of the NPC permeability barrier (Osmanović et al., 2013b; Zahn et al., 2016; Davis et al., 2020). To ensure that the setting of the parameters in our DFT model – outlined in the Computational methods section (below) – is commensurate with the behaviour of FG Nups (Nsp1) and NTRs (NTF2 or Imp-β) as probed in experiments, we make use of experimental data on FG Nup-NTR polymer film assemblies where the macroscopic binding between one type of FG Nup (Nsp1) and one type of NTR (NTF2 or Imp-β) was quantitatively probed (see Figure 1) Zahn et al., 2016. The experiments focused on a polymer film consisting of Nsp1 domains, attached to a flat surface at physiologically relevant densities (≈3.3 polymers/100 nm2), interacting with either NTF2 or Imp-β over a physiologically relevant range of concentrations (0.01–10 µM).
For consistency with the available experimental data we focus on the FG Nup Nsp1, which we treat as a homogeneous, flexible, and cohesive polymer consisting of beads of diameter nm (two amino acids per bead), resulting in the approximately correct persistence length for FG Nups (Lim et al., 2006; Zahn et al., 2016; Hayama et al., 2019; Davis et al., 2020). The intermolecular and intramolecular cohesive properties of FG Nups arise from a combination of hydrophobic motifs, e.g., FG, FxFG, and GLFG, and charged residues along the sequence which, in our model, are subsumed into one single cohesion parameter (equivalent to in the Computational methods section). In addition to the FG Nups, we also include the presence of the NTRs, NTF2 and Imp-, which we model as uniformly cohesive spheres of diameters nm (same as in the Computational methods section) and nm (same as in the Computational methods section), respectively (Zahn et al., 2016). The cohesive properties of the NTRs strictly refer to the attraction between the NTRs and FG Nups, arising at least in part from the hydrophobic grooves and charged regions on the former and the hydrophobic motifs and charged regions on the latter (Kumeta et al., 2012; Kim et al., 2013; Hayama et al., 2018; Frey et al., 2018). We do not consider any attractive interactions between NTRs themselves as there is no empirical evidence that suggests there are cohesive interactions between NTRs, a choice that is consistent with previous work and shown to well replicate experimental data of NTR binding to FG Nup assemblies (Osmanović et al., 2013b; Zahn et al., 2016; Vovk et al., 2016). As with the FG Nup intermolecular and intramolecular cohesive interactions, we subsume all contributions to the respective cohesive interactions FG Nup – NTF2 and FG Nup – Imp- through two more cohesion parameters and (equivalent to the cohesion variables and , respectively).
We begin the parametrisation of our model with the setting of so as to accurately reproduce the experimental thickness of Nsp1 planar film assemblies, at similar anchoring densities, as was done previously (Zahn et al., 2016; Fisher et al., 2018; Davis et al., 2020) (see Figure 1—figure supplement 1). With the here chosen interaction potential, the resulting cohesion strength is (with experiments conducted at ), that yields a film thickness nm, in our model defined as the height above the surface below which 95% of the total polymer density is included. We note that the value of found here is different to that of our previous work (Zahn et al., 2016), mainly due to the different choice of interaction potential. Nonetheless, both the model here and the model in Zahn et al., 2016 are parametrised using the same experimental data and produce the same film thicknesses. To further validate this value of , we checked whether the polymer film would exclude inert molecules, a basic property that has been observed for Nsp1 assemblies (Schmidt and Görlich, 2015; Schmidt and Görlich, 2016). The inert molecules are modelled as non-cohesive spheres of diameter , with i labelling the particle type, and their inclusion/exclusion in the film is quantified through the potential of mean force (PMF) given as
where is the one-body direct correlation function (see equation 19), is the external potential (see equation 11), is the polymer number density, is the one-dimensional (1D), i.e., integrated over the plane, polymer-particle cohesive pair potential (see equation 13), and is the excess chemical potential (here set to 0 ) (Roth et al., 2000; Roth and Kinoshita, 2006). For the inert molecules, the polymer-particle attraction (third term in equation (1)) is nullified. As expected, non-cohesive particles with increasing diameters (1.0, 2.0, 4.0, and 6.0 nm) experience a greater potential barrier upon attempted entry into the polymer film (see Figure 1—figure supplement 2), confirming that our simple model of an Nsp1 film replicates one of the key characteristics of the permeability barrier as seen in the NPC: the degree of exclusion of inert molecules increasing with molecular size (Mohr et al., 2009; Popken et al., 2015; Ghavami et al., 2016). We note that the presence of a maximum, close to the anchoring surface, in the relative density for inert particle diameters and 2 nm is due to the decrease in the polymer density closer to the surface (consistent with a small potential well close to the surface, see Figure 1—figure supplement 2b); the appearance of the maxima in the density depends upon the anchoring density of the FG Nups (not shown here).
We find that particles of diameters ≥4 nm experience PMFs of order 10 . In experiments, the size exclusion limit of the NPC was determined as ∼5 nm, albeit that this limit is rather soft and gradual (Keminer and Peters, 1999; Paine et al., 1975; Mohr et al., 2009). In our model, this empirical size limit corresponds to PMFs at least one order of magnitude greater than the thermal background 1 (Figure 1—figure supplement 2b). This order of magnitude of energetic barrier is therefore indicative of a physiologically relevant barrier, and is consistent with another – independent – modelling work that explicitly takes into account the amino acid sequence of the FG Nups in a pore geometry (Ghavami et al., 2016), as well as with a more recent numerical study that investigated the free energy barriers of nanoparticles of varying size and avidity (Matsuda and Mofrad, 2021).
Having shown that the now parametrised polymer model for Nsp1 films replicates the experimentally observed film thickness and the size selectivity of the NPC, we shift our focus to setting the parameters for the NTRs, NTF2 and Imp-. The cohesion strengths (same as in the Computational methods section) and (same as in the Computational methods section) are set by comparing absorption isotherms as calculated in DFT to those measured in experiment (Zahn et al., 2016) (see Figure 1). Using DFT, we compute the total density of NTRs in the film , , as the total NTR population within the film thickness divided by the area nm2 (converted to units of pmol/cm2).
where () is the number density of the NTRs. With only one free fitting parameter for each NTR (for the NTR-FG Nup interaction strength), the DFT binding isotherms are found to be in excellent agreement with experiment over three orders of magnitude in bulk NTR concentration (Figure 1 [top]), as was previously accomplished by a similar DFT model (where polymers were attached to the base of a cylinder) in Zahn et al., 2016. For the here chosen interaction potential, the resulting parametrised cohesion strengths are and for NTF2 and Imp-, respectively. It turns out that for the two (model) NTRs, despite the Imp- particle having an 1.5-fold larger excluded volume diameter as compared with the NTF2 particle. However, given the differences in diameters, and therefore a difference in the respective ranges of intermolecular interactions (see equation 13), we caution against directly comparing the cohesive properties of the two NTRs based on the cohesion strengths and alone. Of note, the concomitant film thicknesses from DFT – as a function of NTR concentration – are also in good agreement with experiment (Figure 1 [bottom]).
At this point, all the essential interaction parameters , and have been set by quantitative comparisons between DFT and experiment. Next, we investigate the effects of crowding of one type of NTR on the system. We quantify molecular crowding through two observables: (i) the packing fraction , where is the effective 1D number density of a particular NTR (labelled by i), and (ii) the PMF of a particular NTR, in the presence of other NTRs and the FG Nups (see Figure 2 and equation 1). For high crowding, one expects the packing fraction of a particular NTR to increase in magnitude and for the PMF to become more positive (with respect to the same system but with fewer NTRs), which is interpreted as a greater potential barrier (or, somewhat equivalently, a shallower potential well). We observe that both NTF2 and Imp- display higher levels of packing and higher-amplitude density oscillations within the Nsp1 film upon increasing their respective bulk concentrations (0.01, 0.1, 1.0, and 10.0 µM) (Figure 2 [top]). The density oscillations arise from layering/ordering effects mainly caused by packing against a hard planar wall, where particles prefer to pack closer to a flat surface; the layering of hard-spheres next to a planar wall is a well-known phenomenon (Patra, 1999; Roth et al., 2000; Deb et al., 2011). As is expected, in both systems, the maximum observed packing fraction () was located close to the surface (at 10 µM). For the here chosen NTR-particle sizes, it is expected that the packing fraction and PMF will be largely dictated by the interactions with the polymers and the crowding of other NTRs, with less significant effects arising from the particular implementation of the surface hardness. We note that the density oscillations for the Imp- particle show greater amplitudes as compared with the NTF2 particle (for the same concentration above the film), which is expected since the Imp- is larger in size and thus experiences more pronounced layering effects (Padmanabhan et al., 2010).
For both NTF2 and Imp-, the PMFs decrease in magnitude (but remain negative within the bulk of the film) upon an increase in bulk NTR concentration (Figure 2 [bottom]). Specifically, increasing the concentration in solution from 0.01 to 10.0 µM results in an approximate twofold decrease in the absolute value of the PMF (). The implication of this finding is that, at higher levels of packing in the film, it is relatively easier for bound NTRs to unbind from the polymer film, or, equivalently, less favourable for additional NTRs to enter the polymer film from the solution above it. This effect may primarily be attributed to the increased filling of space, i.e., molecular crowding, of the NTRs between the Nsp1 polymers. The results of Figure 2 are particularly relevant to the NPC ‘transport paradox’, where fast transport (∼1000 transport events per second) occurs in conjunction with strong – selective – binding. Whilst there are various explanations of the transport paradox (Hoogenboom et al., 2021), these results show how NTR crowding may facilitate the exit of NTRs from the NPC, noting that a decrease of in PMF well depth would imply a ∼100× faster rate for unbinding if we assume Arrhenius-like kinetics (Figure 2 [bottom]).
Given the constant presence of NTRs in the NPC inner channel (Lowe et al., 2015; Peters, 2009a; Peters, 2009b; Lim et al., 2015; Kim et al., 2018; Hoogenboom et al., 2021), we wondered how varying the available bulk concentration of one type of NTR – in this case NTF2– would affect the entry barrier for inert particles as seen in our modelled FG Nup film (see Figure 1—figure supplement 2). As we increase the bulk concentration of NTF2 from 0.01 to 1.0 µM, we find that the amount of inert particles, with diameters in the range nm, reduces within the FG Nup film, as compared to inert particles in a film with FG Nups only (see Figure 2—figure supplement 1a). The reduced amount of inert particles in the film is also quantified by an increase in the PMF, with the increase in PMF being larger for particles of greater size (see Figure 2—figure supplement 1b). This is in qualitative agreement with experimental observations of an enhancement of the NPC entry barrier to large inert particles through an increased presence of NTRs (Kapinos et al., 2017), and consistent with the notion that NTRs contribute to the integrity of the NPC transport barrier (Peters, 2009a; Peters, 2009b; Lim et al., 2015).
As a next step, we explore how the competition between NTRs may affect the binding, penetration, and distribution of NTRs in FG Nup assemblies. Specifically, we model the mixed crowding effects in a system containing the NTRs, NTF2 and Imp-, in an Nsp1 polymer film (see Figures 3 and 4 and their respective Figure supplements). As in the case with one type of NTR, we probe the packing fractions, PMFs , binding isotherms, and polymer film thickness, but this time keeping the amount of one NTR fixed at a physiologically relevant concentration (1 µM) (Zahn et al., 2016) whilst varying the concentration of the other NTR (Figure 3a). Upon increasing the bulk concentration of NTF2 from 0.01 to 10.0 µM while keeping the bulk concentration (in solution) of Imp- constant at 1 µM, the amount of bound Imp- dramatically drops and the remaining bound Imp- is redistributed towards the surface of the Nsp1 polymer film (Figure 3a [top] and Figure 3—figure supplement 1). Additionally, the density oscillations of Imp- within the film, which are evident at 0.01 µM of NTF2, smooth out upon increasing the amount of NTF2 to 0.1 µM. This shows that the presence of NTF2 directly modulates the distribution of Imp- within the film. Interestingly, upon increasing NTF2 from 0.1 µM while keeping the bulk concentration of Imp- constant, we observe NTR phase separation: an NTF2-rich phase within the FG Nup film and an Imp--rich phase at the film surface. We can attribute this phase separation to the crowding of the NTRs, since by decreasing the affinity between NTF2 and the FG Nups, thereby decreasing the competitive advantage of NTF2 over Imp-, we find an enhanced absorption of Imp- in the FG Nup polymer film (see Figure 3—figure supplement 2).
When considering binary systems of hard-spheres with different diameters subject to packing between planar walls, ignoring any attractive interactions between them, one typically observes the larger particles packing closer to the wall, as compared with the smaller particles (Padmanabhan et al., 2010). This effect, as measured per unit area, is due to the overall system entropy loss being less when the larger particles pack closer to the surface, rather than the smaller ones. Here we observe the opposite effect, with the (smaller) NTF2 particles remaining closer to the grafting surface, which is qualitatively consistent with a theoretical study investigating a binary mixture of attractive particles, where the larger particles were excluded for a distance from a planar surface of up to twice the particle diameter (Padmanabhan et al., 2010). Here, however, we observe the depletion of the larger NTR (Imp-) over much larger distances (in ) for high bulk concentrations of NTF2, apparently dictated by the polymer film thickness.
An intuitive explanation for the crowding-induced phase separation is that the smaller NTF2 competes more readily for binding sites (that are spread uniformly along the polymer in our model) deep within the film, closer to the grafting surface, without paying a substantial entropic penalty for rearranging the polymers. In contrast, closer to the film surface, the larger Imp- binds more readily, because of its overall stronger binding propensity, where the polymers are more diffuse (note , spread over a larger particle surface for Imp-). Indeed, the distribution of NTF2 in the film largely follows the polymer density as a function of distance from the grafting surface, indicating that with its smaller size, NTF2 benefits more from the enhanced concentration of polymer beads (and therewith of binding sites) without having to pay a substantial entropic cost (as for Imp-) for penetrating the polymer film. A similar phenomenon has also been observed in a binary system of small and large positively charged particles soaked in a solution of anions, with a negatively charged surface at the bottom (Fang and Szleifer, 2003).
Throughout the changes in incorporation of NTF2, the density of the polymers did not show noticeable changes. The modulation of Imp- via changes in NTF2 numbers is also articulated in terms of the PMF , where the Imp- PMF is an approximate square well at 0.01 µM of NTF2, but for higher NTF2 concentrations gradually transforms into a pronounced and sharper well located at nm, i.e., at the surface of the film, with the formation of a barrier to enter the rest of the film (Figure 3a [bottom]).
We verified if similar effects resulted when increasing the concentration of Imp- for a given, constant, NTF2 concentration set at 1.0 µM (see Figure 4a and Figure 4—figure supplement 1). We observe no significant change to the Nsp1 or NTF2 packing factions (including the PMF and binding isotherm) upon increasing the concentration of Imp- in solution from 0.01 to 1.0 µM (see also Figure 4b). We have not explored high bulk concentrations (>1 µM) of Imp-, as these yield highly oscillatory packing fractions and therewith are computationally more challenging in our DFT model. However, we expect that further incorporation of Imp- would eventually change the distribution of NTF2 in the film.
To further establish the factors that define NTR crowding and competition in our model Nsp1 system with two types of NTR-like particles, we explored the phase diagram as a function of the respective affinities of the two NTRs to Nsp1, and as a function of their respective sizes . Varying the particle sizes between 2 and 6 nm, and the affinities between ∼0 and ∼10 (Figure 3—figure supplement 2, Figure 5—figure supplement 1, and Figure 5—figure supplement 2), we find mixed and segregated phases of the NTR-like particles in the Nsp1 assembly (Figure 5). Not surprisingly, full mixing between the two types of particles occurs when and , i.e., when the particles are identical. Within the explored parameter range, phase separation occurs when these respective ratios become sufficiently small or large. Physically, the observed behaviour can be readily understood by noting that the entropic cost of penetrating the film increases with particle size (such that, e.g. this cost is larger for Imp- than for NTF2), and that the energy gain of penetrating the film increases with the particle affinity to Nsp1. At the cross-over between these two competing effects, particles will favour binding at the surface of the FG Nup assembly. Part of this is trivial, e.g., small enough particles would readily penetrate the film, as their presence does not significantly affect the entropy related to the various possible polymer configurations; consequently, substantially larger disparities in the respective FG Nup-NTR cohesion strengths are required to phase-separate systems with NTR-like particles sized ≲4 nm (Figure 5—figure supplement 2). The intriguing aspects here are that the polymer-related entropic cost of NTR binding just about kicks in when NTR sizes are a few nanometres (as is true for physiological NTRs), which is at least fourfold larger than the FG Nup persistence length and bead size that are of the order of 1 nm. Additionally, the respective sizes and Nsp1 affinities of NTF2 and Imp- are such that inter-NTR competition readily drives the system from mixed to phase separated states.
Discussion
NTRs are crucial to the function of the NPC and are increasingly recognised to be continuously present in the FG Nup assembly that dictates transport selectivity through the NPC in vivo (Lowe et al., 2015; Peters, 2009a; Peters, 2009b; Lim et al., 2015; Kim et al., 2018; Hoogenboom et al., 2021), yet their configuration inside the NPC remains poorly characterised. In this study, we identified the physical behaviour that can occur when different types of NTRs partition in an FG Nup assembly that mimics the NPC transport channel by its FG Nup density, by being confined to an assembly of only few tens of nm thick, and by binding or excluding NTRs and inert particles in a manner consistent with NPC transport functionality.
Specifically, we have made quantitative and testable predictions regarding the effects of mixed crowding on the binding of different NTRs to FG Nup planar assemblies. Our results are based on a minimal coarse-grained model implemented in a mean-field approach, using classical DFT, similar to our previous model for the binding of a single type of NTR to FG Nup assemblies (Zahn et al., 2016), that treats FG Nups as sticky and flexible homopolymers and NTRs as isotropic and cohesive spheres, with excluded volume interactions between all components based on their known sizes. The model here includes three free interaction parameters, corresponding to the cohesive FG Nup-FG Nup and NTR-FG Nup (for two types of NTR) interactions, i.e., no cohesive interactions between NTRs themselves. These cohesion parameters were parametrised using experimental data for Nsp1 film assemblies and binding thereto of one type of NTR (NTF2 or Imp-).
In Zahn et al., 2016, the FG Nup assembly was modelled via an effectively 2D approach of polymers grafted at the bottom of a wide cylinder with rotational symmetry. Here, we have instead imposed translational invariance across the anchoring surface of the FG Nups therefore reducing the problem to an effectively 1D one, where the distribution of particles depends only on the height above the surface. Apart from changes in the effective interaction parameters that are largely due to the different choice of interaction potential here, the results are fully consistent with (Zahn et al., 2016). Compared with this previous work, however, we have here studied the effects of crowding on the NTR binding to the FG Nup assembly, and more specifically considered the case of simultaneous binding by two different types of NTRs, NTF2 and Imp-, to a planar assembly of Nsp1.
Based on the parametrised model, we have shown that increased crowding of one type of NTR results in a shallower potential well within the film, implying that individual NTRs will have a smaller potential barrier to leaving the film in the presence of more NTRs. The origin of this effect is due to an interplay between the further occupation of volume within the film (entropic) and the increased competition for binding sites. This result has important implications for the NPC: when there is a large influx of material into the channel from either the cytoplasm or nucleoplasm, the exit of said material should be faster since the increased crowding effects will reduce the free energy barrier – thus increasing the likelihood – to leave the pore, with a predicted increase in unbinding rates of two orders of magnitude in the concentration range explored here. While we note that there are multiple factors involved in determining transport speed (Hoogenboom et al., 2021), this scenario highlights the importance of NTRs as an essential component in the NPC transport barrier (Lim et al., 2015) and, specifically, implies that the NPC could perform more efficiently and faster with higher numbers of NTRs present in its inner channel, as has indeed been observed in experiments with Imp- (Yang and Musser, 2006).
In addition, we have found that with increased incorporation of NTF2 within the FG Nup film, the amount of absorbed Imp- was reduced and its distribution within the film was changed. Typically, the smaller NTR binds closer to the bottom of the film (where the polymer packing is higher) and the larger NTR is more likely to bind to the top of the film (where the polymer packing is lower). This observation can be attributed to the smaller entropic penalty incurred when smaller NTRs (here: NTF2) penetrate the FG Nup assembly, compared with larger NTRs (here: Imp-). The intriguing aspect here is that within a physiologically relevant parameter range, NTF2 modulates the absorption and penetration of Imp- in the FG Nup assembly, driving the system to phase separate, resulting in spatially segregated regions that are enriched in the two respective NTRs.
Generally, the here discussed competitive binding phenomenon bears similarity to the observation that NTR binding to FG Nups may be modulated by more weakly, non-specifically binding proteins, provided that the latter are present at sufficiently high concentrations, as in fact validated by experiments on planar FG Nup assemblies (Lennon et al., 2021). Physically, the differential binding of smaller and larger NTRs can be compared with that of smaller positively charged particles packing closer to a negatively charged surface, with the larger positively charged particles sitting on top of the smaller ones, due to the higher entropic costs of packing larger particles close to an attractive surface (Fang and Szleifer, 2003). The case here is different, however, in that the polymer behaviour of the FG Nups further enhances differences in entropic cost for absorption of differently sized particles.
Given the simplified nature of our NTR representation as homogeneous spheres, we cannot draw definite conclusions about how the observed NTR behaviour may depend on any cargo associated with the NTRs. However, given the general trends in the phase diagram (Figure 5), we speculate that the binding of differently sized cargoes may lead to further magnified phase separation between NTR-cargo complexes in FG Nup assemblies. For similar practical considerations, the details of multivalency/avidity at the NTR surface are beyond the scope of this work, but based on our previous work on isolated FG Nups (Davis et al., 2021), we expect that such multivalency/avidity will play a large role in the kinetics of layering and entry/exit of the FG Nup assembly; and note that such multivalency also presents subtle challenges with regard to the thermodynamics of the selective barrier (Matsuda and Mofrad, 2021).
It is important to note that the observed height-dependent phase separation, as predicted in our minimal 1D model (assuming symmetry parallel to the anchoring surface), may not be the only way for NTRs to spatially segregate. Considering the pore geometry of the NPC, it would be interesting to extend the model developed here to two or three dimensions, relaxing the lateral symmetry assumption (Osmanović et al., 2013b), and to consider how radial gradients in polymer/FG Nup density – as, e.g., in Osmanovic et al., 2012; Perez Sirkin et al., 2021 – could facilitate phase separation along the radial axis of the NPC channel. An immediate consequence of this is that transport pathways through the NPC are most likely dependent on the type of NTR, with potentially separate transport pathways mediated and modulated by different NTRs.
Given that there is a stable population of Imp- in the NPC barrier and given that changes to this population affect the selective properties of the NPC (Lowe et al., 2015), our results suggest that NTR crowding plays a substantial role in determining the performance of the NPC barrier. Additionally, the observation of a phase-separated state between two distinct NTRs has implications on how the NPC maintains high-throughput transport despite high NTR densities. Consistent with experimental observations on NPCs (Lowe et al., 2015), Imp- is found to occupy regions of lower FG Nup density (as shown here in planar FG Nup assemblies), where our results here demonstrate that such a distribution of Imp- can at least in part be attributed to competitive binding of other, smaller NTRs to regions of higher FG Nup density. Finally, our results enrich reduction-of-dimensionality and ‘Kap-centric’ perspectives to the NPC transport barrier (Peters, 2005; Wagner et al., 2015), in that small NTRs preferentially occupy regions in the NPC inner channel that are denser in FG motifs whilst larger NTRs preferentially occupy more dilute regions, leading to spatially segregated transport pathways, and/or with the specific transport associated with one type of NTR being modulated by another type of NTR, thus providing additional levels of control for selective transport through the NPC.
Methods
We use classical DFT, an equilibrium mean field theory developed in previous works (Osmanović et al., 2013b; Zahn et al., 2016; Davis et al., 2020), to model the FG Nup-NTR planar film assembly consisting of NTF2, Imp-, and Nsp1.
We first build a rather generic physical model consisting of a ternary mixture, i.e., a -component system with , containing two types of free particles denoted by (which will describe the NTRs) and one type of polymer labelled as (that will describe the anchored FG Nup Nsp1). Here, we do not explicitly describe the small-scale solvent molecules, as the dominant interactions come from the NTRs and the FG Nups. However, the solvent medium is implicitly accounted for when parametrising the intramolecular and intermolecular interactions between FG Nups and the interactions between NTRs and FG Nups (Zahn et al., 2016; Hoogenboom et al., 2021). In this system, the numbers of the two different types of free particles (components ) are given as , diameters are , and chemical potentials are . In addition to the free particles, there are flexible homopolymers each consisting of beads, where each bead has a diameter of 0.76 nm (corresponding to two amino acids per bead). This choice of and produces the approximate contour and persistence length of an Nsp1 FG domain (Lim et al., 2006; Zahn et al., 2016; Hayama et al., 2019; Davis et al., 2020). The polymers are attached uniformly to a flat surface of area nm2, resulting in an attachment/grafting density of ≈3.3 polymers/100 nm2, which is in line with the densities in the native NPC and in in vitro experiments (Zahn et al., 2016; Davis et al., 2020). For simplicity, it is assumed that the system is translationally symmetric along the directions parallel to the grafting surface, therefore, after integrating out the plane the density of component , can be transformed into a function of (height above the anchoring surface) only:
resulting in an effective 1D DFT (Davis et al., 2020).
Furthermore, the grand potential free energy functional , that provides a complete thermodynamic description of the entire system, is approximated as a sum of terms.
The first term is the free energy of a set of two types of ideal gas (for the two types of NTRs) given as
where ( is Boltzmann’s constant) and is the appropriate thermal de Broglie wavelength for component i. The second term describes the ideal polymer free energy (in the presence of a mean field ) and is given as (Fredrickson, 2005)
where the canonical partition function is
where , , is an energy function that imposes a rigid bond length between directly connected beads (pairs for ) on polymer , and is the mean field contribution of the ith bead belonging to polymer (Davis et al., 2020). The third term, , is the additional term that mathematically arises from introducing a mean field description of the affinities between the polymers (Osmanovic et al., 2012) and is given as
The fourth term, , accounts for the external potential imposing the hardness of the anchoring surface and is determined through
where is a repulsive external potential energy function taking a Weeks-Chandler-Anderson form
where is the diameter of the constituent particle for component , is the maximum energy barrier of the wall (chosen sufficiently high so that the number density of all components is zero at and below the surface), and . The fifth term, , imposes an exchange of the NTRs with an external reservoir, leading to a grand canonical ensemble, and is written as
Consistent with our previous work (Davis et al., 2021), the intramolecular and intermolecular cohesive interactions are based upon the Morse potential (in three dimensions)
where is the distance between two particles of type and type , is the cohesion strength, nm-1 is an inverse decay length of the pair potential, and . The potential in equation (13), valid in three spatial dimensions, is then integrated over the plane, making it henceforth only dependent on , and shifted such that for so as to keep the cohesive interactions short ranged. Thus we can now define the sixth term, , as the free energy contribution from intermolecular and intramolecular attractive (‘cohesive’) interactions determined by
The final term, , accounts for the intermolecular and intramolecular excluded volume interactions, including imposing polymer chain connectivity, and is given by
where and are the White Bear (hard-sphere) (Roth et al., 2002) and chain connectivity functionals (Yu and Wu, 2002) given by the equations
where and are, respectively, the1D scalar and vectorial weighted densities and is the radius of a polymer bead (see Roth, 2010). Essentially, the White Bear functional removes – from the free energy – contributions from system configurations with particle overlaps, which is the definition of a hard-sphere system. The hard-sphere functional is built using a standard formalism known as fundamental measure theory, attributed to Rosenfeld, which begins with a decomposition of the Mayer function into weight functions that contain geometrical information of spherical particles (see Rosenfeld, 1989; Roth et al., 2002; Roth, 2010 for details). The chain-connectivity functional is built in a similar manner, but removes from the free energy contributions from system configurations that violate how consecutive beads on a chain should be connected, as specified by the particular polymer model (see Yu and Wu, 2002 for more details).
Thus, the dimensionless grand potential can be written more explicitly as
To find the set of density distributions – for the particles and polymer – and the polymer mean field in the equilibrium state, the following equations must be solved self-consistently:
where the notation represents a functional derivative with respect to and is the one-body direct correlation function given by
For the free particles one can decompose the chemical potential into two terms
where is the bulk density of the free particles of component i and is the excess chemical potential due to the intermolecular and intramolecular interactions. One can solve equations (18) self-consistently by invoking a fictitious time variable , where the solutions are found through an iterative procedure. This is expressed by the following
Finally, discretising space into slices of thickness and discretising fictitious time then yields the resulting discrete update rules which are solved numerically
where is the – midpoint – height above the surface of the spatial slice , labels discrete time, and in the last equation . The simulation parameters used here were , nm (with ), and .
We note that, for the free particles , an excess chemical potential ( for ) is referenced to a zeroed chemical potential that results in a free particle bulk concentration of 10 µM. Changing the excess chemical potential to results in an order of magnitude increase (for +2) or decrease (for –2) of the concentration in solution (see also Osmanović et al., 2013b).
Data availability
The source code used to generate all the simulation data in this manuscript is available on the Github repository: https://github.com/patherlkd/DFT-polymer-colloid (copy archived at swh:1:rev:c3a89c14cd928d37717d5f7a74cc1d1115b4bb7d). Figure 1—source code 1: Simulation parameters for the classical density functional theory code. Figure 3—source code 1: Simulation parameters for the classical density functional theory code.
References
-
Physical modeling of multivalent interactions in the nuclear pore complexBiophysical Journal 120:1565–1577.https://doi.org/10.1016/j.bpj.2021.01.039
-
Hard sphere fluids at a soft repulsive wall: a comparative study using Monte Carlo and density functional methodsThe Journal of Chemical Physics 134:214706.https://doi.org/10.1063/1.3593197
-
Competitive adsorption in model charged protein mixtures: equilibrium isotherms and kinetics behaviorThe Journal of Chemical Physics 119:1053–1065.https://doi.org/10.1063/1.1578992
-
Facilitated transport and diffusion take distinct spatial routes through the nuclear pore complexJournal of Cell Science 123:2773–2780.https://doi.org/10.1242/jcs.070730
-
BookThe Equilibrium Theory of Inhomogeneous PolymersOxford University Press.https://doi.org/10.1093/acprof:oso/9780198567295.001.0001
-
Assembly of Nsp1 nucleoporins provides insight into nuclear pore complex gatingPLOS Computational Biology 10:e1003488.https://doi.org/10.1371/journal.pcbi.1003488
-
Coarse-grained potentials for local interactions in unfolded proteinsJournal of Chemical Theory and Computation 9:432–440.https://doi.org/10.1021/ct300684j
-
Sol–gel transition in solutions of FG-Nups of the nuclear pore complexExtreme Mechanics Letters 22:36–41.https://doi.org/10.1016/j.eml.2018.04.006
-
Thermodynamic characterization of the multivalent interactions underlying rapid and selective translocation through the nuclear pore complexThe Journal of Biological Chemistry 293:4555–4563.https://doi.org/10.1074/jbc.AC117.001649
-
Karyopherins regulate nuclear pore complex barrier and transport functionThe Journal of Cell Biology 216:3609–3624.https://doi.org/10.1083/jcb.201702092
-
Permeability of single nuclear poresBiophysical Journal 77:217–228.https://doi.org/10.1016/S0006-3495(99)76883-9
-
Nanoscale mechanism of molecular transport through the nuclear pore complex as studied by scanning electrochemical microscopyJournal of the American Chemical Society 135:2321–2329.https://doi.org/10.1021/ja311080j
-
Karyopherin-independent spontaneous transport of amphiphilic proteins through the nuclear poreJournal of Cell Science 125:4979–4984.https://doi.org/10.1242/jcs.109520
-
Characterizing binding interactions that are essential for selective transport through the nuclear pore complexInternational Journal of Molecular Sciences 22:10898.https://doi.org/10.3390/ijms221910898
-
How to operate a nuclear pore complex by Kap-centric controlNucleus (Austin, Tex.) 6:366–372.https://doi.org/10.1080/19491034.2015.1090061
-
Super-resolution 3D tomography of interactions and competition in the nuclear pore complexNature Structural & Molecular Biology 23:239–247.https://doi.org/10.1038/nsmb.3174
-
Free energy calculations shed light on the nuclear pore complex’s selective barrier natureBiophysical Journal 120:3628–3640.https://doi.org/10.1016/j.bpj.2021.07.025
-
Transport-related structures and processes of the nuclear pore complex studied through molecular dynamicsStructure (London, England) 17:449–459.https://doi.org/10.1016/j.str.2008.12.021
-
Characterisation of the passive permeability barrier of nuclear pore complexesThe EMBO Journal 28:2541–2553.https://doi.org/10.1038/emboj.2009.200
-
Passive and facilitated transport in nuclear pore complexes is largely uncoupledThe Journal of Biological Chemistry 282:3881–3888.https://doi.org/10.1074/jbc.M608329200
-
Bistable collective behavior of polymers tethered in a nanoporePhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 85:1–8.https://doi.org/10.1103/PhysRevE.85.061917
-
Binary fluid with attractions near a planar wallPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 82:021507.https://doi.org/10.1103/PhysRevE.82.021507
-
Structure of binary hard-sphere mixtures near a hard wall: A simple weighted-density-functional approachThe Journal of Chemical Physics 111:6573–6578.https://doi.org/10.1063/1.479948
-
Translocation through the nuclear pore complex: Selectivity and speed by reduction-of-dimensionalityTraffic (Copenhagen, Denmark) 6:421–427.https://doi.org/10.1111/j.1600-0854.2005.00287.x
-
Functionalization of a nanopore: the nuclear pore complex paradigmBiochimica et Biophysica Acta 1793:1533–1539.https://doi.org/10.1016/j.bbamcr.2009.06.003
-
Translocation through the nuclear pore: Kaps pave the wayBioEssays: News and Reviews in Molecular, Cellular and Developmental Biology 31:466–477.https://doi.org/10.1002/bies.200800159
-
Size-dependent leak of soluble and membrane proteins through the yeast nuclear pore complexMolecular Biology of the Cell 26:1386–1394.https://doi.org/10.1091/mbc.E14-07-1175
-
Depletion potential in hard-sphere mixtures: theory and applicationsPhysical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics 62:5360–5377.https://doi.org/10.1103/physreve.62.5360
-
Fundamental measure theory for hard-sphere mixtures revisited: the White Bear versionJournal of Physics 14:12063–12078.https://doi.org/10.1088/0953-8984/14/46/313
-
Depletion potential between large spheres immersed in a multicomponent mixture of small spheresJournal of Chemical Physics 125:084910.https://doi.org/10.1063/1.2345200
-
Fundamental measure theory for hard-sphere mixtures: a reviewJournal of Physics. Condensed Matter 22:063102.https://doi.org/10.1088/0953-8984/22/6/063102
-
Transport selectivity of nuclear pores, phase separation, and membraneless organellesTrends in Biochemical Sciences 41:46–61.https://doi.org/10.1016/j.tibs.2015.11.001
-
Biomechanics of the transport barrier in the nuclear pore complexSeminars in Cell & Developmental Biology 68:42–51.https://doi.org/10.1016/j.semcdb.2017.05.007
-
Nucleocytoplasmic transport: a role for nonspecific competition in karyopherin-nucleoporin interactionsMolecular & Cellular Proteomics 11:31–46.https://doi.org/10.1074/mcp.M111.013656
-
Simple rules for passive diffusion through the nuclear pore complexThe Journal of Cell Biology 215:57–76.https://doi.org/10.1083/jcb.201601004
-
Gatekeepers of the nucleusScience (New York, N.Y.) 288:1374–1377.https://doi.org/10.1126/science.288.5470.1374
-
Nuclear import time and transport efficiency depend on importin beta concentrationThe Journal of Cell Biology 174:951–961.https://doi.org/10.1083/jcb.200605053
-
Density functional theory for inhomogeneous mixtures of polymeric fluidsThe Journal of Chemical Physics 117:2368–2376.https://doi.org/10.1063/1.1491240
Article and author information
Author details
Funding
Engineering and Physical Sciences Research Council (EP/L504889/1)
- Luke K Davis
- Bart W Hoogenboom
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dino Osmanović, Anton Zilman, and Anđela Šarić for discussions. We thank Ralf Richter for providing feedback on the manuscript. LKD acknowledges the biophysics research computing cluster at UCL that was used to perform the simulations and analysis. This work was funded by the UK Engineering and Physical Sciences Research Council (EP/L504889/1, LKD and BWH).
Copyright
© 2022, Davis et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,257
- views
-
- 240
- downloads
-
- 19
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Physics of Living Systems
The endoplasmic reticulum (ER), the largest cellular compartment, harbours the machinery for the biogenesis of secretory proteins and lipids, calcium storage/mobilisation, and detoxification. It is shaped as layered membranous sheets interconnected with a network of tubules extending throughout the cell. Understanding the influence of the ER morphology dynamics on molecular transport may offer clues to rationalising neuro-pathologies caused by ER morphogen mutations. It remains unclear, however, how the ER facilitates its intra-luminal mobility and homogenises its content. It has been recently proposed that intra-luminal transport may be enabled by active contractions of ER tubules. To surmount the barriers to empirical studies of the minuscule spatial and temporal scales relevant to ER nanofluidics, here we exploit the principles of viscous fluid dynamics to generate a theoretical physical model emulating in silico the content motion in actively contracting nanoscopic tubular networks. The computational model reveals the luminal particle speeds, and their impact in facilitating active transport, of the active contractile behaviour of the different ER components along various time–space parameters. The results of the model indicate that reproducing transport with velocities similar to those reported experimentally in single-particle tracking would require unrealistically high values of tubule contraction site length and rate. Considering further nanofluidic scenarios, we show that width contractions of the ER’s flat domains (perinuclear sheets) generate local flows with only a short-range effect on luminal transport. Only contractions of peripheral sheets can reproduce experimental measurements, provided they are able to contract fast enough.
-
- Biochemistry and Chemical Biology
- Physics of Living Systems
For drugs to be active they have to reach their targets. Within cells this requires crossing the cell membrane, and then free diffusion, distribution, and availability. Here, we explored the in-cell diffusion rates and distribution of a series of small molecular fluorescent drugs, in comparison to proteins, by microscopy and fluorescence recovery after photobleaching (FRAP). While all proteins diffused freely, we found a strong correlation between pKa and the intracellular diffusion and distribution of small molecule drugs. Weakly basic, small-molecule drugs displayed lower fractional recovery after photobleaching and 10- to-20-fold slower diffusion rates in cells than in aqueous solutions. As, more than half of pharmaceutical drugs are weakly basic, they, are protonated in the cell cytoplasm. Protonation, facilitates the formation of membrane impermeable ionic form of the weak base small molecules. This results in ion trapping, further reducing diffusion rates of weakly basic small molecule drugs under macromolecular crowding conditions where other nonspecific interactions become more relevant and dominant. Our imaging studies showed that acidic organelles, particularly the lysosome, captured these molecules. Surprisingly, blocking lysosomal import only slightly increased diffusion rates and fractional recovery. Conversely, blocking protonation by N-acetylated analogues, greatly enhanced their diffusion and fractional recovery after FRAP. Based on these results, N-acetylation of small molecule drugs may improve the intracellular availability and distribution of weakly basic, small molecule drugs within cells.