Quantitative modeling of EGF receptor ligand discrimination via internalization proofreading.

The epidermal growth factor receptor (EGFR) is a central regulator of cell physiology that is stimulated by multiple distinct ligands. Although ligands bind to EGFR while the receptor is exposed on the plasma membrane, EGFR incorporation into endosomes following receptor internalization is an important aspect of EGFR signaling, with EGFR internalization behavior dependent upon the type of ligand bound. We develop quantitative modeling for EGFR recruitment to and internalization from clathrin domains, focusing on how internalization competes with ligand unbinding from EGFR. We develop two model versions: a kinetic model with EGFR behavior described as transitions between discrete states and a spatial model with EGFR diffusion to circular clathrin domains. We find that a combination of spatial and kinetic proofreading leads to enhanced EGFR internalization ratios in comparison to unbinding differences between ligand types. Various stages of the EGFR internalization process, including recruitment to and internalization from clathrin domains, modulate the internalization differences between receptors bound to different ligands. Our results indicate that following ligand binding, EGFR may encounter multiple clathrin domains before successful recruitment and internalization. The quantitative modeling we have developed describes competition between EGFR internalization and ligand unbinding and the resulting proofreading.

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.

Tube geometry controls protein cluster conformation and stability on the endoplasmic reticulum surface.

The endoplasmic reticulum (ER), a cellular organelle that forms a cell-spanning network of tubes and sheets, is an important location of protein synthesis and folding. When the ER experiences sustained unfolded protein stress, IRE1 proteins embedded in the ER membrane activate and assemble into clusters as part of the unfolded protein response (UPR). We use kinetic Monte Carlo simulations to explore IRE1 clustering dynamics on the surface of ER tubes. While initially growing clusters are approximately round, once a cluster is sufficiently large a shorter interface length can be achieved by 'wrapping' around the ER tube. A wrapped cluster can grow without further interface length increases. Relative to wide tubes, narrower tubes enable cluster wrapping at smaller cluster sizes. Our simulations show that wrapped clusters on narrower tubes grow more rapidly, evaporate more slowly, and require a lower protein concentration to grow compared to equal-area round clusters on wider tubes. These results suggest that cluster wrapping, facilitated by narrower tubes, could be an important factor in the growth and stability of IRE1 clusters and thus impact the persistence of the UPR, connecting geometry to signaling behavior. This work is consistent with recent experimental observations of IRE1 clusters wrapped around narrow tubes in the ER network.

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.

Theory of Nonequilibrium Free Energy Transduction by Molecular Machines.

Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared with more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.

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.

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.

Allocating dissipation across a molecular machine cycle to maximize flux.

Biomolecular machines consume free energy to break symmetry and make directed progress. Nonequilibrium ATP concentrations are the typical free energy source, with one cycle of a molecular machine consuming a certain number of ATP, providing a fixed free energy budget. Since evolution is expected to favor rapid-turnover machines that operate efficiently, we investigate how this free energy budget can be allocated to maximize flux. Unconstrained optimization eliminates intermediate metastable states, indicating that flux is enhanced in molecular machines with fewer states. When maintaining a set number of states, we show that-in contrast to previous findings-the flux-maximizing allocation of dissipation is not even. This result is consistent with the coexistence of both "irreversible" and reversible transitions in molecular machine models that successfully describe experimental data, which suggests that, in evolved machines, different transitions differ significantly in their dissipation.

Uniform spatial distribution of collagen fibril radii within tendon implies local activation of pC-collagen at individual fibrils.

Collagen fibril cross-sectional radii show no systematic variation between the interior and the periphery of fibril bundles, indicating an effectively constant rate of collagen incorporation into fibrils throughout the bundle. Such spatially homogeneous incorporation constrains the extracellular diffusion of collagen precursors from sources at the bundle boundary to sinks at the growing fibrils. With a coarse-grained diffusion equation we determine stringent bounds, using parameters extracted from published experimental measurements of tendon development. From the lack of new fibril formation after birth, we further require that the concentration of diffusing precursors stays below the critical concentration for fibril nucleation. We find that the combination of the diffusive bound, which requires larger concentrations to ensure homogeneous fibril radii, and lack of nucleation, which requires lower concentrations, is only marginally consistent with fully processed collagen using conservative bounds. More realistic bounds may leave no consistent concentrations. Therefore, we propose that unprocessed pC-collagen diffuses from the bundle periphery followed by local C-proteinase activity and subsequent collagen incorporation at each fibril. We suggest that C-proteinase is localized within bundles, at fibril surfaces, during radial fibrillar growth. The much greater critical concentration of pC-collagen, as compared to fully processed collagen, then provides broad consistency between homogeneous fibril radii and the lack of fibril nucleation during fibril growth.

