Luminal transport through intact endoplasmic reticulum limits the magnitude of localized Ca2+ signals

The endoplasmic reticulum (ER) forms an interconnected network of tubules stretching throughout the cell. Understanding how ER functionality relies on its structural organization is crucial for elucidating cellular vulnerability to ER perturbations, which have been implicated in several neuronal pathologies. One of the key functions of the ER is enabling Ca[Formula: see text] signaling by storing large quantities of this ion and releasing it into the cytoplasm in a spatiotemporally controlled manner. Through a combination of physical modeling and live-cell imaging, we demonstrate that alterations in ER shape significantly impact its ability to support efficient local Ca[Formula: see text] releases, due to hindered transport of luminal content within the ER. Our model reveals that rapid Ca[Formula: see text] release necessitates mobile luminal buffer proteins with moderate binding strength, moving through a well-connected network of ER tubules. These findings provide insight into the functional advantages of normal ER architecture, emphasizing its importance as a kinetically efficient intracellular Ca[Formula: see text] delivery system.

Endoplasmic reticulum network heterogeneity guides diffusive transport and kinetics

The endoplasmic reticulum (ER) is a dynamic network of interconnected sheets and tubules that orchestrates the distribution of lipids, ions, and proteins throughout the cell. The impact of its complex, dynamic morphology on its function as an intracellular transport hub remains poorly understood. To elucidate the functional consequences of ER network structure and dynamics, we quantify how the heterogeneity of the peripheral ER in COS7 cells affects diffusive protein transport. In vivo imaging of photoactivated ER membrane proteins demonstrates their nonuniform spreading to adjacent regions, in a manner consistent with simulations of diffusing particles on extracted network structures. Using a minimal network model to represent tubule rearrangements, we demonstrate that ER network dynamics are sufficiently slow to have little effect on diffusive protein transport. Furthermore, stochastic simulations reveal a novel consequence of ER network heterogeneity: the existence of "hot spots" where sparse diffusive reactants are more likely to find one another. ER exit sites, specialized domains regulating cargo export from the ER, are shown to be disproportionately located in highly accessible regions, further from the outer boundary of the cell. Combining in vivo experiments with analytic calculations, quantitative image analysis, and computational modeling, we demonstrate how structure guides diffusive protein transport and reactions in the ER.

Diffusive exit rates through pores in membrane-enclosed structures

The function of many membrane-enclosed intracellular structures relies on release of diffusing particles that exit through narrow pores or channels in the membrane. The rate of release varies with pore size, density, and length of the channel. We propose a simple approximate model, validated with stochastic simulations, for estimating the effective release rate from cylinders, and other simple-shaped domains, as a function of channel parameters. The results demonstrate that, for very small pores, a low density of channels scattered over the boundary is sufficient to achieve substantial rates of particle release. Furthermore, we show that increasing the length of passive channels will both reduce release rates and lead to a less steep dependence on channel density. Our results are compared to previously-measured local calcium release rates from tubules of the endoplasmic reticulum, providing an estimate of the relevant channel density responsible for the observed calcium efflux.

Spatiotemporal analysis of axonal autophagosome-lysosome dynamics reveals limited fusion events and slow maturation

Macroautophagy is a homeostatic process required to clear cellular waste. Neuronal autophagosomes form constitutively in the distal tip of the axon and are actively transported toward the soma, with cargo degradation initiated en route. Cargo turnover requires autophagosomes to fuse with lysosomes to acquire degradative enzymes; however, directly imaging these fusion events in the axon is impractical. Here we use a quantitative model, parameterized and validated using data from primary hippocampal neurons, to explore the autophagosome maturation process. We demonstrate that retrograde autophagosome motility is independent of fusion and that most autophagosomes fuse with only a few lysosomes during axonal transport. Our results indicate that breakdown of the inner autophagosomal membrane is much slower in neurons than in nonneuronal cell types, highlighting the importance of this late maturation step. Together, rigorous quantitative measurements and mathematical modeling elucidate the dynamics of autophagosome-lysosome interaction and autophagosomal maturation in the axon.

Mitochondrial mRNA localization is governed by translation kinetics and spatial transport

