Simple biophysics underpins collective conformations of the intrinsically disordered proteins of the Nuclear Pore Complex
Abstract
Nuclear Pore Complexes (NPCs) are key cellular transporter that control nucleocytoplasmic transport in eukaryotic cells, but its transport mechanism is still not understood. The centerpiece of NPC transport is the assembly of intrinsically disordered polypeptides, known as FG nucleoporins, lining its passageway. Their conformations and collective dynamics during transport are difficult to assess in vivo. In vitro investigations provide partially conflicting results, lending support to different models of transport, which invoke various conformational transitions of the FG nucleoporins induced by the cargo-carrying transport proteins. We show that the spatial organization of FG nucleoporin assemblies with the transport proteins can be understood within a first principles biophysical model with a minimal number of key physical variables, such as the average protein interaction strengths and spatial densities. These results address some of the outstanding controversies and suggest how molecularly divergent NPCs in different species can perform essentially the same function.
https://doi.org/10.7554/eLife.10785.001eLife digest
Animal, plant and fungal cells contain a structure called the nucleus, inside which the genetic material of the cell is stored. For the cell to work properly, certain proteins and other molecules need to be able to enter and exit the nucleus. This transport is carried out by pore-like molecular “devices” known as Nuclear Pore Complexes, whose architecture and mode of operation are unique among cellular transporters.
Nuclear Pore Complexes are charged with a daunting task of deciding which of the hundreds of molecules it conducts per second should go through and which should not. Small molecules can pass freely through Nuclear Pore Complexes. However, larger molecules can only pass through the pore efficiently if they are bound to specialized transport proteins that interact with the proteins – called FG nucleoporins – that line the pore. A unique feature of the FG nucleoporins is that, unlike typical proteins, they do not have a defined three-dimensional structure. Instead, they form a soft and pliable lining inside the Nuclear Pore Complex passageway.
Exactly how interacting with transport proteins affects the structure and spatial arrangements of the FG nucleoporins in a way that allows them to control transport is not well understood. This is in part because existing experimental techniques are unable to study the structures of the FG nucleoporins in enough detail to track how they change during transport. The complexity and the diversity of the FG nucleoporins also make them difficult to model in detail.
Vovk, Gu et al. have developed a theoretical model that is based on just three basic physical properties of the FG nucleoporins – their flexibility, their ability to interact with each other, and their binding with the transport proteins. Future work can refine the model by incorporating further molecular details about the interactions between FG nucleoporins and transport proteins.
The predictions made by this simple model agree well with experimental results in a wide range of situations – from single molecules to complex spatial assemblies. They also explain why some of the experimental results appear to contradict each other and suggest how several outstanding controversies in the field can be reconciled. Because the model invokes only fundamental physical principles of FG nucleoporin assemblies, it shows that some of their general properties do not depend on the exact conditions. In particular, this might shed light on why Nuclear Pore Complexes in different organisms perform essentially the same function, although the details of their molecular structure may differ. This also suggests how the FG nucleoporins can be manipulated to build artificial devices based on the same principles.
https://doi.org/10.7554/eLife.10785.002Introduction
Nuclear Pore Complexes (NPCs) are biological 'nanomachines' that conduct all the transport between the nucleus and the cytoplasm in eukaryotic cells. NPCs participate in a vast number of regulatory processes in the cell, as well as pathological conditions such as viral disease and cancer (Dickmanns et al., 2015). Transport through the NPC is fast, highly selective and robust with respect to molecular noise and structural perturbations. Transport of relatively small cargoes up to several nanometers in size, or approximately 30–40 kD, occurs by pure diffusion, without specific interactions with the NPC constituents. By contrast, transport of larger macromolecules, such as import of transcription factors and export of mRNA particles is tightly controlled by the NPC. For efficient transport, macromolecules larger than several nanometers in size must be shuttled through the NPC by soluble nuclear transport proteins from a highly conserved family, known as Karyopherins (Kaps) in yeast or Importins/Transportins in vertebrates. Remarkably, despite its high selectivity and efficiency, the NPC does not consume metabolic energy during transport and does not possess an obvious 'gate' opening or closing during transport (Terry and Wente, 2009; Wente and Rout, 2010; Stewart, 2007; Feldherr and Akin, 1997; Mohr et al., 2009).
The spatial organization of the NPC and its transport mechanism are unique. The passageway through the nuclear envelope of about 35–50 nm in diameter and nm in length is formed by a structural scaffold that comprises multiple proteins of a combined size of MDaltons. This passageway is lined by a set of intrinsically disordered polypeptide chains, collectively known as 'FG nups' due to the large numbers of Phenylalanine-Glycine (FG) repeats in their sequence (Terry and Wente, 2009; Wente and Rout, 2010). Although the actual sequences of the FG nups can vary widely among different species, the overall structure, organization and the transport mechanism of the NPC are conserved (Terry and Wente, 2009; Hülsmann et al., 2012; Schmidt and Görlich, 2015). As the key component of the NPC transport mechanism, the FG nups set up the permeability barrier that prevents free passage of large macromolecules and serve as a template for the transient binding of the cargo-carrying transport proteins. NPCs are also remarkably resilient with respect to structural perturbations: many of the FG nups can be genetically deleted without impairing cell viability and without major effect on transport (Strawn et al., 2004; Popken et al., 2015; Feldherr et al., 2002; Hülsmann et al., 2012). Cargo-carrying transport proteins bind to the FG nups through multiple, yet relatively weak contacts. This binding is crucial for selective transport: interfering with it decreases the transport efficiency, or abolishes the transport altogether. Conversely, particles or molecules that normally cannot penetrate the NPC can be transported after chemical modifications that enables the to interact directly with the FG nups (Wente and Rout, 2010; Bayliss et al., 1999; 2000; Naim et al., 2009; Kim et al., 2013; Kumeta et al., 2012). NPC geometry and architecture are schematically illustrated in Figure 1.
Full understanding of the NPC transport still remains elusive, and various hypotheses remain unsettled. The consensus is that the binding of the transport proteins to the FG nups enables them to overcome the permeability barrier. The strength of this binding controls the transport selectivity and efficiency. Hence, transport proteins can be informally viewed as 'glorified enzymes' that lower the free energy barrier for the translocation through the NPC. Basic models that describe the transport as facilitated diffusion through the FG nup medium, modulated by the interactions with the FG nups, provide a good explanation of the selectivity of the NPC even in the presence of large amounts of molecular noise (Wente and Rout, 2010; Zilman et al., 2007; 2010; Frey and Görlich, 2007; Fernandez-Martinez and Rout, 2012). The overall veracity of these general principles has been demonstrated by creation of artificial nanochannels and nanomaterials that mimic NPC function and recapitulate many of its transport properties (Zilman et al., 2007; 2010; Frey and Görlich, 2007; Schmidt and Görlich, 2015; Zilman, 2009; Jovanovic-Talisman et al., 2009; Kowalczyk et al., 2011; Caspi et al., 2008; Jovanovic-Talisman et al., 2014).
Various models of the mechanistic involvement of the FG nups in transport have been proposed. In the 'virtual gate' model, the permeability barrier arises due to the entropic repulsion from the fluctuating FG nup chains (Zilman et al., 2007; Lim et al., 2007; Rout et al., 2003). In a related idea, an entropically stabilized FG nup 'brush' can be collapsed by the transport proteins thus opening the transport passageway (Lim et al., 2006; 2007; 2008). In another scenario, the permeability barrier arises from a gel-like network, stabilized by the hydrophobic interactions between the FG repeats. Transport proteins disentangle this gel via their binding to the FG repeats thereby allowing their passage through the pore (Hülsmann et al., 2012; Frey and Görlich, 2007; Frey et al., 2006). More complex models have been proposed that take into account the sequence inhomogeneity and local molecular properties of the FG nups, their possible spatial localization and dynamics (Kim et al., 2013; Patel et al., 2007; Yamada et al., 2010; Peters, 2009; Mincer and Simon, 2011; Cardarelli et al., 2012; Ma et al., 2012; Solmaz et al., 2013; Lowe et al., 2015). It is likely that the majority of the effects invoked in all these models contribute to the NPC transport mechanism to some degree. In particular, the FG nups possess various degrees of intra- and inter-chain 'cohesiveness' that can lead to formation of single and multi-chain aggregates (Schmidt and Görlich, 2015; Frey et al., 2006; Patel et al., 2007; Yamada et al., 2010; Hough et al., 2015; Milles and Lemke, 2011).One major contribution to FG nup cohesiveness is believed to arise from the weak binding of the hydrophobic FG repeats to each other. However, intrinsically disordered proteins are notoriously prone to aggregation and the cohesiveness can have multiple sources, including electrostatic, and -charge interactions as well as non-specific interactions between the non-FG parts of the chains (Uversky, 2002; Song et al., 2013; Borg et al., 2007; Milles et al., 2013). Cohesiveness of individual FG nups correlates with the ratio of the numbers of the hydrophobic to charged residues in their sequence (Yamada et al., 2010). However, the relative importance of the cohesiveness and its specific role in the transport mechanism are still under debate.
To add to the complexity of the system, the transport proteins are present in large numbers within the NPC, and can strongly affect FG nup conformations and dynamics (Zilman et al., 2007; 2010; Zilman et al., 2009; Jovanovic-Talisman et al., 2009; Yang and Musser, 2006;Lowe et al., 2015 ; Kapinos et al., 2014; Schoch et al., 2012; Milles et al., 2013). However, several controversies persist with respect to their contributions to the architecture and the function of the NPC. Specifically, it is still under debate whether the transport proteins swell or compact assemblies of FG nups (Kapinos et al., 2014; Wagner et al., 2015; Eisele et al., 2010; 2013).
Understanding the spatial organization and the collective dynamics of the FG nups and the transport proteins during the translocation processes is hindered by the scarcity of experimental methods and insufficient measurement accuracy to assess them in vivo on the relevant time (several milliseconds) and length (several nanometers) scales (Cardarelli et al., 2012; Ma et al., 2012; Yang and Musser, 2006; Yang et al., 2004; Grünwald et al., 2011; Dange et al., 2008). Consequently, computational and theoretical approaches - from atomistic to highly coarse grained - have become increasingly important in the investigations of the possible configurations of the FG nups and the transport dynamics within the NPC (Zilman et al., 2007; 2010; Mincer and Simon, 2011; Opferman et al., 2013; Osmanović et al., 2013a; 2013b; Tagliazucchi et al., 2013; Moussavi-Baygi et al., 2011a; 2011b; Ghavami et al., 2014; Ando et al., 2014; Gamini et al., 2014). However, parameterizations of such models still remain difficult because of the sheer complexity and diversity of the FG nups, exacerbated by conflicting or non-existent measurements of the relevant parameters. The existing parameterizations differ significantly in their physical assumptions and outcomes (Tagliazucchi et al., 2013; Ghavami et al., 2014; Ando et al., 2014; Gamini et al., 2014). As a case in point, even the radii of gyration of the FG nups and their interaction affinities with the transport proteins are still under debate (Yamada et al., 2010; Kapinos et al., 2014; Tetenbaum-Novatt and Rout, 2010; Isgro and Schulten, 2005; Eisele et al., 2010).
For further progress, it is imperative to establish the most pertinent physical features and variables controlling the FG nup conformations induced by the transport proteins. Recent in vitro studies provide the basis for systematic understanding of the conformations of the assemblies of the FG nups with transport proteins in order to relate their molecular properties to their nanomechanical behavior (Hülsmann et al., 2012; Frey and Görlich, 2007; Schmidt and Görlich, 2015; Lim et al., 2007; Frey et al., 2006; Kapinos et al., 2014; Wagner et al., 2015; Eisele et al., 2010; 2013). In one typical experimental setup, an FG nup assembly is grafted to a rigid surface in the presence of varying concentrations of the transport proteins (see Figure 1). The FG nup conformations are inferred from the measurements of layer height. Even this relatively simple experimental setup resulted in conflicting observations. Depending on the experimental conditions and the measurement technique, the transport proteins can either increase or decrease the layer height, or in some cases cause no measurable change. However, despite the variability between different FG nups and transport proteins, general behavior motif is emerging from these experiments. Typically, at low concentrations, the transport proteins do not affect the layer height. This is followed by a partial layer collapse and height decrease, accompanied by accumulation of the transport protein in the layer. Increasing the concentration further reverses the collapse and eventually leads to the swelling of the layer. Virtually all studied mixtures of the FG nups and transport proteins regardless of the species of origin or the natural localization in the NPC exhibit this general pattern of behavior although the degree of collapse and swelling may vary (Lim et al., 2006; 2007; Kapinos et al., 2014; Wagner et al., 2015; Eisele et al., 2010; 2013). Related patterns of behavior are observed in experiments with bulk solutions of FG nups mixed with various transport proteins. At sufficiently high concentrations, the FG nups form a dense phase which either absorbs or excludes the transport proteins, depending on the size of the latter and their interaction strength with the FG nups see Figure 1. The behavior is very general and is observed in a wide range or FG nups from different species (Hülsmann et al., 2012; Frey and Görlich, 2007; Schmidt and Görlich, 2015; Hough et al., 2015).
The generality of these behaviors - despite the large number of molecular factors affecting the FG nup behavior - suggests that it might be understood in terms of a small number of core organizing principles. In this paper we develop a coarse grained theory of the transport proteins-FG nup assemblies that captures only the essential physical features: (i) the flexible nature of the FG nups, (ii) their potential cohesiveness and (iii) the attractive interactions with the transport proteins. The model is investigated using mean field theory supported by coarse-grained simulations (Opferman et al., 2012; Opferman et al., 2013). Our approach is inspired by the successes of the simplified theories in explaining the properties of highly chemically complex and diverse materials in polymer science and soft condensed matter (Doi and Edwards, 1998; Flory, 1953; de Gennes, 1979). Systematic comparison of the model predictions with extensive experimental data shows that the model captures and explains the observed behavior in different regimes. The model suggests a resolution of some of the apparent conflicts in the experimental results and proposes how to reconcile the outstanding controversies regarding the relative importance of the cohesive and the entropic effects for FG nup behavior and NPC selectivity. The model also sheds light on the long standing discrepancies in the measurements of the binding affinities of transport proteins to FG nups.
The model identifies the key physical variables controlling the conformational behavior of FG nup-transport protein assemblies and suggests experimental ways of manipulating them. This essential theoretical framework can be systematically developed in the future with additional molecular and structural details. Beyond the NPC, the results of the model are interesting in a broader context of intrinsically disordered proteins and nanotechnological applications (Stuart et al., 2010; Tagliazucchi and Szleifer, 2015; Coalson et al., 2015). Many aspects of the unfolded protein behavior are still puzzling, and the conceptual and computational frameworks - many of which are built on the foundation of polymer physics - are currently being developed (Uversky, 2002; van der Lee et al., 2014; Tcherkasskaya et al., 2003; Das et al., 2015; Sherman and Haran, 2006). The model described here demonstrates the power of such approaches on a concrete example of an important family of intrinsically disordered proteins.
Mathematical model
Many of the molecular parameters determining the behavior of the FG nup assemblies and their interactions with transport proteins are not fully known. In particular, the number of binding sites on the transport proteins and affinities of their binding to the FG repeats are a matter controversy. The degree of inter- and intra- FG nup cohesion or the relative importance of the electrostatic vs hydrophobic effects are unknown (Bayliss et al., 1999; 2000; Frey et al., 2006; Yamada et al., 2010; Kapinos et al., 2014; Schoch et al., 2012; Tagliazucchi et al., 2013; Tetenbaum-Novatt and Rout, 2010; Isgro and Schulten, 2005; Eisele et al., 2010; Pyhtila and Rexach, 2003; Otsuka et al., 2008). To capture the most salient mechanical and energetic properties responsible for the collective conformational behavior of the FG nups, we model them as flexible polymeric chains. Although not all intrinsically disordered proteins can be described in this way (for instance, due to residual secondary structure), the nanomechanical and the statistical properties of the FG nups can be well described by the polymer worm-like chain model (Lim et al., 2007). The chains can be imagined as consisting of 'monomers', where each monomer roughly corresponds to one or several amino-acids; as will be shown below, the predictions of the model turn out to be insensitive to the exact choice of the monomer size within a realistic range (see the Results section). The transport proteins are modeled as rigid particles of appropriate volume that interact attractively with the FG nup chains. Although the interaction between the transport proteins and FG nups is believed to be predominantly mediated by Phenylalanines of the FG repeats lodging themselves into the hydrophobic grooves on the transport proteins, it is likely that other amino acids and other interactions, such as electrostatic interaction, participate in the binding as well (Bayliss et al., 2000; 2002; Tagliazucchi et al., 2013; Isgro and Schulten, 2005; Colwell et al., 2010; Liu and Stewart, 2005). Thus, the interactions of the FG nups with the transport proteins are taken into account on a coarse-grained level parameterized by an interaction parameter , which is proportional to the average binding energy of a monomer to a transport protein molecule. The main contributions to are enthalpic, although water network re-arrangement entropy contributes to the hydrophobic interaction as well. We discuss the range of experimentally motivated values of below. The inter- and intra-FG nup cohesive interactions are incorporated into the model in a similarly general fashion through the effective interaction parameter , as explained below in Equation (2). This description is appropriate for high salt concentrations conditions (150–300 mM) of the experiments, where the electrostatic interactions are highly screened (Zhulina et al., 2000; Barrat and Joanny, 1997). Similar types of coarse-grained models have been successfully used to describe microtubule-associated unfolded proteins (Leermakers et al., 2010), other bio-polymers (Attili et al., 2012; Akinshina et al., 2013) and DNA-protein interactions (Jung et al., 2012). It was shown in Yamada et al. (2010) that the observed cohesiveness of different FG nup segments correlates with their hydrophobic to charged amino acid content. Our parameter can be thought of as a quantification of this concept and its extension into multi-chain and multi-protein domain.
In addition to the interactions of the FG nups among themselves and with the transport proteins, the main physical factors responsible for the conformations of their assemblies with the transport proteins are the entropy of chain stretching and compression, and the volume constraints. Roughly speaking, the attractive interactions favor more compact structures because they allow formation of higher numbers of energetically favorable contacts, while the entropy of chain configuration favors more diffuse and open conformations. These physical considerations are encapsulated in the free energy that describes an assembly of polymer chains, each of contour length (Opferman et al., 2013; Milner et al., 1988; Halperin et al., 2011; Fredrickson et al., 2002):
The first term in this free energy stands for the entropic elasticity of the polymer chains; describes the trajectory of the -th chain in space and where is the Kuhn length of the chain, roughly corresponding to the bond length between adjacent 'monomers'. In the second term, is the local free energy density that includes all the interactions between the polymers and the transport proteins, written in the 'mean field' (well mixed) approximation in terms of their local volume fractions and , respectively (Opferman et al., 2012; 2013; Flory, 1953; de Gennes, 1979; Lai and Halperin, 1992):
where is the monomer volume and is the ratio of the transport protein and monomer volumes (note that a monomer does not have to be spherical). The first term describes the translational entropy of the transport proteins in the layer. The next two terms describe the excluded volume interaction (steric repulsion) between all the monomers of the chains and the transport proteins. The last two terms describe the attractive interaction between the FG nup chains and the transport proteins and among the FG nups themselves, respectively; negative and correspond to attraction. Other expressions for free energy can be used and provide qualitatively similar results (Osmanović et al., 2013a; Osmanovic et al., 2012).
Surface layer geometry
To further simplify the discussion and establish the key variables controlling the experimentally observed behaviors, we assume that the monomer density is uniform throughout the layer. In reality, the monomer density inside the layer decays away from the grafting surface. We also assume that the entropic elasticity of the chains is described by the Gaussian model (Alexander, 1977). This simple approximation cannot be used to predict the exact layer height, but it is a qualitatively good treatment for the moderately cohesive chains and moderate concentrations of the transport proteins of interest in this paper, because the cohesiveness causes layer compaction, as shown below (Opferman et al., 2012; 2013; Halperin et al., 2011; Alexander, 1977; Zhulina et al., 1991; Lai and Halperin, 1992). The overall conclusions of this paper do not depend on this approximation (see Figure 2—figure supplement 1). The discussion is further simplified by normalizing the average layer height by the chain length , introducing a new variable . With these, the free energy per unit area of a layer of chains grafted at a distance from each other becomes (Alexander, 1977; de Gennes, 1980)
where is the grafting density of the chains normalized by the average monomer cross-section (). The layer height related to the monomer density through the condition which expresses the fact that the total number of the monomers in the layer is constant, . At this level of approximation, the equilibrium layer height is found by the minimization of the free energy over and the transport protein concentration under the constraints that the chemical potential of the transport proteins in the layer and the osmotic pressure in the layer are equal to those in the outside solution of volume fraction . Importantly, because factorizes out of the free energy expression in Equation (3), the resulting equilibrium values of and are independent of the chain length . The chemical potential and the osmotic pressure of the outside dilute solution - assumed to be ideal - are and . This procedure has been described in detail and verified by coarse-grained brownian dynamics simulations in Opferman et al. (2012), Opferman et al. (2013). Finally, it is important to keep in mind that the calculated properties are equilibrium average values. On the molecular scale, the polymers are highly dynamic and their individual conformations fluctuate on the microsecond time scale.
Bulk solutions
In a bulk solution where the chains are not grafted to a surface but are freely floating in solution, the entropic stretching (first term in Equation (3)) is replaced by the translational entropy of the chains, so that the mean field free energy per unit volume is Doi and Edwards (1998), Flory (1953), de Gennes (1979)
This free energy becomes unstable for sufficiently large interaction parameters or , leading to a phase separation where a dense phase of transport proteins mixed with the FG nups coexists with a very dilute solution. The compositions of the dense and the dilute phases are determined from the equality of the osmotic pressures and the chemical potentials of the transport proteins and the FG nups in the coexisting phases (de Gennes, 1979; Zilman and Safran, 2002; Morse, 1969) see Appendix for details.
Dimensions of individual FG nups
The model can be used to calculate the dimensions of individual FG nup molecules in solutions. Due to the thermal motion, each flexible chain dynamically samples multiple spatial conformations that on average occupy a volume of size , which can be found through minimization of the free energy over with the condition ; the first term represents the entropic elasticity of chain conformations in space, while the second term describes the intra-chain interactions where is the free energy of Equation (2) with . Slightly different expressions for the free energy can be used, all leading to qualitatively the same results (de Gennes, 1979; Sherman and Haran, 2006; Sanchez, 1979).
As shown below, the overall qualitative predictions of the theory are robust with respect to the choice of model parameters within a physically feasible range. However, quantitative or semi-quantitative comparison with the experimental data requires a specific choice of the molecular parameters and . The approximate volume of a transport protein molecule can be calculated from its molecular mass and the average protein density g/cm3. For Karyopherin- (molecular mass kD, depending on the attached tag), nm3 and for NTF2, the transport protein specialized for the import of RanGDP into the nucleus (molecular mass kD), nm3. The estimates of the 'monomer' size are somewhat less well defined. The average distance between two adjacent amino acids on a polypeptide chain is nm and the side chain size varies in the range nm (Levitt, 1976; Zamyatnin, 1972; Quillin and Matthews, 2000). However, their effective size can be modulated by the bound ions present within the Debye screening length (of the order of nm at the experimental salt concentrations) or bound denaturant molecules. Thus, for comparison with experiments, the realistic monomer size lies within a range nm and its volume nm3; the upper limit corresponds to a 'monomer' composed of about four amino acids. Finally, the exact structure of the dense phase and its maximal molecular packing fraction are unknown (Frey and Görlich, 2007; Schmidt and Görlich, 2015; Milles et al., 2013). For the conversions between the theoretical volume fractions and the experimentally measured concentrations, we have assumed the maximal packing fraction in the dense phase, typical of the random close packing of dense molecular assemblies which normally lies in the range of (Nolan and Kavanagh, 1992).
Results
Conformations of an FG nup layer in the absence of transport proteins
To establish whether the model captures the basic biophysical characteristics of the FG nup assemblies, we first apply it to the case of an FG nup layer without transport proteins. In this case the transport protein density is , and the free energy per area of Equations (2) and (3) is minimized over to obtain the equilibrium average layer height. We emphasize again that once the layer height has been re-scaled by the polymer contour length , the theoretical predictions become independent of . This will be important in the analysis of the experimental data.
The conclusions of the model are summarized in Figure 2, which shows that the FG nup layer height decreases with the grafting distance , because the steric repulsion between the polymers - that maintains the chains being stretched on average - is higher at lower grafting distances. The theory predicts that the layer height is monotonically decreasing with the cohesion strength. This is expected because the cohesiveness favors more compact conformations with more favorable contacts while the entropic elasticity term favors more diffuse conformations. In polymer physics parlance, increasing the cohesiveness is similar to changing the solvent 'quality' from 'good' to 'bad' (Osmanović et al., 2013a; Eisele et al., 2013; de Gennes, 1979; Milner et al., 1988; Halperin et al., 2011; Lai and Halperin, 1992; Zhulina et al., 1991; Peleg et al., 2011; Moh et al., 2011). The model also shows that sufficiently strong cohesion not only decreases the layer height, but shifts the behavior into a qualitatively different regime. The FG nup internal conformation and cohesiveness can be characterized by the scaling exponent that describes the dependence of the height on the grafting distance , . For an ideal non-cohesive, purely entropically stabilized polymer brush, . Cohesiveness increases the value of the exponent towards , in which regime the layer effectively behaves as a material of constant density with independent of (de Gennes, 1979; Milner et al., 1988; Zhulina et al., 1991; Moh et al., 2011). Nevertheless, even in the highly cohesive regime, the layer has significant 'free space' occupied by the solvent (calculated as ).
Behavior of the FG nup layers reported in Kapinos et al. (2014), Schoch et al. (2012), Wagner et al. (2015) follows the general predictions of the theory. In particular, examination of the data shows that the height of FG nup layers in the absence of transport proteins decays faster than with the grafting distance , but slower that , as shown in Figure 3. The accuracy of the experimental measurements does not allow one to quantitatively differentiate between different FG nups based on the scaling exponent . However, interpreted in light of the theoretical model, the data strongly indicate the presence of significant cohesion in all FG nups. Importantly, Nsp1 segments of different lengths fall within the same family of curves, once normalized by their length, in accord with the theoretical predictions. These results indicate that the salient physical mechanisms responsible for the behavior of grafted FG nup layers are adequately captured by the model.
Effects of the transport proteins on the layer height
We now compare the theoretical predictions with the experimental data in the presence of transport proteins. Although at a first glance experimentally observed responses of different FG nups to the addition of transport proteins appear rather different (Kapinos et al., 2014; Eisele et al., 2010), closer inspection shows that most FG nups exhibit the same general pattern of behavior. With progressive addition of the transport proteins, the layer height decreases to some extent followed by recovery and eventual swelling. The initiation of the collapse is correlated with the penetration of the transport proteins into the layer. We have recently shown that this behavior is expected on very general grounds for polymer layers infiltrated with nanoparticles (Opferman et al., 2012; 2013). This common behavior is shown in Figure 4, which renders the experimental data from Kapinos et al. (2014), Wagner et al. (2015).
Intuitively, the penetration of an individual transport protein into the layer is determined by the balance between the energetic (enthalpic) gain of creating more contacts with the FG nups and the entropic cost of displacing and crowding the FG nup chains. If the overall free energy change upon insertion of one transport protein into the layer is negative, it will typically penetrate the layer. Otherwise, the penetration is exponentially suppressed, although there still will be some particles in the layer. For a single particle of radius , at low particle concentration, the entropic cost of penetrating an ideal polymer brush layer can be estimated as , where the prefactor depends on the grafting density and the degree of cohesiveness (Halperin et al., 2011; Egorov, 2012; Milchev et al., 2008). On the other hand, a rough estimate for the energetic/enthalpic gain contacts is , where is the number of the interaction sites on the protein and is the energy per contact. Thus, for the entropic repulsion dominates, and one does not expect significant penetration into the layer. In principle, this entropic repulsion from the flexible polymer layer is sufficient for creating the permeability barrier for non-binding molecules. It is crucial to bear in mind that the polymer chains are not static but highly fluctuating entities. It is the entropy of these molecular motions that is responsible for the penetration barrier; thinking of a polymer layer as a static entity with some amount of free space can lead to erroneous conclusions. The barrier could be enhanced by other effects, such as, for instance, inter-chain cohesion.
Thus, at low concentrations the transport proteins penetrate the layer if their attractive interaction with the FG nups is strong enough. However, they do not cause significant conformational changes - essentially occupying the available empty space inside the layer. At higher densities of the transport proteins, or higher interaction strengths, the number of transport proteins in the layer increases, and collective effects start to play a role (Opferman et al., 2013; Halperin and Kröger, 2011; Kim and O’Shaughnessy, 2006). Further addition of the transport proteins causes a cooperative conformational transition of the FG nups leading to either collapse or swelling of the layer, correlated with the accumulation of the transport proteins inside the layer. The magnitude of the collapse and swelling depend on the transport protein size, concentration, interaction strength with the FG nups, the grafting density, and the cohesion strength. Typical height responses to transport protein concentration are shown in Figure 5. One important conclusion of the theory is that the FG nup response to the addition of the transport proteins (for instance, 'swelling' vs. 'collapse') is not an intrinsic property of an FG nup, but can be modulated by the grafting density, transport protein size and the interaction strength. The overall repertoire of predicted behavior as a function of parameters is shown in the 'phase diagrams' in Figures 6 and 7. The model captures the general salient features of the experimental observations. In particular, at the same grafting distance and the cohesion strength, the collapse is more pronounced for a smaller protein such as NTF2, in accord with the experimental observations (Wagner et al., 2015). Similarly, higher grafting distance results in more collapse - because the entropic repulsion of the transport proteins by the chains is lower at lower chain density. Different experimental results can be placed in the different parts of the 'phase diagram', potentially explaining the observed discrepancies. Notably, the theory predicts a high degree of collapse at large grafting distances, in agreement with the experiments of (Lim et al., 2007). Another interesting prediction of the model is that changing the grafting distance affects small and large transport proteins in a different way. For large ones (e.g. Karyopherin), increasing the grafting distance enhances the penetration into the layer because it reduces the repulsive barrier due to lower monomer density in the layer. By contrast, for small transport factors (such as NTF2), for which the repulsion is less important, increasing the grafting distance decreases the penetration into the layer because it reduced the density of available binding sites in the layer.
Can this simplified theory be qualitatively or semi-quantitatively related to the experimental data for realistic values of parameters? Due to the experimental accuracy limitations and large uncertainly in the known values of all the parameters, the fitting of the parameters to the data is not unique and has limited information content. Rather, we focus on clearly distinguishable trends, such as the difference between the behavior of Karyopherin- vs. NTF2 on the Nsp1 layer, reported in Wagner et al. (2015). These measurements, performed on the same FG nup in approximately same range of grafting distances, allow us to examine the effects of the transport protein size and the binding strength, unconfounded by other factors. Karyopherin- and NTF2 have significantly different sizes and binding strengths. NTF2 is ∼ 4 times smaller in volume and has only two binding sites, while Karyopherin- can have up to ten specific FG-binding sites and a significantly larger surface area with potentially much larger number of non-FG interactions with the FG nups (Bayliss et al., 1999; 2000; 2002; Isgro and Schulten, 2005; Liu and Stewart, 2005). We focus on one clearly discernible difference in the experimental behavior: significant penetration of NTF2 starts at higher concentrations but causes stronger compaction of the layer compared to Karyopherin-. Comparison of the theoretical predictions with the experimental results is shown in Figure 8 in the range of concentrations where direct comparison is possible. Because of the relatively large uncertainty in the measurements of the absolute layer height, the grafting distance (see Figure 4) and the binding strengths, the theoretical predications are shown for a range of values approximately corresponding to the experimental ones. The model reproduces the observed differences in the behavior of Karyopherin- and NTF2 on Nsp1 at physically plausible values of the parameters in the regime of its validity. It might also explain why no significant change in the layer height (or very limited swelling) was observed upon addition of transport proteins by other experimental groups (Eisele et al., 2010; 2012): the behavior of both NTF2 and Karyopherin- can be easily shifted into the swelling regime by relatively small changes in the grafting distance, cohesion or interaction strength, the latter of which can be modulated by small changes in the pH or salt and denaturant concentrations. Other patterns revealed by the experimental data shown in Figure 4, such as the dependence of the maximal compaction concentration on the grafting distance, are also qualitatively explained by the theory. Quantitative comparison of these features requires more analysis of the data and more detailed approximations, including the sparse and high transport protein density regimes, and will be presented elsewhere.
The chosen values of and correspond to the 'monomer' size of approximately four amino acids - roughly the size of one FG patch, reflecting the fact that FG nup interaction with the transport proteins requires this particular local sequence of amino acids (Bayliss et al., 2000; 2002); we emphasize that our model does not correspond to a randomized amino acid sequence of the FG nups (Tagliazucchi et al., 2013; Ghavami et al., 2014). The interaction parameter roughly corresponds to the attractive part of the second virial coefficient of the interaction between the transport proteins and the monomers, so that , where is the number of available binding sites on the protein, and is the average binding energy of one site (Doi and Edwards, 1998; Pathria, 1996). Taking into account the possible non-FG interactions, for Karyopherin and for NTF2. This translates to for the chosen parameter values, which is a reasonable estimate for the average energies of the weak hydrophobic and electrostatic interactions in question. These numbers also agree with recent estimates using other methods (Kapinos et al., 2014; Tu et al., 2013). Similarly, assuming that the main contribution to cohesiveness comes from the FG-FG interactions, the value of corresponds to . These values will guide our analysis of bulk solutions of FG nups and transport proteins in the next section. Importantly, the general agreement between the theory and experiment is robust with respect to parameterization choices (see Figure 6—figure supplement 1).
The model also sheds light on another long-standing controversy in the field - the high variability among different experimental groups of the measured affinities of the transport protein binding to the FG nups, with values of the measured dissociation constant ranging from several nanomolars to several micromolars (Bayliss et al., 1999; Kapinos et al., 2014; Eisele et al., 2010; Pyhtila and Rexach, 2003; Tetenbaum-Novatt et al., 2012). As can be seen in Figure 8, due to the collective effects and the conformational changes during binding, the adsorption curves are not well described by the standard Langmiur curve typically used to quantify the interaction in binding assays (Kapinos et al., 2014). An attempt to fit it with one or a combination of Langmiur isotherms would lead to different results depending on the concentration range and the grafting density (Schmidt and Görlich, 2015). The model developed here provides fundamental physical reasons for the discrepancies in experimental measurements.
Phase separation in bulk solutions of FG nups mixed with transport proteins
The fundamental physical considerations underlying the behavior of the surface assemblies of FG nups manifest themselves also in the behavior of mixtures of the FG nups and the transport proteins in bulk solutions. Such mixtures were systematically studied in Schmidt and Görlich (2015) over a wide range of FG nups from different species (recombinantly expressed in bacteria), systematically varying the size of transport protein-cargo complexes. It was found that even in the absence of transport proteins, solutions of Nup98 FG nucleoporin phase separate into a dense phase with a very high protein concentration, in equilibrium with a very dilute solution. Importantly, unlike the previously reported 'gels' of Nup98 and other FG nucleoporins (Hülsmann et al., 2012; Frey and Görlich, 2007), these phases form via an equilibrium phase separation mechanism, allowing comparison with our model. Upon addition of cargo-carrying transport proteins bound to such FG nup solutions, the transport protein-cargo complexes either penetrate the dense phase, or stay predominantly in solution, depending on their size and the interaction strength with the FG nups. These results are summarized in Figure 9. In this section, we show that the observed patterns naturally follow from the minimal model of this paper.
Binary FG nup/buffer solutions
Phase separations of polymer solutions with significantly strong inter-chain cohesiveness are well understood. Phase separation is expected when the cohesiveness is strong enough to make the formation of the dense phase energetically favorable to overcome the loss of the chain entropy in the compact aggregate, (in the mean field approximation for ); note that this condition is independent of any other parameters such as the monomer volume (de Gennes, 1979).FG nup densities in the co-existing phases can be found from the equality of the FG nup chemical potentials and the osmotic pressures, as explained in the Methods and the Appendix.
The model predictions for the phase separation in FG nup solutions agree very well with the experimental data, as shown in Figure 10. The comparison is made for the for the same values of the parameters and as in Figure 3 that describes the surface layers FG nups. The inferred value of the cohesiveness strength for Nup98, known to be more cohesive than Nsp1, is also consistent with the behavior of the grafted layers shown in Figure 3. Taken together, this indicates that the model captures the essential biophysical properties of the FG nup assemblies both in surface grafted geometry and in the bulk solutions.
The model also captures biophysical properties of different FG nups reported in Yamada et al. (2010), where dimensions of individual FG nups were shown to correlate with their hydrophobic to charged amino acid content ratio. The differences between the cohesive collapsed coils - such as Nsp1n (N-terminal domain of Nsp1), Nup100n, Nup116m - and the non-cohesive extended coils, such as Nsp1m (C-terminal domain of Nsp1), are captured by different values of , as shown in Figure 11 and Figure 11—figure supplement 1, in accord with the interpretation of Yamada et al (Yamada et al., 2010). Importantly, the values of arising from the analysis of the data from Yamada et al. (2010) are entirely consistent with those obtained from the analysis of the surface assemblies of Nsp1 in experiments by Lim and coworkers (Kapinos et al., 2014; Wagner et al., 2015) (which used an overlapping segment of the C-terminal domain) and of bulk assemblies of Nup98, Nup110 and 116 in the work by Gorlich and collaborators (Schmidt and Görlich, 2015). Transition from the 'extended coil' regime to the 'collapsed coil' regime upon increase in the chain cohesiveness strength is known in polymer physics as the 'coil-globule' transition and stems from the same physical factors that govern the changes in height of the grafted layers described in Figure 2.
Ternary FG nup/transport proteins/buffer solutions
Phase separation in bulk solutions of FG nups and the transport proteins can be understood through the same theoretical prism. From the theoretical standpoint, at sufficiently high concentrations, such mixtures phase separate into a dilute phase in equilibrium with the dense phase, driven by the inter-chain cohesiveness and the interactions of the FG nups with the transport proteins. The fundamental physics behind this process is the same as in the pure FG nup phase separation: in the dense phase the magnitude of the energetic/enthalpic interactions is maximized at the expense of the entropy.
The overall predictions and conclusions of the theory are presented in the phase diagram in Figure 12. In a range of concentrations, the system phase separates into a dense phase of FG nups mixed with transport proteins, in equilibrium with a dilute solution. The corresponding phases lying on the boundary of the phase separation region (shown in black line) are connected with straight lines. The inset shows the effect of the transport protein-cargo complex size and their interaction strength with the FG nups on the phase separation. Increasing the interaction strength or decreasing the complex size enhances the penetration of the transport proteins and into the dense phase and increases its overall density (green line). On the other hand, sufficiently big or too weakly interacting particles are excluded from the dense phase, as shown in the red phase separation line, where the density of the transport proteins is lower inside the dense FG nup phase than in the outside solution.
The model semi-quantitatively reproduces experimentally observed patterns of phase separation in bulk FG nup solutions mixed with different transport protein-cargo complexes as shown in Figure 13 for the mixtures of Nup98 with Importin--cargo complexes studied in Schmidt and Görlich (2015) (Importin- is the vertebrate homologue of the yeast Karyopherin-). Specifically, both NTF2 and Importin- penetrate the dense FG nup phase, but the complexes of Importin- with either medium (IBB-GFP) or large (IBB-MBP-GFP) cargo are excluded. By contrast, a very large complex of four Importin- molecules with four ZsGreen proteins does partition into the dense FG nup phase, showing that the large size exclusion is not absolute but can be overcome with stronger interactions. Other observations, such as the lack of penetration of the dense phase by the small TEV-mCherry (approx kDa) complex, which does not interact with the FG nups (), are also captured by the model. The parameters for this analysis were based on those inferred from the analysis of phase separation of Nup98 in the absence of the transport proteins: N=167 nm3, . The approximate volumes of different constructs were estimated based on their molecular weights as explained in the Model section and Figure 9. The interaction parameters of Imp- with the FG nups are , . To incorporate the potential direct repulsive interactions between the cargo and the FG nups described in Schmidt and Görlich (2015), we assumed the interaction of the Imp--IBB-GFP and Imp--IBB-MBP-GFP complexes with the FG nups to be slightly weaker (). As a parsimonious estimate, we assumed that the interaction strength of the Imp--ZsGreen tetramer with the FG nups is four times that of an individual Importin-. The correct penetration pattern is produced by the FG nup-Imp interaction strengths lying in the range , in a good agreement with the parameter values inferred from the analysis of the surface layers of Nup98 from Refs. Wagner et al. (2015) and Kapinos et al. (2014); see Figure 13—figure supplement 1. Small differences in the values of between different experimental systems are expected due to the differences in the buffers and the intrinsic sequence differences between the proteins fragments used in each study.
The model also predicts how the penetration of the transport proteins into the dense phase is affected by the cohesiveness , whose role in the NPC selectivity and permeability has been the subject of a vigorous debate (Frey and Görlich, 2007; Schmidt and Görlich, 2015; Lim et al., 2007; 2008; Wagner et al., 2015; Eisele et al., 2010). At low cohesion, insufficient to cause the phase separation of the FG nups alone, addition of the transport proteins can trigger the phase separation and the formation of the mixed dense phase. Further increase in the cohesiveness promotes formation of an even denser phase and enhances Importin partitioning into it, contrary to the naive expectation that cohesiveness always inhibits particle penetration in to the FG Nup medium. Nevertheless, further increase in cohesiveness prevents penetration of the transport proteins into the dense phase, thereby enhancing the permeability barrier as has been suggested before (Hülsmann et al., 2012; Frey and Görlich, 2007). These results will be explored in detail elsewhere.
The model also might explain a recent report that the aggregation of cohesive FG nup segments can be reversed by addition of highly concentrated cell lysate (up to 50 mM) (Hough et al., 2015). Cell lysate has been independently shown to interact with the FG nups non-specifically but strongly enough to compete with the transport proteins for binding sites (Tetenbaum-Novatt et al., 2012). In our model this is captured through medium value of , the interaction strength of the transport proteins with the FG nups. As can be seen from Figure 12 and Figure 12—figure supplement 1, addition of a high concentration of attractive proteins to the phase separated FG nup aggregate can take the system out of the phase separation region, in agreement with the experimental observations.
Discussion
Mechanistic understanding of transport through the Nuclear Pore Complex is hindered by the complexity of the NPC organization and the absence of experimental methods for directly probing FG nup conformations and dynamics during transport. One has to rely on the interpretation of indirect measurements of the FG nup properties in vitro, which are partially incomplete or conflicting. On the theoretical side, the progress is hampered by the lack of a universally accepted comprehensive theory of intrinsically disordered proteins that is able to incorporate all the factors dictating the physico-chemical and nanomechanical properties of the FG nups.
The theoretical model developed here provides a rigorous physical framework for organizing the known phenomenology of the observed FG nup behaviors and suggests ways to reconcile the apparently contradictory experimental findings and models of transport. The model relies on the minimal number of key physical concepts and variables to describe the FG nups and their interactions with the transport proteins, integrating the multitude of molecular details into a small number of variables, thus avoiding the over-fitting problem, inherent to models with multiple unknown parameters. It has been already suggested that the biophysical properties and dimensions of individual FG nups can be, to a large extent, captured in one phenomenological parameter, the hydrophobic to charged amino acid content ratio, paralleling our cohesiveness strength (Yamada et al., 2010). Our model extends these concepts into the regime of multiple chain assemblies interacting with transport proteins. We have found that this limited set of concepts is sufficient to qualitatively, and even semi-quantitatively, explain the experimentally observed patterns and trends. We expect that explicit inclusion of several important factors such as the spatial distributions of ions, discrete nature of the binding sites on the transport proteins, non-linear elasticity of the chains, potential heterogeneity of the FG nup sequence, and more careful modeling of the microscopic structure of the FG nups assemblies will lead to better quantitative agreement. Additional sources of discrepancy are the possible length polydispersity of the FG nups, transport protein aggregation, and inherent biases of the experimental techniques (for instance, in the SPR method, for technical reasons the measurement of the layer height lags behind in time after the measurement of the protein adsorption (Wagner et al., 2015). The theory provides clear experimental predictions for the variation of the properties of FG nup-transport protein assemblies with the experimental conditions such as the grafting density. Such measurements would allow further systematic refinement of the model.
The analysis of this paper shows that both entropic ('brush-like') and entalpic/cohesive ('gel-like') effects naturally cooperate in determining the spatial structures of the assemblies of FG nucleoporins with the transport proteins, and suggests how the 'brush' and the 'gel' concepts can be reconciled. Although the overall qualitative predicted behavior of the surface layers is similar for cohesive and non-cohesive chains alike, quantitative comparison with experiments was only possible by assuming a certain amount of cohesiveness for the studied FG nups. Our analysis indicates that all FG nups likely possess some degree of intra- and inter-chain cohesiveness, although only for some of them it is sufficiently strong to cause aggregation in bulk solutions. The model also shows that the different classes of 'extended' and 'collapsed' FG nups (Yamada et al., 2010) can be accounted for by different values of the cohesiveness . Another puzzling observation of the resilience of the NPC transport with respect to the deletion of large numbers of FG nups might be attributed to the fact that the permeability and the selectivity of the FG nup assemblies are relatively insensitive to the grafting density in a significant range - at lower grafting densities, the neighboring chains simply expand, maintaining selective permeability properties of the layer (Popken et al., 2015). Finally, relatively weak interaction energies of the FG nups among themselves and with the transport proteins inferred from our analysis are consistent with the high local molecular mobility inside the dense aggregates observed experimentally in Schmidt and Görlich (2015), Hough et al. (2015).
The success of the theory relies on the very robust physical mechanisms underlying it. Attractive interactions between long flexible filaments and compact objects, such as folded proteins, cause their surface assemblies to attain more compact conformations at low concentrations of the transport proteins and swollen conformations at higher concentrations. Similarly, in bulk solutions of ungrafted chains, inter-chain interactions and the interactions with transport proteins lead to phase separation and formation of a dense phase at sufficiently high concentrations. These behavior motifs are always expected irrespective of the nature of the interactions on the molecular scale. In this sense, all such systems lie in the same 'universality class' (de Gennes, 1979). The ability of the model to capture the behavior of different FG nups in different geometries under a variety of experimental conditions makes it a useful tool for the development of further, more refined, models. Finally, the analysis of this paper underscores the importance of always considering the presence of the transport factors when thinking about NPC architecture.
The theory also sheds light on the discrepancies in the experimental measurements of the binding affinities of the FG nups to transport proteins. Depending on the experimental procedure, the values of the measured dissociation constants range from several nanomolars to several micromolars. Moreover, some of these measured affinities appear to be inconsistent with the observed transport times in the millisecond range (Yang and Musser, 2006; Yang et al., 2004; Tetenbaum-Novatt and Rout, 2010; Tetenbaum-Novatt et al., 2012; Tu et al., 2013; Denning et al., 2003; Ma, 2010). Our theory shows how these discrepancies might stem from the fundamental statistical thermodynamics of the transport protein-FG nup interaction. First, as shown in this paper, penetration of the transport proteins into an FG nup layer is a cooperative process, not described well by a single Langmiur isotherm typically used in the interpretation of binding assays. Second, penetration of the transport proteins into the layer is determined not only by enthalpic but also entropic effects and therefore the measured effective affinity does not directly reflect the interaction energies. Finally, even the purely enthalpic part of the interaction can vary with the experimental conditions, because the average number of monomers available to bind to a transport protein depends on the monomer concentration, which in turn depends on the layer height and grafting distance (Kapinos et al., 2014; Schoch et al., 2012; Tu et al., 2013; Sethi et al., 2011). All this highlights the fact that the classical characterization of inter-molecular interactions by a single affinity value is not informative for complex multivalent interactions of spatially extended objects such as the FG nups. It is also worth bearing in mind that a -fold difference in the dissociation constant translates into difference in the effective interaction strength - a relatively small difference that can be easily influenced by many factors.
The results of this paper have important implications for the behavior of the Nuclear Pore Complex and design of bio-molecular sorters based on the same principles. One has to be careful in making inference about the NPC properties based on the in vitro results because the detailed features of the actual spatial morphologies of FG nup assemblies in the channel-like geometry of the NPC are likely to differ from flat and bulk geometries. Analysis of this paper establishes the pertinent parameters that constrain possible scenarios and guide future model building. Recent work on coarse grained models of flexible chains in channel geometries allows us to gauge the implications of the findings in the flat and bulk geometries for FG nup morphologies within the NPC (Osmanović et al., 2013a; Coalson et al., 2015; Peleg et al., 2011). Behavior of non-cohesive or weakly cohesive chains in relatively wide pores - wider than the natural height of the planar layer - is expected to be qualitatively similar to that of flat layers. In this case, the chains form a relatively dense layer along the inner surface of the channel, and increase in cohesiveness compacts the layer towards the walls. Addition of transport proteins causes either collapse or swelling of this surface layer, analogously to the behavior in the flat geometry. In the other limit of strongly cohesive chains in relatively narrow channels, the monomers of the chains accumulate at the pore center in a plug-like shape. It is still unknown in what regime the NPC lies, and further computational and experimental work is required. These findings are also interesting in the more general context of intrinsically disordered proteins, where physics concepts are often invoked to organize and explain the experimental observations (Uversky, 2002; van der Lee et al., 2014; Das et al., 2015).
Although the qualitative behavior motifs predicted by the theory and borne out by the experiments are very general, the specific quantitative features such as the exact height or the degree of compaction are rather sensitive to the parameter values, such as the grafting distance, the density and the interaction strength. This has been already noted in other computational theories (Popken et al., 2015; Osmanović et al., 2013a; Tagliazucchi et al., 2013; Gamini et al., 2014). This raises a question - how the Nuclear Pore Complex functioning remains so invariant across species despite the large variations in sizes and spatial organization. Similarly puzzling in this light is the ability of the NPC to maintain its function despite large structural perturbations (Strawn et al., 2004; Hülsmann et al., 2012). One possible solution to this puzzle is that the NPC is exquisitely ”fine tuned” in a sense that it works correctly only when all the molecular details are right - the FG nup sequence and localization, local pH, ionic strength and concentrations of other molecules - and these conditions are maintained by the cellular homeostasis. On the other hand, the results of this model suggest another possibility - that the NPC is 'robust' in a sense that any structure with the approximately right physical properties will function nearly optimally, which could explain the NPC resilience with respect to structural damage and re-arrangements. This provides an interesting example of functional conservation in the absence of sequence conservation. The importance of this question goes beyond the NPC and arises in the discussion of many cellular machines and networks (Alon et al., 1999; Bialek, 2012).
Appendix
Details of the phase separation calculations
Two-component (FG nup/buffer) solutions
In a system where the only solute are the FG Nups without the transport proteins, phase separation occurs for sufficiently high . The volume fractions of the monomers and in the co-existing phases are determined from the equality of the chemical potentials and the osmotic pressures in the separated phases:
where the free energy is given by Equation (4) in the main text with . Numerical solution of these equations produces the phase diagram of Figure 10 in the main text.
Three-component (FG nups/transport proteins/buffer) solutions
In mixed solutions of the FG nups and transport proteins, the chemical potentials of each species in the coexisting phases are equal to each other, as is the osmotic pressure in both phases, similar to the two component case above (de Gennes, 1979):
However, unlike in the two component solution, three Equation (2) in themselves are not sufficient to determine the four unknowns and . In the ternary solutions FG nups/transport proteins/buffer, the densities of both components in the coexisting phases depend also on the initial concentrations and that satisfy the following conservation laws:
where is the fraction of volume occupied by the phase with .
Together, Equation (2-3) provide five equations for five variables , and . Their solution results in the phase diagrams of Figures 11–12 in the main text.
References
-
Adsorption of chain molecules with a polar head a scaling descriptionJournal De Physique 38:983–987.https://doi.org/10.1051/jphys:01977003808098300
-
GLFG and FxFG nucleoporins bind to overlapping sites on importin-betaJournal of Biological Chemistry 277:50597–50606.https://doi.org/10.1074/jbc.M209037200
-
Interaction between NTF2 and xFxFG-containing nucleoporins is required to mediate nuclear import of RanGDPJournal of Molecular Biology 293:579–593.https://doi.org/10.1006/jmbi.1999.3166
-
Polyelectrostatic interactions of disordered ligands suggest a physical basis for ultrasensitivityProceedings of the National Academy of Sciences of the United States of America 104:9650–9655.https://doi.org/10.1073/pnas.0702580104
-
Capturing directed molecular motion in the nuclear pore complex of live cellsProceedings of the National Academy of Sciences of the United States of America 109:9863–9868.https://doi.org/10.1073/pnas.1200486109
-
Domain topology of nucleoporin Nup98 within the nuclear pore complexJournal of Structural Biology 177:81–89.https://doi.org/10.1016/j.jsb.2011.11.004
-
A polymer-brush-based nanovalve controlled by nanoparticle additives: design principlesThe Journal of Physical Chemistry. B 119:11858–11866.https://doi.org/10.1021/acs.jpcb.5b02623
-
Charge as a selection criterion for translocation through the nuclear pore complexPLoS Computational Biology 6:e1000747.https://doi.org/10.1371/journal.pcbi.1000747
-
Autonomy and robustness of translocation through the nuclear pore complex: a single-molecule studyThe Journal of Cell Biology 183:77–86.https://doi.org/10.1083/jcb.200806173
-
Relating sequence encoded information to form and function of intrinsically disordered proteinsCurrent Opinion in Structural Biology 32:102–112.https://doi.org/10.1016/j.sbi.2015.03.008
-
Conformations of polymers attached to an interfaceMacromolecules 13:1069–1075.https://doi.org/10.1021/ma60077a009
-
Disorder in the nuclear pore complex: The FG repeat regions of nucleoporins are natively unfoldedProceedings of the National Academy of Sciences of the United States of America 100:2450–2455.https://doi.org/10.1073/pnas.0437902100
-
Nuclear pore complexes and nucleocytoplasmic transport: From structure to function to diseaseInternational Review of Cell and Molecular Biology 320:171–233.https://doi.org/10.1016/bs.ircmb.2015.07.010
-
Insertion of nanoparticles into polymer brush under variable solvent conditionsThe Journal of Chemical Physics 137:134905.https://doi.org/10.1063/1.4757017
-
The molecular mechanism of translocation through the nuclear pore complex is highly conservedJournal of Cell Science 115:2997–3005.
-
The location of the transport gate in the nuclear pore complexJournal of Cell Science 110:3065–3070.
-
A jumbo problem: mapping the structure and functions of the nuclear pore complexCurrent Opinion in Cell Biology 24:92–99.https://doi.org/10.1016/j.ceb.2011.12.013
-
Assembly of Nsp1 nucleoporins provides insight into nuclear pore complex gatingPLoS Computational Biology 10:e1003488.https://doi.org/10.1371/journal.pcbi.1003488
-
Colloid-brush interactions: The effect of solvent qualityMacromolecules 44:3622–3638.https://doi.org/10.1021/ma200068d
-
Collapse of thermoresponsive brushes and the tuning of protein adsorptionMacromolecules 44:6986–7005.https://doi.org/10.1021/ma201006h
-
NPC mimics: probing the mechanism of nucleocytoplasmic transportMethods in Cell Biology 122:379.https://doi.org/10.1016/B978-0-12-417160-2.00017-5
-
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
-
Single-molecule transport across an individual biomimetic nuclear pore complexNature Nanotechnology 6:433–438.https://doi.org/10.1038/nnano.2011.88
-
Nucleoporins as components of the nuclear pore complex core structure and Tpr as the architectural element of the nuclear basketMolecular Biology of the Cell 15:4261–4277.https://doi.org/10.1091/mbc.E04-03-0165
-
Karyopherin-independent spontaneous transport of amphiphilic proteins through the nuclear poreJournal of Cell Science 125:4979–4984.https://doi.org/10.1242/jcs.109520
-
Polymer brushes in mixed solvents: chromatography and collapseMacromolecules 25:6693–6695.https://doi.org/10.1021/ma00050a045
-
A simplified representation of protein conformations for rapid simulation of protein foldingJournal of Molecular Biology 104:59–107.https://doi.org/10.1016/0022-2836(76)90004-8
-
Towards reconciling structure and function in the nuclear pore complexHistochemistry and Cell Biology 129:105–116.https://doi.org/10.1007/s00418-007-0371-x
-
Flexible phenylalanine-glycine nucleoporins as entropic barriers to nucleocytoplasmic transportProceedings of the National Academy of Sciences of the United States of America 103:9512–9517.https://doi.org/10.1073/pnas.0603521103
-
Self-regulated viscous channel in the nuclear pore complexProceedings of the National Academy of Sciences of the United States of America 109:7326–7331.https://doi.org/10.1073/pnas.1201724109
-
Three-dimensional distribution of transient interactions in the nuclear pore complex obtained from single-molecule snapshotsProceedings of the National Academy of Sciences of the United States of America 107:7305–7310.https://doi.org/10.1073/pnas.0908269107
-
Single molecule study of the intrinsically disordered FG-repeat nucleoporin 153Biophysical Journal 101:1710–1719.https://doi.org/10.1016/j.bpj.2011.08.025
-
Theory of the grafted polymer brushMacromolecules 21:2610–2619.https://doi.org/10.1021/ma00186a051
-
Simulations of nuclear pore transport yield mechanistic insights and quantitative predictionsProceedings of the National Academy of Sciences of the United States of America 108:E351.https://doi.org/10.1073/pnas.1104521108
-
Solvent quality effects on scaling behavior of poly(methyl methacrylate) brushes in the moderate- and high-density regimesLangmuir : The ACS Journal of Surfaces and Colloids 27:3698–3702.https://doi.org/10.1021/la2002139
-
Characterisation of the passive permeability barrier of nuclear pore complexesThe EMBO Journal 28:2541–2553.https://doi.org/10.1038/emboj.2009.200
-
Computer simulation of random packing of hard spheresPowder Technology 72:149–155.https://doi.org/10.1016/0032-5910(92)88021-9
-
Morphology of polymer brushes infiltrated by attractive nanoinclusions of various sizesLangmuir : The ACS Journal of Surfaces and Colloids 29:8584–8591.https://doi.org/10.1021/la4013922
-
Bistable collective behavior of polymers tethered in a nanoporePhysical Review E 85:061917.https://doi.org/10.1103/PhysRevE.85.061917
-
Physical modelling of the nuclear pore complexSoft Matter 9:10442–10451.https://doi.org/10.1039/c3sm50722j
-
Individual binding pockets of importin- for FG-nucleoporins have different binding properties and different sensitivities to RanGTPProceedings of the National Academy of Sciences of the United States of America 105:16101–16106.https://doi.org/10.1073/pnas.0802647105
-
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
-
A gradient of affinity for the karyopherin Kap95p along the yeast nuclear pore complexThe Journal of Biological Chemistry 278:42699–42709.https://doi.org/10.1074/jbc.M307135200
-
Accurate calculation of the density of proteinsActa Crystallographica Section D Biological Crystallography 56:791–794.https://doi.org/10.1107/S090744490000679X
-
Virtual gating and nuclear transport: the hole pictureTrends in Cell Biology 13:622–628.https://doi.org/10.1016/j.tcb.2003.10.007
-
Phase transition behavior of the isolated polymer chainMacromolecules 12:980–988.https://doi.org/10.1021/ma60071a040
-
Nuclear transport receptor binding avidity triggers a self-healing collapse transition in FG-nucleoporin molecular brushesProceedings of the National Academy of Sciences of the United States of America 109:16911–16916.https://doi.org/10.1073/pnas.1208440109
-
Coil-globule transition in the denatured state of a small proteinProceedings of the National Academy of Sciences of the United States of America 103:11539–11543.https://doi.org/10.1073/pnas.0601395103
-
Ring cycle for dilating and constricting the nuclear poreProceedings of the National Academy of Sciences of the United states of America 110:5858–5863.https://doi.org/10.1073/pnas.1302655110
-
Molecular mechanism of the nuclear protein import cycleNature Reviews. Molecular Cell Biology 8:195–208.https://doi.org/10.1038/nrm2114
-
Minimal nuclear pore complexes define FG repeat domains essential for transportNature Cell Biology 6:197–206.https://doi.org/10.1038/ncb1097
-
Emerging applications of stimuli-responsive polymer materialsNature Materials 9:101–113.https://doi.org/10.1038/nmat2614
-
Effect of charge, hydrophobicity, and sequence of nucleoporins on the translocation of model particles through the nuclear pore complexProceedings of the National Academy of Sciences of the United States of America 110:3363–3368.https://doi.org/10.1073/pnas.1212909110
-
Biophysical constraints for protein structure predictionJournal of Proteome Research 2:37–42.https://doi.org/10.1021/pr025552q
-
Nucleocytoplasmic transport: a role for nonspecific competition in karyopherin-nucleoporin interactionsMolecular & Cellular Proteomics 11:31–46.https://doi.org/10.1074/mcp.M111.013656
-
The mechanism of nucleocytoplasmic transport through the nuclear pore complexCold Spring Harbor Symposia on Quantitative Biology 75:567–584.https://doi.org/10.1101/sqb.2010.75.033
-
Natively unfolded proteins: a point where biology waits for physicsProtein Science : A Publication of the Protein Society 11:739–756.https://doi.org/10.1110/ps.4210102
-
Classification of intrinsically disordered regions and proteinsChemical Reviews 114:6589–6631.https://doi.org/10.1021/cr400525m
-
The nuclear pore complex and nuclear transportCold Spring Harbor Perspectives in Biology 2:a000562.https://doi.org/10.1101/cshperspect.a000562
-
A bimodal distribution of two distinct categories of intrinsically disordered structures with separate functions in FG nucleoporinsMolecular & Cellular Proteomics 9:2205–2224.https://doi.org/10.1074/mcp.M000035-MCP201
-
Imaging of single-molecule translocation through nuclear pore complexesProceedings of the National Academy of Sciences of the United States of America 101:12887–12892.https://doi.org/10.1073/pnas.0403675101
-
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
-
Protein volume in solutionProgress in Biophysics and Molecular Biology 24:107–123.https://doi.org/10.1016/0079-6107(72)90005-3
-
Screening effects in a polyelectrolyte brush: Self-consistent-field theoryMacromolecules 33:4945–4953.https://doi.org/10.1021/ma990187i
-
Efficiency, selectivity, and robustness of nucleocytoplasmic transportPLoS Computational Biology 3:e125.https://doi.org/10.1371/journal.pcbi.0030125
-
Enhancement of transport selectivity through nano-channels by non-specific competitionPLoS Computational Biology 6:e1000804.https://doi.org/10.1371/journal.pcbi.1000804
-
Thermodynamics and structure of self-assembled networksPhysical Review E 66:051107.https://doi.org/10.1103/PhysRevE.66.051107
Article and author information
Author details
Funding
National Science and Engineering Research Council of Canada (Discovery Grant RGPIN 4025901)
- Andrei Vovk
- Anton Zilman
- Chad Gu
Swiss National Science Foundation (Grant 31003A_146614)
- Larisa E Kapinos
- Roderick YH Lim
National Science Foundation (CHE-1464551)
- Rob D Coalson
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
AZ acknowledges the support of the National Science and Engineering Research Council of Canada (NSERC) and Compute Canada. RDC acknowledges support form NSF grant CHE-1464551. RYHL is supported by the Swiss National Science foundation through Grant 31003A-146614. We thank Dr. B Hoogenboom, Dr. S Musser, Dr. P Onck, Dr. R Richter and Dr. R Schoch for helpful discussions.
Copyright
© 2016, Vovk 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
-
- 2,705
- views
-
- 580
- downloads
-
- 71
- 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.