Single file diffusion into a semi-infinite tube.

We investigate single file diffusion (SFD) of large particles entering a semi-infinite tube, such as luminal diffusion of proteins into microtubules or flagella. While single-file effects have no impact on the evolution of particle density, we report significant single-file effects for individually tracked tracer particle motion. Both exact and approximate ordering statistics of particles entering semi-infinite tubes agree well with our stochastic simulations. Considering initially empty semi-infinite tubes, with particles entering at one end starting from an initial time t = 0, tracked particles are initially super-diffusive after entering the system, but asymptotically diffusive at later times. For finite time intervals, the ratio of the net displacement of individual single-file particles to the average displacement of untracked particles is reduced at early times and enhanced at later times. When each particle is numbered, from the first to enter (n = 1) to the most recent (n = N), we find good scaling collapse of this distance ratio for all n. Experimental techniques that track individual particles, or local groups of particles, such as photo-activation or photobleaching of fluorescently tagged proteins, should be able to observe these single-file effects. However, biological phenomena that depend on local concentration, such as flagellar extension or luminal enzymatic activity, should not exhibit single-file effects.

PEX5 and ubiquitin dynamics on mammalian peroxisome membranes.

Peroxisomes are membrane-bound organelles within eukaryotic cells that post-translationally import folded proteins into their matrix. Matrix protein import requires a shuttle receptor protein, usually PEX5, that cycles through docking with the peroxisomal membrane, ubiquitination, and export back into the cytosol followed by deubiquitination. Matrix proteins associate with PEX5 in the cytosol and are translocated into the peroxisome lumen during the PEX5 cycle. This cargo translocation step is not well understood, and its energetics remain controversial. We use stochastic computational models to explore different ways the AAA ATPase driven removal of PEX5 may couple with cargo translocation in peroxisomal importers of mammalian cells. The first model considered is uncoupled, in which translocation is spontaneous, and does not immediately depend on PEX5 removal. The second is directly coupled, in which cargo translocation only occurs when its PEX5 is removed from the peroxisomal membrane. The third, novel, model is cooperatively coupled and requires two PEX5 on a given importomer for cargo translocation--one PEX5 with associated cargo and one with ubiquitin. We measure both the PEX5 and the ubiquitin levels on the peroxisomes as we vary the matrix protein cargo addition rate into the cytosol. We find that both uncoupled and directly coupled translocation behave identically with respect to PEX5 and ubiquitin, and the peroxisomal ubiquitin signal increases as the matrix protein traffic increases. In contrast, cooperatively coupled translocation behaves dramatically differently, with a ubiquitin signal that decreases with increasing matrix protein traffic. Recent work has shown that ubiquitin on mammalian peroxisome membranes can lead to selective degradation by autophagy, or 'pexophagy.' Therefore, the high ubiquitin level for low matrix cargo traffic with cooperatively coupled protein translocation could be used as a disuse signal to mediate pexophagy. This mechanism may be one way that cells could regulate peroxisome numbers.

Cluster coarsening on drops exhibits strong and sudden size-selectivity.

Autophagy, an important process for degradation of cellular components, requires the targeting of autophagy receptor proteins to potential substrates. Receptor proteins have been observed to form clusters on membranes. To understand how receptor clusters might affect autophagy selectivity, we model cluster coarsening on a polydisperse collection of spherical drop-like substrates. Our model receptor corresponds to NBR1, which supports peroxisome autophagy. We recover dynamical scaling of cluster sizes, but find that changing the drop size distribution changes the cluster-size scaling distribution. The magnitude of this effect is similar to how changing the spatial-dimension affects scaling in bulk systems. We also observe a sudden onset of size-selection of the remaining drops with clusters, due to clusters evaporating from smaller drops and growing on larger drops. This coarsening-driven size selection provides a physical mechanism for autophagy selectivity, and may explain reports of size selection during peroxisome degradation.

A storage-based model of heterocyst commitment and patterning in cyanobacteria.