For many nuclear-encoded mitochondrial genes, mRNA localizes to the mitochondrial surface co-translationally, aided by the association of a mitochondrial targeting sequence (MTS) on the nascent peptide with the mitochondrial import complex. For a subset of these co-translationally localized mRNAs, their localization is dependent on the metabolic state of the cell, while others are constitutively localized. To explore the differences between these two mRNA types we developed a stochastic, quantitative model for MTS-mediated mRNA localization to mitochondria in yeast cells. This model includes translation, applying gene-specific kinetics derived from experimental data; and diffusion in the cytosol. Even though both mRNA types are co-translationally localized we found that the steady state number, or density, of ribosomes along an mRNA was insufficient to differentiate the two mRNA types. Instead, conditionally-localized mRNAs have faster translation kinetics which modulate localization in combination with changes to diffusive search kinetics across metabolic states. Our model also suggests that the MTS requires a maturation time to become competent to bind mitochondria. Our work indicates that yeast cells can regulate mRNA localization to mitochondria by controlling mitochondrial volume fraction (influencing diffusive search times) and gene translation kinetics (adjusting mRNA binding competence) without the need for mRNA-specific binding proteins. These results shed light on both global and gene-specific mechanisms that enable cells to alter mRNA localization in response to changing metabolic conditions.

Morphology and Transport in Eukaryotic Cells

Transport of intracellular components relies on a variety of active and passive mechanisms, ranging from the diffusive spreading of small molecules over short distances to motor-driven motion across long distances. The cell-scale behavior of these mechanisms is fundamentally dependent on the morphology of the underlying cellular structures. Diffusion-limited reaction times can be qualitatively altered by the presence of occluding barriers or by confinement in complex architectures, such as those of reticulated organelles. Motor-driven transport is modulated by the architecture of cytoskeletal filaments that serve as transport highways. In this review, we discuss the impact of geometry on intracellular transport processes that fulfill a broad range of functional objectives, including delivery, distribution, and sorting of cellular components. By unraveling the interplay between morphology and transport efficiency, we aim to elucidate key structure-function relationships that govern the architecture of transport systems at the cellular scale. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Optimizing microtubule arrangements for rapid cargo capture

Cellular functions such as autophagy, cell signaling, and vesicular trafficking involve the retrograde transport of motor-driven cargo along microtubules. Typically, newly formed cargo engages in slow undirected movement from its point of origin before attaching to a microtubule. In some cell types, cargo destined for delivery to the perinuclear region relies on capture at dynein-enriched loading zones located near microtubule plus ends. Such systems include extended cell regions of neurites and fungal hyphae, where the efficiency of the initial diffusive loading process depends on the axial distribution of microtubule plus ends relative to the initial cargo position. We use analytic mean first-passage time calculations and numerical simulations to model diffusive capture processes in tubular cells, exploring how the spatial arrangement of microtubule plus ends affects the efficiency of retrograde cargo transport. Our model delineates the key features of optimal microtubule arrangements that minimize mean cargo capture times. Namely, we show that configurations with a single microtubule plus end abutting the distal tip and broadly distributed other plus ends allow for efficient capture in a variety of different scenarios for retrograde transport. Live-cell imaging of microtubule plus ends in Aspergillus nidulans hyphae indicates that their distributions exhibit these optimal qualitative features. Our results highlight important coupling effects between the distribution of microtubule tips and retrograde cargo transport, providing guiding principles for the spatial arrangement of microtubules within tubular cell regions.

Diffusive search and trajectories on tubular networks: a propagator approach

Several organelles in eukaryotic cells, including mitochondria and the endoplasmic reticulum, form interconnected tubule networks extending throughout the cell. These tubular networks host many biochemical pathways that rely on proteins diffusively searching through the network to encounter binding partners or localized target regions. Predicting the behavior of such pathways requires a quantitative understanding of how confinement to a reticulated structure modulates reaction kinetics. In this work, we develop both exact analytical methods to compute mean first passage times and efficient kinetic Monte Carlo algorithms to simulate trajectories of particles diffusing in a tubular network. Our approach leverages exact propagator functions for the distribution of transition times between network nodes and allows large simulation time steps determined by the network structure. The methodology is applied to both synthetic planar networks and organelle network structures, demonstrating key general features such as the heterogeneity of search times in different network regions and the functional advantage of broadly distributing target sites throughout the network. The proposed algorithms pave the way for future exploration of the interrelationship between tubular network structure and biomolecular reaction kinetics.

Optimizing mitochondrial maintenance in extended neuronal projections

