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.001
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.002
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.
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).
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.
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.
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).
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.
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.
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.
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.
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.
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).
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:
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 .
Adsorption of chain molecules with a polar head a scaling descriptionJournal De Physique 38:983–987.https://doi.org/10.1051/jphys:01977003808098300
Theory of polyelectrolyte solutionsAdvances in Chemical Physics 94:66.
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
Nanoscale stiffness topography reveals structure and mechanics of the transport barrier in intact nuclear pore complexesNature Nanotechnology, 10, 10.1038/nnano.2014.262.
Biophysics: Searching for PrinciplesPrinceton University Press.
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
Scaling Concepts in Polymer ScienceCornell University Press.
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
The Theory of Polymer DynamicsClarendon Press.
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
Principles of Polymer ChemistryCornell University Press.
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
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
Thermal PhysicsNew York: WA Benjamin.
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
Ph.D. ThesisUniversity of Pittsburgh.
Bistable collective behavior of polymers tethered in a nanoporePhysical Review E 85:061917.https://doi.org/10.1103/PhysRevE.85.061917
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
Statistical Mechanics (2nd Edition)Elsevier.
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
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
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
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
Arup K ChakrabortyReviewing Editor; Massachusetts Institute of Technology, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your work entitled "Simple biophysical model explains the conformational transitions of the unfolded proteins of the Nuclear Pore Complex" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and John Kuriyan as the Senior Editor.
The reviewers have discussed the reviews with one another and the Reviewing editor has drafted this decision to help you prepare a revised submission.
This paper develops a self-consistent-field theory for a model of transport through nuclear pores. The model is properly formulated based on well-established and long-standing concepts in polymer physics. Though several theoretical models have been previously published attempting to describe the biophysical nature of nuclear transport, the model presented here has the considerable attraction that it is less inclined than these other models to try to pitch a particular mechanism – rather, it sets up a computational framework that can be adapted and added to upon additional data, giving researchers a computational framework for a more rigorous analysis of data. To demonstrate the effectiveness of the model, it is used to perform a new set of analyses of previously published work from two of the co-authors, Kapinos and Lim. The model re-capitulates certain previously puzzling behaviors seen for FG nups and provides a biophysical hypothesis for these behaviors. This model could therefore be a valuable new tool. However, our enthusiasm for the paper is considerably diluted by two issues.
Firstly, the manuscript oversells itself as having resolved the controversies, and furthermore, the comparisons with experiment are qualitative. Qualitative comparisons are fine, but rephrasing some of the text so that it accurately reflects the accomplishments is essential.
Secondly, and more importantly, to establish itself as a useful computational platform, additional work is necessary. In the paper, the model is fit to existing data from some of the co-authors. To make a compelling case that the model indeed captures the essential features of nuclear pore transport, it is necessary to show that the model is predictive – that is, take published data from several groups (not the co-authors') that was explicitly not included in building the model, and without re-adjusting the parameters used to fit the authors' data, show that the model describes these additional data. Without demonstrating this capability, the model cannot be considered predictive or explanatory. We suggest that you consider describing the data presented in the following references: 1) Strawn et al., 2004: can the model describe the observation here that so much of the FG mass can be deleted without completely disrupting transport? 2) Yamada et al.: can the model explain why might there are two classes of cohesiveness as suggested in this paper? 3) Popken et al., 2015: Compare the coarse grained simulations and data in this paper with those described in the manuscript. 4) Patel et al., 2007: which also contains interesting data. Of course, you could choose other data sets that you believe are illustrative.
Some additional detailed points that need to be addressed are listed below.
1) The model completely ignores sequence information and electrostatic interactions. The problem with this is that there is clear evidence from bioinformatics studies (see ref. 68) and theory (see ref. 50) that clearly show the need to consider the role of electrostatics. These studies seem to consider similar screening conditions as that in the present manuscript. Moreover, ref. 50 (a theoretical study) shows that smearing the sequence leads to qualitatively different behavior. Therefore, these points need to be mentioned explicitly as refs. 68 and 50 suggest that it may be difficult to describe the NPC using minimal models adequately, as is claimed in the last sentence of the manuscript.
2) If the minimal model provides an excellent description of nuclear pores, why can synthetic systems not recapitulate the behavior of nuclear pores? The manuscript should emphasize that the model describes some universal features, and more details will have to be added upon considering additional data – this point might become more vivid or refuted by addressing the second major point made above.
3) Figure 2 presents well known results, as is noted in the manuscript, but the original work of Alexander is not cited. Also use of the term "cross-link" here is misleading since this is just a model van der Waals interaction.
4) Figure 3 shows the two limits that a brush can have from maximal repulsion (stretching) to completely compact. The fact that the experimental observations fall between these two asymptotic limits is expected. But, why is the same nup at different grafting densities suggesting very different values of χ?
5) It appears that the model assumes macroscopic uniform planar layers grafted with one type of FG repeat. Does the model also make accurate predictions for geometries that more accurately resemble the NPC i.e. nanoscopic pores, or mixtures of different FG repeat types?
[Editors' note: further revisions were requested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled "Simple biophysics underpins collective conformations of intrinsically disordered proteins of the Nuclear Pore Complex" for further consideration at eLife. Your revised article has been favorably evaluated by John Kuriyan as Senior editor, a Reviewing editor, and one reviewer. The revised manuscript is much improved, and most of the points we raised previously have been fully addressed. If the remaining points detailed below can be addressed, we think that there is a very good chance that the paper will be accepted for publication.
1) In the Abstract, we still feel "These results reconcile some of the outstanding controversies…" is strong, given that this work does not concretely reconcile, but rather suggest solutions, to the current controversies. Perhaps "These results address some of the outstanding controversies…" would be better.
2) In the Discussion, the sentence "The suggests how the ‘brush’ and… effects" is grammatically incorrect and unclear – can you please re-write?
3) On use of data by others: You addressed our comments by including a work that focuses only on Nup98, which represents only one type of FG, one that can form gels in vitro, but are not necessarily representative of others. Even if an experiment is not possible, could you still replicate computational findings by others, if not at a quantitative, at a qualitative level? In the absence of an independent test, again, we would recommend a more cautious approach on claims that the current work resolves different findings. Instead, perhaps emphasize rather that this work presents a tool of utility to future studies.
4) On the omitting of some aspects of FGs, etc.: A) Sequence: the explanation about the randomized AA sequence is much improved. However, it would be nice to have a demonstration of the thesis in a case with a simple heterogeneity in sequence, such as in Nsp1 (i.e. difference in property in N-terminal and C-terminal regions). Such a demonstration would also reinforce the paper by possibly recapitulating the computational finding in Yamada et al. (difference in expected structure on N and C-termini). B) On the geometry of FG grafts: Although further computational investigation of this topic may be beyond the scope of this paper, it is worth discussing the limitations of the current presentation of the simple monolayer geometry and the possible differences/changes in FG behavior/conformation (and/or lack thereof) that could arise from change in geometry (e.g. grafting them in a pore lumen, etc.).
5) In the Abstract, there is a sentence: "NPC transport relies on the conformational transitions of assembly…". It is true that the conformation is known to change in in vitro experimental setups, but is there any evidence in NPC-context that supports this statement (i.e. conformational transitions)?
6) Figure 1 is rather crudely rendered, and does not accurately represent the positions of the FG nups indicated. Krull et al., 2004 can be used for Nups98 and 153, while recent crystallographic work (Chug et al., 2015) indicates the position of Nup62 to be symmetric and immediately adjacent to Nup93 (Krull et al., 2004). Also, the lower two panels as they stand don't really help readers understand what you call the "experimental situation." This figure should be polished further.
7) For Figure 2, it can be guessed from the context that different curves correspond to different χs, but it's better to specify which curve correspond to what value (or add some direction like an arrow with "increasing χ," etc.)
8) For Figure 3, please consider omitting "experimental evidence of cohesiveness" from the figure title. Strictly speaking, it should go something like "polymer model based on cohesion theory fits well with the data" and it's not really evidence for cohesiveness.https://doi.org/10.7554/eLife.10785.024
- Andrei Vovk
- Anton Zilman
- Chad Gu
- Larisa E Kapinos
- Roderick YH Lim
- Rob D Coalson
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
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.
- Arup K Chakraborty, Reviewing Editor, Massachusetts Institute of Technology, United States
- Received: August 11, 2015
- Accepted: March 24, 2016
- Version of Record published: May 20, 2016 (version 1)
© 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.