When deprived of fixed nitrogen (fN), certain filamentous cyanobacteria differentiate nitrogen-fixing heterocysts. There is a large and dynamic fraction of stored fN in cyanobacterial cells, but its role in directing heterocyst commitment has not been identified. We present an integrated computational model of fN transport, cellular growth, and heterocyst commitment for filamentous cyanobacteria. By including fN storage proportional to cell length, but without any explicit cell-cycle effect, we are able to recover a broad and late range of heterocyst commitment times and we observe a strong indirect cell-cycle effect. We propose that fN storage is an important component of heterocyst commitment and patterning in filamentous cyanobacteria. The model allows us to explore both initial and steady-state heterocyst patterns. The developmental model is hierarchical after initial commitment: our only source of stochasticity is observed growth rate variability. Explicit lateral inhibition allows us to examine ΔpatS, ΔhetN, and ΔpatN phenotypes. We find that ΔpatS leads to adjacent heterocysts of the same generation, while ΔhetN leads to adjacent heterocysts only of different generations. With a shortened inhibition range, heterocyst spacing distributions are similar to those in experimental ΔpatN systems. Step-down to non-zero external fN concentrations is also investigated.

Protein translocation without specific quality control in a computational model of the Tat system.

The twin-arginine translocation (Tat) system transports folded proteins of various sizes across both bacterial and plant thylakoid membranes. The membrane-associated TatA protein is an essential component of the Tat translocon, and a broad distribution of different sized TatA-clusters is observed in bacterial membranes. We assume that the size dynamics of TatA clusters are affected by substrate binding, unbinding, and translocation to associated TatBC clusters, where clusters with bound translocation substrates favour growth and those without associated substrates favour shrinkage. With a stochastic model of substrate binding and cluster dynamics, we numerically determine the TatA cluster size distribution. We include a proportion of targeted but non-translocatable (NT) substrates, with the simplifying hypothesis that the substrate translocatability does not directly affect cluster dynamical rate constants or substrate binding or unbinding rates. This amounts to a translocation model without specific quality control. Nevertheless, NT substrates will remain associated with TatA clusters until unbound and so will affect cluster sizes and translocation rates. We find that the number of larger TatA clusters depends on the NT fraction f. The translocation rate can be optimized by tuning the rate of spontaneous substrate unbinding, [Formula: see text]. We present an analytically solvable three-state model of substrate translocation without cluster size dynamics that follows our computed translocation rates, and that is consistent with in vitro Tat-translocation data in the presence of NT substrates.

An equilibrium double-twist model for the radial structure of collagen fibrils.

Mammalian tissues contain networks and ordered arrays of collagen fibrils originating from the periodic self-assembly of helical 300 nm long tropocollagen complexes. The fibril radius is typically between 25 to 250 nm, and tropocollagen at the surface appears to exhibit a characteristic twist-angle with respect to the fibril axis. Similar fibril radii and twist-angles at the surface are observed in vitro, suggesting that these features are controlled by a similar self-assembly process. In this work, we propose a physical mechanism of equilibrium radius control for collagen fibrils based on a radially varying double-twist alignment of tropocollagen within a collagen fibril. The free-energy of alignment is similar to that of liquid crystalline blue phases, and we employ an analytic Euler-Lagrange and numerical free energy minimization to determine the twist-angle between the molecular axis and the fibril axis along the radial direction. Competition between the different elastic energy components, together with a surface energy, determines the equilibrium radius and twist-angle at the fibril surface. A simplified model with a twist-angle that is linear with radius is a reasonable approximation in some parameter regimes, and explains a power-law dependence of radius and twist-angle at the surface as parameters are varied. Fibril radius and twist-angle at the surface corresponding to an equilibrium free-energy minimum are consistent with existing experimental measurements of collagen fibrils. Remarkably, in the experimental regime, all of our model parameters are important for controlling equilibrium structural parameters of collagen fibrils.

Endoplasmic reticulum morphology regulation by RTN4 modulates neuronal regeneration by curbing luminal transport.

Cell functions rely on intracellular transport systems distributing bioactive molecules with high spatiotemporal accuracy. The endoplasmic reticulum (ER) tubular network constitutes a system for delivering luminal solutes, including Ca2+, across the cell periphery. How the ER structure enables this nanofluidic transport system is unclear. Here, we show that ER membrane-localized reticulon 4 (RTN4/Nogo) is sufficient to impose neurite outgrowth inhibition in human cortical neurons while acting as an ER morphoregulator. Improving ER transport visualization methodologies combined with optogenetic Ca2+ dynamics imaging and in silico modeling, we observed that ER luminal transport is modulated by ER tubule narrowing and dilation, proportional to the amount of RTN4. Excess RTN4 limited ER luminal transport and Ca2+ release, while RTN4 elimination reversed the effects. The described morphoregulatory effect of RTN4 defines the capacity of the ER for peripheral Ca2+ delivery for physiological releases and thus may constitute a mechanism for controlling the (re)generation of neurites.