Neurons rely on localized mitochondria to fulfill spatially heterogeneous metabolic demands. Mitochondrial aging occurs on timescales shorter than the neuronal lifespan, necessitating transport of fresh material from the soma. Maintaining an optimal distribution of healthy mitochondria requires an interplay between a stationary pool localized to sites of high metabolic demand and a motile pool capable of delivering new material. Interchange between these pools can occur via transient fusion / fission events or by halting and restarting entire mitochondria. Our quantitative model of neuronal mitostasis identifies key parameters that govern steady-state mitochondrial health at discrete locations. Very infrequent exchange between stationary and motile pools optimizes this system. Exchange via transient fusion allows for robust maintenance, which can be further improved by selective recycling through mitophagy. These results provide a framework for quantifying how perturbations in organelle transport and interactions affect mitochondrial homeostasis in neurons, a key aspect underlying many neurodegenerative disorders.

Elastic wrinkling of keratocyte lamellipodia driven by myosin-induced contractile stress

During actin-based cell migration, the actin cytoskeleton in the lamellipodium both generates and responds to force, which has functional consequences for the ability of the cell to extend protrusions. However, the material properties of the lamellipodial actin network and its response to stress on the timescale of motility are incompletely understood. Here, we describe a dynamic wrinkling phenotype in the lamellipodium of fish keratocytes, in which the actin sheet buckles upward away from the ventral membrane of the cell, forming a periodic pattern of wrinkles perpendicular to the cell's leading edge. Cells maintain an approximately constant wrinkle wavelength over time despite new wrinkle formation and the lateral movement of wrinkles in the cell frame of reference, suggesting that cells have a preferred or characteristic wrinkle wavelength. Generation of wrinkles is dependent upon myosin contractility, and their wavelength scales directly with the density of the actin network and inversely with cell adhesion. These results are consistent with a simple physical model for wrinkling in an elastic sheet under compression and suggest that the lamellipodial cytoskeleton behaves as an elastic material on the timescale of cell migration despite rapid actin turnover.

Design principles for the glycoprotein quality control pathway

Newly-translated glycoproteins in the endoplasmic reticulum (ER) often undergo cycles of chaperone binding and release in order to assist in folding. Quality control is required to distinguish between proteins that have completed native folding, those that have yet to fold, and those that have misfolded. Using quantitative modeling, we explore how the design of the quality-control pathway modulates its efficiency. Our results show that an energy-consuming cyclic quality-control process, similar to the observed physiological system, outperforms alternative designs. The kinetic parameters that optimize the performance of this system drastically change with protein production levels, while remaining relatively insensitive to the protein folding rate. Adjusting only the degradation rate, while fixing other parameters, allows the pathway to adapt across a range of protein production levels, aligning with in vivo measurements that implicate the release of degradation-associated enzymes as a rapid-response system for perturbations in protein homeostasis. The quantitative models developed here elucidate design principles for effective glycoprotein quality control in the ER, improving our mechanistic understanding of a system crucial to maintaining cellular health.

We explore the architecture and limitations of the quality-control pathway responsible for efficient folding of secretory proteins. Newly-synthesized proteins are tagged by the attachment of a ‘glycan’ sugar chain which facilitates their binding to a chaperone that assists protein folding. Removal of a specific sugar group on the glycan ends the interaction with the chaperone, and not-yet-folded proteins can be re-tagged for another round of chaperone binding. A degradation pathway acts in parallel with the folding cycle, to remove those proteins that have remained unfolded for a sufficiently long time. We develop and solve a mathematical model of this quality-control system, showing that the cyclical design found in living cells is uniquely able to maximize folded protein throughput while avoiding accumulation of unfolded proteins. Although this physiological model provides the best performance, its parameters must be adjusted to perform optimally under different protein production loads, and any single fixed set of parameters leads to poor performance when production rate is altered. We find that a single adjustable parameter, the protein degradation rate, is sufficient to allow optimal performance across a range of conditions. Interestingly, observations of living cells suggest that the degradation speed is indeed rapidly adjusted.

Getting around the cell: physical transport in the intracellular world

