Jun 28, 2023
Visualisierung der ungeordneten nuklearen Transportmaschinerie vor Ort
Nature Band 617, Seiten 162–169 (2023) Diesen Artikel zitieren 18.000 Zugriffe 3 Zitate 214 altmetrische Metrikdetails Der etwa 120 MDa große Säugetier-Kernporenkomplex (NPC) fungiert als Gatekeeper für
Nature volume 617, pages 162–169 (2023)Cite this article
18k Accesses
3 Citations
214 Altmetric
Metrics details
The approximately 120 MDa mammalian nuclear pore complex (NPC) acts as a gatekeeper for the transport between the nucleus and cytosol1. The central channel of the NPC is filled with hundreds of intrinsically disordered proteins (IDPs) called FG-nucleoporins (FG-NUPs)2,3. Although the structure of the NPC scaffold has been resolved in remarkable detail, the actual transport machinery built up by FG-NUPs—about 50 MDa—is depicted as an approximately 60-nm hole in even highly resolved tomograms and/or structures computed with artificial intelligence4,5,6,7,8,9,10,11. Here we directly probed conformations of the vital FG-NUP98 inside NPCs in live cells and in permeabilized cells with an intact transport machinery by using a synthetic biology-enabled site-specific small-molecule labelling approach paired with highly time-resolved fluorescence microscopy. Single permeabilized cell measurements of the distance distribution of FG-NUP98 segments combined with coarse-grained molecular simulations of the NPC allowed us to map the uncharted molecular environment inside the nanosized transport channel. We determined that the channel provides—in the terminology of the Flory polymer theory12—a ‘good solvent’ environment. This enables the FG domain to adopt expanded conformations and thus control transport between the nucleus and cytoplasm. With more than 30% of the proteome being formed from IDPs, our study opens a window into resolving disorder–function relationships of IDPs in situ, which are important in various processes, such as cellular signalling, phase separation, ageing and viral entry.
IDPs are flexible, dynamic macromolecules that lack a fixed tertiary structure and can adopt a range of conformations to perform various functions across the cell. IDPs are highly relevant for human physiology and have central roles, among others, in neurodegenerative ageing diseases and cancer. IDPs are also key players in phase separation and are involved in the formation of biomolecular condensates13,14,15,16,17,18,19,20,21. In the nanosized NPC, which has a total molecular weight of approximately 120 MDa in mammals, there are hundreds of IDPs enriched in phenylalanine (F) and glycine (G) residues, known as FG-NUPs1. The FG-NUPs form a permeability barrier in the central channel of the NPC, which regulates nucleocytoplasmic transport by restricting the passage of large cargo unless it presents a nuclear localization sequence or a nuclear export sequence2,3. Nuclear transport receptors can specifically recognize these sequences and efficiently shuttle the cargo through the barrier. With recent advances in cryo-electron tomography, crystallography, proteomics and artificial intelligence (AI)-based structure prediction, approximately 70 MDa of the NPC scaffold enclosing the central channel has been resolved with near-atomic resolution4,5,6,7,8,9,10,11. However, signals from the highly dynamic FG-NUPs are by and large not accessible to those structural biology techniques, and the actual transport machinery inside the central channel—another approximately 50 MDa—is not captured, leaving an approximately 60-nm hole in the centre of the scaffold structure. Consequently, the protein conformational state inside the NPC remains elusive, which has led to several partially conflicting hypotheses for the morphologies of the FG domains in their functional state22,23,24,25,26,27,28. With approximately 30% of the entire eukaryotic proteome being intrinsically disordered, the problem that the conformational state is not easily studied in cells extends far beyond NPC biology. Besides magnetic resonance and scattering techniques13,14, single-molecule fluorescence of purified and labelled proteins has become a powerful tool for probing the conformations of proteins in solution; advanced studies have even shown that this is possible in cells if such probes are microinjected29,30,31. However, the NPC is assembled only in late mitosis and during nuclear growth in interphase32, and its labelling thus requires genetic encoding. Established fluorescent protein-based technologies such as GFP or self-labelling protein tags such as SNAP-tag33, however, do not readily enable the extraction of multiple distance distributions for the same protein, owing to the sheer size of the fluorescent label and the inherently limited freedom of labelling.
In this study, we developed a method to probe distance distributions of FG-NUPs inside the NPCs by combining fluorescence lifetime imaging of fluorescence resonance energy transfer (FLIM–FRET) with a site-specific synthetic biology approach. We show that the methods deliver quantitative results when using permeabilized cells with functional transport machinery, and offer sound agreement with qualitative measurements from live cells. We focused on NUP98 because it is the essential constituent of the NPC permeability barrier and is accessible to this technology34,35,36. By measuring the distance distribution for 18 labelled chain segments of NUP98 in the NPC using FLIM–FRET, we showed that the FG domain is exposed to—in the terminology of the Flory polymer model12—‘good solvent’ conditions inside the NPC. This enables the protein to adopt much more extended conformations in the functional state than in the highly collapsed state of single chains in solution at ‘poor solvent’ conditions. We combined our residue-specific measurements with coarse-grained molecular dynamics (MD) at residue resolution. This enabled us to integrate the distance distribution data and the recently solved scaffold structure7 into a molecular picture of FG-NUP distribution and motion in the central channel of a functional NPC.
High-precision fluorescence measurements of the conformation of FG-NUPs in their functional state require the introduction of labelling tags with minimal linkage errors and minimal disruption of the structures and functions of the labelled proteins. To this end, we performed site-specific labelling with a non-canonical amino acid (ncAA) using genetic code expansion (GCE)37. We used the pyrrolysine orthogonal tRNA–synthetase suppressor pair to reassign the amber stop codon (TAG) to incorporate the ncAA trans-cyclooct-2-en-l-lysine (TCO*A) at that site. The chemical functionality of this ncAA residue can then be reacted with an organic fluorophore containing a tetrazine moiety to undergo an inverse-electron-demand Diels–Alder reaction (click chemistry38). Thus, the dye is stably attached to the protein via a small chemical linker, causing minimal disruption to the protein structure and its function. One potential downside of this technique is that it is not mRNA-specific, leading to background labelling of untargeted proteins with their naturally occurring stop codons being suppressed. To circumvent this problem, we utilized our recently developed film-like, synthetic orthogonally translating organelles (OTOs) to form a distinct protein translational machinery on the outer mitochondrial membrane surface39. These organelles exclusively reassigned two amber codons for the target FG-NUP and incorporated TCO*A at the two specified sites with high selectivity, ensuring minimal interference with endogenous protein translation and negligible background staining (see Fig. 1a and Extended Data Fig. 1 for the improved contrast when comparing OTO technology with conventional GCE). The incorporated ncAAs were reacted with a mixture of donor and acceptor dyes for FRET measurements. This results in FRET species mixed with donor-only and acceptor-only species. Therefore, the chosen measurement method must be able to distinguish the FRET species from the other two. Furthermore, quantitative high-precision FRET measurements have high requirements on the properties of the FRET dye pair, such as photostability, Förster radius, monoexponential decay of fluorescence lifetime and fluorescence emission clearly distinguishable from background.
a, Schematic of site-specific labelling of target NUP98 inside the NPC. The genetic code was exclusively expanded for target NUP98. ncAAs were introduced into NUP98 at specific sites by synthetic orthogonal translating organelle-enabled genetic code expansion (OTO-GCE). Tetrazine-modified dye molecules were added to the cells and reacted with ncAAs (click chemistry). b, Passive exclusion (70-kDa dextran was excluded) and facilitated assay, or active transport assay (IBB–MBP–GFP supplied with transport mixture was imported) showing that the NPCs were functional with site-specific labelled NUP98 in permeabilized COS-7 cells (with the plasma membrane selectively permeabilized with low-dosage digitonin). Scale bars, 20 μm. c, Purified NUP98 FG domain was phase-separated in vitro by rapidly diluting a denatured highly concentrated stock solution into physiological buffer. The permeability of the droplet-like condensates formed was measured rapidly where they still obeyed liquid-like characteristics26. The droplets recapitulated the function of the permeability barrier, as shown in the passive exclusion and facilitated transport assays. Scale bars, 5 μm. For b,c, n = 3 experiments were repeated independently with the same conclusion.
We first performed live-cell labelling of NUP98. We co-transfected COS-7 cells with plasmids encoding NUP98 with an amber codon (NUP98A221TAG) and the OTO-GCE system, and treated the cells with various cell-permeable dyes, including JF549-tetrazine, JF646-tetrazine, silicon rhodamine-tetrazine and TAMRA-tetrazine (see Methods). We could hardly identify a clear nuclear envelope due to the low signal-to-noise ratio (from nonspecific dye sticking), unless we highly overexpressed NUP98 (Extended Data Fig. 2).
We also performed labelling with cell-impermeable dyes, which are more hydrophilic and show less nonspecific sticking. To deliver the dyes to the nuclear envelope, we permeabilized the cells with low-dosage digitonin, in which the plasma membrane was selectively permeabilized while leaving the nuclear membrane and endoplasmic reticulum intact40. We chose a FRET dye pair of AZDye594-tetrazine (orange) and LD655-tetrazine (red) for their exceptional photostability and suitable Förster radius (R0 ~ 7.7 nm; see Methods), and because they are spectrally distinct from typically green cellular autofluorescence. To verify that our procedure did not perturb the functionality of the NPC, we performed a transport assay on permeabilized COS-7 cells labelled with LD655-tetrazine41 (Fig. 1b and Extended Data Fig. 3). A large inert cargo (70-kDa FITC-labelled dextran) was excluded from the nucleus, suggesting that both the permeability barrier and the nuclear envelope were intact, whereas IBB–MBP–GFP (a triple fusion of the importin-β-binding domain (which is recognized by the import receptor importin-β) with maltose-binding protein (to make the construct bigger) and green fluorescent protein) supplied with transport mixtures (which contained the nuclear transport receptor importin-β; see Methods) was actively imported into the nucleus, demonstrating that NPCs containing labelled NUP98 in the permeabilized cells were fully functional. Therefore, we performed the following FRET measurements of the NUP98 FG domain with AZDye594 and LD655, unless otherwise stated.
FLIM–FRET is exquisitely sensitive to the spatial distance between pairs of donor and acceptor fluorescent dyes. As one of the few FRET methods that provide quantitative FRET information without a priori knowledge about the actual labelling stoichiometry, FLIM–FRET is also independent of fluorophore concentration and excitation intensity, providing an ideal tool for probing the dimension of FG-NUPs in the cellular milieu42. By quantifying the decrease in donor fluorescence lifetime when it undergoes FRET coupling with an acceptor molecule, we determined the average spatial distances between the dyes. The robustness of the method was further enhanced by combining FLIM–FRET measurements with acceptor photobleaching, from which we could measure the fluorescence lifetime of the donor-only population and the cellular background with high precision. In the Supplementary Text (equations (15–17)) and Supplementary Figs. 1 and 2, we further detail how the combination of FLIM–FRET and acceptor photobleaching increases robustness of the analysis.
To make sure that the spectral properties of the chosen dye pair are independent of the labelling site, we labelled 19 different positions of the NUP98 FG domain with only the donor dye. Extended Data Fig. 4 shows that the donor fluorescence lifetimes and anisotropies were similar for each labelling site, demonstrating that the fluorophores experience a similar microenvironment within the NPC and that the dyes are mobile enough to allow for high-precision FRET distance measurements. Another concern is that intermolecular FRET could occur, that is, FRET between different NUP98 molecules in the same NPC, owing to the high FG-NUP density in the nanosized central channel. To verify that no intermolecular FRET was measured, we transfected COS-7 cells with a plasmid encoding NUP98 with only a single amber codon (NUP98A221TAG) and treated them with a mixture of donor and acceptor dyes, such that each modified copy of NUP98 could only be singly labelled with either donor or acceptor dye. We selected nuclear envelopes with an acceptor intensity per pixel (excited by a 660-nm laser) below a determined threshold at which to measure FRET to eliminate risks of intermolecular FRET (see Methods; Fig. 2 and Supplementary Fig. 3). We then measured the average fluorescence lifetime of the donor dye before and after acceptor photobleaching. If FRET occurs, the donor intensity and the fluorescence lifetime will increase when the acceptor is bleached selectively by a high-power laser43. We found the donor signal to be unchanged in the permeabilized cells, thus excluding the existence of intermolecular FRET (Fig. 2b). However, when we performed the same assay on the living cells labelled with JF549 and JF646, where we could identify a clear nuclear rim, we noticed an increase of donor lifetime after acceptor photobleaching (Fig. 2c). This indicates the presence of intermolecular FRET under live-cell labelling conditions, where due to the poorer signal to noise, nuclear rims could only be identified when NUP98 was highly overexpressed. Intermolecular FRET does not permit to measure at the same time quantitative intramolecular FRET, and thus we continued with permeabilized cell measurements for the rest of the work.
a, Schematic of the FLIM–FRET analysis pipeline. Different chain segments of the NUP98 FG domain were labelled with a FRET dye pair, and the donor fluorescence intensity was measured on a cell-by-cell basis. Each nuclear rim was selected as a region of interest, and the measured donor fluorescence intensity profiles before and after acceptor photobleaching were extracted and analysed. b,c, Acceptor photobleaching assays were performed for a single-amber-mutated sample (NUP98A221TAG) in permeabilized cells labelled with a AZDye594–LD655 mixture (b) and living cells labelled with a JF549–JF646 mixture (c). In b, the average fluorescence lifetime of the donor dye did not change before and after acceptor photobleaching, indicating the absence of intermolecular FRET. In c, the average fluorescence lifetime of the donor dye changed before and after acceptor photobleaching changed, indicating that intermolecular FRET was detected in highly overexpressing living cells with a bright nuclear rim. d, Acceptor photobleaching assay was performed for a double-amber-mutated sample (NUP98A221TAG–A312TAG) labelled with the AZDye594–LD655 mixture. The average fluorescence lifetime of the donor dye changed before and after acceptor photobleaching, validating the presence of intramolecular FRET. For b–d, n = 5 experiments were repeated independently with the same conclusion; scale bars, 5 μm. e, Fluorescence decay profile before photobleaching was subtracted from the one after photobleaching for the 18 different chain segments of the NUP98 FG domain. Each profile represents an averaged result of approximately 100 cells. The higher peak shows a greater difference in the intensity profiles, indicating higher FRET efficiency and smaller inter-residue distance. f, Phasor plot showing donor lifetimes of the measured 18 chain segments on a single-cell basis (here approximately 2,000 cells in total), in which each point represents the fluorescence decay of one nuclear rim. The left-shifted points represent longer lifetimes.
Source Data
Next, we conducted the acceptor photobleaching assay on the permeabilized cells expressing double-amber-mutated NUP98 (NUP98A221TAG–A312TAG; Fig. 2d). We observed an increase in both the donor intensity on the nuclear rim and its average fluorescence lifetime when the acceptor was photobleached, confirming the occurrence of FRET. With these results, we validated the experimental setup as sufficiently sensitive for measuring FRET between two labelled sites of NUP98 in the permeabilized cells with a functional transport machinery.
We then created a series of NUP98 mutants to form a set of chain segments of different length between the labels for FLIM–FRET measurements. We chose A221TAG as the reference site and kept it constant while varying the second site along the FG domain (Fig. 2e). We define Nres as the number of amino acid residues between the two labelled sites and RE as the root-mean-square inter-residue distance between the fluorophores at these sites. We measured 18 chain segments in COS-7 cells using our developed pipeline (Fig. 2a and Supplementary Figs. 3–6; see Methods for details). In brief, our pipeline involved auto-segmentation to select the nuclear rim from each cell as a region of interest and extraction of measured donor fluorescence intensity profiles before and after acceptor photobleaching, defined as I(t) and I′(t), respectively. By subtracting I(t) from I′(t), the signals from the donor-only population and the background were eliminated, and the difference was taken as the pure FRET population (here we refer to the FRET population as the NUP98 chains specifically labelled with a pair of donor and acceptor dyes). We directly observed the differences in FRET efficiencies among all mutants in the lifetime fluorescence decay profiles (Fig. 2e). We detected a clear trend by which a smaller Nres showed a larger difference before and after acceptor photobleaching, indicating a higher FRET efficiency and smaller RE. To further visualize and compare the fluorescence lifetimes across cells, we converted individual fluorescence decay curves from the time domain into a phasor plot—a graphical view enabling a cell-by-cell analysis of the complicated lifetime curves (Fig. 2f and Extended Data Fig. 5; see Methods)44. Each point in the phasor plot represents the fluorescence decay of the nuclear envelope of a single cell. The phasor plot revealed not only some heterogeneity across cells but also an overall trend, in which greater Nres exhibited more left-shifted phasor values, corresponding to longer fluorescence lifetime and greater RE. These results suggested that FLIM–FRET could be successfully used to spatially distinguish and map the chain dimensions of the NUP98 FG domain.
We determined the chain dimensions of the NUP98 FG domain in situ according to the polymer scaling law, RE ~ Nresν, which relates the root-mean-square inter-residue distance RE to the chain length Nres, where ν is the scaling exponent12. In Flory’s homopolymer theory12, ν ~ 0.3 indicates that the polymer is very compact, as self-interactions dominate over interactions with the poor solvent; at ν = 0.5 those interactions are balanced, and at ν ~ 0.6, the interactions between the polymer segments and the solvent are maximized. Disordered protein chains thus expand in good solvent conditions (ν > 0.5) and compactify in poor solvents (ν < 0.5)45. Flory’s theory can be further extended to describe densely grafted polymer brushes where ν ~ 1 (ref. 46). Note that the scaling law is derived for infinitely long homopolymers. Despite its limiting definition, the law has been applied to calculate an apparent scaling exponent for finite-length proteins30,47. In brief, the apparent scaling exponent captures a complex distance distribution in one number and thus provides excellent economy in describing how protein conformational changes are tuned by their environment.
Here, we measured lifetime decay curves for the 18 chain segments of NUP98, probing the FG domain (Extended Data Fig. 5a), and globally fitted them to the Gaussian chain model to extract the scaling exponent ν (see Supplementary Text). Gratifyingly, the lifetime decays showed qualitative agreement between living cells and permeabilized cells (Extended Data Fig. 6). However, as concluded before, we refrained from quantitively analysing the scaling law for living cells due to the co-existence of intermolecular and intramolecular FRET. In the permeabilized cells, we obtained a scaling exponent ν = 0.56 ± 0.03 (more than 0.5) using a global fit that makes use of the data obtained from all measured cells (approximately 2,000). This indicates that the probed NUP98 FG domain adopted a rather extended conformation inside the NPC in cells (Fig. 3). We also performed a mutant-by-mutant analysis, in which no scaling law model is assumed a priori, but the error is higher, as only approximately 100 cells go into the average of each mutant. Reassuringly, we obtained a scaling exponent ν = 0.55 ± 0.05, which agrees well with the global fitting (see Supplementary Text and Supplementary Figs. 7 and 8 for the details of the analysis as well as the error estimate). Encouraged by these results, we further performed a global fitting of our data to a self-avoiding walk-ν model, which provides a better description of protein chains under good solvent conditions48. We obtained a scaling exponent ν = 0.56 ± 0.001 (see Extended Data Fig. 7 and Supplementary Text), showing the robustness of our scaling exponent for different polymer models.
On the left, the NUP98 FG domain showed more extended conformations inside the functional NPC (scaling exponent ν = 0.56 ± 0.03) and in the phase-separated condensate (NUP98 FG droplets on coverslips) (ν = 0.56 ± 0.04). The NUP98 FG domain at picomolar concentration showed a more collapsed conformation in physiological buffer on a single-molecule level (ν = 0.29 ± 0.01). Three reference lines at the theta-solvent state (ν = 0.5; dashed grey lines) are plotted, because the prefactor of the scaling law is unknown for a heteropolymer under different solvent conditions. The NUP98 FG domain in the NPC was labelled with AZDye594 tetrazine and LD655 tetrazine. The purified NUP98 FG domain in the phase-separated condensate and in the solution was labelled with Alexa Fluor 594 maleimide and LD655 maleimide. On the right, schematics show the conformations of FG-NUPs inside the NPC associated with extreme scaling exponents. If ν ~ 1, FG domains behave like densely grafted polymer brushes. If ν ~ 0.3, FG domains adopt a collapsed conformation.
Source Data
The result inside the NPC contrasts with our in vitro single-molecule FRET measurements on the purified NUP98 FG domain at picomolar concentration (Extended Data Fig. 8), for which the scaling exponent was ν = 0.29 ± 0.01 (less than 0.5). With water acting as a poor solvent for the hydrophobic side chains, the single NUP98 FG chain tends to become buried in a globule-like protein conformation in vitro. By contrast, the central channel of the NPC is enriched with other FG-NUPs, whose mutual attractive interactions establish conditions of a good solvent with expanded chains (ν > 0.5), reminiscent of the chain conformations in polymer melts12. On top of that, the presence of nuclear transport receptors, which exist in large quantities in the NPC49 (Extended Data Fig. 3e,f), and post-translational modifications (for example, glycosylation50 and phosphorylation), as well as transport cargos such as proteins and RNAs, could also contribute to good solvent conditions.
In reconstituted condensates, NUP98 FG domains tend to gelate over time, potentially by forming β-structures26,51,52. During the required FLIM–FRET measurement time of approximately 5 min, the resulting distance could thus be affected by condensate ageing despite using freshly prepared FG droplets. With this caveat, we obtained a scaling exponent of ν = 0.56 ± 0.04 for FG condensates, consistent with good solvent conditions (Fig. 3) and with the ability to function as a permeability barrier (Fig. 1c and Extended Data Figs. 3b and 9). However, the large prefactor of 7.8 Å in the scaling fit, compared with 5.5 Å in the self-avoiding walk-ν model, indicates heterogeneities in the ageing condensate.
The FLIM–FRET experiments assist MD simulations in providing us with a 3D view of the organization of FG-NUPs in functional NPCs. The disordered FG domains are grafted onto the NPC scaffold via folded domains. The recently published structure of the human NPC scaffold7 enabled us to anchor FG-NUPs with the correct positions, orientations and grafting densities. This allowed us to build a coarse-grained bead–spring polymer model53 and perform MD simulations of the so-far elusive FG domains attached to the NPC scaffold (see Methods). We first parameterized the effective NUP–NUP interaction strength as \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44 (defined in equation (14) in the Methods) by matching the phase behaviour of in vitro-reconstituted NUP98 FG condensates that mimic the permeability barrier27,54 (Extended Data Fig. 10 and Supplementary Video 1). However, if the same interaction strength of \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44 was applied to the whole-NPC simulations, the inner-ring FG-NUPs collapsed onto the scaffold to form surface condensates, despite only weak direct interactions between FG-NUPs and the scaffold (\(\mathop{{\epsilon }}\limits^{ \sim }\)scaffold = 0.1). The surface condensates left a void at the centre of the pore with a diameter of approximately 20 nm, which would not seem consistent with the function of the NPC to block the unaided passage of large cargo. After a slight adjustment of the interaction strength to \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42, the FG-NUPs recovered the RE matching our FLIM–FRET data (Fig. 4a, Supplementary Figs. 9–11 and 18 and Supplementary Videos 2–7). We observed that under such a condition, the FG-NUPs formed extended coil configurations and the inner-ring FG-NUPs fluctuated extensively to form a dynamic barrier across the central channel (Fig. 4b and Supplementary Video 3). Remarkably, the optimal FG-NUP interaction strength determined in this way nearly coincides with the critical strength \(\mathop{{\epsilon }}\limits^{ \sim }\)c ≈ 0.42 at which condensation occurs. The large FG-NUP motions seen in the MD simulations thus amount to critical fluctuations that create a highly dynamic polymer network well suited for the rapid but selective molecular transport in the permeability barrier55,56,57,58. In an NPC model with explicit solvent (Supplementary Table 4), the extension of NUP98 follows the experimental FRET measurement at an interaction strength of \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42, again close to the respective critical value of \(\mathop{{\epsilon }}\limits^{ \sim }\)c ≈ 0.41 (see Supplementary Text and Supplementary Figs. 12 and 13). We note that \(\mathop{{\epsilon }}\limits^{ \sim }\) should be considered effective, as it accounts, for example, for the presence of cargo and nuclear transport receptors unresolved in our in situ experiments. Phosphorylation, glycosylation50 and other post-translational modifications of the FG-NUPs in the NPC will alter \(\mathop{{\epsilon }}\limits^{ \sim }\) and may in turn impact their state, given the sensitivity to \(\mathop{{\epsilon }}\limits^{ \sim }\) for states close to the critical point (Supplementary Fig. 12).
a, Root-mean-square inter-residue distances RE of beads on the same NUP98 FG domain in the NPC as a function of the number of amino acid residues Nres with different effective NUP–NUP interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) on linear and (inset) log-scale. The symbols and error bars represent the average and standard error of the mean, respectively, as estimated from four blocks of size 104τ. Lines are guides to the eye. The distances obtained from FLIM–FRET experiments (Exp.) are represented by filled circles, which overlap with the simulations for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42. b, Top and side views of the NPC at the end of the MD simulations with \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.35 (left), \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42 (middle) and \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44 (right) (scaffold is shown in blue, NUP98 is in yellow and other FG-NUPs are in red). Scale bar, 20 nm. See Supplementary Table 3 for details on NPC model II and Supplementary Fig. 18 for a view of the 48 NUP98 chains alone.
Source Data
We probed for the possible effects of sequence heterogeneity by adapting the stickers-and-spacers concept59 for our MD simulations. In the FYW model, we treated aromatic residues (F, Y and W) as mutually attractive stickers separated by weakly interacting spacer regions (see Supplementary Text and Supplementary Fig. 14). As for the homopolymer models, we calibrated the effective sticker interaction strength \({\mathop{{\epsilon }}\limits^{ \sim }}_{\text{FYW}}\) against the FLIM–FRET measurements of the NUP98 FG domain and applied it to all FG-NUPs in the NPC. We obtained very good agreement of the calculated NUP98 extension with the FLIM–FRET data for \({\mathop{{\epsilon }}\limits^{ \sim }}_{\text{FYW}}\) = 3.5, which again lies just close to the critical value (Supplementary Figs. 15–17), as in the homopolymer models. We conclude that a more detailed model accounting for sequence heterogeneity does not substantially alter the conclusions.
We further explored the possibility of heterogeneity in the chain ensembles by comparing the calculated extensions of the 8 × 6 = 48 NUP98 chains grafted to the 8-fold symmetric NPC scaffold at six non-symmetry-equivalent positions (Supplementary Figs. 19 and 20). Consistently across the homopolymer models I and II (Supplementary Tables 2 and 3) and the heteropolymer FYW model, we found that the NUP98 chains at the UCc positions in the cytoplasmic ring are somewhat more extended than at the Uno positions in the inner ring on the nuclear side (Supplementary Fig. 20). Our MD simulations thus indicate some heterogeneity in NUP98 chain extension depending on their precise ‘grafting’ point within the NPC, including in the scaling behaviour (Supplementary Fig. 20d). A detailed account of electrostatic interactions could further enhance the heterogeneity60. The unimodal distributions of the inter-residue distances provide no indication of co-existing collapsed and extended populations (Supplementary Fig. 21).
Here we developed an experimental approach using site-specific fluorescent labelling of IDPs in mammalian cells and FLIM to directly decipher their plasticity by FRET measurements in the functional state. We showed that this approach works for the sub-resolution permeability barrier of the NPC, a nanocavity with a diameter of approximately 60 nm, filled with approximately 50 MDa of highly concentrated FG-NUPs. By measuring the inter-residue distances of different segments of the labelled FG-NUP98 using FLIM–FRET in permeabilized cells with functional transport machinery, we obtained the distribution from which we could estimate an apparent scaling exponent and revealed that the intact NPC environment provides a good solvent (ν > 0.5) in which the FG domains adopt extended conformations compared with the collapsed solution state in vitro.
Owing to recent advances in determining the structure of the NPC scaffold, the actual grafting points for the FG domains are now known with high confidence4,5,6,7,8,9,10,11. We synergistically combined our residue-specific FLIM–FRET measurements with computational modelling and investigated the conformational behaviour and interaction mechanisms of FG-NUPs at the molecular level with residue resolution. We found that a simple polymer model could capture the motions of FG-NUPs. The model reproduced the FLIM–FRET results at a near-critical interaction strength, which defines the energetic threshold for forming protein condensates and is associated with large fluctuations in the dynamic polymer network. For weaker interactions, the FG-NUPs were too loose to generate a polymer network. For stronger interactions, the FG-NUPs formed a surface condensate, the collapse of which onto the scaffold was driven by the geometry of the grafting points, not by the weak direct interactions with the scaffold. Indeed, the chain extension is exquisitely sensitive to the NUP–NUP interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) in the NPC (Fig. 4a), whereas in bulk condensates, this dependence is negligible (Extended Data Fig. 10b).
A key finding is that the ensembles of FG-NUPs within functional nuclear pores inside cells are distinct from the chain ensembles sampled in solution, untethered from the pore. The intriguing inference presented is that the apparent solvent quality of the nuclear pore is better than that of a bulk, aqueous solvent, with conformational statistics measured in the pore resembling what one would observe in at least a theta solvent (ν = 0.5) and more closely, a good solvent (ν > 0.5). Our study removes much speculation about the conformational state of FG-NUPs in the NPC and provides a sound coarse-grained model with amino acid precision to explain the function of the permeability barrier. The measured apparent scaling exponent ν = 0.56 ± 0.03 disaccords with, for example, the polymer brush model46 (ν ~ 1) and the forest model24 (ν ~ 0.3), which describes the low charge-content, cohesive FG domain as a globular structure. Remarkably, we also showed that the parameterization based on in vitro reconstitution studies failed to reproduce a functional pore. Despite having similar permeability barrier properties as the intact NPC, the bulk condensate formed from phase separating NUP98 is an incomplete approximation of the actual permeability barrier, the materials properties of which are modulated by the anchoring of a distinct number of FG-NUPs with 3D precision on a half-toroidal NPC scaffold. In terms of nuclear transport selectivity, there are consequences: whereas a surface condensate would leave a substantial hole at the centre, we found the hole to be filled by FG-NUPs at near-critical conditions (Fig. 4b). These results emphasize the importance of interrogating the permeability barrier in situ to reconcile different transport models and understand the molecular basis for nuclear transport.
In this work, a range of technologies have been used, which finally enabled us to measure the distance distribution of a key FG-NUP inside the functional transport machinery of the NPC. The combination of fluorescent lifetime technologies and chemical synthetic biology tools compensates for the weaknesses or ambiguities that can originate from either of the methods alone. We show that our labelling technology is also live-cell compatible at the cost of a lower signal-to-noise ratio, which originates from stronger background sticking of the tested membrane-permeable dyes (Extended Data Fig. 2). Although this problem might not be an issue for studying biological systems that are more abundant, for the NPC we could only solve it by increasing expression levels. However, owing to the fact that 48 NUP98 chains are cramped into the nanosized cavity of the NPC, with approximately 250 FG-NUPs in total (Fig. 4), intermolecular FRET can become a concern. We show a procedure to detect intermolecular FRET; with this, the method should be generally extendable to other systems. To avoid intermolecular FRET within the nanosized NPC, we had to work at very low expression levels. Sufficient signal-to-noise was then ensured by mild permeabilization of the plasma membrane, which enabled us to use hydrophilic dyes that stick less. At the same time, we validated that the nuclear transport machinery is still intact in situ. Owing to the potential protein truncation problem arising from translation termination at the reassigned codon in GCE39,61, we refrained from probing the FG-NUPs with N-terminal anchoring domains (for example, FG-NUP214), where the truncated proteins could dock on the nuclear pore and interfere with labelling or even transport functionality. We currently could only infer the molecular behaviours of other FG-NUPs by MD simulations. It is highly desired to experimentally probe all FG-NUPs in living cells when more advanced GCE and highly fluorogenic cell-permeable dyes compatible with GCE are available in the future. However, as NUP98 has been shown to be the only NUP that is vital to reconstitute functional transport activity of intact NPCs36, our study captures fundamental properties of NPC function.
The most advanced computational modelling tools, such as the AI-based AlphaFold, can very precisely predict the structures of folded proteins. However, when applied to IDPs, they still require experimental constraints and validation62. It is necessary to probe IDPs with in situ experimental measurements. As we showed here, a protein chain with globular-like character in vitro can adopt a much more expanded conformation in the functional state due to the cellular environment. Our work is an example of how the structural knowledge from cryo-electron tomography—which revealed the anchoring sites of IDPs—paired with in situ FLIM–FRET—which revealed the conformational state of the proteins inside the nanosized NPC—can yield a complete coarse-grained picture of a cellular machine enriched in IDPs inside a cell. The tools developed here could be applied generally to study the plasticity and functions of many other IDPs in the cell, filling a major technology gap in the field.
COS-7 cells (87021302, Sigma) were maintained in Dulbecco’s modified eagle medium (41965-039, Thermo Fisher) supplemented with 10% v/v fetal bovine serum (F7524, Sigma), 1% penicillin–streptomycin (15140-122, Thermo Fisher), 1% l-glutamine (25030-081, Thermo Fisher) and 1% sodium pyruvate (11360-070, Thermo Fisher) at 37 °C and 5% CO2, and passaged every 2–3 days up to 15–20 passages. The COS-7 cells were authenticated by the manufacturer, validated by morphology and regularly tested for mycoplasma contamination, with negative results.
Cells were trypsinized (trypsin-EDTA; 25300-054, Thermo Fisher) and seeded into a 35-mm imaging dish (81158, ibidi) 24 h before transfection. The cells were transfected at a confluency of 60–70% with plasmids of interest listed in Supplementary Table 1 in Supplementary Information (for amber suppression with synthetic organelles, for example, pcDNA3.1-hsNUP98TAG-BoxB and pcDNA3.1-TOM20-FUS-λN22-PylRS(AF)-tRNA at a mass ratio of 1:1) using jetPRIME according to the manufacturer’s protocol. After 4–5 h, the medium was changed, and 10 mM HEPES with 50 μM trans-cyclooct-2-en-l-lysine (TCO*A; SciChem; for in-cell labelling) or t-butyloxycarbonyl-l-lysine (BOC; Iris Biotech; for control measurements) were added, respectively.
At 20–24 h post-transfection, COS-7 cells were washed with fresh culture medium supplemented with 50 μM BOC and 10 mM HEPES and incubated at 37 °C for 2 h to get rid of the excessive TCO*A. Then, the cells were incubated with 250 nM cell-permeable dye (Janelia Fluor 549-tetrazine (Tocris), Janelia Fluor 646-tetrazine (Tocris), silicon rhodamine-tetrazine (Spirochrome) or TAMRA-tetrazine (click chemistry tools)) in the serum-free culture medium at 37 °C for 45 min. After the labelling, the cells were washed with fresh culture medium four times in 2 h at 37 °C. Cells were imaged immediately at room temperature in the phenol red-free culture medium supplemented with 10 mM HEPES. The lid of the imaging dish was kept closed for the duration of cell imaging for 2 h.
At 20–24 h post-transfection, COS-7 cells were washed twice with transport buffer40 (TB; 20 mM HEPES, 110 mM KOAc, 5 mM NaOAc, 2 mM Mg(OAc)2, 1 mM EGTA and 2 mM DTT, pH 7.3, supplemented with PEG6000 (5 mg ml−1) to avoid osmotic shock41), and permeabilized for 2 min with 20 μg ml−1 digitonin (A1905, AppliChem) in TB. After the permeabilization, cells were washed twice with TB to remove digitonin and labelled with a dye solution containing 33.3 nM AZDye594-tetrazine (Click Chemistry Tools) and 66.6 nM LD655-tetrazine (Lumidyne Technologies) in TB for 2 min. Here we optimized the molar ratio between donor and acceptor dyes as 1:2 to reduce the donor-only population. To remove residual dyes, the cells were washed in TB and incubated at 37 °C for 30 min. Cells were imaged immediately at room temperature. The lid of the imaging dish was kept closed for the duration of cell imaging for 2 h.
COS-7 cells were transfected, permeabilized and labelled with 100 nM LD655-tetrazine as described above. After the last wash, the imaging dish was mounted on the custom-built confocal microscope. After gently removing the washing buffer, 0.5 μM of 70-kDa FITC–dextran (53471, Sigma) in TB was added to the dish, incubated for 10 min and the cells were imaged.
To allow the import complex to form, 0.5 μM IBB–MBP–GFP cargo was first pre-incubated with 1 μM importin-β on ice for 10 min, and then combined with the rest of the transport mix (5 μM RanGDP, 4 μM NTF2 and 2 mM GTP in TB, where importin-β, Ran and NTF2 were purified as previously described55,63). After that, the complete import mixture was added to the permeabilized and labelled COS-7 cells as described above.
To ensure that the NPC was functional during the 2-h time window of FLIM–FRET measurements for each imaging dish, passive exclusion assay and active transport reactions were performed for both the freshly labelled cells and the cells incubated in TB at room temperature for 2 h after the labelling steps (Extended Data Fig. 3c,d).
After labelling with 100 nM LD655-tetrazine in TB as described above, the COS-7 cells were immediately fixed with 2% PFA in PBS for 10 min or kept at room temperature for 2 h and then fixed. After washing with PBS twice, the cells were permeabilized with 0.5% Triton X-100 in PBS for 15 min. Then, the cells were washed with PBS twice and incubated with blocking buffer (3% BSA in PBS) for 90 min. The cells were subsequently immunolabelled with anti-KPNB1 antibody (ab2811, Abcam; at 1:1,000) and anti-mouse Alexa Fluor 488 secondary antibody (A-11001, Thermo Fisher; at 1:1,000). Immunofluorescence demonstrated that endogenous importin-β colocalized at the nuclear envelope after permeabilization and was retained during the 2-h time window of FLIM–FRET measurements (Extended Data Fig. 3e,f).
Homo sapiens NUP98 FG domain (1–505 amino acids, with or without the Gle2-binding domain (GLEBS; 157–213 amino acids), a structured domain in between the two FG domains) was cloned into pQE-14his-TEV vectors. We also designed a construct without GLEBS, because the structured GLEBS domain might misfold and trigger aggregation when rapidly changing the buffer condition from denaturing to native in the droplet assay (see ‘In vitro droplet assay and FLIM–FRET measurements’). The GLEBS domain is outside the residue segments probed by the in-cell experiments. Unless otherwise stated, we used the NUP98 construct without GLEBS for in vitro experiments.
Escherichia coli BL21 AI cells containing NUP98 FG construct with or without GLEBS were grown in terrific broth medium containing 50 µg ml−1 of kanamycin at 37 °C and 200 rpm, and protein expression was induced with 0.02% (w/v) arabinose and 1 mM IPTG at OD600 = 0.6. After 16 h of expression at 18 °C, cells were harvested by centrifugation and lysed in a cell disruptor (CF1, Constant Systems) in the lysis buffer containing 6 M GdmCl, 0.2 mM Tris(2-carboxyethyl)phosphine (TCEP), 20 mM imidazole and 50 mM Tris-HCl, pH 8. The lysate was centrifuged for 1 h at 12,000g at 4 °C to remove cell debris, and the supernatant was incubated with Ni-beads for 2 h at 4 °C. The Ni-beads with lysate were loaded in polypropylene tubes (Qiagen), washed twice with 2 M GdmCl and 20 mM imidazole, pH 8, and eluted with buffer containing 2 M GdmCl and 500 mM imidazole, pH 8. To remove the His-tag, the elution was dialysed in 0.5 M GdmCl and 50 mM Tris-HCl, pH 8, and cleaved overnight with TEV protease at room temperature. Proteins were incubated again with Ni-beads in 2 M GdmCl and 50 mM Tris-HCl, pH 8, to remove cleaved His-tags, TEV protease and nonspecific proteins with Ni-bead affinity. The flow-through containing the NUP98 FG domain was collected and further purified by size-exclusion chromatography (Superdex 200, Akta pure protein purification system, Cytiva) in 2 M GdmCl, 0.2 mM TCEP and 50 mM Tris-HCl, pH 8. Fractions were analysed by SDS–PAGE and stained with Coomassie blue. Pure fractions were pooled and concentrated to around 15 mg ml−1 in 4 M GdmCl using 3-kDa MWCO centrifugal filters (Merck Millipore), with the concentration measured by a BCA protein assay kit (Thermo Fisher). The proteins were flash-frozen and stored at −80 °C.
For labelling, the purified NUP98 FG domain with single or double cysteine mutations was exchanged to 4 M GdmCl, 1× PBS, 0.1 mM EDTA and 0.2 mM TCEP, pH 7. Labelling with Alexa Fluor 594 maleimide (A10256, Thermo Fisher) and LD655-maleimide (Lumidyne Technologies) was done at the molar ratio of 1:2 (dye:protein) overnight at 4 °C. The reaction was quenched with 10 mM DTT in 4 M GdmCl and 1× PBS, pH 7. Unreacted dye was washed off using a 3-kDa MWCO centrifugal filter, and the labelled protein was further purified with Superdex 200. Pure fractions were chosen, pooled and concentrated as described above, and the final concentration was measured by the absorbance spectrometer Duetta (Horiba). The proteins were flash-frozen and stored at −80 °C.
The custom-built FLIM–FRET imaging setup (Supplementary Fig. 6) was equipped with picosecond pulsed laser diode heads including the wavelengths of 485 nm (LDH-D-C-485, PicoQuant), 560 nm (LDH-D-TA-560, PicoQuant) and 660 nm (LDH-D-C-660, PicoQuant). The laser heads were controlled through a multichannel picosecond diode laser driver (Sepia II PDL 828, PicoQuant). The beams were coupled into a single-mode polarization-maintaining optical fibre (KineFLEX-P-2-S-405/640-2.5-2.5-p2) and fibre coupler (60FC-4-RGBV11-47, Schäfter + Kirchhoff). The beam travelled through a Glan-laser polarizer (Thorlabs) and was directed into a laser scanning system (FLIMbee, PicoQuant). The three galvo mirrors in the scanning system were imaged onto the backfocal plane of the objective (×60 SR Plan Apo IR, 1.27 NA; Nikon) with 200-mm tube lens. The fluorescence emission was focused onto a pinhole (100 μm for cell measurements, 50 μm for droplet measurements and 100 μm for single-molecule FRET measurements), and then separated into parallel and perpendicular components using a 50/50 polarizing beam splitter (Thorlabs). Each component was further separated by two sets of beamsplitters (ZT561 RDC and T647 LPXR, Chroma), passed through three sets of bandpass filters (green channels: 525/50 BrightLine HC; orange channels: 609/57 BrightLine HC, Semrock; red channels: ET700/75m, Chroma), and focused onto the single-photon counting detectors (green channels and orange channels: PMA Hybrid 40, PicoQuant; red channels: τ-SPAD, PicoQuant). The signals from the photon detectors were recorded by a TCSPC system (HydraHarp 400, PicoQuant) at a time resolution of 16 ps. Data acquisition was carried out with SymPhoTime 64 software v2.6 (PicoQuant).
The Förster radius R0 is the distance between a pair of fluorophores at which the FRET efficiency is 50%, which was calculated as43,
where \({\kappa }^{2}\) is the orientation factor, n is the refractive index (n = 1.375 for cell measurements64, and n = 1.426 for in vitro condensates65), \({\varphi }_{D}\) is the quantum yield of the donor without energy transfer, and \(J(\lambda )\) is the overlap integral between the donor emission and acceptor absorption spectra at wavelength \(\lambda \), given by43,
where \({F}_{D}\left(\lambda \right)\) is the radiation emission intensity of the donor at wavelength \(\lambda \), and \({\varepsilon }_{A}\left(\lambda \right)\) is the extinction coefficient of the acceptor. The emission spectra of the donor and the excitation spectra scan of the acceptor were measured on a Leica SP8 STELLARIS microscope using \(xy\lambda \) or \(xy\Lambda \) acquisition mode, for single-labelled cells with the donor or the acceptor, respectively. The acquired spectra were then compared with the spectra measured for free dyes in TB using an absorbance spectrometer (Duetta) and no spectral shift was detected. \({\kappa }^{2}\) can be assumed to be approximately two-thirds when the dyes can freely rotate, which is expected to be the case as all our dyes have a C5 flexible linker between the conjugating group and the chromophore. Here we assumed near-parallel excitation and emission dipoles. We also measured the fluorescence anisotropies for the sample labelled with the donor dye inside the NPC (0.27 ± 0.015), in the FG condensates in vitro (0.28 ± 0.003), and in the solution on a single-molecule level (0.15 ± 0.003). The measured anisotropies were found to be smaller than 0.3 in all cases, and the error in the distance measurements was shown to be below 10% in such a regime43. The R0 for in cell measurements was determined as 77.1 Å, and for in vitro condensate measurements was 77.3 Å. Note that the apparent scaling exponent is robust to the error in R0 measurements (see Supplementary Fig. 7).
The average power of laser excitation was optimized to collect enough photons from the cell within a reasonable time but avoid photon pile-up and other artefacts in the fluorescence lifetime measurements. The instrument response function was measured on a daily basis using a freshly prepared saturated solution of KI and erythrosine B66. The temporal offsets of the parallel and perpendicular detectors were pre-aligned by the measured instrument response functions. The cell measurements were performed using the following imaging settings: pixel size of 100 nm, image size of 256 × 256 pixels, pixel dwell time of 150 μs and the time resolution of 16 ps. The fluorescence photons were detected in T3 mode and collected from the perpendicular and parallel detectors for each colour individually.
To measure FLIM–FRET, the morphology of the labelled cell was first checked with 660-nm laser excitation to ensure that no GLFG bodies (normally a sign of overexpressed NUP98 in the cell67) existed in the nucleus. The acceptor intensity per pixel of the nuclear rim was then checked with 660-nm laser excitation to further assess the expression level of the mutant NUP98, ensuring that cells with similar expression levels and not highly overexpressed mutant NUP98 were chosen. Specifically, we determined such an acceptor intensity range based on the criteria that, on the nuclear rim, the relative FRET efficiency (that is, the proximity ratio, \({E}_{{\rm{rel}}.}=\frac{{I}_{{\rm{A}}}}{{I}_{{\rm{A}}}+{I}_{{\rm{D}}}}\,\); IA and ID are the total acceptor and donor fluorescence intensities, respectively, excited by a 560-nm laser) did not correlate with the acceptor intensity per pixel excited by a 660-nm laser (Supplementary Fig. 3). This intensity threshold was further verified by acceptor photobleaching assays on cells that expressed single-amber-mutant NUP98 and labelled with donor and acceptor dye mixture. The average fluorescence lifetime of the donor dye did not change before and after acceptor photobleaching as shown in Fig. 2b, indicating that no intermolecular FRET could be detected.
After checking the expression level of the mutant NUP98 with 660-nm laser excitation, we imaged the selected cell for 5 min using 560-nm laser excitation with an average power of 40 μW at 40 MHz, and then for 30 s using 660-nm laser excitation with an average power of 35 μW at 40 MHz. Next, the acceptor labelling was photobleached using 660-nm laser excitation with an average power of 300 μW at 40 MHz for 2 min. The donor signal was measured again post-photobleaching for 5 min using 560-nm laser excitation with an average power of 40 μW at 40 MHz.
The recorded images were processed using an automatic segmentation pipeline developed based on the software package PAM in MATLAB68. The nuclear rim was selected as a region of interest (ROI) using a thresholding algorithm according to the intensity and average lifetime of each pixel. The time-resolved donor fluorescent intensity profiles before and after acceptor photobleaching were extracted from the selected ROI, respectively. The total fluorescence decay was calculated by combining the parallel and perpendicular fluorescence decays (\({I}_{\parallel }\) and \({I}_{\perp }\), respectively)69,
where \({L}_{1}\) and \({L}_{2}\) are factors accounting for polarization mixing caused by the high numerical aperture objective lens, and G is the factor accounting for the difference in the detection efficiencies \(\eta \) between parallel and perpendicular polarization, given by,
The detected fluorescence decays in parallel and perpendicular channels, Dpar and Dper, respectively, for a reference sample are expressed as64
To determine \({L}_{1}\) and \({L}_{2}\), the fluorescence decays of YFP with known fluorescence rotational relaxation time (ρ = 16 ns)70 was measured and globally fitted with equations (5 and 6).
G factor was determined as a ratio of the average intensities of the perpendicular and parallel donor channels when measuring a solution of Tris(2,2′-bipyridyl)dichlororuthenium(II) chloride ([Ru(bpy)3]Cl2). For each mutant, the total fluorescence decays for approximately 100 cells were added up for further fitting analysis (see the sections about FLIM analysis methods and fitting the scaling law for NUP98 FG in the NPC in the Supplementary Text).
The purified and labelled NUP98 FG domain was mixed with unlabelled protein at a molar ratio of 1:5,000 in 2 M GdmCl and 1× PBS, pH 7. Of such a mixture, 1 µl was then quickly mixed with 24 µl TB supplemented with 5 mg ml−1 PEG6000 in a chambered coverslip (81507, ibidi) and imaged immediately with the custom-built confocal microscope. The final total concentration of NUP98 in the system was 10 µM for unlabelled and 2 nM for labelled. The trace amount of GdmCl left in solution was 80 mM (a negligible level).
For the passive exclusion assay and the facilitated/active transport assay on the phase-separated condensates, unlabelled NUP98 and the LD655-labelled condensates were mixed with the ratio described above, to avoid the crosstalk with the fluorescently labelled cargos (that is, FITC or GFP). After mixing the protein with TB for 5 min, the buffer was carefully replaced by 0.5 μM 70-kDa FITC–dextran in TB or 0.5 μM IBB–MBP–GFP in the transport mixture as described for the labelled cells avoiding disturbing the condensates on the coverslip surfaces and imaged immediately.
To measure FLIM–FRET on the NUP98 FG condensates, we focused on the first 5 min as the condensates behaved more liquid-like during this time26 (for example, FG-droplets merged quickly as shown in Extended Data Fig. 9a). We applied the same conformation distribution model (that is, the Gaussian chain model) to fit the lifetime curves of the FG condensates as we used for in cell measurements (see Supplementary Text).
The formed condensates were measured for 5 min using 560-nm laser excitation with an average power of 70 μW at 40 MHz, and the imaging settings were as follows: pixel size of 200 nm, image size of 256 × 256 pixels, pixel dwell time of 100 μs and time resolution of 16 ps. The fluorescence photons were detected in T3 mode and collected from the perpendicular and parallel detectors for each colour individually. The FG-droplets were selected as the ROI based on the intensity and average lifetime in the SymPhoTime 64 software. The time-resolved donor fluorescent intensity profiles were extracted from the ROI pixels for analysis.
To ensure that no intermolecular FRET was detected, single-cysteine-mutated NUP98 was labelled with either the donor dye or with the mixture of donor and acceptor dyes. The labelled protein was then mixed with the unlabelled protein with the same concentration as described above to perform the droplet assay. The donor fluorescence lifetimes of the donor-only and the donor–acceptor mixture showed no difference, indicating that no intermolecular FRET could be detected (Extended Data Fig. 9b,c).
To analyse the lifetime decays on a single-cell basis, the raw intensity profile I(t) of each selected nuclear rim was plotted as a single point in a phasor plot by applying the Fourier transform to the measured decay data, given by71
where \(g\,(\omega )\) and \(s\,(\omega )\) are the x and y coordinates of the phasor plot, respectively, m and φ are the modulation and the phase delay of emission with respect to the laser excitation that has the repetition rate of f, that is, 40 MHz, and T is the repeat frequency of the acquisition, that is, 25 ns. To establish the correct scale for the plotted phasor points, the coordinates of the phasor plot were first calibrated by applying Fourier transform to the measured instrument response function trace and setting it as the zero lifetime71. Each phasor point from the acquired FLIM data was then calibrated accordingly using the same calibration parameters so that the final phase plot was referenced relative to the calibration standard. The above procedure was performed with a self-written code in MATLAB.
Single-molecule FRET measurements were performed with the custom-built confocal microscope mentioned above. The sample was illuminated with 560-nm and 660-nm lasers in pulsed interleaved excitation mode at a repetition rate of 32 MHz. Photon counts were recorded with a resolution of 16 ps. Purified and double-labelled NUP98 was diluted to the final protein concentration of 50 pM in 1× PBS supplemented with 10 mM fresh DTT. The donor dye was excited with a 560-nm laser at an average power of 70 μW and a 660-nm laser at an average power of 20 μW. Data acquisition was carried out with the SymPhoTime 64 software (Picoquant).
Acquired data were analysed with the PAM software package68 for burst search, and the identified bursts were further analysed in the BurstBrowser. The FRET efficiency E in a burst is defined as72,
and the stoichiometry S is defined as67,
where \({I}_{A}^{D}\) is the acceptor fluorescence detected upon donor excitation, \({I}_{D}^{D}\) is the donor fluorescence upon donor excitation, \({I}_{A}^{A}\) is the acceptor fluorescence detected upon acceptor excitation, and γcor is the factor accounting for the detection efficiency of the acceptor and donor channels. After performing the correction for the leakage of donor fluorescence into the acceptor channel (α = 1.00) and direct acceptor excitation (δ = 0.30) and confirming minimal variation of quantum yields among different mutants, the γcor parameter was extracted from the apparent FRET efficiency, Eapp, and the apparent stoichiometry, Sapp (ref. 72). A linear fit to a plot of 1/Sapp versus Eapp yields intercept a and slope b, which relates to γcor in the following way,
We determined the γcor parameter of labelled NUP98 in 1× PBS to be 2.51. Single-molecule transfer efficiency histograms are shown in Extended Data Fig. 8. Each dataset shows the donor-only species at zero FRET efficiency (owing to incomplete labelling or photophysically inactive acceptor such as due to bleaching) and a FRET population. The FRET population was selected and fitted with an asymmetric Gaussian function to determine the centre of the population while accounting for non-linear effects at high FRET efficiency. The fitted values were taken as the average FRET efficiency, which is defined as73,
where ρ(r) describes the distance distribution of the inter-residue distance. A typical single-molecule FRET (smFRET) experiment does not contain enough information to retrieve a model-free distance distribution, and one must therefore choose a model to fit. Here we adopted the Gaussian chain model, which assumes that the monomers occupy zero volume and excluded volume effects are not considered. Despite its simplicity, it is commonly used for the analysis of IDPs30,74. The distance distribution function takes the form given by the Gaussian chain model as73
where the root-mean-square inter-residue distance \(\sqrt{\langle {r}^{2}\rangle }\equiv {R}_{{\rm{E}}}\). By comparing the fitted \(\left\langle E\right\rangle \) with the calculated average FRET efficiency derived from the model, \({R}_{E}\) was obtained for each mutant.
We performed MD simulations to model the large-scale conformations and motions of FG-NUPs in the NPC and in in vitro-reconstituted NUP98 FG condensates. To ensure efficient sampling, we used a coarse-grained polymer model, in which each amino acid is represented by a single bead without local structure.
The model complements the recently resolved scaffold structure7 of the human NPC with a nearly complete set of FG-NUPs. In our model I, we did not include NUP153 and POM121, whose anchor positions are not yet resolved with high confidence (see Supplementary Table 2). In our model II, we included NUP153 and POM121 (see Supplementary Table 3). The FG-NUPs were grafted onto the NPC scaffold in the constricted state (PDB ID: 7R5K) at established positions as indicated7. The NUP98 simulation models included the GLEBS domain, which is outside the residue segments whose distances were probed in simulations or experiments.
We used a FENE potential UFENE53 to describe the disordered FG-NUPs, Lennard–Jones (LJ) interactions ULJ between the CG beads, and a potential confining NUP98 chains. The total potential energy is
The length and interaction energy scales of the LJ potential between a pair of beads i and j are \(\sigma \) and \({\widetilde{{\epsilon }}}_{{ij}}\), respectively. The latter is normalized by the thermal energy kBT, where kB is the Boltzmann constant and T is the system temperature. We fixed the interaction strength between scaffold and FG beads at \({\widetilde{{\epsilon }}}_{{ij}}={\widetilde{{\epsilon }}}_{\text{scaffold}}=0.1\), which ensures that chains do not stick to the scaffold. The interaction strength between FG beads is given uniformly by \({\widetilde{{\epsilon }}}_{{ij}}=\widetilde{{\epsilon }}\). The sums over ⟨i,j = i + 1⟩ extend over the bonded beads of the FG-NUPs. The bond contraction and stretching are controlled by the repulsive LJ and logarithmic terms, whose maximum bond limit is RFE = 1.5σ with kFE = 80 kBT. To mimic the membrane envelope underneath the scaffold, we applied an axial confinement on NUP98-FG and POM121-FG beyond a radius of Rc = 81σ ≡ 48.6 nm, where \({r}_{{xy}}=\sqrt{{x}^{2}+{y}^{2}}\) and kc = 80 kBT.
We used the LAMMPS software package75 to simulate the polymeric systems. The systems were thermalized with a Langevin thermostat76 (at kBT = 1) with a damping coefficient of 10τ. We used a uniform mass \(m\) for all monomers and a characteristic time scale \(\tau =\sqrt{m{\sigma }^{2}/{k}_{{\rm{B}}}T}\,\). The time step of all simulations was set to 0.01 τ for the single-chain and condensate simulations and to 0.001 τ for NPC simulations. We used block averaging with four non-overlapping blocks to estimate uncertainties. In the MD simulations of the FG-NUPs attached to the NPC scaffold, the residues of the FG-NUPs grafted to the scaffold (see Supplementary Tables 2 and 3 for details) and the scaffold itself were kept frozen.
We determined the bond length σ by matching the root-mean-square inter-residue distance distribution of the CG model of a single NUP98 FG chain (1–499) to that in the Martini 2.2 level of coarse graining77,78. With a focus on the geometric extensibility, we worked in the limit of weak non-bonded interactions between distant amino acids. For the Martini simulations of the single chain, we used α-scaling (α = 0.1)79, and for the CG model, a weak cohesive strength between FG beads (\(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.1). The Martini MD simulations were performed with GROMACS 2020.6 (refs. 80,81). The first 499 N-terminal amino acids of NUP98 were converted into a Martini coarse-grained model using the martinize.py python script77,78. A cubic simulation box of length 30 nm was built and solvated with coarse-grained Martini water and 10% anti-freezing water. Ions were added to neutralize the system. The system was initially energy minimized and then equilibrated for 500 ns in NVT and NPT ensembles, respectively, using the velocity scaling thermostat82 and the Berendsen barostat83 at temperature T = 300 K and pressure P = 1 atm. The time constant of the thermostat was 1 ps, and that of the barostat was 5 ps with a compressibility of 3 × 10−4 bar−1. For a 10-μs long production run in the NPT ensemble, we used a velocity scaling thermostat82 and a Parrinello–Rahman barostat84 with time constants of 1 ps and 12 ps, respectively. As shown in Supplementary Fig. 9b, the probability distribution of the end-to-end distance of the Martini model peaked at the same position as that of the CG model with \(\sigma \equiv 0.6\,{\rm{nm}}\), and the averaged end-to-end distance of both models were equal. In the following, we fixed the FENE bond length at \(\sigma \equiv 0.6\,{\rm{nm}}\).
In a second step, we adjusted the interaction energy between the FG beads (\(\mathop{{\epsilon }}\limits^{ \sim }\)) by matching the model to the measured thermodynamic properties of the NUP98 FG chain (1–499) condensate formation27,54. As shown in Extended data Fig. 10a, 500 chains were simulated inside a simulation box size of 360 × 90 × 90 σ3. We adjusted the interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) to match the experimentally measured concentration of a NUP98 FG condensate27,54 (see Extended data Fig. 10b). We note that in this way, we also got a good fit to the measured concentration of the dilute phase and thus the transfer free energy. The simulations of the condensate and NPC systems were performed for durations exceeding 1.8 × 106 τ and 6 × 104 τ, respectively. We used the VMD software to visualize all systems85.
We calculated the root-mean-square distances between the labelled sites across NUP98 FG chains and simulation trajectories to compare to the smFRET and FLIM–FRET experiments. As in the experiments, we kept one site fixed (amino acid position 221 for NUP98 in the NPC) and swept across the C-terminal residues. For an isolated NUP98 FG (1–499) chain, the results are shown in Supplementary Fig. 9c. We found that the experiments are best explained with a NUP–NUP interaction strength of \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.5, which is mostly in the collapsed state (Supplementary Video 8). According to the radius of gyration of single chains (Supplementary Fig. 9a), the midpoint of the coil–globule transition is at \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44. We conclude that isolated NUP98 FG chains in aqueous solution are mostly collapsed.
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
The NPC scaffold used in the MD simulation is PDB ID: 7R5K. The source data for the main and extended data figures as well as the coordinates for MD simulations are provided in Supplementary Information. All other data are available in the main text or Supplementary Information. Source data are provided with this paper.
Initial configurations and trajectories of the MD simulations are available at https://doi.org/10.5281/zenodo.7648957 under CC-BY license.
Lin, D. H. & Hoelz, A. The structure of the nuclear pore complex (an update). Annu. Rev. Biochem. 88, 725–783 (2019).
Article CAS PubMed PubMed Central Google Scholar
Ribbeck, K. & Görlich, D. Kinetic analysis of translocation through nuclear pore complexes. EMBO J. 20, 1320–1330 (2001).
Article CAS PubMed PubMed Central Google Scholar
Timney, B. L. et al. Simple rules for passive diffusion through the nuclear pore complex. J. Cell Biol. 215, 57–76 (2016).
Article CAS PubMed PubMed Central Google Scholar
Zimmerli, C. E. et al. Nuclear pores dilate and constrict in cellulo. Science 374, eabd9776 (2021).
Article CAS PubMed Google Scholar
Schuller, A. P. et al. The cellular environment shapes the nuclear pore complex architecture. Nature 598, 667–671 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Akey, C. W. et al. Comprehensive structure and functional adaptations of the yeast nuclear pore complex. Cell 185, 361–378 (2022).
Article CAS PubMed PubMed Central Google Scholar
Mosalaganti, S. et al. AI-based structure prediction empowers integrative structural analysis of human nuclear pores. Science 376, eabm9506 (2022).
Article CAS PubMed Google Scholar
Petrovic, S. et al. Architecture of the linker-scaffold in the nuclear pore. Science 376, eabm9798 (2022).
Article CAS PubMed PubMed Central Google Scholar
Zhu, X. et al. Structure of the cytoplasmic ring of the Xenopus laevis nuclear pore complex. Science 376, eabl8280 (2022).
Article CAS PubMed Google Scholar
Bley, C. J. et al. Architecture of the cytoplasmic face of the nuclear pore. Science 376, eabm9129 (2022).
Article CAS PubMed PubMed Central Google Scholar
Fontana, P. et al. Structure of cytoplasmic ring of nuclear pore complex by integrative cryo-EM and AlphaFold. Science 376, abm9326 (2022).
Article Google Scholar
Flory, P. J. The configuration of real polymer chains. J. Chem. Phys. 17, 303–310 (1949).
Article ADS CAS Google Scholar
Dyson, H. J. & Wright, P. E. Intrinsically unstructured proteins and their functions. Nat. Rev. Mol. Cell Biol. 6, 197–208 (2005).
Article CAS PubMed Google Scholar
Oldfield, C. J. & Dunker, A. K. Intrinsically disordered proteins and intrinsically disordered protein regions. Annu. Rev. Biochem. 83, 553–584 (2014).
Article CAS PubMed Google Scholar
Das, R. K., Ruff, K. M. & Pappu, R. V. Relating sequence encoded information to form and function of intrinsically disordered proteins. Curr. Opin. Struct. Biol. 32, 102–112 (2015).
Article CAS PubMed PubMed Central Google Scholar
Hyman, A. A., Weber, C. A. & Julicher, F. Liquid–liquid phase separation in biology. Annu. Rev. Cell Dev. Biol. 30, 39–58 (2014).
Article CAS PubMed Google Scholar
Elbaum-Garfinkle, S. et al. The disordered P granule protein LAF-1 drives phase separation into droplets with tunable viscosity and dynamics. Proc. Natl Acad. Sci. USA 112, 7189–7194 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Folkmann, A. W., Putnam, A., Lee, C. F. & Seydoux, G. Regulation of biomolecular condensates by interfacial protein clusters. Science 373, 1218–1224 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Bremer, A. et al. Deciphering how naturally occurring sequence features impact the phase behaviours of disordered prion-like domains. Nat. Chem. 14, 196–207 (2022).
Article CAS PubMed Google Scholar
Patel, A. et al. A liquid-to-solid phase transition of the ALS protein FUS accelerated by disease mutation. Cell 162, 1066–1077 (2015).
Article CAS PubMed Google Scholar
Bergeron-Sandoval, L.-P. et al. Endocytic proteins with prion-like domains form viscoelastic condensates that enable membrane remodeling. Proc. Natl Acad. Sci. USA 118, e2113789118 (2021).
Article PubMed PubMed Central Google Scholar
Frey, S., Richter, R. P. & Görlich, D. FG-rich repeats of nuclear pore proteins form a three-dimensional meshwork with hydrogel-like properties. Science 314, 815–817 (2006).
Article ADS CAS PubMed Google Scholar
Rout, M. P., Aitchison, J. D., Magnasco, M. O. & Chait, B. T. Virtual gating and nuclear transport: the hole picture. Trends Cell Biol. 13, 622–628 (2003).
Article CAS PubMed Google Scholar
Yamada, J. et al. A bimodal distribution of two distinct categories of intrinsically disordered structures with separate functions in FG nucleoporins. Mol. Cell Proteomics 9, 2205–2224 (2010).
Article PubMed PubMed Central Google Scholar
Kapinos, L. E., Schoch, R. L., Wagner, R. S., Schleicher, K. D. & Lim, R. Y. H. Karyopherin-centric control of nuclear pores based on molecular occupancy and kinetic analysis of multivalent binding with FG nucleoporins. Biophys. J. 106, 1751–1762 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Celetti, G. et al. The liquid state of FG-nucleoporins mimics permeability barrier properties of nuclear pore complexes. J. Cell Biol. 219, e201907157 (2020).
Article PubMed Google Scholar
Najbauer, E. E., Ng, S. C., Griesinger, C., Görlich, D. & Andreas, L. B. Atomic resolution dynamics of cohesive interactions in phase-separated Nup98 FG domains. Nat. Commun. 13, 1494 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Lim, R. Y. H. et al. Nanomechanical basis of selective gating by the nuclear pore complex. Science 318, 640–643 (2007).
Article ADS CAS PubMed Google Scholar
Lerner, E. et al. FRET-based dynamic structural biology: challenges, perspectives and an appeal for open-science practices. eLife 10, e60416 (2021).
Article CAS PubMed PubMed Central Google Scholar
Fuertes, G. et al. Decoupling of size and shape fluctuations in heteropolymeric sequences reconciles discrepancies in SAXS vs. FRET measurements. Proc. Natl Acad. Sci. USA 114, E6342–E6351 (2017).
Article CAS PubMed PubMed Central Google Scholar
König, I. et al. Single-molecule spectroscopy of protein conformational dynamics in live eukaryotic cells. Nat. Methods 12, 773–779 (2015).
Article PubMed Google Scholar
Otsuka, S. et al. A quantitative map of nuclear pore assembly reveals two distinct mechanisms. Nature 613, 575–581 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Liu, S., Hoess, P. & Ries, J. Super-resolution microscopy for structural cell biology. Annu. Rev. Biophys. 51, 301–326 (2022).
Article PubMed Google Scholar
Wu, X. et al. Disruption of the FG nucleoporin NUP98 causes selective changes in nuclear pore complex stoichiometry and function. Proc. Natl Acad. Sci. USA 98, 3191–3196 (2001).
Article ADS CAS PubMed PubMed Central Google Scholar
Ratner, G. A., Hodel, A. E. & Powers, M. A. Molecular determinants of binding between Gly-Leu-Phe-Gly nucleoporins and the nuclear pore complex. J. Biol. Chem. 282, 33968–33976 (2007).
Article CAS PubMed Google Scholar
Hülsmann, B. B., Labokha, A. A. & Görlich, D. The permeability of reconstituted nuclear pores provides direct evidence for the selective phase model. Cell 150, 738–751 (2012).
Article PubMed Google Scholar
Nikić, I. et al. Debugging eukaryotic genetic code expansion for site-specific click-PAINT super-resolution microscopy. Angew. Chem. Int. Ed. 55, 16172–16176 (2016).
Article Google Scholar
Blackman, M. L., Royzen, M. & Fox, J. M. Tetrazine ligation: fast bioconjugation based on inverse-electron-demand Diels–Alder reactivity. J. Am. Chem. Soc. 130, 13518–13519 (2008).
Article CAS PubMed PubMed Central Google Scholar
Reinkemeier, C. D. & Lemke, E. A. Dual film-like organelles enable spatial separation of orthogonal eukaryotic translation. Cell 184, 4886–4903.e21 (2021).
Article CAS PubMed PubMed Central Google Scholar
Adam, S. A., Marr, R. S. & Gerace, L. Nuclear protein import in permeabilized mammalian cells requires soluble cytoplasmic factors. J. Cell Biol. 111, 807–816 (1990).
Article CAS PubMed Google Scholar
Paci, G., Zheng, T., Caria, J., Zilman, A. & Lemke, E. A. Molecular determinants of large cargo transport into the nucleus. eLife 9, e55963 (2020).
Article CAS PubMed PubMed Central Google Scholar
Periasamy, A., Mazumder, N., Sun, Y., Christopher, K. G. & Day, R. N. in Advanced Time-Correlated Single Photon Counting Applications (ed. Becker, W.) 249–276 (Springer, 2015).
Lakowicz, J. R. Principles of Fluorescence Spectroscopy (Springer, 2006).
Redford, G. I. & Clegg, R. M. Polar plot representation for frequency-domain analysis of fluorescence lifetimes. J. Fluoresc. 15, 805–815 (2005).
Article CAS PubMed Google Scholar
Hofmann, H. et al. Polymer scaling laws of unfolded and intrinsically disordered proteins quantified with single-molecule spectroscopy. Proc. Natl Acad. Sci. USA 109, 16155–16160 (2012).
Article ADS CAS PubMed PubMed Central Google Scholar
de Gennes, P.-G. Conformations of polymers attached to an interface. Macromolecules 13, 1069–1075 (1980).
Article ADS Google Scholar
Kohn, J. E. et al. Random-coil behavior and the dimensions of chemically unfolded proteins. Proc. Natl Acad. Sci. USA 101, 12491–12496 (2004).
Article ADS CAS PubMed PubMed Central Google Scholar
Zheng, W. et al. Inferring properties of disordered chains from FRET transfer efficiencies. J. Chem. Phys. 148, 123329 (2018).
Article ADS PubMed PubMed Central Google Scholar
Kalita, J. et al. Karyopherin enrichment and compensation fortifies the nuclear pore complex against nucleocytoplasmic leakage. J. Cell Biol. 221, e202108107 (2022).
Article CAS PubMed PubMed Central Google Scholar
Tan, P. S. et al. Two differential binding mechanisms of FG-nucleoporins and nuclear transport receptors. Cell Rep. 22, 3660–3671 (2018).
Article CAS PubMed PubMed Central Google Scholar
Ibáñez de Opakua, A. et al. Molecular interactions of FG nucleoporin repeats at high resolution. Nat. Chem. 14, 1278–1285 (2022).
Article PubMed PubMed Central Google Scholar
Ng, S. C., Güttler, T. & Görlich, D. Recapitulation of selective nuclear import and export with a perfectly repeated 12mer GLFG peptide. Nat. Commun. 12, 4047 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Kremer, K. & Grest, G. S. Dynamics of entangled linear polymer melts: a molecular‐dynamics simulation. J. Chem. Phys. 92, 5057–5086 (1990).
Article ADS CAS Google Scholar
Schmidt, H. B. & Görlich, D. Nup98 FG domains from diverse species spontaneously phase-separate into particles with nuclear pore-like permselectivity. eLife 4, e04251 (2015).
Article PubMed PubMed Central Google Scholar
Milles, S. et al. Plasticity of an ultrafast interaction between nucleoporins and nuclear transport receptors. Cell 163, 734–745 (2015).
Article CAS PubMed PubMed Central Google Scholar
Raveh, B. et al. Slide-and-exchange mechanism for rapid and selective transport through the nuclear pore complex. Proc. Natl Acad. Sci. USA 113, E2489–E2497 (2016).
Article CAS PubMed PubMed Central Google Scholar
Hough, L. E. et al. The molecular mechanism of nuclear transport revealed by atomic-scale measurements. eLife 4, e10027 (2015).
Article PubMed PubMed Central Google Scholar
Sakiyama, Y., Mazur, A., Kapinos, L. E. & Lim, R. Y. H. Spatiotemporal dynamics of the nuclear pore complex transport barrier resolved by high-speed atomic force microscopy. Nat. Nanotechnol. 11, 719–723 (2016).
Article ADS CAS PubMed Google Scholar
Wang, J. et al. A molecular grammar governing the driving forces for phase separation of prion-like RNA binding proteins. Cell 174, 688–699.e16 (2018).
Article CAS PubMed PubMed Central Google Scholar
Zeng, X., Ruff, K. M. & Pappu, R. V. Competing interactions give rise to two-state behavior and switch-like transitions in charge-rich intrinsically disordered proteins. Proc. Natl Acad. Sci. USA 119, e2200559119 (2022).
Article CAS PubMed PubMed Central Google Scholar
Reinkemeier, C. D., Girona, G. E. & Lemke, E. A. Designer membraneless organelles enable codon reassignment of selected mRNAs in eukaryotes. Science 363, eaaw2644 (2019).
Tesei, G., Schulze, T. K., Crehuet, R. & Lindorff-Larsen, K. Accurate model of liquid–liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties. Proc. Natl Acad. Sci. USA 118, e2111696118 (2021).
Article CAS PubMed PubMed Central Google Scholar
Kapinos, L. E., Huang, B., Rencurel, C. & Lim, R. Y. H. Karyopherins regulate nuclear pore complex barrier and transport function. J. Cell Biol. 216, 3609–3624 (2017).
Article CAS PubMed PubMed Central Google Scholar
Rappaz, B. et al. Measurement of the integral refractive index and dynamic cell morphometry of living cells with digital holographic microscopy. Opt. Express 13, 9361 (2005).
Article ADS PubMed Google Scholar
Ng, S. C. et al. Barrier properties of Nup98 FG phases ruled by FG motif identity and inter-FG spacer length. Nat. Commun. 14, 747 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Szabelski, M. et al. Collisional quenching of erythrosine B as a potential reference dye for impulse response function evaluation. Appl. Spectrosc. 63, 363–368 (2009).
Article ADS CAS PubMed Google Scholar
Griffis, E. R., Altan, N., Lippincott-Schwartz, J. & Powers, M. A. Nup98 is a mobile nucleoporin with transcription-dependent dynamics. Mol. Biol. Cell 13, 16 (2002).
Article Google Scholar
Schrimpf, W., Barth, A., Hendrix, J. & Lamb, D. C. PAM: a framework for integrated analysis of imaging, single-molecule, and ensemble fluorescence data. Biophys. J. 114, 1518–1528 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Koshioka, M., Sasaki, K. & Masuhara, H. Time-dependent fluorescence depolarization analysis in three-dimensional microspectroscopy. Appl. Spectrosc. 49, 224–228 (1995).
Article ADS CAS Google Scholar
Schaffer, J. et al. Identification of single molecules in aqueous solution by time-resolved fluorescence anisotropy. J. Phys. Chem. A 103, 331–336 (1999).
Article CAS Google Scholar
Ranjit, S., Malacrida, L., Jameson, D. M. & Gratton, E. Fit-free analysis of fluorescence lifetime imaging data using the phasor approach. Nat. Protoc. 13, 1979–2004 (2018).
Article CAS PubMed Google Scholar
Hellenkamp, B. et al. Precision and accuracy of single-molecule FRET measurements—a multi-laboratory benchmark study. Nat. Methods 15, 669–676 (2018).
Article CAS PubMed PubMed Central Google Scholar
Hoffmann, A. et al. Mapping protein collapse with single-molecule fluorescence and kinetic synchrotron radiation circular dichroism spectroscopy. Proc. Natl Acad. Sci. USA 104, 105–110 (2007).
Article ADS CAS PubMed Google Scholar
Choi, U. B., McCann, J. J., Weninger, K. R. & Bowen, M. E. Beyond the random coil: stochastic conformational switching in intrinsically disordered proteins. Structure 19, 566–576 (2011).
Article CAS PubMed PubMed Central Google Scholar
Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).
Article ADS CAS MATH Google Scholar
Dünweg, B. & Paul, W. Brownian dynamics simulations without Gaussian random numbers. Int. J. Mod. Phys. C 2, 817–827 (1991).
Article ADS Google Scholar
Monticelli, L. et al. The MARTINI coarse-grained force field: extension to proteins. J. Chem. Theory Comput. 4, 819–834 (2008).
Article CAS PubMed Google Scholar
de Jong, D. H. et al. Improved parameters for the martini coarse-grained protein force field. J. Chem. Theory Comput. 9, 687–697 (2013).
Article PubMed Google Scholar
Benayad, Z., von Bülow, S., Stelzl, L. S. & Hummer, G. Simulation of FUS protein condensates with an adapted coarse-grained model. J. Chem. Theory Comput. 17, 525–537 (2020).
Article PubMed PubMed Central Google Scholar
Berendsen, H. J. C., van der Spoel, D. & van Drunen, R. GROMACS: a message-passing parallel molecular dynamics implementation. Comput. Phys. Commun. 91, 43–56 (1995).
Article ADS CAS Google Scholar
Abraham, M. J. et al. GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1, 19–25 (2015).
Article ADS Google Scholar
Bussi, G., Donadio, D. & Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 126, 14101 (2007).
Article Google Scholar
Berendsen, H. J. C., Postma, J. P. M., Van Gunsteren, W. F., DiNola, A. & Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 3684–3690 (1984).
Article ADS CAS Google Scholar
Parrinello, M. & Rahman, A. Polymorphic transitions in single crystals: a new molecular dynamics method. J. Appl. Phys. 52, 7182–7190 (1981).
Article ADS CAS Google Scholar
Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics. J. Mol. Graphics 14, 33–38 (1996).
Article CAS Google Scholar
Download references
We thank all of the members of the Lemke laboratory for helpful discussions; the core facilities of the Faculty of Biology at Johannes Gutenberg University Mainz; and the Protein Production Core Facility at the Institute of Molecular Biology Mainz for expert assistance. M.Y. was funded by the MSCA Individual Fellowship (TFNUP 89410) and a Humboldt Research Fellowship for Postdoctoral Researchers. E.A.L. acknowledges funding from the ERC-ADG grant ‘MultiOrganelleDesign’, as well as CRC1551 ‘Polymer concepts in cellular function’ of the Deutsche Forschungsgemeinschaft (DFG project number 464588647) and SPP2191 (DFG project number 419070619). E.A.L., G.H. and M.B. acknowledge funding from CRC 1507 ‘Membrane-associated protein assemblies, machineries, and supercomplexes’ of the Deutsche Forschungsgemeinschaft (DFG project number 450648163). M.H., A.O.-K., M.S., M.B. and G.H. thank the Max Planck Society for support.
Marc Siggel
Present address: Centre for Structural Systems Biology, Hamburg, Germany
Marc Siggel
Present address: European Molecular Biology Laboratory Hamburg, Hamburg, Germany
These authors contributed equally: Miao Yu, Maziar Heidari, Sofya Mikhaleva, Piau Siong Tan
Biocenter, Johannes Gutenberg University Mainz, Mainz, Germany
Miao Yu, Sofya Mikhaleva, Sara Mingu, Hao Ruan, Christopher D. Reinkemeier & Edward A. Lemke
Institute of Molecular Biology Mainz, Mainz, Germany
Miao Yu, Sofya Mikhaleva, Sara Mingu, Hao Ruan, Christopher D. Reinkemeier & Edward A. Lemke
Structural and Computational Biology, European Molecular Biology Laboratory, Heidelberg, Germany
Miao Yu, Sofya Mikhaleva, Piau Siong Tan & Christopher D. Reinkemeier
Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Frankfurt am Main, Germany
Maziar Heidari, Marc Siggel & Gerhard Hummer
Department of Molecular Sociology, Max Planck Institute of Biophysics, Frankfurt am Main, Germany
Agnieszka Obarska-Kosinska & Martin Beck
Institute of Biophysics, Goethe University Frankfurt, Frankfurt am Main, Germany
Gerhard Hummer
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
E.A.L. and G.H. conceived the project. M.Y., S. Mikhaleva, P.S.T. and S. Mingu designed and performed the in cell experiments. M.Y., S. Mingu and H.R. designed and performed the in vitro-reconstituted condensates and smFRET experiments. M.H. performed the MD simulations with help from M.S. M.Y. and M.H. analysed the experimental and simulation data, and co-wrote the original draft together with E.A.L. and G.H. C.D.R., A.O.-K. and M.B. contributed analytical tools or reagents. All authors contributed to manuscript editing and approved the final manuscript.
Correspondence to Gerhard Hummer or Edward A. Lemke.
The authors declare no competing interests.
Nature thanks Abhishek Chatterjee and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
(a) Schematic of the OTOs that form a distinct protein translational machinery on the outer mitochondrial membrane surface. The OTOs preferably expanded the genetic codon of the target NUP98, ensuring minimal interference with the endogenous protein translation in the cytoplasm. Non-canonical amino acids (ncAAs) were introduced into NUP98 at the sites specified by amber mutation. The genetically encoded NUP98 with ncAAs incorporated into the nuclear pore complex. The ncAAs were labelled with tetrazine-modified organic dyes via click chemistry. (b) Schematic showing sample preparation pipeline including cell seeding, cell transfection with NUP98 and OTO-GCE plasmids, labelling the incorporated ncAAs with clickable dyes, and fluorescent lifetime imaging of labelled NUP98 on the nuclear rims using the custom-built FLIM-FRET optical setup. COS-7 cells were co-transfected with NUP98A221TAG-boxB plasmid and (c) OTO-GCE system or (d) cytoplasmic GCE system in the presence of trans-cyclooct-2-en-l-lysine (TCO*A, an ncAA that can be labelled) or t-butyloxycarbonyl-l-lysine (BOC, a control ncAA that cannot be labelled), and labelled with LD655-tetrazine. Clearer nuclear rims and less off-target-amber-suppression-induced background labelling were observed in OTO-GCE system compared to the cytoplasmic GCE system. The BOC control groups prove that the observed off-target-labelling was not only due to non-specific sticking of the dye molecules. (Scale bars: 50 μm).
Source Data
The reporter protein was expressed in the presence of trans-cyclooct-2-en-L-lysine (TCO*A, an ncAA that can be labelled) or t-butyloxycarbonyl-L-lysine (BOC, a control ncAA that cannot be labelled). Cells are labelled with the cell-permeable dyes (JF549, JF646, SiR, and TAMRA) under live cell conditions. For the cell-impermeable dyes (AZDye594 and LD655), we permeabilized the plasma membrane with low-dosage digitonin in order to deliver the dyes across the plasma membrane. (Scale bars: 20 μm).
IBB-MPB-GFP cargo in the transport mixture without importin-β and RanGTP was not imported into (a) the labelled COS-7 cells (Scale bars: 20 μm) or (b) the in vitro NUP98 FG-condensates. (Scale bars: 5 μm.) 2 hrs after labelling, the permeabilized COS-7 cells maintained their transport functions in (c) passive exclusion assay (70 kDa dextran was excluded) and (d) facilitated/active transport assay (IBB-MBP-GFP supplied with transport mixture was imported). Immunofluorescence reveals that endogenous importin-β colocalized at the nuclear envelope directly after (e) permeabilization and labelling, and (f) was retained after 2 hrs. For each condition, n = 3 experiments were repeated independently with same conclusion.
Different positions along the FG-domain of NUP98 labelled with the donor dye showed similar (a) fluorescence lifetime and (b) fluorescence anisotropy of the donor dye inside the NPC (n = 9 cells for each mutant), indicating the fluorophores at different labelling sites shared similar fluorescent properties.
Source Data
(a) Normalized total fluorescence intensity profiles of the donor channel I(t) for the eighteen chain segments of NUP98 FG domain for the FLIM-FRET measurements in functional NPCs. Each profile represents an averaged result of ~100 cells. (b) Normalized FRET intensity profiles IFRET(t) by subtracting the fitted donor-only signal IDonly (t) and fitted background signal Ibg (t) from the total signal I(t) in (a) by using Supplementary equation 15. (c) Extended phasor plot for Fig. 2f where the lifetime decays of NUP98 encoded with a nonreactive ncAA BOC are also included. The lifetime decays of the FRET samples could be well separated from the BOC controls. The averaged lifetime decay of the BOC controls was used to determine the cellular background signal for the lifetime fitting procedure.
Source Data
(a) Fluorescent lifetime curves of double-amber-mutated sample NUP98221TAG-251TAG (Nres = 30 aa) and NUP98221TAG-482TAG (Nres = 261 aa) labelled with donor (JF549) and acceptor dye (JF646) mixture. Each profile represents an averaged result of 6 cells. (b) Fluorescent lifetime curves of double-amber-mutated samples same as in (a) but labelled with donor (AZDye594) and acceptor dye (LD655) in permeabilized cells (same dataset as in Extended Data Fig. 5a for the selected mutants).
Source Data
Global fitting was performed by fitting the lifetime decays of all mutants over ~2000 cells with Supplementary equations 15, 16 and 18 for the Gaussian chain model and Supplementary equations 18 and 19 for the SAW-ν model. We obtained ν = 0.56 ± 0.03 (Gaussian chain model) and ν = 0.56 ± 0.001 (SAW-ν model).
Source Data
FRET efficiency histograms for different mutants of purified NUP98 FG-domain in 1x PBS. The peaks in the grey shadow indicate the donor-only population. The FRET population was fitted with an asymmetric-Gaussian function to retrieve the FRET efficiency.
Source Data
(a) Time-lapse images showing phase separation of the purified NUP98 FG domain in vitro. The early stage of the formed condensates demonstrated liquid-like behavior where fast merging events were observed in line with other observations of the early state of NUP98 FG droplets before they harden out to gels26. Scale bar 5 μm. (b) Normalized donor fluorescence intensity profiles of six NUP98 double-mutants in phase-separated condensates that were allowed to settle on the coverslip (n = 5 independent experimental repeats). Nres refers to the number of amino acid residues between the two labelled sites. The normalized donor fluorescence intensity profiles of the single-mutant NUP98A221C labelled with donor only or donor and acceptor mixture showed no difference, indicating no intermolecular FRET was detected and the lifetime differences detected in the double-mutants were due to intramolecular FRET. (c) Phasor plot showing the donor lifetimes of the six double-mutants for individual experimental repeats (n = 5). On the same phasor plot, the lifetimes for the single-mutant labelled with donor-only and a mixture of donor and acceptor were also compared, which showed no difference, ensuring no intermolecular FRET was involved during the measurements.
Source Data
(a) Condensate formation of homopolymer model of NUP98 FG chain at different cohesive interaction strengths \(\mathop{{\epsilon }}\limits^{ \sim }\). Shown are side-views of an elongated simulation box (blue outline) with 500 NUP98-FG (1-499) chains (red). (b) Concentration profiles of the NUP98 FG system at different cohesive interaction strengths. The symbols and error bars represent the average and standard deviation, respectively, of time-averaged density profiles over 3×105τ. For \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42, at the critical interaction strength, the condensate is not stable and the chains are rather uniformly distributed inside the system. (c) Phase diagram of the NUP98 FG condensate in the plane of cohesive interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) and the concentration c of coexisting phases. The solid black line is the result of Flory-Huggins theory for a homopolymer of a length N = 499. The symbols and error bars represent the average and standard deviation, respectively, of time-averaged densities of dense and dilute phases over 3×105τ. The critical cohesive strength is obtained as \(\mathop{{\epsilon }}\limits^{ \sim }\)c ≈ 0.417. For \(\widetilde{{\epsilon }}=0.45\), no chain escaped the condensate during the MD simulation. Based on capillary wave theory, we estimated the dilute-phase density by fitting a double error function to the averaged concentration profile, c(x) = A{erf[(x + B)/w]-erf[(x - B)/w]}+ cdilute, where A = (cdense-cdilute)/2erf(B/w). For \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.45, we obtained cdense = 227.7 mg/mL, cdilute = 0.13 mg/mL, B = 30.69 nm and w = 6.41 nm. The experimental values of the dense and dilute phases are marked in green bars27,54. (d) Root-mean-square inter-residue-distance RE of beads on the same NUP98 FG chain inside the condensate as a function of residue separation Nres. The symbols and error bars represent the average and standard error of the mean, respectively, as estimated from four blocks of size 8×104τ. Issues with condensate aging and gelation during the FLIM-FRET experiments (empty circles: with GLEBS; filled circles: without GLEBS) likely explain the larger distances compared to the simulations.
Source Data
This file contains Supplementary Text, Supplementary Figs. 1–21, Supplementary Tables 1–4 and Supplementary References.
Coordinates for simulation NPC model I
Coordinates for simulation NPC model II
MD simulation of NUP98 FG domain forming condensates. NUP98 FG domain phase-separates and forms condensates for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44.
MD simulation of the FG-NUPs in model II of the nuclear pore complex with interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.35. Upper and lower panels show top view and side views. The scaffold structure known from cryo-ET is rendered in blue. The simulation shows that FG-NUPs (red) are too loose to generate a polymer network in the NPC for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.35. The NUP98 chains are shown in yellow.
MD simulation of the FG-NUPs in model II of the nuclear pore complex with interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42. Upper and lower panels show top view and side views. The scaffold structure known from cryo-ET is rendered in blue. The simulations show that FG-NUPs (red) form extended coil configurations and the inner-ring FG-NUPs fluctuate extensively to form a dynamic barrier across the central channel of the NPC for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42\(.\) The NUP98 chains are shown in yellow.
MD simulation of the FG-NUPs in model II of the nuclear pore complex with interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44. Upper and lower panels show top view and side views. The scaffold structure known from cryo-ET is rendered in blue. The simulations show that FG-NUPs (red) form a surface condensate and collapse onto the NPC scaffold for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44. The NUP98 chains are shown in yellow.
MD simulation of the FG-NUPs in model I of the nuclear pore complex with interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.35. Upper and lower panels show top view and side views. The scaffold structure known from cryo-ET is rendered in blue. The simulation shows that FG-NUPs (red) are too loose to generate a polymer network in the NPC for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.35.
MD simulation of the FG-NUPs in model I of the nuclear pore complex with interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42. Upper and lower panels show top view and side views. The scaffold structure known from cryo-ET is rendered in blue. The simulations show that FG-NUPs (red) form extended coil configurations and the inner-ring FG-NUPs fluctuate extensively to form a dynamic barrier across the central channel of the NPC for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.42.
MD simulation of the FG-NUPs in model I of the nuclear pore complex with interaction strength \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44. Upper and lower panels show top view and side views. The scaffold structure known from cryo-ET is rendered in blue. The simulations show that FG-NUPs (red) form a surface condensate and collapse onto the NPC scaffold for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.44.
MD simulation of a single NUP98 FG-chain. A single NUP98 FG-chain adopts mostly globular-like structure in aqueous solution for \(\mathop{{\epsilon }}\limits^{ \sim }\) = 0.5.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Yu, M., Heidari, M., Mikhaleva, S. et al. Visualizing the disordered nuclear transport machinery in situ. Nature 617, 162–169 (2023). https://doi.org/10.1038/s41586-023-05990-0
Download citation
Received: 28 July 2022
Accepted: 21 March 2023
Published: 26 April 2023
Issue Date: 04 May 2023
DOI: https://doi.org/10.1038/s41586-023-05990-0
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.