Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport

Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport Seok Joo Chae , Dae Wook Kim , View ORCID Profile Oleg A. Igoshin , Seunggyu Lee , View ORCID Profile Jae Kyoung Kim doi: https://doi.org/10.1101/2024.01.12.575351 Seok Joo Chae 1 Department of Mathematical Sciences, KAIST , Daejeon 34141, Republic of Korea 2 Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science , Daejeon 34126, Republic of Korea Find this author on Google Scholar Find this author on PubMed Search for this author on this site Dae Wook Kim 3 Department of Mathematics, University of Michigan , Ann Arbor, MI 48109, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Oleg A. Igoshin 4 Department of Bioengineering, Rice University , Houston, TX 77005, USA 5 Center for Theoretical Biological Physics, Rice University , Houston, TX 77005, USA 6 Department of Chemistry, Rice University , Houston, TX 77005, USA 7 Department of Biosciences, Rice University , Houston, TX 77005, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Oleg A. Igoshin Seunggyu Lee 2 Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science , Daejeon 34126, Republic of Korea 8 Division of Applied Mathematical Sciences, Korea University , Sejong 30019, Republic of Korea Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jae Kyoung Kim 1 Department of Mathematical Sciences, KAIST , Daejeon 34141, Republic of Korea 2 Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science , Daejeon 34126, Republic of Korea Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jae Kyoung Kim For correspondence: jaekkim{at}kaist.ac.kr Abstract Full Text Info/History Metrics Preview PDF Summary All proteins are translated in the cytoplasm, yet many, including transcription factors, play vital roles in the nucleus. While previous research has concentrated on molecular motors for the transport of these proteins to the nucleus, recent observations reveal perinuclear accumulation even in the absence of an energy source, hinting at alternative mechanisms. Here, we propose that structural properties of the cellular environment, specifically the endoplasmic reticulum (ER), can promote molecular transport to the perinucleus without requiring additional energy expenditure. Specifically, physical interaction between proteins and the ER impedes their diffusion and leads to their accumulation near the nucleus. This result explains why larger proteins, more frequently interacting with the ER membrane, tend to accumulate at the perinucleus. Interestingly, such diffusion in a heterogeneous environment follows Chapman’s law rather than the popular Fick’s law. Our findings suggest a novel protein transport mechanism arising solely from characteristics of the intracellular environment. Highlights The interaction of proteins with ER slows down their diffusion at the perinucleus. This leads proteins to migrate toward the perinucleus without ATP consumption. Frequent ER interaction of larger proteins promotes perinuclear accumulation. Diffusion with the physical interaction can be described by Chapman’s law. Download figure Open in new tab Introduction In eukaryotic cells, transcription factors and other proteins that function in the nucleus are translated throughout the cytoplasm but must be translocated across the nuclear membrane 1 - 3 . The translocation of proteins is regulated by the nuclear pore complexes, which are localized in the nuclear envelope 4 , 5 . To reach the vicinity of the envelope, nuclear proteins must traverse the cytoplasm. This mainly occurs via intracellular transport mechanisms facilitated by cytoskeletal filament systems, including microtubules and actin, with the assistance of cytoskeletal motors 6 - 10 . Moreover, actomyosin contraction can also contribute to protein transport by generating cytoplasmic flow 11 . These mechanisms consume adenosine triphosphate (ATP) to fuel motor proteins or contract myosin fibers 6 - 11 . Recent studies have shed light on an as-yet-unidentified alternative mechanism that operates passively (i.e., without consuming ATP). Specifically, exogenous particles like nanoparticles and particulate matter, accumulate in the perinuclear region upon entering cells 12 - 14 . This perinuclear accumulation was also still observed when the energy source (ATP) for the motor proteins was depleted ( Figure 1(i) ) 14 . Moreover, recombinant EGFP-dynactin can migrate toward the nucleus and accumulate in the perinuclear region even without its microtubule binding motif ( Figure 1(ii) ) 15 . Altogether, these studies support the existence of the transport mechanism independent of microtubule-based transport. Interestingly, the accumulation depends on temperature, i.e., decreasing as the temperature decreases ( Figure 1(iii) ) 14 . These characteristics suggest that the unknown mechanism leading to perinuclear accumulation is likely based on diffusion. Notably, diffusion is not homogenous across the cell. In particular, the endoplasmic reticulum (ER), a cellular organelle surrounding the nucleus, is known to impede diffusion 16 . This observation raises the question of whether the heterogeneous diffusive environment in the cytoplasm leads to the migration of proteins toward the nucleus. Download figure Open in new tab Figure 1. Perinuclear accumulation of protein occurs without consuming ATP. Proteins that need to be transported to the nucleus typically rely on motor proteins powered by ATP for their transport. However, surprisingly, exogenous particle can accumulate in the perinuclear region even when the ATP is depleted (i) 14 , suggesting an alternative transport mechanism. Furthermore, disruption of microtubule binding does not impede this accumulation. Specifically, removing the binding motif of dynactin, a protein assisting dynein in transportation, does not hinder its movement toward the nucleus (ii) 15 . Notably, such accumulation was weakened when the temperature was lowered (iii) 14 . This indicates that the accumulation may have a diffusive character. Reprinted by permission from Rockefeller University Press: Journal of Cell Biology 2007 (doi: 10.1083/jcb.200608128) 15 . Our study uncovers that the structural characteristics of the intracellular environment that lead to heterogeneous diffusion of proteins also enable proteins to migrate toward the nucleus. Specifically, our agent-based model (ABM) simulation reveals that such migration occurs due to the ER, which impedes protein diffusion. Notably, our ABM also predicts that larger proteins interact with the ER more frequently, leading to their increased accumulation at the perinucleus compared to smaller ones. Indeed, when we analyzed all the proteins whose sizes and subcellular localization are known simultaneously, we found that perinuclear and nuclear proteins tend to have larger radii than cytoplasmic proteins. This passive mechanism offers an energy-efficient alternative for transporting proteins. Furthermore, we found that the diffusion of proteins in a heterogeneous cellular environment can be effectively described using a partial differential equation (PDE) based on Chapman’s law, rather than the conventional Fick’s law that is typically employed to describe diffusion in the cell. This result suggests that employing Fick’s law to explain diffusion in the presence of physical interactions within a cell might not provide accurate results. Taken together, these findings shift our understanding of protein transport, demonstrating that passive mechanisms can be remarkably efficient and that the cell’s own architecture plays a crucial role in guiding protein movement. Results The cellular environment at the ER leads to perinuclear accumulation The diffusion process within the cytoplasm is not uniform due to the presence of cellular organelles and macromolecules. Among these organelles, the ER surrounds the nucleus with a maze-like structure ( Figure 2A ). A previous study suggested that due to the ER, the movement of particles around the nucleus is hindered, and the diffusion significantly slowed down around the nucleus 16 . Furthermore, the diffusion coefficient was decreased as the ER density increased 16 . Intriguingly, when the ER was disassembled, the diffusion coefficient near the nucleus was comparable to that in the peripheral cytoplasm ( Figure 2B ) 16 . Download figure Open in new tab Figure 2. The endoplasmic reticulum (ER) hampers the diffusion of large-sized protein near the nucleus, resulting in migration toward the nucleus. ( A ) The complex structure of ER can impede the diffusion of protein molecules near the nucleus. ( B ) The diffusion coefficient of a quantum dot near the nucleus was much lower compared to the peripheral cytoplasm. When the ER was degraded, the diffusion coefficients became similar throughout the entire cell. The figure is obtained by analyzing Figures 3 and 5 in a previous study 16 . ( C and D ) To mimic diffusion in the presence of the ER, an agent-based model (ABM) describing intracellular diffusion was utilized. In particular, we located obstacles with narrow intervals (dark grey region) near the nucleus to describe ER (C). When a protein collided with an obstacle, the protein could bind to the obstacle and become immobile, thereby restricting its movement. This restriction in movement affected larger proteins more than smaller ones since larger proteins more frequently collide with the ER (Figure S1) and thus larger proteins more likely to be accumulate near the nucleus (D). In (C), proteins and obstacles are magnified 20 times for better visualization while biologically realistic sizes are used in the simulation for (D) (see STAR Methods for details). ( E ) The perinuclear diffusion coefficient relative to cytoplasm decreases with increasing protein radius. Protein diffusion coefficients were estimated from their mean-squared displacement ( MSD ) using the Einstein relation ( MSD = 4 Dt ). ( F ) Consistent with the ABM simulation, diffusion of larger nanoparticles is more likely to be disrupted by the ER compared to that of smaller particles 21 . ( G ) Aligning with the ABM simulation, proteins which are localized in the perinucleus and the nucleus (blue) have larger sizes, while proteins that are localized in the cytoplasm (red) have smaller sizes. The protein size was determined using the Stokes radii (i.e., hydrodynamic radii) from previous studies 27 , 30 - 33 . To explore how proteins move in a heterogeneous diffusive cytoplasm due to the ER, we utilized an ABM, which has previously been used to track molecule movement in various biological systems 17 - 20 (see STAR Methods for details). Protein molecules were introduced into the cytoplasm and were set to move randomly in different directions at each time step, simulating diffusive movement. Additionally, we incorporated the ER structure surrounding the nucleus, which impeded protein diffusion. We modeled that this hindrance of protein diffusion stems from protein binding to the ER membrane upon collision with it (see STAR Methods for details). Interestingly, ABM simulation revealed that as the protein size increases, more proteins accumulate in the perinuclear region ( Figure 2C ). Specifically, as the radius of protein increases, the portion of protein accumulation at the perinucleus increases ( Figure 2D ). This is because larger proteins collide more frequently with the ER membrane than smaller proteins despite lower diffusion coefficient of larger proteins (Figure S1). Notably, proteins with a radius of ∼10nm exhibit 2.2 times more collisions in the perinuclear region, compared to proteins with a radius of ∼ 2nm (Figure S1). The higher collision frequency leads to a steeper decline in the diffusion coefficient of larger proteins at the perinucleus compared to their smaller counterparts ( Figure 2E ), consistent with a previously observed experiment by Reisch et al.: ER impedes the diffusion of nanoparticles more as their size increases ( Figure 2F ) 21 . Our findings indicate that the ER selectively impedes the diffusion of proteins depending on their sizes, which leads to size-dependent perinuclear accumulation. To further validate the size-dependent perinuclear accumulation with experimental data, we comprehensively analyzed all the proteins with known sizes and intracellular distribution (Table S1). Specifically, we compared the Stokes radius of proteins 22 - 34 based on their subcellular localization, as annotated in previous experiments 35 , 36 . The proteins localized in the perinucleus or the nucleus have Stokes radii larger than proteins localized in the cytoplasm ( Figure 2G ), consistent with the model prediction. An effective PDE can recapitulate the perinuclear accumulation of proteins Even though our ABM explains size-dependent accumulation of proteins in perinuclear region ( Figure 2 ), it is challenging to mathematically analyze it. To enable mathematical analysis, we aimed to convert the ABM to a PDE that accounts for the perinuclear accumulation of proteins without explicitly modeling ER geometry. For this, we utilized a simple microscopic model where proteins freely diffuse in the cytoplasm but become immobile upon binding to cellular components, such as the ER (see STAR Methods for details). Under the assumption that the binding and unbinding is faster than the diffusion (see STAR Methods for details), we can derive that this effective diffusion of proteins impeded by the binding can be described by Chapman’s law, i.e., the probability of finding a protein p ( x, t ) at a given position x and time t is governed by the following equation ( Figure 3A ; see STAR Methods for details): where D ( x ) is the effective diffusion coefficient at x . Note that Chapman’s law is distinct from the commonly used Fick’s law ( Figure 3A ). Download figure Open in new tab Figure 3. A partial differential equation (PDE) based on Chapman’s law captures the movement of proteins toward the nucleus by heterogeneous diffusion. ( A ) The hindered diffusion due to the physical interaction with the ER near the nucleus can be described by the PDE based on Chapman’s law. This is different from the PDE based on Fick’s law, which has been widely utilized to describe diffusion. Here, D ( x ) is the diffusion coefficient at the position x and p ( x, t ) is the probability of finding the protein at the position x and the time t . ( B ) Three different environments of diffusion coefficient. ( C-E ) While the PDE based on Chapman’s law resulted in uniform concentration distribution in the homogeneous environment (C), proteins gathered near the nucleus in heterogeneous environments (D and E). ( F ) Fick’s law, which has been widely used to describe diffusion 60 - 62 , did not exhibit perinuclear accumulation even in the highly heterogeneous environment. ( G ) The flux of Chapman’s law (i) has two components: a drift term ( p ∇ D ) resulting from the gradient of the diffusion coefficient (ii) and a diffusive flux ( D ∇ p ) resulting from the concentration gradient (iii). In the heterogeneous intracellular environment, the drift term (ii) directs proteins toward the nucleus. This leads the diffusive flux (iii) to push proteins outward from the nucleus. While the drift term (ii) remains stronger than the diffusive flux (iii), the protein continues to accumulate near the nucleus (t=0, 50 s ). Eventually, when the drift term (ii) and the diffusive flux (iii) become similar in strength, the net flux reaches the equilibrium, meaning there is no net movement of proteins. (t=100 s ). The effective diffusion coefficient, D ( x ), accounts for how protein diffusion is impeded by physical interactions. In the absence of such interactions, D ( x ) remains constant throughout the cytoplasm, resulting in a homogeneous distribution of D ( x ) ( Figure 3B ). Conversely, when protein interact with the ER near the nucleus, D ( x ) in this region becomes smaller than that of the peripheral cytoplasm, leading to a heterogeneous distribution of D ( x ) across the cytoplasm. Notably, proteins with larger radii more frequently interact with the ER (Figure S1), causing increased hindrance of diffusion compared to proteins with smaller radii ( Figure 2E ). Thus, D ( x ) for larger proteins near the nucleus decreases more than that of smaller proteins. To explore the impact of size-dependent D ( x ) on protein movement and accumulation near the nucleus, we investigated moderately and highly heterogeneous D ( x ) ( Figure 3B ). In the homogeneous environment, perinuclear accumulation did not occur ( Figure 3C ). On the other hand, in heterogeneous environments, proteins gather near the nucleus ( Figures 3D and 3E ), resulting in a higher concentration (i.e., higher p ( x, t ) near the nucleus). Furthermore, p ( x, t ) near the nucleus is higher in the highly heterogeneous environment compared to the moderate heterogeneous environment. Notably, unlike Chapman’s law, the PDE based on conventional Fick’s law cannot capture perinuclear accumulation even in a highly heterogeneous environment ( Figure 3F ). To understand why only Chapman’s law, not Fick’s law, leads to perinuclear accumulation, we analyzed the flux (i.e., the amount of protein flows) of PDEs in the highly heterogeneous environment ( Figure 3G ). In Fick’s law, the flux ( D ∇ p ) is directly proportional to the concentration gradient, causing proteins to migrate from regions of high concentration to regions of low concentration. This flux leads to a steady state with a uniform concentration distribution (Figure S2). In contrast, Chapman’s law incorporates an additional drift term ( p ∇ D ) that arises from the gradient of the effective diffusion coefficient in its flux ( D ∇ p + p ∇ D ) ( Figure 3G(i) ) 37 . Since the drift term heads toward the nucleus, proteins begin to gather near the nucleus from their initial homogeneous concentration distribution ( Figure 3G(ii) ). As proteins gather near the nucleus, the diffusive flux ( D ∇ p ) directs proteins outward from the nucleus ( Figure 3G(iii) ), opposing the drift term. Despite this opposition, the drift term ( p ∇ D ) remains stronger than the diffusive flux ( D ∇ p ), causing the total flux ( D ∇ p + p ∇ D ) to still be directed toward the nucleus ( Figure 3G(i) ). When the diffusive flux cancels out the drift term, the total flux becomes zero, resulting in a steady state with protein accumulation near the nucleus ( Figure 3G(i) ). At the steady state (i.e., D ∇ p + p ∇ D = ∇( Dp ) = 0), p ( x, t ) is inversely proportional to D ( x ), highlighting the crucial role of physical interactions in perinuclear accumulation. Taken together, our findings suggest that Chapman’s law is more suitable than Fick’s law for describing the existing experimental results on diffusion in cellular environments 36 . An ABM, which corresponds to PDE, accurately captures the perinuclear accumulation Although the PDE based on Chapman’s law provides a simple way to account for the perinuclear accumulation, it can become cumbersome for complex cell shapes and intracellular environments. In such cases, using ABM equivalent to PDE based on Chapman’s law can be a versatile alternative. Specifically, incorporating heterogenous effective diffusion in ABM without explicitly modeling binding and unbinding reactions unlike the ABM in Figure 2 can be an efficient way to simulate perinuclear accumulation of proteins within a heterogeneous diffusive environment. Thus, we aimed to explore the conditions under which ABM becomes equivalent to the PDE based on Chapman’s law, enabling the reproduction of perinuclear accumulation of proteins within a heterogeneous diffusive environment ( Figure 3A ). In ABMs, the diffusion of proteins is characterized by a random walk 18 , 20 . The step size of each step of random walk is proportional to the square root of the effective diffusion coefficient 38 . Since the effective diffusion coefficient varies within the cytoplasm ( Figure 3B ), the effective diffusion coefficient can change even during a single step. Thus, different choices of step size can be made based on the position where the effective diffusion coefficient is determined. For example, determination of the effective diffusion coefficient can be based on the current position before the step or the position after the step. The choice of where to evaluate the diffusion coefficient results in distinct types of random walk, leading to distinct dynamics. Notably, when the effective diffusion coefficient is evaluated at the current position, the random walk becomes equivalent to the discretization of the Ito integral, and the dynamics of ABM follows Chapman’s law 39 ( Figure 4A ). Therefore, when this random walk scheme was incorporated into the ABM, perinuclear accumulation was observed over time in a highly heterogeneous diffusion environment ( Figure 4B ). Taken together, ABM utilizing the effective diffusion coefficient of current position of the protein can effectively simulate protein diffusion in heterogeneous intracellular environments to capture perinuclear accumulation of proteins. Download figure Open in new tab Figure 4. ABM scheme that corresponds to Chapman’s law results in the perinuclear accumulation of proteins. ( A ) For the ABM simulation, the step size of the protein movement is determined by the diffusion coefficient. When the diffusion coefficient varies over the cytoplasm, the step size can vary depending on the position at which the diffusion coefficient is obtained. When the step size is determined by the molecule’s current position, the dynamics of proteins is equivalent to Chapman’s law. ( B ) This ABM scheme was simulated in the highly heterogeneous environment, which is used in Figure 3E . Initially (t=0 s ), proteins were uniformly dispersed throughout the cytoplasm. Over time, the proteins gradually accumulated near the nucleus (t=50, 100 s ). ( C ) To simulate the ER partially surrounding the nucleus, we selectively reduced the diffusion coefficient in specific perinuclear regions (blue). Even when the ER is partially surrounding the nucleus, perinuclear accumulation still occurred. ( D ) Besides, ABM provided a convenient approach for describing diffusion in cells of various shapes, including elliptical cells. Proteins again accumulated near the nucleus in elliptical cells. Such ABM framework can depict the diffusion in various intracellular environments. Specifically, the ABM can be easily adapted to describe more realistic but complex environments, such as cells with the ER structure partially surrounding the nucleus ( Figure 4C ), or elliptical-shaped cells ( Figure 4D ) (see STAR Methods for details). Remarkably, even in these complex cell geometries, the ABM simulations replicate perinuclear accumulation ( Figures 4C and 4D ) observed in idealized cells ( Figure 4B ). Discussion To maintain normal physiological processes within cells, it is crucial to effectively and selectively transport to the nucleus proteins translated in the cytoplasm. Our research has shown that proteins can passively diffuse to accumulate near the nucleus without requiring active transport mechanisms, such as transport by motor proteins along microtubules. This accumulation occurs because physical interactions between proteins and the ER slow the diffusion at the perinucleus relative to the peripheral cytoplasm ( Figures 2E and 3B ). Larger proteins are more likely to accumulate at the perinucleus ( Figures 2C and 2D ) because physical interactions occur more frequently as the size of the molecules increases (Figure S1). Our findings reveal that diffusion within a heterogeneous cellular environment enables the ATP-free transport of large-sized proteins, contrasting with the conventional microtubule-based transport, which requires ATP as an energy source 6 , 7 . In addition to the ER, the perinucleus is more crowded than the cell periphery because the perinucleus contains more organelles and denser microtubule networks 40 . This molecular crowding can also impede protein diffusion if proteins can bind to organelles or microtubules. Since physical interactions occur more frequently as the size of the molecules increases (Figure S1), larger proteins are more likely to accumulate at the perinucleus ( Figures 2C and 2D ). This provides insight for previous experimental observations: 1) proteins destined to reach the nucleus often form complexes before translocating to the nucleus 41 , 42 , 2) attaching specific moieties to anticancer drug particles, thereby increasing their size, facilitates targeting to the ER 43 , 44 , and 3) different types of particles, such as lipid nanoparticles and CaCO3, tend to accumulate in the perinuclear region based on their sizes 12 - 14 . While we focused on proteins with sizes comparable to the ER width (∼20 nm), cells also contain significantly larger particles (>>20nm) like vesicles and lysosomes 45 , 46 . These extremely large particles might interact with actin filaments, whose mesh size is approximately 50nm 40 . Investigating the diffusion of these extremely large particles and their potential interactions with the cytoskeleton would be an interesting avenue for future work. The ER consists of two main structural components: a tubular network structure in the cell periphery and a sheet structure near the nucleus. While both structures contribute to the heterogeneity in diffusion coefficients of proteins 16 , their effects differ significantly. The tubular network results in the compartmentalization of the diffusion in the cell periphery, and the diffusion coefficient changes minimally due to a sparse ER density 16 . In contrast, the sheet structure is highly dense near the nucleus, and thus large proteins slow down as they approach the nucleus ( Figure 2E ) 16 . This suggests that the sheet structure is important for helping proteins move towards the nucleus, while the tubular network likely has little effect on translocation. Among various types of PDEs describing diffusion 47 - 55 , it is challenging to determine a priori the most appropriate diffusion law (e.g., Fick’s law or Chapman’s law) for modeling diffusion in heterogeneous environments 56 , 57 . This challenge arises from the diverse physical causes giving rise to heterogeneous diffusion, each manifested in distinct forms of diffusion equations 57 . Consequently, a thorough understanding of the underlying physical phenomena governing heterogeneity is essential. Nevertheless, intracellular diffusion has been primarily described by Fick’s law 58 - 62 . Our study suggests that describing protein diffusion with binding to cellular components like the ER ( Figure 3B ) is equivalent to the “trap model” in Sokolov, resulting in Chapman’s law 57 . When such cellular components impede the diffusion near the nucleus, the PDE based on Chapman’s law but not Fick’s Law yields a non-uniform concentration distribution in a steady state, and accurately describes the perinuclear accumulation ( Figures 3C-3F ). However, mathematical models describing circadian clocks have adopted Fick’s law to describe the diffusion of clock proteins 18 , 59 , although a recent study showed that a clock protein is accumulated at the perinucleus 18 . Thus, utilizing Chapman’s law to develop the circadian clock model would be more appropriate to investigate the spatiotemporal dynamics of molecular circadian clock 18 , 20 , 59 , 63 - 65 . Furthermore, it would be an interesting future work to investigate dynamics of various biological oscillators based on transcriptional feedback loops requiring nuclear transport of transcriptional factors (e.g., p53 oscillator 66 - 68 , NF- k b oscillator 69 , and synthetic oscillator 70 - 73 ) when Chapman’s law is used rather than Fick’s law. Previous studies of diffusion in heterogeneous environments, particularly focusing on highly crowded cellular environments, have highlighted the prevalence of subdiffusion 74 - 76 for various particles, including small labeled proteins, lipid and insulin granules, and virus particle 77 - 80 . This subdiffusive behavior, characterized by slower-than-linear growth in displacement variance, suggests hindered transport and potential inefficiencies in cellular processes. On the other hand, our study demonstrates how heterogeneous diffusion, specifically under the influence of physical interactions with the ER, can facilitate directional protein transport toward the nucleus, as implied by Chapman’s law. Consistently, heterogeneous diffusion was known to be beneficial for reducing the first-passage time of reaction events, even though it slows down the overall dynamics of proteins 81 . Taken together, diffusion in heterogeneous environments can also contribute to efficient cellular functioning. Limitations of the study We have comprehensively examined all the proteins for which both their sizes and subcellular localizations (i.e., locations of proteins) are known, and this analysis supports that larger proteins tend to accumulate near the nucleus ( Figure 2G ). However, the scarcity of proteins with known sizes and subcellular localizations restricts the generalizability of our conclusions. It would be an important topic to explore how subcellular localization and proteins sizes are correlated across a broader range of proteins. Moreover, it is crucial to verify whether a protein does not contain NLS when examining its subcellular localization, since proteins and particles with NLS are reported to induce microtubule-based transport 82 , 83 . To further strengthen our conclusion, it would also be helpful to investigate whether smaller proteins can accumulate near the nucleus after fusing with fluorescent proteins, as these modifications can significantly increase the size of protein 84 . Author Contributions S.J.C., D.W.K., S.L. and J.K.K. designed the study. S.J.C. performed computational modeling and simulation and all authors contributed to analysis. J.K.K. supervised the project. S.J.C. and J.K.K. wrote the manuscript draft, and all authors revised the manuscript. Declaration of interests The authors declare no competing interests. STAR★Methods RESOURCE AVAILABILITY Lead contact Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Jae Kyoung Kim ( jaekkim{at}kaist.ac.kr ). Materials availability This study did not generate new unique reagents. Data and code availability - This paper analyzes publicly available existing data. These accession numbers for the datasets are listed in the key resources table. - The MATLAB and Netlogo codes of the computational package are available in the following Database: The link will be available upon acceptance of the manuscript. - Any additional information required to analyze the data is available from the lead contact upon request. METHOD DETAILS ABM describing diffusion in the presence of the ER To describe the heterogeneous diffusion in a cell, we utilized an ABM which describes the intracellular diffusion of proteins. Following previous studies 18 , 20 , we simplified the cell geometry as a two-dimensional circular shape with the radius of 10µm 85 . Within the cell, a circular nucleus resided, sharing the same center and possessing a radius that is one-third of the cell radius (i.e., 3.33µm), so that the nucleus occupies about 10% of the cell area 18 . Furthermore, we incorporated the ER by introducing obstacles near the nucleus ( Figure 2C ). These obstacles were distributed across 80 concentric layers. Starting from 3.58µm, each layer increased in radius by 0.25µm, reflecting the width of the ER sheet structure 85 . Within each layer, 1,000 uniformly distributed obstacles were uniformly distributed to mimic increased ER density as proteins approached the nucleus. Each obstacle was modeled as a circle with a 0.0025µm radius. A protein molecule was represented as a circle with the radius r . In our simulation, we varied r between 2 nm, corresponding to the radius of carbonic anhydrase 33 , to 10 nm, corresponding to the PER2 complex 86 . The diffusion of proteins was described by using a random walk. During each time step of length Δt , each molecule was moved in a randomly chosen direction with a step size δ . The angles of movement were randomly sampled from the uniform distribution, ranging from 0° to 360°. The step size δ was chosen to be sufficiently smaller than the protein size to ensure precise and fine-scale movements. For example, for a protein with 10nm radius, δ was set to be 1.6nm. For proteins with smaller radii, the step size was adjusted by using the Stokes-Einstein equation (i.e., and thus ). The time step length ( Δt ) was determined by using the Einstein relation δ 2 = 4DΔt 38 , where D is the diffusion coefficient. In this study, we utilized the known diffusion coefficient of PER2 protein ( D= 0.2µ m 2 /s ) 87 and its radius of 10nm 86 , (i.e., δ = 1.6nm), resulting in Δt = 0.0000032 s . When proteins collided to an obstacle, they were repositioned to the point of collision. Subsequently, they bound to the obstacle with a binding probability of 1 − exp(− k bind Δt ), where k bind is the binding rate. Once bound, at each time step, proteins could unbind with a probability of 1 − exp(− k unbind Δt ), where k unbind is the unbinding rate. In our simulation, k bind and k unbind were set to 1,562,500 /s and 625 /s , respectively. We implemented reflecting boundary conditions to prevent the protein from leaving the cytoplasm. That is, we repositioned proteins as if they were elastically reflected by the boundary (cell membrane or the nuclear membrane). Specifically, proteins exceeding the cytoplasmic boundary (cytoplasmic membrane or nuclear membrane) are mirrored across the tangent line of the boundary at the collision point with the boundary. Additionally, for the initial condition, protein began diffusion from a uniform distribution of the protein throughout the cytoplasm. Physical interpretation of Chapman’s law Since our mechanism of protein transport operates without the need for energy sources, such as ATP hydrolysis, we expect that steady state distribution will correspond to the thermodynamic equilibrium. This steady state with non-uniform concentration distribution should result from the energy difference between spatial locations. We postulate that the difference in energy level comes from the physical interaction between proteins and the cellular organelles, including the ER. To investigate this, we considered a one-dimensional diffusion of protein in the cytoplasm. Specifically, we considered a protein moving along a linear chain of N sites ( i = 1 ∼N ), where the chain represents the cytoplasm and each site is separated by an interval h . The protein located at the site i can interact with cellular organelles, such as the ER, with a binding rate of and it can also detach from the organelles with an unbinding rate of . If the protein at the site i is not bound, it can diffuse to the adjacent site (e.g., i + 1 or i − 1) with a diffusion coefficient of D . Denoting the probability of finding the bounded and unbounded protein at the site i at time t by and , respectively, we obtain the following ordinary differential equation (ODE): Here, we use the convention of and by imposing Neumann boundary conditions. Since we are primarily interested in the probability of finding proteins at the i-th site, regardless of whether they are bound or not, we introduce a new variable the probability of finding the protein at the site i, i.e., . Then, Eq. (1) can be written as shown below: We assume that binding and unbinding reactions occur significantly faster than diffusion. Then is affected by only slow diffusion, whereas is affected by fast binding and unbinding reactions, leading to a complete time scale separation. This allows us to assume that reaches the quasi-equilibrium (the fast binding and unbinding equilibrium) at each site before changes appreciably (i.e., for each i ) 88 - 92 . Then, can be expressed by Consequently, follows a reduced ODE shown below: where Taking the limit of infinitesimal interval ( h → 0 ), we obtain the following Chapman’s law: where To validate the model reduction, we compared whether the original model ( Eq. (1) ) and the reduced model exhibit similar behavior. As long as the binding and unbinding occur significantly faster than the diffusion, the two models exhibit similar dynamics (Figure S3). Numerical simulation for the PDE After deriving the PDE based on Chapman’s law in the previous section, we solved the PDE numerically. Depicting diffusion within a cell accurately using the standard finite difference method on a rectangular domain presents challenges because cells usually have a non-rectangular shape and contain a nucleus where diffusion does not occur. To address this issue, we adopted an approach from a previous study which described the diffusion under a domain with a complex geometry 93 . Specifically, we solved the PDE based on Chapman’s law in a rectangular domain Ω = [− L, L ] × [− L, L ], where, L was set to be 1.1 ⋅ R c , where R c is the radius of the cell. By letting we obtained a uniform mesh with a set of grid edges as Ω h ={( x i , y j ): 1 ≤ i, j ≤ N }, where x i = − L + ( i – 1) h and y j = − L + ( j − 1) h . Here, we used N = 100. At the edge ( x i , y j ), the boundary control function G -i,j was defined as where R N is the radius of the nucleus. In the grid edge, G is defined as and Utilizing G allowed us to utilize the finite difference method in a complex domain mimicking the cell by simply computing PDE in a rectangular domain 93 . We denote the discrete approximation p ( x i , y j , nΔt ) as . Then, we can define the discrete differential operator ∇ d as follows: Furthermore, to calculate ∇ d at the boundary, the value of at ghost points located outside of the domain also required ( i, j = 0 or N + 1). To determine values out of the domain, we utilized the Neumann boundary condition. Specifically, we set for i = 0, N + 1 and for j = 0, N + 1. The boundary condition and boundary control function G ensures the conservation of p ( x, t ) over time, i.e., ∫ p ( x, t ) dx remains constant over time 94 . Finally, the PDE based on Chapman’s law can be written as the following: for 1 ≤ i, j ≤ N , where h =0.2µm and Δt = 0.00004s were used. To determine G , the same value of R N and R C as in the ABM are used. Using the PDE ( Eq. (6) ), we obtained the solution of Chapman’s law in Figure 3 . To describe the diffusion with the existence of physical interactions, we have utilized various effective diffusion coefficients at ( x, y ), D ( x, y ). A homogeneous case described by using D ( x ) = 0 .2µm 2 /s . For moderately and highly heterogeneous cases, we employed functions that decrease as the protein approaches the nucleus. Specifically, we used D ( x, y ) = 0.2 · and to describe moderately heterogeneous and highly heterogeneous cases, respectively. Here, r N is the radius of the nucleus, and ‖( x, y )‖ is the distance between the protein and the cell center. Simulating ABM based on Chapman’s law To simulate an ABM equivalent to Chapman’s law, we introduced a variation in the protein’s movement ( Figures 4B-4D ). Specifically, the protein molecule at position ( x, y ) moves by a step size δ ( x, y ) in a random direction. The step size of the protein was determined by the diffusion coefficient of the current position. We modeled the decrease in the step size by using a Hill function: µm, where r N is the radius of the nucleus, and ‖( x, y )‖ is the distance between the protein and the cell center. This choice of step size of the same form of Hill function to PDE ( δ 2 α D ) ensures that the ABM aligns with the highly heterogeneous environment. As the protein approaches the nucleus, the step size decreases to be 10% of the maximum value (=0.02µm). The length of time of step Δ t was obtained by using the relation , where δ M is the maximum value of the step size and D is the diffusion coefficient of PER2 protein ( D= 0.2µm 2 /s ) 87 . In Figures 4C and 4D , we have extended the ABM to describe cells with diverse geometries. To model the ER structure partially surrounding the nucleus in Figure 4C , we implemented the following step size function: . where O = { θ : 6 0 ° ≤ θ< 21 0 ° or 27 0 ° ≤ θ< 36 0 °} is and 1 0 ( x) is the indicator function of O . In Figure 4D , a cell was represented by an ellipse with a semi-major axis of 15µm and a semi-minor axis of 10µm, centered at the origin ( Figure 4D ). The nucleus was depicted as a circle with a radius of 3.33µm, centered at (−5µm, 0µm). Acknowledgements We thank Yong Jung Kim and Ho-Youn Kim for the helpful discussion. We thank Life Science Editors for editorial assistance. This work was supported by the Institute for Basic Science (grant no. IBS-R029-C3) (J.K.K.) and Welch Foundation grant no. C-1995 (O.A.I.). Some parts of figure 2 were retrieved from Biorender. Footnotes ↵ 9 Lead Contact Error in figure is revised References 1. ↵ Tsai , Y.J. , Lee , H.I. , and Lin , A. ( 2012 ). Ribosome distribution in HeLa cells during the cell cycle . PLoS One 7 , e32820 . doi: 10.1371/journal.pone.0032820 . OpenUrl CrossRef PubMed 2. Kim , Y.H. , Han , M.E. , and Oh , S.O. ( 2017 ). The molecular mechanism for nuclear transport and its application . Anat Cell Biol 50 , 77 – 85 . doi: 10.5115/acb.2017.50.2.77 . OpenUrl CrossRef PubMed 3. ↵ Lu , J. , Wu , T. , Zhang , B. , Liu , S. , Song , W. , Qiao , J. , and Ruan , H. ( 2021 ). Types of nuclear localization signals and mechanisms of protein import into the nucleus . Cell Communication and Signaling 19 , 60 . doi: 10.1186/s12964-021-00741-y . OpenUrl CrossRef 4. ↵ Beck , M. , and Hurt , E. ( 2017 ). The nuclear pore complex: understanding its function through structural insight . Nat Rev Mol Cell Biol 18 , 73 – 89 . doi: 10.1038/nrm.2016.147 . OpenUrl CrossRef PubMed 5. ↵ Lin , D.H. , and Hoelz , A. ( 2019 ). The Structure of the Nuclear Pore Complex (An Update) . Annu Rev Biochem 88 , 725 – 783 . doi: 10.1146/annurev-biochem-062917-011901 . OpenUrl CrossRef PubMed 6. ↵ Reck-Peterson , S.L. , Redwine , W.B. , Vale , R.D. , and Carter , A.P. ( 2018 ). The cytoplasmic dynein transport machinery and its many cargoes . Nat Rev Mol Cell Biol 19 , 382 – 398 . doi: 10.1038/s41580-018-0004-3 . OpenUrl CrossRef PubMed 7. ↵ Sweeney , H.L. , and Holzbaur , E.L.F. ( 2018 ). Motor Proteins . Cold Spring Harb Perspect Biol 10 . doi: 10.1101/cshperspect.a021931 . OpenUrl Abstract / FREE Full Text 8. Olenick , M.A. , and Holzbaur , E.L.F. ( 2019 ). Dynein activators and adaptors at a glance . Journal of Cell Science 132 . doi: 10.1242/jcs.227132 . OpenUrl Abstract / FREE Full Text 9. Xiang , X. , and Qiu , R. ( 2020 ). Cargo-Mediated Activation of Cytoplasmic Dynein in vivo . Front Cell Dev Biol 8 , 598952 . doi: 10.3389/fcell.2020.598952 . OpenUrl CrossRef 10. ↵ Egan , Martin J. , McClintock Mark A. , Hollyer Ian H.L. , Elliott Hunter L. , and Reck-Peterson Samara L. ( 2015 ). Cytoplasmic Dynein Is Required for the Spatial Organization of Protein Aggregates in Filamentous Fungi . Cell Reports 11 , 201 – 209 . doi: 10.1016/j.celrep.2015.03.028 . OpenUrl CrossRef PubMed 11. ↵ Murrell , M. , Oakes , P.W. , Lenz , M. , and Gardel , M.L. ( 2015 ). Forcing cells into shape: the mechanics of actomyosin contractility . Nature Reviews Molecular Cell Biology 16 , 486 – 498 . doi: 10.1038/nrm4012 . OpenUrl CrossRef PubMed 12. ↵ Chakraborty , S.K. , Fitzpatrick , J.A.J. , Phillippi , J.A. , Andreko , S. , Waggoner , A.S. , Bruchez , M.P. , and Ballou , B. ( 2007 ). Cholera Toxin B Conjugated Quantum Dots for Live Cell Labeling . Nano Letters 7 , 2618 – 2626 . doi: 10.1021/nl0709930 . OpenUrl CrossRef PubMed 13. Wei , W. , Ma , G.-H. , Hu , G. , Yu , D. , McLeish , T. , Su , Z.-G. , and Shen , Z.-Y. ( 2008 ). Preparation of Hierarchical Hollow CaCO3 Particles and the Application as Anticancer Drug Carrier . Journal of the American Chemical Society 130 , 15808 – 15810 . doi: 10.1021/ja8039585 . OpenUrl CrossRef PubMed 14. ↵ Rivolta , I. , Panariti , A. , Collini , M. , Lettiero , B. , D’Alfonso , L. , Sironi , L. , Miserocchi , G. , and Chirico , G. ( 2011 ). A biophysical model of intracellular distribution and perinuclear accumulation of particulate matter . Biophysical Chemistry 158 , 134 – 140 . doi: 10.1016/j.bpc.2011.06.009 . OpenUrl CrossRef PubMed 15. ↵ Kim , H. , Ling , S.-C. , Rogers , G.C. , Kural , C. , Selvin , P.R. , Rogers , S.L. , and Gelfand , V.I. ( 2007 ). Microtubule binding by dynactin is required for microtubule organization but not cargo transport . Journal of Cell Biology 176 , 641 – 651 . doi: 10.1083/jcb.200608128 . OpenUrl Abstract / FREE Full Text 16. ↵ Li , H. , Dou , S.-X. , Liu , Y.-R. , Li , W. , Xie , P. , Wang , W.-C. , and Wang , P.-Y. ( 2015 ). Mapping Intracellular Diffusion Distribution Using Single Quantum Dot Tracking: Compartmentalized Diffusion Defined by Endoplasmic Reticulum . Journal of the American Chemical Society 137 , 436 – 444 . doi: 10.1021/ja511273c . OpenUrl CrossRef 17. ↵ Longo , D. , Peirce , S.M. , Skalak , T.C. , Davidson , L. , Marsden , M. , Dzamba , B. , and DeSimone , D.W. ( 2004 ). Multicellular computer simulation of morphogenesis: blastocoel roof thinning and matrix assembly in Xenopus laevis . Dev Biol 271 , 210 – 222 . doi: 10.1016/j.ydbio.2004.03.021 . OpenUrl CrossRef PubMed Web of Science 18. ↵ Beesley , S. , Kim , D.W. , D’Alessandro , M. , Jin , Y. , Lee , K. , Joo , H. , Young , Y. , Tomko , R.J. , Faulkner , J. , Gamsby , J. , et al. ( 2020 ). Wake-sleep cycles are severely disrupted by diseases affecting cytoplasmic homeostasis . Proceedings of the National Academy of Sciences 117 , 28402 – 28411 . doi: 10.1073/pnas.2003524117 . OpenUrl Abstract / FREE Full Text 19. K.N. Das , J.C. Bansal , K. Deep , A.K. Nagar , P. Pathipooranam , and R.C. Naidu Shinde , S.B. , and Kurhekar , M.P. ( 2020 ). Agent-Based Modeling of the Adaptive Immune System Using Netlogo Simulation Tool . held in Singapore, 2020// . K.N. Das , J.C. Bansal , K. Deep , A.K. Nagar , P. Pathipooranam , and R.C. Naidu , eds. ( Springer Singapore ), pp. 463 – 474 . 20. ↵ Chae , S.J. , Kim , D.W. , Lee , S. , and Kim , J.K. ( 2023 ). Spatially coordinated collective phosphorylation filters spatiotemporal noises for precise circadian timekeeping . iScience 26 , 106554 . doi: 10.1016/j.isci.2023.106554 . OpenUrl CrossRef 21. ↵ Reisch , A. , Heimburger , D. , Ernst , P. , Runser , A. , Didier , P. , Dujardin , D. , and Klymchenko , A.S. ( 2018 ). Protein-Sized Dye-Loaded Polymer Nanoparticles for Free Particle Diffusion in Cytosol . Advanced Functional Materials 28 , 1805157 . doi: 10.1002/adfm.201805157 . OpenUrl CrossRef 22. ↵ Kushner , I. , and Somerville , J.A. ( 1970 ). Estimation of the molecular size of C-reactive protein and Cx-reactive protein in serum . Biochimica et Biophysica Acta (BBA)-Protein Structure 207 , 105 – 114 . doi: 10.1016/0005-2795(70)90140-6 . OpenUrl CrossRef PubMed 23. Morgenstern , R. , Guthenberg , C. , and Depierre , J.W. ( 1982 ). Microsomal Glutathione S-Transferase . European Journal of Biochemistry 128 , 243 – 248 . doi: 10.1111/j.1432-1033.1982.tb06958.x . OpenUrl CrossRef PubMed Web of Science 24. Randahl , H. , Slabý , I. , and Holmgren , A. ( 1982 ). An Improved Purification Method and a Physical Characterization of Phage T7 DNA Polymerase . European Journal of Biochemistry 128 , 445 – 449 . doi: 10.1111/j.1432-1033.1982.tb06984.x . OpenUrl CrossRef PubMed Web of Science 25. Triebel , H. , Bär , H. , Walter , A. , Osipova , T.N. , Ramm , E.I. , Kostyleva , E.I. , and Vorob’ev , V.I. ( 1989 ). Structural differences between histone H1 molecules from sea urchin (Strongylocentrotus intermedius) sperm and calf thymus: hydrodynamic and c.d . studies. International Journal of Biological Macromolecules 11 , 153 – 158 . doi: 10.1016/0141-8130(89)90060-3 . OpenUrl CrossRef 26. Tomkinson , A.E. , Lasko , D.D. , Daly , G. , and Lindahl , T. ( 1990 ). Mammalian DNA ligases . Catalytic domain and size of DNA ligase I. Journal of Biological Chemistry 265 , 12611 – 12617 . doi: 10.1016/S0021-9258(19)38387-5 . OpenUrl CrossRef 27. ↵ Forsberg , E.C. , Lam , L.T. , Yang , X.J. , Nakatani , Y. , and Bresnick , E.H. ( 1997 ). Human histone acetyltransferase GCN5 exists in a stable macromolecular complex lacking the adapter ADA2 . Biochemistry 36 , 15918 – 15924 . doi: 10.1021/bi971664x . OpenUrl CrossRef PubMed Web of Science 28. Raab-Graham , K.F. , and Vandenberg , C.A. ( 1998 ). Tetrameric Subunit Structure of the Native Brain Inwardly Rectifying Potassium Channel Kir 2.2 * . Journal of Biological Chemistry 273 , 19699 – 19707 . doi: 10.1074/jbc.273.31.19699 . OpenUrl Abstract / FREE Full Text 29. Adamczyk , Z. , Cichocki , B. , Ekiel-Jeżewska , M.L. , Słowicka , A. , Wajnryb , E. , and Wasilewska , M. ( 2012 ). Fibrinogen conformations and charge in electrolyte solutions derived from DLS and dynamic viscosity measurements . Journal of Colloid and Interface Science 385 , 244 – 257 . doi: 10.1016/j.jcis.2012.07.010 . OpenUrl CrossRef PubMed 30. ↵ Rajani , K.R. , Pettit Kneller , E.L. , McKenzie , M.O. , Horita , D.A. , Chou , J.W. , and Lyles , D.S. ( 2012 ). Complexes of vesicular stomatitis virus matrix protein with host Rae1 and Nup98 involved in inhibition of host transcription . PLoS Pathog 8 , e1002929 . doi: 10.1371/journal.ppat.1002929 . OpenUrl CrossRef PubMed 31. Moor , N.A. , Vasil’eva , I.A. , Anarbaev , R.O. , Antson , A.A. , and Lavrik , O.I. ( 2015 ). Quantitative characterization of protein-protein complexes involved in base excision DNA repair . Nucleic Acids Res 43 , 6009 – 6022 . doi: 10.1093/nar/gkv569 . OpenUrl CrossRef PubMed 32. Ritterhoff , T. , Das , H. , Hofhaus , G. , Schröder , R.R. , Flotho , A. , and Melchior , F. ( 2016 ). The RanBP2/RanGAP1 * SUMO1/Ubc9 SUMO E3 ligase is a disassembly machine for Crm1-dependent nuclear export complexes . Nature Communications 7 , 11482 . doi: 10.1038/ncomms11482 . OpenUrl CrossRef PubMed 33. ↵ La Verde , V. , Dominici , P. , and Astegno , A. ( 2017 ). Determination of Hydrodynamic Radius of Proteins by Size Exclusion Chromatography . Bio Protoc 7 , e2230 . doi: 10.21769/BioProtoc.2230 . OpenUrl CrossRef 34. ↵ Kilgore , H.R. , Latham , A.P. , Ressler , V.T. , Zhang , B. , and Raines , R.T. ( 2020 ). Structure and Dynamics of N-Glycosylated Human Ribonuclease 1 . Biochemistry 59 , 3148 – 3156 . doi: 10.1021/acs.biochem.0c00191 . OpenUrl CrossRef 35. ↵ Rastogi , S. , and Rost , B. ( 2011 ). LocDB: experimental annotations of localization for Homo sapiens and Arabidopsis thaliana . Nucleic Acids Res 39 , D230 – 234 . doi: 10.1093/nar/gkq927 . OpenUrl CrossRef PubMed Web of Science 36. ↵ Shaiken , T.E. , and Opekun , A.R. ( 2014 ). Dissecting the cell to nucleus, perinucleus and cytosol . Scientific Reports 4 , 4923 . doi: 10.1038/srep04923 . OpenUrl CrossRef 37. ↵ Mehrer , H. ( 2007 ). Diffusion in solids: fundamentals, methods, materials, diffusioncontrolled processes ( Springer Science & Business Media ). 38. ↵ Berg , H.C. ( 1993 ). Random Walks in Biology ( Princeton University Press ). 39. ↵ Bressloff , P.C. ( 2014 ). Stochastic Processes in Cell Biology . Interdisciplinary Applied Mathematics ,. 1st ed. Springer International Publishing : Imprint: Springer ,. 40. ↵ Jiang , C. , Yang , M. , Li , W. , Dou , S.X. , Wang , P.Y. , and Li , H. ( 2022 ). Spatiotemporal three-dimensional transport dynamics of endocytic cargos and their physical regulations in cells . iScience 25 , 104210 . doi: 10.1016/j.isci.2022.104210 . OpenUrl CrossRef 41. ↵ Takahashi , J.S. ( 2017 ). Transcriptional architecture of the mammalian circadian clock . Nat Rev Genet 18 , 164 – 179 . doi: 10.1038/nrg.2016.150 . OpenUrl CrossRef PubMed 42. ↵ Wintrebert , M. , Nguyen , M.C. , and Smith , G.R. ( 2021 ). Activation of meiotic recombination by nuclear import of the DNA break hotspot-determining complex in fission yeast . J Cell Sci 134 . doi: 10.1242/jcs.253518 . OpenUrl Abstract / FREE Full Text 43. ↵ Ghosh , C. , Nandi , A. , and Basu , S. ( 2019 ). Lipid Nanoparticle-Mediated Induction of Endoplasmic Reticulum Stress in Cancer Cells . ACS Applied Bio Materials 2 , 3992 – 4001 . doi: 10.1021/acsabm.9b00532 . OpenUrl CrossRef 44. ↵ Liu , Y. , Jia , H.-R. , Han , X. , and Wu , F.-G. ( 2021 ). Endoplasmic reticulum-targeting nanomedicines for cancer therapy . Smart Materials in Medicine 2 , 334 – 349 . doi: 10.1016/j.smaim.2021.09.001 . OpenUrl CrossRef 45. ↵ Huang , C. , Quinn , D. , Sadovsky , Y. , Suresh , S. , and Hsia , K.J. ( 2017 ). Formation and size distribution of self-assembled vesicles . Proceedings of the National Academy of Sciences 114 , 2910 – 2915 . doi: 10.1073/pnas.1702065114 . OpenUrl Abstract / FREE Full Text 46. ↵ de Araujo , M.E.G. , Liebscher , G. , Hess , M.W. , and Huber , L.A. ( 2020 ). Lysosomal size matters . Traffic 21 , 60 – 75 . doi: 10.1111/tra.12714 . OpenUrl CrossRef PubMed 47. ↵ Landsberg , P.T. ( 1984 ). D grad v or grad(Dv)? Journal of Applied Physics 56 , 1119 – 1122 . doi: 10.1063/1.334083 . OpenUrl CrossRef 48. Van Kampen , N.G. ( 1988 ). Diffusion in inhomogeneous media . Journal of Physics and Chemistry of Solids 49 , 673 – 677 . doi: 10.1016/0022-3697(88)90199-0 . OpenUrl CrossRef 49. Lançon , P. , Batrouni , G. , Lobry , L. , and Ostrowsky , N. ( 2002 ). Brownian walker in a confined geometry leading to a space-dependent diffusion coefficient . Physica A: Statistical Mechanics and its Applications 304 , 65 – 76 . doi: 10.1016/S0378-4371(01)00510-6 . OpenUrl CrossRef 50. van Milligen , B.P. , Bons , P.D. , Carreras , B.A. , and Sánchez , R. ( 2005 ). On the applicability of Fick’s law to diffusion in inhomogeneous systems . European Journal of Physics 26 , 913 . doi: 10.1088/0143-0807/26/5/023 . OpenUrl CrossRef 51. Sattin , F. ( 2008 ). Fick’s law and Fokker–Planck equation in inhomogeneous environments . Physics Letters A 372 , 3941 – 3945 . doi: 10.1016/j.physleta.2008.03.014 . OpenUrl CrossRef 52. Huang , F. , Chakraborty , P. , Lundstrom , C.C. , Holmden , C. , Glessner , J.J. , Kieffer , S.W. , and Lesher , C.E. ( 2010 ). Isotope fractionation in silicate melts by thermal diffusion . Nature 464 , 396 – 400 . doi: 10.1038/nature08840 . OpenUrl CrossRef GeoRef PubMed Web of Science 53. Choi , B. , and Kim , Y.-J. ( 2019 ). Diffusion of Biological Organisms: Fickian and Fokker--Planck Type Diffusions . SIAM Journal on Applied Mathematics 79 , 1501 – 1527 . doi: 10.1137/18m1163944 . OpenUrl CrossRef 54. Won , Y.-Y. , and Ramkrishna , D. ( 2019 ). Revised Formulation of Fick’s, Fourier’s, and Newton’s Laws for Spatially Varying Linear Transport Coefficients . ACS Omega 4 , 11215 – 11222 . doi: 10.1021/acsomega.9b00736 . OpenUrl CrossRef 55. ↵ Kim , Y.-J. , and Seo , H. ( 2021 ). Model for Heterogeneous Diffusion . SIAM Journal on Applied Mathematics 81 , 335 – 354 . doi: 10.1137/19M130087X . OpenUrl CrossRef 56. ↵ van Kampen , N.G. ( 1981 ). Itô versus Stratonovich . Journal of Statistical Physics 24 , 175 – 187 . doi: 10.1007/BF01007642 . OpenUrl CrossRef 57. ↵ Sokolov , I.M. ( 2010 ). Ito, Stratonovich, Hänggi and all the rest: The thermodynamics of interpretation . Chemical Physics 375 , 359 – 363 . doi: 10.1016/j.chemphys.2010.07.024 . OpenUrl CrossRef 58. ↵ Bressloff , P.C. , and Newby , J.M. ( 2013 ). Stochastic models of intracellular transport . Reviews of Modern Physics 85 , 135 – 196 . doi: 10.1103/RevModPhys.85.135 . OpenUrl CrossRef 59. ↵ Eliaš , J. , and Clairambault , J. ( 2014 ). Reaction–diffusion systems for spatio-temporal intracellular protein networks: A beginner’s guide with two examples . Computational and Structural Biotechnology Journal 10 , 12 – 22 . doi: 10.1016/j.csbj.2014.05.007 . OpenUrl CrossRef 60. ↵ Yuan , J. , and Sundén , B. ( 2014 ). On mechanisms and models of multi-component gas diffusion in porous structures of fuel cell electrodes . International Journal of Heat and Mass Transfer 69 , 358 – 374 . doi: 10.1016/j.ijheatmasstransfer.2013.10.032 . OpenUrl CrossRef 61. Soh , S. , Byrska , M. , Kandere-Grzybowska , K. , and Grzybowski , B.A. ( 2010 ). Reactiondiffusion systems in intracellular molecular transport and control . Angew Chem Int Ed Engl 49 , 4170 – 4198 . doi: 10.1002/anie.200905513 . OpenUrl CrossRef PubMed Web of Science 62. ↵ Aly , S. , Kim , I. , and Sheen , D. ( 2011 ). Turing instability for a ratio-dependent predator– prey model with diffusion . Applied Mathematics and Computation 217 , 7265 – 7281 . doi: 10.1016/j.amc.2011.02.018 . OpenUrl CrossRef 63. ↵ Gotoh , T. , Kim , J.K. , Liu , J. , Vila-Caballer , M. , Stauffer , P.E. , Tyson , J.J. , and Finkielstein , C.V. ( 2016 ). Model-driven experimental approach reveals the complex regulatory distribution of p53 by the circadian factor Period 2 . Proceedings of the National Academy of Sciences 113 , 13516 – 13521 . doi: 10.1073/pnas.1607984113 . OpenUrl Abstract / FREE Full Text 64. Asgari-Targhi , A. , and Klerman , E.B. ( 2019 ). Mathematical modeling of circadian rhythms . WIREs Systems Biology and Medicine 11 , e1439 . doi: 10.1002/wsbm.1439 . OpenUrl CrossRef 65. ↵ Rojas , P. , Piro , O. , and Garcia , M. ( 2023 ). Biological rhythms generated by a single activator-repressor loop with heterogeneity and diffusion . 66. ↵ Purvis , J.E. , Karhohs , K.W. , Mock , C. , Batchelor , E. , Loewer , A. , and Lahav , G. ( 2012 ). p53 Dynamics Control Cell Fate . Science 336 , 1440 – 1444 . doi: 10.1126/science.1218351 . OpenUrl Abstract / FREE Full Text 67. Kim , J.K. , and Jackson , T.L. ( 2013 ). Mechanisms That Enhance Sustainability of p53 Pulses . PLOS ONE 8 , e65242 . doi: 10.1371/journal.pone.0065242 . OpenUrl CrossRef 68. ↵ Eliaš , J. , Dimitrio , L. , Clairambault , J. , and Natalini , R. ( 2014 ). The dynamics of p53 in single cells: physiologically based ODE and reaction–diffusion PDE models . Physical Biology 11 , 045001 . doi: 10.1088/1478-3975/11/4/045001 . OpenUrl CrossRef 69. ↵ Son , M. , Frank , T. , Holst-Hansen , T. , Wang , A.G. , Junkin , M. , Kashaf , S.S. , Trusina , A. , and Tay , S. ( 2022 ). Spatiotemporal NF-κB dynamics encodes the position, amplitude, and duration of local immune inputs . Science Advances 8 , eabn6240 . doi: 10.1126/sciadv.abn6240 . OpenUrl CrossRef 70. ↵ Tigges , M. , Marquez-Lago , T.T. , Stelling , J. , and Fussenegger , M. ( 2009 ). A tunable synthetic mammalian oscillator . Nature 457 , 309 – 312 . doi: 10.1038/nature07616 . OpenUrl CrossRef PubMed Web of Science 71. Purcell , O. , Savery , N.J. , Grierson , C.S. , and di Bernardo , M. ( 2010 ). A comparative analysis of synthetic genetic oscillators . Journal of The Royal Society Interface 7 , 1503 – 1524 . doi: 10.1098/rsif.2010.0183 . OpenUrl CrossRef PubMed 72. Kim , J.K. , Chen , Y. , Hirning , A.J. , Alnahhas , R.N. , Josić , K. , and Bennett , M.R. ( 2019 ). Long-range temporal coordination of gene expression in synthetic microbial consortia . Nature Chemical Biology 15 , 1102 – 1109 . doi: 10.1038/s41589-019-0372-9 . OpenUrl CrossRef 73. ↵ Macnamara , C.K. , Mitchell , E.I. , and Chaplain , M.A.J. ( 2019 ). Spatial-Stochastic modelling of synthetic gene regulatory networks . Journal of Theoretical Biology 468 , 27 – 44 . doi: 10.1016/j.jtbi.2019.02.003 . OpenUrl CrossRef 74. ↵ Cherstvy , A.G. , Chechkin , A.V. , and Metzler , R. ( 2013 ). Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes . New Journal of Physics 15 , 083039 . doi: 10.1088/1367-2630/15/8/083039 . OpenUrl CrossRef 75. Cherstvy , A.G. , and Metzler , R. ( 2014 ). Nonergodicity, fluctuations, and criticality in heterogeneous diffusion processes . Physical Review E 90 , 012134 . doi: 10.1103/PhysRevE.90.012134 . OpenUrl CrossRef 76. ↵ Ghosh , S.K. , Cherstvy , A.G. , Grebenkov , D.S. , and Metzler , R. ( 2016 ). Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments . New Journal of Physics 18 , 013027 . doi: 10.1088/1367-2630/18/1/013027 . OpenUrl CrossRef 77. ↵ Wachsmuth , M. , Waldeck , W. , and Langowski , J. ( 2000 ). Anomalous diffusion of fluorescent probes inside living cell nuclei investigated by spatially-resolved fluorescence correlation spectroscopy11Edited by W . Baumeister. Journal of Molecular Biology 298 , 677 – 689 . doi: 10.1006/jmbi.2000.3692 . OpenUrl CrossRef PubMed Web of Science 78. Seisenberger , G. , Ried , M.U. , Endreß , T. , Büning , H. , Hallek , M. , and Bräuchle , C. ( 2001 ). Real-Time Single-Molecule Imaging of the Infection Pathway of an Adeno-Associated Virus . Science 294 , 1929 – 1932 . doi: 10.1126/science.1064103 . OpenUrl Abstract / FREE Full Text 79. Banks , D.S. , and Fradin , C. ( 2005 ). Anomalous Diffusion of Proteins Due to Molecular Crowding . Biophysical Journal 89 , 2960 – 2971 . doi: 10.1529/biophysj.104.051078 . OpenUrl CrossRef PubMed Web of Science 80. ↵ Jeon , J.-H. , Tejedor , V. , Burov , S. , Barkai , E. , Selhuber-Unkel , C. , Berg-Sørensen , K. , Oddershede , L. , and Metzler , R. ( 2011 ). In Vivo Anomalous Diffusion and Weak Ergodicity Breaking of Lipid Granules . Physical Review Letters 106 , 048103 . doi: 10.1103/PhysRevLett.106.048103 . OpenUrl CrossRef PubMed 81. ↵ Lanoiselée , Y. , Moutal , N. , and Grebenkov , D.S. ( 2018 ). Diffusion-limited reactions in dynamic heterogeneous media . Nature Communications 9 , 4398 . doi: 10.1038/s41467-018-06610-6 . OpenUrl CrossRef 82. ↵ Salman , H. , Abu-Arish , A. , Oliel , S. , Loyter , A. , Klafter , J. , Granek , R. , and Elbaum , M. ( 2005 ). Nuclear localization signal peptides induce molecular delivery along microtubules . Biophys J 89 , 2134 – 2145 . doi: 10.1529/biophysj.105.060160 . OpenUrl CrossRef PubMed Web of Science 83. ↵ Halbi , G. , Fayer , I. , Aranovich , D. , Gat , S. , Bar , S. , Erukhimovitch , V. , Granek , R. , and Bernheim-Groswasser , A. ( 2021 ). Nano-Particles Carried by Multiple Dynein Motors Self-Regulate Their Number of Actively Participating Motors . International Journal of Molecular Sciences 22 , 8893 . OpenUrl 84. ↵ Snapp , E. ( 2005 ). Design and use of fluorescent fusion proteins in cell biology . Curr Protoc Cell Biol Chapter 21 , 21 24 21 – 21 24 13 . doi: 10.1002/0471143030.cb2104s27 . OpenUrl CrossRef 85. ↵ Milo , R. , and Phillips , R. ( 2015 ). Cell Biology by the Numbers ( CRC Press ). 86. ↵ Aryal , R.P. , Kwak , P.B. , Tamayo , A.G. , Gebert , M. , Chiu , P.L. , Walz , T. , and Weitz , C.J. ( 2017 ). Macromolecular Assemblies of the Mammalian Circadian Clock . Mol Cell 67 , 770 – 782 e776. doi: 10.1016/j.molcel.2017.07.017 . OpenUrl CrossRef PubMed 87. ↵ Koch , A.A. , Bagnall , J.S. , Smyllie , N.J. , Begley , N. , Adamson , A.D. , Fribourgh , J.L. , Spiller , D.G. , Meng , Q.J. , Partch , C.L. , Strimmer , K. , et al. ( 2022 ). Quantification of protein abundance and interaction defines a mechanism for operation of the circadian clock . Elife 11 . doi: 10.7554/eLife.73976 . OpenUrl CrossRef 88. ↵ Tzafriri , A.R. ( 2003 ). Michaelis–Menten kinetics at high enzyme concentrations . Bulletin of Mathematical Biology 65 , 1111 – 1129 . doi: 10.1016/S0092-8240(03)00059-4 . OpenUrl CrossRef PubMed 89. Kim , J.K. , and Tyson , J.J. ( 2020 ). Misuse of the Michaelis-Menten rate law for protein interaction networks and its remedy . PLoS Comput Biol 16 , e1008258 . doi: 10.1371/journal.pcbi.1008258 . OpenUrl CrossRef 90. Song , Y.M. , Hong , H. , and Kim , J.K. ( 2021 ). Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities . PLOS Computational Biology 17 , e1008952 . doi: 10.1371/journal.pcbi.1008952 . OpenUrl CrossRef 91. Eilertsen , J. , Schnell , S. , and Walcher , S. ( 2022 ). On the anti-quasi-steady-state conditions of enzyme kinetics . Mathematical Biosciences 350 , 108870 . doi: 10.1016/j.mbs.2022.108870 . OpenUrl CrossRef 92. ↵ Jeong , E.M. , Song , Y.M. , and Kim , J.K. ( 2022 ). Combined multiple transcriptional repression mechanisms generate ultrasensitivity and oscillations . Interface Focus 12 , 20210084 . doi: 10.1098/rsfs.2021.0084 . OpenUrl CrossRef PubMed 93. ↵ Shin , J. , Choi , Y. , and Kim , J. ( 2019 ). The Cahn–Hilliard Equation with Generalized Mobilities in Complex Geometries . Mathematical Problems in Engineering 2019 , 1710270 . doi: 10.1155/2019/1710270 . OpenUrl CrossRef 94. ↵ Shin , J. , Jeong , D. , and Kim , J. ( 2011 ). A conservative numerical method for the Cahn– Hilliard equation in complex domains . Journal of Computational Physics 230 , 7441 – 7455 . doi: 10.1016/j.jcp.2011.06.009 . OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted January 24, 2024. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport Seok Joo Chae , Dae Wook Kim , Oleg A. Igoshin , Seunggyu Lee , Jae Kyoung Kim bioRxiv 2024.01.12.575351; doi: https://doi.org/10.1101/2024.01.12.575351 Share This Article: Copy Citation Tools Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport Seok Joo Chae , Dae Wook Kim , Oleg A. Igoshin , Seunggyu Lee , Jae Kyoung Kim bioRxiv 2024.01.12.575351; doi: https://doi.org/10.1101/2024.01.12.575351 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Systems Biology Subject Areas All Articles Animal Behavior and Cognition (7651) Biochemistry (17746) Bioengineering (13928) Bioinformatics (42066) Biophysics (21499) Cancer Biology (18650) Cell Biology (25579) Clinical Trials (138) Developmental Biology (13409) Ecology (19947) Epidemiology (2067) Evolutionary Biology (24374) Genetics (15633) Genomics (22557) Immunology (17775) Microbiology (40505) Molecular Biology (17217) Neuroscience (88796) Paleontology (667) Pathology (2845) Pharmacology and Toxicology (4836) Physiology (7664) Plant Biology (15179) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9839) Zoology (2272)