Eukaryotic cells face the challenging task of transporting a variety of particles through the complex intracellular milieu in order to deliver, distribute, and mix the many components that support cell function. In this review, we explore the biological objectives and physical mechanisms of intracellular transport. Our focus is on cytoplasmic and intra-organelle transport at the whole-cell scale. We outline several key biological functions that depend on physically transporting components across the cell, including the delivery of secreted proteins, support of cell growth and repair, propagation of intracellular signals, establishment of organelle contacts, and spatial organization of metabolic gradients. We then review the three primary physical modes of transport in eukaryotic cells: diffusive motion, motor-driven transport, and advection by cytoplasmic flow. For each mechanism, we identify the main factors that determine speed and directionality. We also highlight the efficiency of each transport mode in fulfilling various key objectives of transport, such as particle mixing, directed delivery, and rapid target search. Taken together, the interplay of diffusion, molecular motors, and flows supports the intracellular transport needs that underlie a broad variety of biological phenomena.

Thermal fracture kinetics of heterogeneous semiflexible polymers

The fracture and severing of polymer chains plays a critical role in the failure of fibrous materials and the regulated turnover of intracellular filaments. Using continuum wormlike chain models, we investigate the fracture of semiflexible polymers via thermal bending fluctuations, focusing on the role of filament flexibility and dynamics. Our results highlight a previously unappreciated consequence of mechanical heterogeneity in the filament, which enhances the rate of thermal fragmentation particularly in cases where constraints hinder the movement of the chain ends. Although generally applicable to semiflexible chains with regions of different bending stiffness, the model is motivated by a specific biophysical system: the enhanced severing of actin filaments at the boundary between stiff bare regions and mechanically softened regions that are coated with cofilin regulatory proteins. The results presented here point to a potential mechanism for disassembly of polymeric materials in general and cytoskeletal actin networks in particular by the introduction of locally softened chain regions, as occurs with cofilin binding.

Drive, filter, and stick: A protein sorting conspiracy in photoreceptors