Confinement of unliganded EGFR by tetraspanin nanodomains gates EGFR ligand binding and signaling.

The epidermal growth factor receptor (EGFR) is a central regulator of cell physiology. EGFR is activated by ligand binding, triggering receptor dimerization, activation of kinase activity, and intracellular signaling. EGFR is transiently confined within various plasma membrane nanodomains, yet how this may contribute to regulation of EGFR ligand binding is poorly understood. To resolve how EGFR nanoscale compartmentalization gates ligand binding, we developed single-particle tracking methods to track the mobility of ligand-bound and total EGFR, in combination with modeling of EGFR ligand binding. In comparison to unliganded EGFR, ligand-bound EGFR is more confined and distinctly regulated by clathrin and tetraspanin nanodomains. Ligand binding to unliganded EGFR occurs preferentially in tetraspanin nanodomains, and disruption of tetraspanin nanodomains impairs EGFR ligand binding and alters the conformation of the receptor's ectodomain. We thus reveal a mechanism by which EGFR confinement within tetraspanin nanodomains regulates receptor signaling at the level of ligand binding.

Reconciling cyanobacterial fixed-nitrogen distributions and transport experiments with quantitative modelling.

Filamentous cyanobacteria growing in media with insufficient fixed nitrogen (fN) differentiate some cells into heterocysts, which fix nitrogen for the remaining vegetative cells. Transport studies have shown both periplasmic and cytoplasmic connections between cells that could transport fN along the filament. Two experiments have imaged fN distributions along filaments. In 1974, Wolk et al found a peaked concentration of fN at heterocysts using autoradiographic techniques. In contrast, in 2007, Popa et al used nanoSIMS to show large dips at the location of heterocysts, with a variable but approximately level distribution between them. With an integrated model of fN transport and cell growth, we recover the results of both Wolk et al and of Popa et al using the same model parameters. To do this, we account for immobile incorporated fN and for the differing durations of labelled nitrogen fixation that occurred in the two experiments. The variations seen by Popa et al are consistent with the effects of cell-by-cell variations of growth rates, and mask diffusive gradients. We are unable to rule out a significant amount of periplasmic fN transport.

Heterocyst placement strategies to maximize the growth of cyanobacterial filaments.

Under conditions of limited fixed-nitrogen, some filamentous cyanobacteria develop a regular pattern of heterocyst cells that fix nitrogen for the remaining vegetative cells. We examine three different heterocyst placement strategies by quantitatively modelling filament growth while varying both external fixed-nitrogen and leakage from the filament. We find that there is an optimum heterocyst frequency which maximizes the growth rate of the filament; the optimum frequency decreases as the external fixed-nitrogen concentration increases but increases as the leakage increases. In the presence of leakage, filaments implementing a local heterocyst placement strategy grow significantly faster than filaments implementing random heterocyst placement strategies. With no extracellular fixed-nitrogen, consistent with recent experimental studies of Anabaena sp. PCC 7120, the modelled heterocyst spacing distribution using our local heterocyst placement strategy is qualitatively similar to experimentally observed patterns. As external fixed-nitrogen is increased, the spacing distribution for our local placement strategy retains the same shape, while the average spacing between heterocysts continuously increases.

Impact of global structure on diffusive exploration of organelle networks.

We investigate diffusive search on planar networks, motivated by tubular organelle networks in cell biology that contain molecules searching for reaction partners and binding sites. Exact calculation of the diffusive mean first-passage time on a spatial network is used to characterize the typical search time as a function of network connectivity. We find that global structural properties - the total edge length and number of loops - are sufficient to largely determine network exploration times for a variety of both synthetic planar networks and organelle morphologies extracted from living cells. For synthetic networks on a lattice, we predict the search time dependence on these global structural parameters by connecting with percolation theory, providing a bridge from irregular real-world networks to a simpler physical model. The dependence of search time on global network structural properties suggests that network architecture can be designed for efficient search without controlling the precise arrangement of connections. Specifically, increasing the number of loops substantially decreases search times, pointing to a potential physical mechanism for regulating reaction rates within organelle network structures.