The sorting of proteins into different functional compartments is a fundamental cellular task. In this issue, Maza et al. (2019. J. Cell Biol https://doi.org/10.1083/jcb.201906024) demonstrate that distinct protein populations are dynamically generated in specialized regions of photoreceptors via an interplay of protein-membrane affinity, impeded diffusion, and driven transport.

High-resolution and high-accuracy topographic and transcriptional maps of the nucleosome barrier

Nucleosomes represent mechanical and energetic barriers that RNA Polymerase II (Pol II) must overcome during transcription. A high-resolution description of the barrier topography, its modulation by epigenetic modifications, and their effects on Pol II nucleosome crossing dynamics, is still missing. Here, we obtain topographic and transcriptional (Pol II residence time) maps of canonical, H2A.Z, and monoubiquitinated H2B (uH2B) nucleosomes at near base-pair resolution and accuracy. Pol II crossing dynamics are complex, displaying pauses at specific loci, backtracking, and nucleosome hopping between wrapped states. While H2A.Z widens the barrier, uH2B heightens it, and both modifications greatly lengthen Pol II crossing time. Using the dwell times of Pol II at each nucleosomal position we extract the energetics of the barrier. The orthogonal barrier modifications of H2A.Z and uH2B, and their effects on Pol II dynamics rationalize their observed enrichment in +1 nucleosomes and suggest a mechanism for selective control of gene expression.

Listeria monocytogenes cell-to-cell spread in epithelia is heterogeneous and dominated by rare pioneer bacteria

Listeria monocytogenes hijacks host actin to promote its intracellular motility and intercellular spread. While L. monocytogenes virulence hinges on cell-to-cell spread, little is known about the dynamics of bacterial spread in epithelia at a population level. Here, we use live microscopy and statistical modeling to demonstrate that L. monocytogenes cell-to-cell spread proceeds anisotropically in an epithelial monolayer in culture. We show that boundaries of infection foci are irregular and dominated by rare pioneer bacteria that spread farther than the rest. We extend our quantitative model for bacterial spread to show that heterogeneous spreading behavior can improve the chances of creating a persistent L. monocytogenes infection in an actively extruding epithelium. Thus, our results indicate that L. monocytogenes cell-to-cell spread is heterogeneous, and that rare pioneer bacteria determine the frontier of infection foci and may promote bacterial infection persistence in dynamic epithelia.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Eating food that has been contaminated with bacteria called Listeria monocytogenes can result in life-threatening infections. The bacteria first invade the epithelial cells that line the small intestine. After this, L. monocytogenes can move from one host cell to another, which allows the infection to reach other organs.

Most studies into how L. monocytogenes infections spread have focused either on how single bacterial cells move from one host cell to the next, or on how millions of bacteria damage host tissues. Little was known about the intermediate steps of an infection, where the bacteria start to colonize the small intestine.

To investigate, Ortega et al. recorded videos of L. monocytogenes spreading between epithelial cells grown on a glass coverslip, and developed computer simulations to try to reproduce how the bacteria spread. This revealed that the bacteria do not all move in the same way. Instead, less than 1% of the bacteria move in ‘steps’ that are up to 10 times longer than those taken by the others. Ortega et al. named these bacteria ‘pioneers’.

Ortega et al. propose that the pioneers form long protrusions that allow them to spread directly from an infected cell to a non-neighboring cell. By taking these large steps, the pioneers may increase the chances that the bacteria will cause a long-lasting infection.

Future research will be needed to answer further questions about the pioneers. For example, how do the pioneer bacteria differ from the majority of bacterial cells? Would targeting anti-bacterial treatments at pioneers make it easier to treat infections? It also remains to be seen if other types of bacteria also show this pioneer behavior.

Spatial control of neuronal metabolism through glucose-mediated mitochondrial transport regulation

Eukaryotic cells modulate their metabolism by organizing metabolic components in response to varying nutrient availability and energy demands. In rat axons, mitochondria respond to glucose levels by halting active transport in high glucose regions. We employ quantitative modeling to explore physical limits on spatial organization of mitochondria and localized metabolic enhancement through regulated stopping of processive motion. We delineate the role of key parameters, including cellular glucose uptake and consumption rates, that are expected to modulate mitochondrial distribution and metabolic response in spatially varying glucose conditions. Our estimates indicate that physiological brain glucose levels fall within the limited range necessary for metabolic enhancement. Hence mitochondrial localization is shown to be a plausible regulatory mechanism for neuronal metabolic flexibility in the presence of spatially heterogeneous glucose, as may occur in long processes of projection neurons. These findings provide a framework for the control of cellular bioenergetics through organelle trafficking.

Cells are equipped with power factories called mitochondria that turn nutrients into chemical energy to fuel processes in the cell. Hundreds of mitochondria move throughout the cell, shifting their positions in response to energy demands. This happens via molecular motors that pick the mitochondria up and carry them to new locations. Such movements enable the mitochondria to accumulate in parts of the cell with the greatest energy needs.

Mitochondria of nerve cells or neurons have a particular challenging job, as neurons can be very long and different parts within the cells can have different energy needs. It has been shown that mitochondria stop in regions where nutrients such as sugar are most concentrated. So far, it has been unclear whether this regulated stopping helps control energy balance in neurons.

Here, Agrawal et al. used a computational model of rat neurons to find out whether sugar levels are sufficient in guiding mitochondria. The results showed that the mitochondria only accumulated in high-nutrient regions when the sugar concentrations were moderate – not too low and not too high. A specific range of sugar levels was necessary to make this mechanism useful for increasing the efficiency of energy production. Such concentrations match the ones observed in healthy rat brains.

When neurons are unable to meet their energy demands, they stop working and sometimes even die. This is the case in many diseases, including diabetes, dementia, and Alzheimer’s disease. Computer models allow us to explore the complex energy regulation in detail. A better understanding of how neurons regulate their energy production and demand may help us discover how they become faulty in these diseases.

Multimodal transport and dispersion of organelles in narrow tubular cells

Intracellular components explore the cytoplasm via active motor-driven transport in conjunction with passive diffusion. We model the motion of organelles in narrow tubular cells using analytical techniques and numerical simulations to study the efficiency of different transport modes in achieving various cellular objectives. Our model describes length and time scales over which each transport mode dominates organelle motion, along with various metrics to quantify exploration of intracellular space. For organelles that search for a specific target, we obtain the average capture time for given transport parameters and show that diffusion and active motion contribute to target capture in the biologically relevant regime. Because many organelles have been found to tether to microtubules when not engaged in active motion, we study the interplay between immobilization due to tethering and increased probability of active transport. We derive parameter-dependent conditions under which tethering enhances long-range transport and improves the target capture time. These results shed light on the optimization of intracellular transport machinery and provide experimentally testable predictions for the effects of transport regulation mechanisms such as tethering.

Cytoplasmic Flow and Mixing Due to Deformation of Motile Cells

The cytoplasm of a living cell is a dynamic environment through which intracellular components must move and mix. In motile, rapidly deforming cells such as human neutrophils, bulk cytoplasmic flow couples cell deformation to the transport and dispersion of cytoplasmic particles. Using particle-tracking measurements in live neutrophil-like cells, we demonstrate that fluid flow associated with the cell deformation contributes to the motion of small acidic organelles, dominating over diffusion on timescales above a few seconds. We then use a general physical model of particle dispersion in a deforming fluid domain to show that transport of organelle-sized particles between the cell periphery and the bulk can be enhanced by dynamic deformation comparable to that observed in neutrophils. Our results implicate an important mechanism contributing to organelle transport in these motile cells: cytoplasmic flow driven by cell shape deformation.

Disentangling Random Motion and Flow in a Complex Medium

We describe a technique for deconvolving the stochastic motion of particles from large-scale fluid flow in a dynamic environment such as that found in living cells. The method leverages the separation of timescales to subtract out the persistent component of motion from single-particle trajectories. The mean-squared displacement of the resulting trajectories is rescaled so as to enable robust extraction of the diffusion coefficient and subdiffusive scaling exponent of the stochastic motion. We demonstrate the applicability of the method for characterizing both diffusive and fractional Brownian motion overlaid by flow and analytically calculate the accuracy of the method in different parameter regimes. This technique is employed to analyze the motion of lysosomes in motile neutrophil-like cells, showing that the cytoplasm of these cells behaves as a viscous fluid at the timescales examined.

Bacterial division. Mechanical crack propagation drives millisecond daughter cell separation in Staphylococcus aureus

When Staphylococcus aureus undergoes cytokinesis, it builds a septum, generating two hemispherical daughters whose cell walls are only connected via a narrow peripheral ring. We found that resolution of this ring occurred within milliseconds ("popping"), without detectable changes in cell volume. The likelihood of popping depended on cell-wall stress, and the separating cells split open asymmetrically, leaving the daughters connected by a hinge. An elastostatic model of the wall indicated high circumferential stress in the peripheral ring before popping. Last, we observed small perforations in the peripheral ring that are likely initial points of mechanical failure. Thus, the ultrafast daughter cell separation in S. aureus appears to be driven by accumulation of stress in the peripheral ring and exhibits hallmarks of mechanical crack propagation.

Thermodynamic model of heterochromatin formation through epigenetic regulation

Gene regulation in eukaryotes requires the segregation of silenced genomic regions into densely packed heterochromatin, leaving the active genes in euchromatin regions more accessible. We introduce a model that connects the presence of epigenetically inherited histone marks, methylation at histone 3 lysine-9, to the physical compaction of chromatin fibers via the binding of heterochromatin protein 1 (HP1). Our model demonstrates some of the key physical features that are necessary to explain experimental observations. In particular, we demonstrate that strong cooperative interactions among the HP1 proteins are necessary to see the phase segregation of heterochromatin and euchromatin regions. We also explore how the cell can use the concentration of HP1 to control condensation and under what circumstances there is a threshold of methylation over which the fibers will compact. Finally, we consider how different potential in vivo fiber structures as well as the flexibility of the histone 3 tail can affect the bridging of HP1. Many of the observations that we make about the HP1 system are guided by general thermodynamics principles and thus could play a role in other DNA organizational processes such as the binding of linker histones.

Multiscale dynamics of semiflexible polymers from a universal coarse-graining procedure

Simulating the dynamics of a semiflexible polymer across time and length scales that bridge the rigid and flexible regimes requires a physically sound method for generating coarse-grained representations of the polymer. Here, we study the dynamic behavior of the discrete stretchable, shearable wormlike chain model, which can be used to coarse-grain a continuous semi-elastic chain at an arbitrary discretization. We show that the dynamics of this universal model match those of the wormlike chain at length scales above the discretization length. The evolution of the stress correlation, as probed through Brownian dynamics simulations, is found to reproduce the predicted behavior in both the rigid and flexible regimes, spanning over six orders of magnitude in time scales. The coarse-graining approach employed here thus enables dynamic simulation of semiflexible polymers at lengths and times that were previously inaccessible with conventional methods.

An in vitro assay for entry into cilia reveals unique properties of the soluble diffusion barrier

The ciliary permeability barrier is mechanistically distinct from other cellular diffusion barriers and allows soluble proteins under ∼100 kD in size to enter cilia in the absence of active transport.

Specific proteins are concentrated within primary cilia, whereas others remain excluded. To understand the mechanistic basis of entry into cilia, we developed an in vitro assay using cells in which the plasma membrane was permeabilized, but the ciliary membrane was left intact. Using a diffusion-to-capture system and quantitative analysis, we find that proteins >9 nm in diameter (∼100 kD) are restricted from entering cilia, and we confirm these findings in vivo. Interference with the nuclear pore complex (NPC) or the actin cytoskeleton in permeabilized cells demonstrated that the ciliary diffusion barrier is mechanistically distinct from those of the NPC or the axon initial segment. Moreover, applying a mass transport model to this system revealed diffusion coefficients for soluble and membrane proteins within cilia that are compatible with rapid exploration of the ciliary space in the absence of active transport. Our results indicate that large proteins require active transport for entry into cilia but not necessarily for movement inside cilia.

Single molecule imaging reveals a major role for diffusion in the exploration of ciliary space by signaling receptors

The dynamic organization of signaling cascades inside primary cilia is key to signal propagation. Yet little is known about the dynamics of ciliary membrane proteins besides a possible role for motor-driven Intraflagellar Transport (IFT). To characterize these dynamics, we imaged single molecules of Somatostatin Receptor 3 (SSTR3, a GPCR) and Smoothened (Smo, a Hedgehog signal transducer) in the ciliary membrane. While IFT trains moved processively from one end of the cilium to the other, single SSTR3 and Smo underwent mostly diffusive behavior interspersed with short periods of directional movements. Statistical subtraction of instant velocities revealed that SSTR3 and Smo spent less than a third of their time undergoing active transport. Finally, SSTR3 and IFT movements could be uncoupled by perturbing either membrane protein diffusion or active transport. Thus ciliary membrane proteins move predominantly by diffusion, and attachment to IFT trains is transient and stochastic rather than processive or spatially determined.

DOI:http://dx.doi.org/10.7554/eLife.00654.001

Primary cilia are tiny protrusions from the cell surface, which have a central role in processing sensory stimuli, such as light or odorants. Cilia are also involved in mediating the response to developmental signaling molecules, including Sonic Hedgehog, and may help to convert mechanical signals into electrical or chemical ones. Primary cilia are made up of an axoneme—a core structure that consists of microtubules extending along the length of the cilium—ensheathed by a membrane that contains a number of receptor proteins.

These receptor proteins travel up and down the cilium, and it is generally assumed that an active process known as intraflagellar transport is responsible for their movement. This process is mediated by motor proteins called kinesins and dyneins, which carry cargo proteins along axonemal microtubules. However, it has been difficult to study the transport of individual receptor proteins directly because they are uniformly distributed over the membranes of the cilia.

Now, Ye et al. have shown that intraflagellar transport is not the most important mode of transport for membrane proteins within primary cilia. By labelling individual receptors with a fluorescent dye and then filming their movements under a microscope, Ye et al. found that the receptors generally did not show the directed, linear motion that would be expected from intraflagellar transport. Instead, much of their movement occurred through passive diffusion, with occasional short bursts of directed motion.

To investigate how rapidly receptor molecules could move through the cilium in this way, Ye et al. used a technique called fluorescence recovery after photobleaching (FRAP). This involves using light to bleach the fluorescent dye attached to receptor molecules in part of the cilium, and then measuring how long it takes for the fluorescence to return as a result of other labelled molecules moving into the bleached area: the shorter this time, the faster the movement of the molecules. It took less than a minute for fluorescence to be restored within a primary cilium, indicating that passive diffusion with occasional active transport can move proteins rapidly through the structure.

By using drugs to inhibit intraflagellar transport, Ye et al. confirmed that the majority of membrane protein transport within primary cilia occurs via diffusion. Further studies are now required to determine whether this is also the case for other molecules that travel along cilia, and whether intraflagellar transport may have a more important role in the assembly of these structures.

DOI:http://dx.doi.org/10.7554/eLife.00654.002

Force fluctuations impact kinetics of biomolecular systems

A wide array of biological processes occur at rates that vary significantly with force. Instantaneous molecular forces fluctuate due to thermal noise and active processes, leading to concomitant fluctuations in biomolecular rate constants. We demonstrate that such fluctuations have a dramatic effect on the transition kinetics of force-dependent processes. As an illustrative, biologically relevant example, we model the pausing of eukaryotic RNA polymerase as it transcribes nucleosomal DNA. Incorporating force fluctuations in the model yields qualitatively different predictions for the pausing time scales when compared to behavior under the average force alone. We use our model to illustrate the broad range of behaviors that can arise in biomolecular processes that are susceptible to force fluctuations. The fluctuation time scale, which varies significantly for in vivo biomolecular processes, yields very different results for overall rates and dramatically alters the force regime of relevance to the transition. Our results emphasize the importance of transient high-force behavior for determining kinetics in the fluctuating environment of a living cell.

Local geometry and elasticity in compact chromatin structure

The hierarchical packaging of DNA into chromatin within a eukaryotic nucleus plays a pivotal role in both the accessibility of genomic information and the dynamics of replication. Our work addresses the role of nanoscale physical and geometric properties in determining the structure of chromatin at the mesoscale level. We study the packaging of DNA in chromatin fibers by optimization of regular helical morphologies, considering the elasticity of the linker DNA as well as steric packing of the nucleosomes and linkers. Our model predicts a broad range of preferred helix structures for a fixed linker length of DNA; changing the linker length alters the predicted ensemble. Specifically, we find that the twist registry of the nucleosomes, as set by the internucleosome repeat length, determines the preferred angle between the nucleosomes and the fiber axis. For moderate to long linker lengths, we find a number of energetically comparable configurations with different nucleosome-nucleosome interaction patterns, indicating a potential role for kinetic trapping in chromatin fiber formation. Our results highlight the key role played by DNA elasticity and local geometry in regulating the hierarchical packaging of the genome.

Dynamic strategies for target-site search by DNA-binding proteins

Gene regulatory proteins find their target sites on DNA remarkably quickly; the experimental binding rate for lac repressor is orders-of-magnitude higher than predicted by free diffusion alone. It has been proposed that nonspecific binding aids the search by allowing proteins to slide and hop along DNA. We develop a reaction-diffusion theory of protein translocation that accounts for transport both on and off the strand and incorporates the physical conformation of DNA. For linear DNA modeled as a wormlike chain, the distribution of hops available to a protein exhibits long, power-law tails that make the long-time displacement along the strand superdiffusive. Our analysis predicts effective superdiffusion coefficients for given nonspecific binding and unbinding rate parameters. Translocation rate exhibits a maximum at intermediate values of the binding rate constant, while search efficiency is optimized at larger binding rate constant values. Thus, our theory predicts a region of values of the nonspecific binding and unbinding rate parameters that balance the protein translocation rate and the efficiency of the search. Published data for several proteins falls within this predicted region of parameter values.

Twist- and tension-mediated elastic coupling between DNA-binding proteins

We study the effective interaction between DNA-binding proteins that arises from elastic stresses in the DNA when tension is applied. Using the wormlike chain model, we calculate the free energy cost of introducing multiple nearby bends in the DNA. We find that the bend deformation energy promotes aggregation to straighten the linker DNA, while twist resistance of the linker leads to damped oscillations in the coupling free energy between two proteins. We calculate the mean first encounter time for proteins sliding along DNA, indicating, in some cases, an optimal applied tension for protein assembly. Our results highlight the need to consider DNA twist even when no torsion is applied and the DNA ends are free to rotate. The variable-range oscillatory coupling between DNA-binding proteins may provide a versatile mechanism for tension-mediated gene regulation.

End-to-end distribution for a wormlike chain in arbitrary dimensions

We construct an efficient methodology for calculating wormlike chain statistics in arbitrary D dimensions over all chain rigidities, from fully rigid to completely flexible. The structure of our exact analytical solution for the end-to-end distribution function for a wormlike chain in arbitrary D dimensions in Fourier-Laplace space (i.e., Fourier-transformed end position and Laplace-transformed chain length) adopts the form of an infinite continued fraction, which is advantageous for its compact structure and stability for numerical implementation. We then proceed to present a step-by-step methodology for performing the Fourier-Laplace inversion in order to make full use of our results in general applications. Asymptotic methods for evaluating the Laplace inversion (power-law expansion and Rayleigh-Schrödinger perturbation theory) are employed in order to improve the accuracy of the numerical inversions of the end-to-end distribution function in real space. We adapt our results to the evaluation of the single-chain structure factor, rendering simple, closed-form expressions that facilitate comparison with scattering experiments. Using our techniques, the accuracy of the end-to-end distribution function is enhanced up to the limit of the machine precision. We demonstrate the utility of our methodology with realizations of the chain statistics, giving a general methodology that can be applied to a wide range of biophysical problems.