Abstract
MreB is an actin homolog that is essential for coordinating the cell wall synthesis required for the rod shape of many bacteria. Previously we have shown that filaments of MreB bind to the curved membranes of bacteria and translocate in directions determined by principal membrane curvatures to create and reinforce the rod shape (Hussain et al., 2018). Here, in order to understand how MreB filament dynamics affects their cellular distribution, we model how MreB filaments bind and translocate on membranes with different geometries. We find that it is both energetically favorable and robust for filaments to bind and orient along directions of largest membrane curvature. Furthermore, significant localization to different membrane regions results from processive MreB motion in various geometries. These results demonstrate that the in vivo localization of MreB observed in many different experiments, including those examining negative Gaussian curvature, can arise from translocation dynamics alone.
https://doi.org/10.7554/eLife.40472.001Introduction
The role of membrane curvature in influencing the cellular location and function of proteins has been increasingly appreciated (McMahon and Gallop, 2005; Zimmerberg and Kozlov, 2006; Baumgart et al., 2011). In cells, membrane curvature can both be induced by proteins that bind to membranes as well as recruit proteins that bind to this curvature (Peter et al., 2004; Renner and Weibel, 2011; Wang and Wingreen, 2013; Ursell et al., 2014; Hussain et al., 2018; Wu et al., 2018). Here, we focus on MreB, an actin homolog essential for coordinating the cell wall synthesis required for the rod shape of many bacteria (Jones et al., 2001). We examine how MreB orientation and motion along directions of principal curvature affect its localization in different cellular geometries. MreB filaments polymerize onto membranes (Figure 1A) (Salje et al., 2011; van den Ent et al., 2014; Hussain et al., 2018), creating short filaments that move along their lengths in live cells. Viewed dynamically, MreB filaments are seen to rotate around the rod width, a motion powered by the activity of associated cell wall synthesis enzymes (Figure 1B and Figure 1—video 1) (Garner et al., 2011; DomínguezEscobar et al., 2011; van Teeffelen et al., 2011; Reimold et al., 2013; Hussain et al., 2018). The orientation of MreB filaments coincides with the direction of their motion, as filaments move along the direction in which they point (Figure 1C) (Hussain et al., 2018; Olshausen et al., 2013). However, studies examining the localization of MreB in kymographs (Ursell et al., 2014) or at single time points (Bratton et al., 2018) have found that, in Escherichia coli, MreB filaments are enriched at regions of negative Gaussian curvature or small mean curvature (Figure 1D). This prompts the question of how the collective motion of MreB filaments could affect, or give rise to, their enrichment. As filaments are constantly moving, their instantaneous localization to any one point is transient: a typical filament (~250 nm long) in Bacillus subtilis moves through a 1 µm^{2} region in ~30 s. To understand how the previously observed enrichment at negative Gaussian curvatures could arise from the dynamics of filaments moving around the cell, we sought to relate the binding and motion of MreB filaments to their distribution in different geometries.
In previous work (Hussain et al., 2018), we demonstrated that MreB filaments orient and translocate along directions of largest principal curvature inside differently shaped B. subtilis cells and liposomes. To understand how MreB filaments orient along different membrane curvatures and how their motion along this orientation affects their cellular localization, we first model the mechanics of MreBmembrane binding and provide a quantitative description of how MreB filaments bind both in vivo and in vitro. Next, we model the curvaturedependent motion of MreB filaments in different geometries and examine how this motion affects their distribution to different membrane regions. Strikingly, we find that the dynamics of MreB translocation alone, without requiring any intrinsic preference of MreB filaments for curved regions of the cell, results in differential enrichment of MreB filaments at regions of negative Gaussian curvature or small mean curvature similar to those observed in cells.
Results
Mechanics of binding
We first model how inwardly curved (Salje et al., 2011; van den Ent et al., 2014; Hussain et al., 2018) MreB filaments bind and orient on membranes. Previous theoretical studies have modeled the binding of protein filaments to membranes and demonstrated that binding conformations can be influenced by both filament thickness and twist. In a seminal theoretical work, Wang and Wingreen modeled twisted bundles of MreB approximately sixfold thicker than typical filaments (Wang and Wingreen, 2013). This study nicely demonstrated that the mechanics of binding alone could orient bundles and suggested that bundle length could be limited by twist. Another elegant theoretical study by Quint et al. demonstrated that twisted filaments of varying rigidities could bind to regions of negative Gaussian curvature in manners that are particularly energetically favorable (Quint et al., 2016). While it is intriguing to examine the effects of filament rigidity and twist on general filament systems (discussed below), here we focus on modeling thin, inwardly curved MreB filaments with no twist. We focus on these parameters as they reflect the observations of all available in vitro studies of membraneassociated MreB: three different cryoelectron microscopy studies have shown that membranebound MreB filaments are flat and untwisted, binding to the membrane on one filament face (Figure 1A) (Salje et al., 2011; van den Ent et al., 2014; Hussain et al., 2018). These studies also suggest that inward curvature, and not twist, limits filament length: MreB filaments are short when polymerized onto nondeforming planarsupported lipid bilayers, but become extremely long when polymerized inside deformable liposomes (Salje et al., 2011; van den Ent et al., 2014; Hussain et al., 2018). The model we present here extends our previous work (Hussain et al., 2018) and demonstrates that thin, untwisted filaments orient robustly along directions of largest principal curvature, thus providing a generic mechanism for orienting their motion.
We model an MreB filament as a polymer which binds linearly along its length to a membrane in an energetically favorable manner, namely via burial of hydrophobic residues on one face (Figure 2A; see Figure 2—figure supplement 1 for additional details) (Salje et al., 2011). We assume the filament to be a curved, cylindrical, linearelastic rod which is free to bend to maximize membrane interaction, so that its elastic energy of deformation is ${E}_{\mathrm{b}\mathrm{e}\mathrm{n}\mathrm{d}}=(\pi Y{r}_{f}^{4}/8)\times \int (\kappa {\kappa}_{s}{)}^{2}d\ell$, where ${r}_{f}$ is the filament crosssectional radius, $Y$ is the filament Young’s modulus, $\kappa $ is the curvature of the deformed state, ${\kappa}_{s}$ is the intrinsic filament curvature, and the integration is over the filament length (Landau and Lifshitz, 1970). We will minimize the free energy change due to binding with respect to the deformed curvature, $\kappa $.
Next, we assume an isotropic, fluid, bilayer membrane, where there are no inplane shears and the only inplane deformations are compressions and expansions (Safran, 2003). The membrane free energy is given by the Helfrich form,
where ${k}_{b}$ is the bending rigidity of the membrane, ${k}_{t}$ is the saddlesplay modulus of the membrane, ${H}_{s}$ is the spontaneous curvature, $\gamma $ is the membrane surface tension, $p$ is the pressure difference across the membrane, $H$ and $K$ are the mean and Gaussian curvatures of the membrane surface, $S$, respectively, and $dA$ and $dV$ denote area and volume elements, respectively (Helfrich, 1973; Safran, 2003; Zhongcan and Helfrich, 1987; Zhongcan and Helfrich, 1989). For the small membrane deformations considered in this work, we assume that excess phospholipids can be freely added to the membrane to compensate for stretching, so that the membrane surface tension $\gamma =0$ (Safran, 2003). Note that a nonzero surface tension would make the preferred binding orientation determined below more energetically favorable. For simplicity, we also assume ${H}_{s}=0$ and note that the case of a nonzero ${H}_{s}$ can be considered similarly. Finally, we assume that the membrane surface can be parameterized in the Monge gauge by a function $h=h(x,y)$, where $x$ and $y$ are real numbers and terms of quadratic order or higher in the gradient of $h$, $\nabla h$, are neglected. As shown in Appendix 1, the mechanical energy associated with membrane deformations is determined by the solution of the shape equation
which is similar to the equilibrium equation of a thin plate (Ventsel and Krauthammer, 2001; Timoshenko and WoinowskyKrieger, 1959). Here ${\mathrm{\Delta}}^{2}$ is the biharmonic operator and Equation (2) is subject to Dirichlet boundary conditions enforcing continuity of mean curvature and surface height in a manner compatible with the deformed filament. Equation (2) can then be decomposed as two Poisson equations, each with Dirichlet boundary conditions, and solved numerically using the finite element method (Appendix 1).
The free energy change due to filament binding is determined by ${E}_{\mathrm{bend}}$, $\mathcal{F}$, and the solution of Equation (2). For characteristic parameter values relevant to the binding of MreB filaments to bacterial membranes, as summarized in Supplementary file 1, our model predicts a preferred, circumferential orientation of MreB binding in a rodshaped cell. This result arises because the intrinsic curvature of filaments is smaller than characteristic cell radii. While the filament bends to conform to the membrane for physiological values of $p$—as shown in previous work (Hussain et al., 2018)—and grossly deforms the membrane for small values of $p$, the preferred binding orientation is robust to changes in $p$ (Figure 2B). In fact, across a wide range of parameter values including the filament bending rigidity ($B=\pi Y{r}_{f}^{4}/4$), the filament intrinsic curvature (${\kappa}_{s}$), and the membrane pressure difference ($p$), a preferred binding orientation exists and coincides with the direction of largest principal curvature for any membrane which is less curved than the filament. In the case of $p=0$, as discussed in Appendix 1, the prediction that MreB binding induces large membrane deformations is consistent with cryoelectron microscopy images of MreB binding to vesicles (Figure 1A) (Salje et al., 2011; van den Ent et al., 2014). Importantly, the energetic penalties for deviatory binding conformations are larger than the energy of thermal fluctuations across a large range of $p$, suggesting the empirically observed variation in binding orientation (Figure 1C) to be caused by other sources of stochasticity. The energetic penalties are also decreased in wider membranes for both small and physiological values of $p$, the latter of which is consistent with the gradual widening of the distributions of MreB trajectory angles in wider B. subtilis protoplasts (Hussain et al., 2018).
In general, the mechanics of filament binding are well described by the pressure difference across the membrane and the filament bending energy, which can be viewed as order parameters that largely dictate whether the filament predominantly bends the membrane, bends to conform to the membrane, or both. An approximate phase diagram for MreB binding to any membrane which is less curved than the filament can be determined (Figure 2C and Appendix 1). Below, we suppose the membrane surface to be less curved than the filament—so that it is always energetically favorable to orient along directions of largest membrane curvature—and model filament translocation along these directions.
Dynamics of translocation
We next examine how the translocation of MreB filaments, once bound to the membrane, affects their distribution in different geometries. Inside cells, MreB filaments move along the membrane in the direction of their orientation (Hussain et al., 2018; Olshausen et al., 2013) (Figure 1C). This directional and processive motion is driven by cell wall synthesis (Sliusarenko et al., 2010; van Teeffelen et al., 2011; Olshausen et al., 2013), and filaments may reorient according to different membrane geometries (Hussain et al., 2018). Thus, highlybent MreB filaments translocate along the direction of largest curvature, a direction that minimizes the energetic cost of binding. As discussed below, this hypothesis is supported by observations that (1) MreB moves circumferentially in live B. subtilis cells, but this motion becomes disoriented if cells become round, (2) circumferential motion is reestablished when round cells are confined into rods, (3) MreB moves directionally in bulges protruding from round cells, and (4) MreB filaments rapidly translocate out of poles in rods, reorienting when filaments reach the cylindrical bulks (Hussain et al., 2018). Assuming translocation on a static surface, we may model the trajectories of filaments as biased random walks as follows (with more details provided in Appendix 1). The case of a dynamical surface, as expected for MreBdirected growth, can be considered similarly. Note that a ‘biased random walk’ refers to a succession of random steps which may be processive: while the meansquared displacement of a filament will be approximately quadratic, and not linear, in time, the processive motion we consider is random only because the translocation direction can deviate from directions of largest membrane curvature due to sources of stochasticity (Figure 1C). We will show that our biased random walk model of filament trajectories leads to predictions of MreB localization.
We consider the membrane as a parametric surface, $\mathbf{\mathbf{r}}=r(u,v)$, embedded in threedimensional space (${\mathbb{R}}^{3}$) with surface coordinates $u$ and $v$ and a filament as a point on this surface which, at any moment in time, translocates along the largest principal direction $\mathbf{\mathbf{d}}$—that is, the direction of largest curvature of the surface. As $\mathbf{\mathbf{d}}$ is a vector in ${\mathbb{R}}^{3}$, arbitrarily moving in the direction of $\mathbf{\mathbf{d}}$ may move the filament off of the surface. To define the translocation consistently, we set $\eta ={\mathrm{cos}}^{1}\frac{\mathbf{d}\cdot {\mathbf{r}}_{\theta}}{\mathbf{d}\cdot {\mathbf{r}}_{\theta}}$, where $\eta$ is an angular deviation from the largest principal direction on the surface introduced by possible sources of stochasticity, the modified direction corresponds to an angle $\theta $ relative to the $u$axis in parametric coordinates, $\mathbf{r}}_{\theta}\in {\mathbb{R}}^{3$ is the derivative of $\mathbf{r}$ in the direction of $\theta $, and distances are defined by the surface metric. Translocating along an angle $\theta $ with respect to the $u$axis in $(u,v)$coordinates then ensures that the filament remains on the surface, and the direction of translocation corresponds to that on a patch of $\mathbf{\mathbf{r}}$.
As a discretetime flow in $(u,v)$coordinates, and with suitable units of time so that the filament may reorient at every timestep, the 2D equation of filament motion is
where ${X}_{n}$, ${\mathrm{\ell}}_{n}$, and $\theta}_{n$ are the position, step size, and translocation angle, respectively, of the filament at a timestep $n$. Here ${\theta}_{n}$ is the value of $\theta $ computed at the surface point corresponding to ${X}_{n}$ and assuming $\eta \sim \mathcal{\mathcal{N}}(0,{\sigma}^{2})$—that is, the angular noise is normallydistributed, with mean zero and a variance, ${\sigma}^{2}$, to be inferred from data—and note that the translocation noise, $\sigma $, may depend on quantities such as the principal curvatures, as discussed later. ${\chi}_{n}$ is a random sign, which accounts for the possibility of both lefthanded and righthanded translocation, and may not substantially vary in $n$ if the filament does not backtrack, as is assumed for the remainder of this work. We assume that ${\mathrm{\ell}}_{n}$ satisfies an integral equation which relates it to a constant and finite filament step size, $L$, on the surface (Appendix 1), and note that inertia in MreB motion, as measured previously by the velocity autocorrelation (Kim et al., 2006; van Teeffelen et al., 2011), is modeled by the finite step size. Note that, when both principal curvatures are equal, so that $\mathbf{\mathbf{d}}$ is not well defined, we will assume steps of uniformly random angles and size $L$ on the surface. Furthermore, for finite step sizes, Equation (3) extrapolates the translocation direction in surface coordinates: while motion along the largest principal direction can be maintained irrespective of parameterization by parallel transport, we anticipate the distinction from Equation (3) to be insignificant for many geometries due to the small step sizes considered in this work. Determining the probability distribution of $X$ then determines the probability distribution of the filament on $\mathbf{\mathbf{r}}$, and this can be done analytically and numerically for several geometries as discussed below.
We next consider the dynamics of activating and deactivating filaments as follows. We suppose that a filament may be activated at a position $X$ at any timestep at a constant rate $k\ge 0$ with probability proportional to the membrane surface area, $dA(X)$, and deactivated at a constant rate $\lambda \ge 0$ which determines the filament’s processivity—that is, the mean number of steps that a filament takes on the membrane surface before becoming inactive (Figure 3A). The case of $k$ being dependent on fields, such as mechanical strains (Wong et al., 2017), can be considered similarly but is not necessary for the results below. An ensemble of filaments produced by such dynamics will exhibit filament numbers, ${N}_{F}$, that vary in space and time, and likewise for the filament concentration ${C}_{F}(X,n)={N}_{F}(X,n)/dA(X)$. Below, we discuss characteristic parameter values relevant to MreB and show that the dynamics of Equation (3) gives rise to localization. We then examine the model in detail and describe how localization depends on different parameters of the model.
Implications to MreB localization
Previous fluorescence microscopy measurements provide estimates for the step size ($L$), deactivation rate ($\lambda$), and translocation noise ($\sigma$) of MreB filaments in cells. We assume $L$ to be $200\phantom{\rule{thinmathspace}{0ex}}\mathrm{n}\mathrm{m}$ as a modeling choice, but show in Appendix 1 that the results discussed below are qualitatively similar for significantly larger $L$ ($\sim 2\phantom{\rule{thinmathspace}{0ex}}\mu \mathrm{m}$). Similar experiments have estimated the persistence time of MreB filaments in E. coli as ~5 min (Ursell et al., 2014), while a characteristic translocation noise of $\sigma \approx 0.3\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}\mathrm{r}\mathrm{a}\mathrm{d}$ in B. subtilis has been found separately by (1) measuring filament trajectory angles relative to the midline and (2) measuring binding angles in confined protoplasts (Figure 1C) (Hussain et al., 2018). While we assume the values of λ and σ to be based on these measurements, we examine the effects of varying λ and σ in the following section. Furthermore, in recent studies, rodshaped cells have been perturbed to be in geometries other than a spherocylinder (Ursell et al., 2014; Wong et al., 2017; Hussain et al., 2018; Renner et al., 2013; Amir et al., 2014). As the distribution of MreB filament angles gradually becomes broader as B. subtilis cells become wider (Hussain et al., 2018), it may also be reasonable to suppose that $\sigma $ depends on the difference, $\mathrm{\Delta}c$, of principal curvatures at the location of any MreB filament: $\sigma =\alpha {(\mathrm{\Delta}c)}^{1}$, where $\alpha \approx 0.6\phantom{\rule{thinmathspace}{0ex}}\mathrm{r}\mathrm{a}\mathrm{d}\cdot \mu {\mathrm{m}}^{1}$ is a constant of proportionality determined by experimental data. While all our results pertaining to MreB below assume this dependence as to be consistent with data, we show in Appendix 1 that our results are similar for different dependencies of σ on $\mathrm{\Delta}c$.
Given the aforementioned parameters, Equation (3) leads to predictions for the statistics of the filament position ($X$) and the ensuing filament concentration (${C}_{F}$) across different membrane geometries. For a cylindrical cell, analytical expressions for the statistics of $X$ show that translocation noise does not significantly affect the mean or variance of the circumferential displacement of a filament (Appendix 1). In contrast, the value of $\sigma \approx 0.3\phantom{\rule{thinmathspace}{0ex}}\mathrm{r}\mathrm{a}\mathrm{d}$ corresponds to a standard deviation of approximately $0.2\mu \text{m}$ for the axial displacement of a filament per hoop of wall material inserted. This value is consistent with experimental measurements (Figure 3B), showing that deviations from a circumferential translocation direction can significantly contribute to wall insertions in the axial direction and disordered wall architecture.
MreB filaments have been observed to be depleted from the hemispherical poles of spherocylindrical cells compared to the cylindrical bulks (Kawazura et al., 2017; Ursell et al., 2014). Observations of filament dynamics revealed a possible explanation: MreB filaments reorient rapidly in, and translocate out of, the poles and into the bulks, where motion then becomes aligned (Hussain et al., 2018). Consistent with this observation, simulations of Equation (3) on a spherocylindrical surface show that the concentration of filaments in the bulk is enhanced (Figure 3C). The average filament concentration is predicted to be approximately twofold higher in the bulk than the poles, in agreement with experimental measurements in E. coli (Ursell et al., 2014). Simulations of Equation (3) on a toroidal surface are also quantitatively consistent with prior measurements of MreB fluorescence in E. coli cells confined to donutshaped microchambers, which have shown that MreB intensity is increased at the inner edges by a factor of ∼1.1 relative to the midlines (Figure 3D) (Wong et al., 2017). For a spherocylinder, filament enrichment arises because the cylindrical bulk retains filaments: oriented motion is preserved in the bulk, while disordered motion at the poles eventually becomes ordered. In contrast, filament enrichment arises in a curved cell because filaments become uniformly distributed along circumferential hoops. The smaller arclength along the inner edge then results in a greater density of filaments.
In our previous study, we found that MreB rotation and localization at small protrusions in B. subtilis protoplasts preceded rod shape generation from these protrusions (Hussain et al., 2018). To model the geometry observed in these experiments, we consider a cylindrical body with a protruding bulge in which filament trajectories become parallel to the cylinder long axis. Simulations of Equation (3) on this geometry reveal that the filament concentration is larger in the bulge and that the predicted enrichment is quantitatively consistent with the MreB enhancement observed in bulged cells, without any fitting parameters (Figure 3E and Figure 3—figure supplement 1). Similar to the case of a spherocylinder, localization arises due to the bulge attracting filaments. The dynamics of Equation (3) therefore results in localization which contributes to de novo generation of rod shape.
Finally, previous work has examined MreB localization in E. coli cells (1) with submicronscale shape fluctuations or (2) confined in sinusoidal chambers (Ursell et al., 2014; Shi et al., 2017; Bratton et al., 2018). The empirically observed magnitudes of MreB enrichment at regions of negative Gaussian curvature or small mean curvature in these studies are consistent with our modeling. To model the cell shapes observed in these experiments, we consider filament translocation on a geometry with both negative and positive Gaussian curvatures and undulations of smaller wavelengths than the surface size (Figure 3F and Figure 3—figure supplement 2). As discussed in Appendix 1, the Gaussian and mean curvatures in this geometry are positively correlated and consistent with experimental observations (Ursell et al., 2014). For this geometry, filament translocation results in increased values of concentration (${C}_{F}$) at regions of negative Gaussian curvature or small mean curvature (Figure 3F). This effect arises because the principal curvatures away from these regions reorient filaments axially, instead of circumferentially, so that regions of negative Gaussian curvature or small mean curvature attract filaments. Furthermore, the magnitude of this enhancement is consistent with the amount of MreB enrichment observed (Figure 3G and Figure 3—figure supplement 3), demonstrating that translocation dynamics alone can negatively correlate filament concentration with Gaussian or mean curvature in cells with similar short wavelength undulations.
Dependence of localization on processivity and Gaussian curvature
As we anticipate our model to be applicable to general filament systems, we now explore the response of the filament concentration (${C}_{F}$) to (1) different parameter values and (2) other geometries. We fix the filament step size ($L$) and suppose the translocation noise (σ) and the deactivation rate (λ) to be constants which are varied within a broad range. We show in Appendix 1 that, for any value of processivity and zero translocation noise, ${C}_{F}$ is uniform over the surface of an ellipsoid, as is generally the case for any surface when the processivity is small (or, equivalently, λ is large). In contrast, in the case of small λ corresponding to large processivity—a limiting case that is relevant to MreB—and over a range of σ, ${C}_{F}$ is larger at the inner edge of a torus, at the inner edge of a helix, and at the tips of an ellipsoid (Figure 4A and Figure 4—figure supplements 1 and 2). As discussed above, localization occurs geometrically in these cases due to the filament number (${N}_{F}$) becoming uniform over the surface and spatial variations in the surface area element. The magnitude of the localization can be quantitatively predicted by geometric parameters alone (Appendix 1). The mechanism underlying localization is different for a spherocylinder or a bulged cylinder, for which surface regions attract filaments. Nevertheless, a nonzero processivity is required for localization even in geometries which attract filaments (Figure 4B and Figure 4—figure supplements 1 and 2).
Since filament enrichment depends on both processivity and geometry, we wondered if the localization of processive filaments always correlates with the Gaussian or mean curvatures, regardless of overall geometry. Although Figure 3G demonstrates that filament enrichment correlates with negative Gaussian curvature or small mean curvature in a specific, undulating geometry, this correlation is reversed in bulged cylinders (Figure 3E). Furthermore, Figure 4C illustrates a surface of zero Gaussian curvature exhibiting regions which attract filaments, as filaments change from moving circumferentially to moving axially in such regions (see Figure 4—figure supplement 3 for additional details). Examining Equation (3) on different surfaces therefore shows that ${C}_{F}$ need not depend on Gaussian curvature at all, and the dynamics modeled in this work cannot act as a generic mechanism for sensing Gaussian curvature.
Finally, while large filament bundles or twist have not been observed in MreB filaments reconstituted in vitro (Salje et al., 2011; van den Ent et al., 2014; Hussain et al., 2018), it is possible that general filament systems could exhibit these properties (Wang and Wingreen, 2013; Quint et al., 2016). The binding and activation of twisted filaments may also depend on membrane Gaussian curvature, as previously demonstrated (Quint et al., 2016). We systematically explore the effects of varying filament bending rigidity, filament twist, and Gaussian curvaturedependent activation in Appendix 1, where we show that our model predictions remain largely robust across a broad range of these parameters (Figure 4—figure supplement 4). Thus, we expect filament dynamics to contribute to localization in different filament systems, regardless of the details of filament rigidity, twist, and other parameters of our model.
Discussion
An outstanding problem in bacterial physiology has been to understand how short and disconnected filaments distribute themselves within cells to conduct different cellular functions (Eun et al., 2015). In this work, we have examined an aspect of this problem by modeling the direct binding of protein filaments to membranes and the curvaturebased translocation of an ensemble of such filaments. Our results provide a theoretical framework for prior work examining MreB dynamics and localization (Hussain et al., 2018; Salje et al., 2011; Wong et al., 2017; Ursell et al., 2014; Shi et al., 2017; Bratton et al., 2018; Renner et al., 2013). Furthermore, our results are consistent with the cellular localization observed in all these works and demonstrate that filament motion alone can correlate enrichment with Gaussian curvature in specific geometries. Our work may be extended by modeling an evolving membrane surface, as expected for MreBdirected growth, and it would be intriguing to explore whether and how principal curvaturebased translocation contributes to determining cell width.
The main contribution of this work is to show that the biological results of MreB localization, as observed in many different experiments involving a range of cell shapes (Hussain et al., 2018; Wong et al., 2017; Ursell et al., 2014; Shi et al., 2017; Bratton et al., 2018), can arise from processivity and principal curvaturedependent motion alone. Our study therefore helps to unravel how rod shape formation may be achieved through subcellularscale mechanisms (Amir and van Teeffelen, 2014; Shi et al., 2018; Surovtsev and JacobsWagner, 2018). More broadly, our work shows that the localization of translocating protein filaments can vary significantly depending on membrane geometry. This paves the way for exploring similar behavior in other contexts, such as bacterial cytokinesis and eukaryotic membrane trafficking and transport. For example, in bacterial cytokinesis, filaments of the tubulin homolog FtsZ assemble at, and treadmill around, the septum, a process which directs the insertion of new PG and constricts the cell (BissonFilho et al., 2017; Yang et al., 2017). Like MreB, FtsZ filaments are curved and could orient along the largest principal direction on membranes through bending alone (Osawa et al., 2009; Erickson et al., 2010). Treadmilling along such directions would then allow filaments to drive PG synthesis circumferentially at the septum.
Aside from MreB and FtsZ, septins, BARdomaincontaining proteins, dynamins, and endoproteins are known to exhibit similar, curvaturedependent membrane binding behaviors important for membrane trafficking, growth, and movement in both prokaryotes and eukaryotes (Baumgart et al., 2011; Zimmerberg and Kozlov, 2006; McMahon and Gallop, 2005; Peter et al., 2004; Low and Löwe, 2006; Raiborg and Stenmark, 2009; Teo et al., 2006; Kostelansky et al., 2007). Like MreB filaments, many such proteins sense membrane curvature through mechanical deformations of either the membrane or the protein itself. Unlike MreB or FtsZ, these proteins do not translocate; rather, they often induce membrane curvature to facilitate downstream processes. One example is BARdomaincontaining proteins, which scaffold higherorder assemblies of dynamin that actively constrict for vesicle scission (McMahon and Gallop, 2005). It would be interesting to apply the methods introduced here to this and other biological systems where molecules are known to bind to membranes or sense membrane curvature. These systems are widespread and involved in pathogenesis (Baumgart et al., 2011; Frost et al., 2009), cell division (Renner et al., 2013; Ramamurthi and Losick, 2009; Ramamurthi et al., 2009; Frost et al., 2009), intracellular trafficking (Zimmerberg and Kozlov, 2006; McMahon and Gallop, 2005; Raiborg and Stenmark, 2009; Ford et al., 2002; Frost et al., 2009; Römer et al., 2007), and cell migration (Frost et al., 2009; Zhao et al., 2013). The mathematical model introduced in this work, which requires minimal assumptions as to how filaments bind to and translocate on membranes, should be widely applicable to these and other broader contexts.
Appendix 1
1. Mechanics of binding
1.1. Model of a protein filament binding to a membrane
We consider the protein as a filament with monomeric subunits that bind to a membrane in an energetically favorable manner, such as burial of hydrophobic residues (Hussain et al., 2018). When a filament binds to a membrane, an energetic cost ${E}_{\mathrm{def}}({\mathrm{\ell}}_{b})$ is associated to deformations which deviate from a position at mechanical equilibrium, while the free energy may be lowered by an amount ${E}_{\mathrm{int}}({\mathrm{\ell}}_{b})$ due to interaction (Figure 2—figure supplement 1a). Both the deformation and interaction energies are expressed as functions of the bound filament length, ${\mathrm{\ell}}_{b}$, which is less than or equal to the total filament length, ${L}_{f}$. We wish to minimize the free energy due to filament binding, $\mathrm{\Delta}E={E}_{\mathrm{def}}{E}_{\mathrm{int}}$. If $\mathrm{\Delta}E({\mathrm{\ell}}_{b})$ is negative, then it is energetically favorable for the filament to bind to the membrane along a length ${\mathrm{\ell}}_{b}$. We estimate ${E}_{\mathrm{int}}$ as
where ${N}_{\mathrm{int}}$ denotes the total number of membrane binding sites of the filament and ${\epsilon}_{0}$ denotes an independent and additive single binding site energy, which is given for MreB along with other parameter values in Supplementary file 1.
We assume that the binding sites are arranged linearly along the filament, and in particular, that the filament is not twisted. In this case, it suffices to account only for filament bending, and we may decompose the deformation energy ${E}_{\mathrm{def}}$ into the bending energy of the filament, ${E}_{\mathrm{bend}}$, and the deformation energy of the membrane, $E}_{\mathrm{m}\mathrm{e}\mathrm{m}}:{E}_{\mathrm{d}\mathrm{e}\mathrm{f}}={E}_{\mathrm{b}\mathrm{e}\mathrm{n}\mathrm{d}}\phantom{\rule{thinmathspace}{0ex}}+\phantom{\rule{thinmathspace}{0ex}}{E}_{\mathrm{m}\mathrm{e}\mathrm{m}$. With notation as in the main text, we model the filament as a curved, cylindrical elastic rod with a circular crosssection of radius ${r}_{f}$ and curvature $1/{R}_{s}$, so that the elastic energy density per unit length of bending the filament from a curvature of $1/{R}_{s}$ to a curvature of $1/R$ is
where $Y$ is the elastic modulus of the filament and $B=\pi Y{r}_{f}^{4}/4$ is its flexural rigidity (Landau and Lifshitz, 1970). The resulting filament bending energy is ${E}_{\mathrm{bend}}={\epsilon}_{\mathrm{bend}}{\mathrm{\ell}}_{b}$. For simplicity, we have assumed the filament to be bent uniformly, but the case of a curvature which varies with position along the filament length can be considered similarly.
As stated in the main text, we assume an isotropic, fluid bilayer membrane, where there is no inplane shear modulus and the only inplane deformations are compressions and expansions. The membrane free energy assumes the form of Equation (1) in the main text. The mechanical energy required to bend a membrane from a surface ${S}_{0}$, with a mean curvature ${H}_{0}$, to a surface $S$ is then the difference of the corresponding free energies:
where the surface integrals of the Gaussian curvature are topological invariants by the GaussBonnet theorem and therefore cancel in the difference, and the volume integral is understood to be a difference of the volumes in the deformed and undeformed states.
Minimizing $\mathrm{\Delta}E$ requires the minimization of ${E}_{\mathrm{mem}}$ given some value $1/R$ of the deformed filament curvature. To minimize ${E}_{\mathrm{mem}}$ in Equation (S3), we assume that the surface $S$ can be parameterized in the Monge gauge $h=h(x,y)$, where $(x,y)\in {\mathbb{R}}^{2}$, and furthermore that $\nabla h\ll 1$: this means that the membrane surface is not excessively curved or kinked. We assume the same for the undeformed surface ${S}_{0}$, which is parameterized by a function ${h}_{0}$ in the Monge gauge. In the case of binding to a cylindrical membrane, we may, for instance, take the undeformed surface to correspond to a cylinder with radius ${R}_{\mathrm{cell}}$:
In the Monge gauge, the mean curvature can be expanded as $H=\frac{1}{2}{\nabla}^{2}h+O[{(\nabla h)}^{2}]$, where the big$O$ notation signifies $H\frac{1}{2}{\nabla}^{2}h\le M{(\nabla h)}^{2}$ when $0\phantom{\rule{thinmathspace}{0ex}}<(\mathrm{\nabla}h{)}^{2}\phantom{\rule{thinmathspace}{0ex}}<\delta$ for some positive numbers $\delta $ and $M$. The membrane bending energy in Equation (S3) can then be rewritten with the Laplacian, $\mathrm{\Delta}$, as
for some domain $\mathrm{\Omega}\subset {\mathbb{R}}^{2}$ of $h$ and ${h}_{0}$ not containing the domain $U$ of the filament surface (Figure 2—figure supplement 1a). Setting the first variation of $\mathcal{F}[h]$ to zero, we find that the equilibrium membrane shape is given by the solution of the shape equation
where ${\mathrm{\Delta}}^{2}$ is the biharmonic operator. Equation (S6) is subject to the Dirichlet boundary conditions
Here $\varphi $ and $\psi $ are indicator functions defined by their values on the boundary of $\mathrm{\Omega}$, $\mathrm{\partial}\mathrm{\Omega}$. In the case of a cylindrical membrane, for instance,
where, as above, ${h}_{0}$ parameterizes the undeformed surface ${S}_{0}$, ${p}_{0}(y)$ is a quadratic function describing the values of the filament height at the curve $\partial U$ parameterizing the binding region, and ${C}_{0}$ is the mean curvature of the filament along $\partial U$. Thus, the first condition in Equation (S7) comes from imposing continuity of membrane height with respect to the filament surface, while the second condition comes from imposing continuity of mean curvature. In live cells, we treat the periplasm as a rigid, undeformable body (Hussain et al., 2018), so that, for instance, ${p}_{0}(y)\ge {h}_{0}(x,y)$ for all $(x,y)\in \mathrm{\Omega}$ when ${h}_{0}$ assumes the form of Equation (S4) above. With the boundary conditions of Equation (S7), Equation (S6) can be conveniently decoupled as two Poisson equations, each with Dirichlet boundary conditions:
Since any solution to the Poisson equation with Dirichlet boundary conditions is unique, the decomposition above yields a unique solution for $h$.
As the foregoing considerations assume a fixed $\mathrm{\Omega}$ in determining h, the size of $\mathrm{\Omega}$ is an additional variable that must be considered. Since shape space is infinitedimensional, determining $h$ for an arbitrarily large $\mathrm{\Omega}$ does not necessarily imply that the global minimum of $\mathrm{\Delta}E$ is achieved; neither does it necessarily determine the appropriate decay length of the indentation, since $h(x,y)={h}_{0}(x,y)$ is generally not a solution to Equation (S9). It is possible that $\mathrm{\Delta}E$ may be minimized at a finite $\mathrm{\Omega}$. Due to the boundary conditions (Equation (S7) and (S8)), solutions of Equation (S9) over finite $\mathrm{\Omega}$ are continuous, with continuous mean curvatures, and could in fact be physically plausible. This subtlety can be addressed by choosing $\mathrm{\Omega}$ so that the numerically computed value of $\mathrm{\Delta}E$ is minimal among differently sized $\mathrm{\Omega}$, as discussed below.
1.2. Finite element solutions of the shape equation
Given values of ${k}_{b},p,{R}_{\mathrm{cell}},{C}_{0}$, and the filament height along $\partial U$, we numerically solved Equation (S6) by individually solving the decoupled equations, Equation (S9), with a twodimensional finite element Poisson equation solver (Figure 2—figure supplement 1b–c) (Burkardt, 2011). The value of $\mathrm{\Delta}E$ was then calculated numerically from the height function, $h$, by triangulating $h$ and extracting the mean curvature and enclosed volume of the resulting mesh using preexisting MATLAB (Mathworks, Natick, MA) software (Mecklai, 2004; Kroon, 2014; Suresh, 2010). To find the energyminimizing filamentmembrane conformation, $\mathrm{\Delta}E$ was numerically computed while varying the deformed filament curvature, $1/R$, the size of $\mathrm{\Omega}$, and the size of, and the filament height in, $U$ (the size of $U$ corresponds to the bound length of the filament; consistent with the discussion below, we find that ${\mathrm{\ell}}_{b}={L}_{f}$ in all cases of interest). As the simulation details above assume a perpendicular binding orientation relative to the long axis of a cylinder, for simplicity we model binding to deviatory angles by perpendicular binding to a different cell radius, $R$, where $\theta ={\mathrm{cos}}^{1}(R/{R}_{\mathrm{cell}})$. These implementation details were used to generate Figure 2—figure supplement 1b–c and Figure 4—figure supplement 4.
1.3. Preferred orientation of filament binding and binding phase diagram
In this section, following Hussain et al. (2018), we provide details of analytic calculations complementing the numerical calculations discussed above. For a cylindrical membrane whose radius is larger than the radius of curvature of a filament, it is energetically favorable for the filament to bind at an angle of $\theta ={90}^{\circ}$ relative to the long axis of the cell, as this orientation requires minimal bending of both the protein and the membrane. For deviatory angles, $\theta {90}^{\circ}\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}{0}^{\circ}$, an effective correction to the cell radius ${R}_{\mathrm{cell}}$ is a multiplicative factor of $1/\mathrm{cos}\theta $, which makes binding less energetically favorable. The energetic penalty for MreB filaments binding at deviatory angles is approximately $40kT$ for a broad range of membrane pressures and depending on the membrane radius (below and Figure 2B of the main text). In general, higher membrane pressures make it more energetically favorable for the protein filament to bend, which minimizes the amount of volume displaced by the proteinmembrane interaction (see below) (Hussain et al., 2018), while smaller membrane pressures make it more energetically favorable for the membrane to bend. Since a filament can always bend to conform to the membrane curvature, we see that large pressure differences across the membrane may enhance the energetic preference of an MreB filament for the perpendicular orientation.
We may summarize our results over a wide range of parameter values with an approximate phase diagram, assuming a cylindrical membrane (Figure 2C in the main text). We use the membrane pressure $p$ and the filament Young modulus $Y$, which varies with the bending energy ${E}_{\mathrm{bend}}$ of the filament, as order parameters. By considering only the volume and surface height displaced directly underneath the binding region of the membrane, the deformation energy of the binding region, $U$, is (Hussain et al., 2018)
where $b$ is the fraction of the filament crosssection needed to adhere to the membrane. The critical pressure at which the deformation is dominated by filament bending can be estimated as the value that sets ${E}_{\mathrm{def},U}$ to be minimal at $R={R}_{\mathrm{cell}}$:
which, for the parameter values summarized in Supplementary file 1, estimates $p}^{\ast}\approx 20\text{kPa$ (Hussain et al., 2018). This value of ${p}^{*}$ is smaller than estimates of the turgor pressures of both B. subtilis and E. coli (Supplementary file 1), suggesting that in vivo, MreB filaments always bend to adhere to the inner membrane. Equation (S10) is used under this assumption to generate the curves in Figure 2B of the main text. Furthermore, the linear dependence on ${\mathrm{\ell}}_{b}$ of the first term of Equation (S10) implies that MreB filaments bind fully along their lengths. If $p\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}{p}^{\ast}$, as is the case for vesicles, then both the membrane and the MreB filament can deform each other in a manner that minimizes the total energy, with the membrane shape determined by Equation (S6). For a large range of membrane pressures $0\phantom{\rule{thinmathspace}{0ex}}\le \phantom{\rule{thinmathspace}{0ex}}p\phantom{\rule{thinmathspace}{0ex}}\lesssim \phantom{\rule{thinmathspace}{0ex}}{p}^{\ast}$, we find that bound MreB filaments induce membrane curvature, and for vesicles where the pressure difference across the membrane is vanishingly small, the shape equation also predicts that MreB filaments can grossly deform the membrane (Figure 2—figure supplement 1c), a prediction consistent with experimental observations (Salje et al., 2011; van den Ent et al., 2014). For a filament with fixed dimensions, both ${E}_{\mathrm{bend}}$ and the expression for ${p}^{*}$ in Equation (S11) are proportional to $Y$. Hence, ${E}_{\mathrm{bend}}\propto {p}^{*}$ as $Y$ increases, and this relation is shown as the diagonal line in Figure 2C of the main text.
Similarly, considering only the membrane curvature induced by the binding region $U$, the filament does not bend if $E}_{\mathrm{b}\mathrm{e}\mathrm{n}\mathrm{d}}+E\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}{E}_{\mathrm{i}\mathrm{n}\mathrm{t}$, where $E=\pi b{k}_{b}{L}_{f}/{r}_{f}$. If this inequality were satisfied, then the interaction energy ${E}_{\mathrm{int}}$ may be too small to justify membrane binding, which requires a combination of polymer and membrane bending. Note that, while we assume $b=1/6$ in this work, our results do not significantly change for different $b$, as shown in Figure 4 of (Hussain et al., 2018), and the regimes delineated in this section are summarized in Figure 2C of the main text.
Finally, we note that, for the parameter values summarized in Supplementary file 1, the prediction that it is energetically favorable for MreB filaments to align along the circumferential direction of a rodlike cell is robust in the case where the intrinsic curvature varies with position along the filament. In particular, a calculation based on Equation (S2) shows that this conclusion follows given that the filaments are, on average, more curved than the membrane. To see this, let ${\kappa}_{s}(\mathrm{\ell})$ denote the intrinsic filament curvature as a function of position along the filament length, $\mathrm{\ell}$, and $\kappa $ denote the deformed curvature. In live cells, $\kappa $ does not vary as a function of $\mathrm{\ell}$ because it is most energetically favorable for the filament to bend completely to match the ambient membrane curvature, as we have shown above. The total bending energy of the filament is then
When binding to an angle that deviates from the circumferential direction in a cylinder, the deformed curvature will be smaller: let $\kappa ={\kappa}_{0}{\kappa}^{\prime}$, where ${\kappa}_{0}$ is the curvature along the circumferential direction of a cylinder and ${\kappa}^{\prime}\ge 0$ is a constant correction to ${\kappa}_{0}$ depending on the binding angle. Then, the difference in bending energies between binding in the direction of $\kappa $ as opposed the circumferential ${\kappa}_{0}$ direction is
which is larger than zero provided the filament is, on average, more curved than the cell (that is, the second term above is nonnegative). Hence, because the binding orientation is robust, our model predictions will remain the same. Nevertheless, we note that cryoEM experiments (Hussain et al., 2018; Salje et al., 2011; van den Ent et al., 2014) support the assumption of a uniformly bent MreB filament, and we have therefore focused on this case in the main text.
In the foregoing argument, we have assumed that $\kappa $ does not vary as a function of $\mathrm{\ell}$. The case in which $\kappa $ varies significantly as a function of $\mathrm{\ell}$, which is relevant to the geometry considered in Figure 3F of the main text, is more delicate. Note, however, that in this geometry $\kappa (\mathrm{\ell})$ can be explicitly calculated to show that MreB still orients in the largest principal direction. We consider, in particular,
where, at an axial coordinate $v$, ${E}_{1}$ is the bending energy of a (constantly curved) filament aligning in the axial direction (with principal curvature ${\kappa}_{1}$) and ${E}_{2}$ is the bending energy of a filament aligning in the circumferential direction (with principal curvature ${\kappa}_{2}$). Here the integration is over the axial coordinates corresponding to filament length, $[vd,v+d]$. Figure 3—figure supplement 2c–d shows that, in this case, binding along the greatest principal direction generally still incurs the least bending energy.
2. Curvaturebased translocation of filaments
Consider a surface parameterized by $\mathbf{\mathbf{r}}=r(u,v)\subset {\mathbb{R}}^{3}$, with $({x}_{1},{x}_{2})=(u,v)\in {\mathbb{R}}^{2}$, which, for simplicity, we assume to be smooth almost everywhere so that the following quantities are well defined. The considerations below can be readily extended to the case of a piecewise smooth surface, as discussed in the following section (§2.1). We model a curved filament as a point $\mathbf{\mathbf{p}}$ on this surface which translocates in a direction $\mathbf{\mathbf{d}}\in {\mathbb{R}}^{3}$. At $\mathbf{\mathbf{p}}$, the principal direction, $\mathbf{\mathbf{w}}$, represented in the basis of the tangent space satisfies
where $S$ is the shape operator at $\mathbf{\mathbf{p}}$ and $c$ is one of the principal curvatures of $\mathbf{\mathbf{r}}$, hereafter taken to be the largest. Throughout this work, we assume the sign convention that the curvature is positive when any normal vector at $\mathbf{\mathbf{p}}$ points towards the interior of $\mathbf{\mathbf{r}}$, so that, for instance, the largest principal curvature is always positive for a cylinder (§2.2). When the filament translocates in the direction of the largest principal curvature, $\mathbf{\mathbf{d}}$ satisfies $\mathbf{\mathbf{d}}=\mathbf{\mathbf{w}}\cdot ({\mathbf{\mathbf{r}}}_{u},{\mathbf{\mathbf{r}}}_{v})$, where ${\mathbf{\mathbf{r}}}_{u}={\partial}_{u}\mathbf{\mathbf{r}}$ and ${\mathbf{\mathbf{r}}}_{v}={\partial}_{v}\mathbf{\mathbf{r}}$.
For convenience, we recall the following statements from the main text. We set $\eta ={\mathrm{cos}}^{1}\frac{\mathbf{d}\cdot {\mathbf{r}}_{\theta}}{\mathbf{d}\cdot {\mathbf{r}}_{\theta}}$, where $\eta$ is an angular deviation from the largest principal direction on the surface introduced by possible sources of stochasticity, the modified direction corresponds to an angle $\theta $ relative to the $u$axis in parametric coordinates, $\mathbf{r}}_{\theta}\in {\mathbb{R}}^{3$ is the derivative of $\mathbf{r}$ in the direction of $\theta $, and distances are defined by the surface metric. Note that, when the largest principal direction is well defined, it is independent of the parameterization up to a sign (Lee, 2009). After the filament orients, it may translocate in the direction of $\mathbf{\mathbf{d}}$ for a certain distance, dependent on the filament speed, before reorienting its direction of motion. As $\mathbf{\mathbf{d}}$ is a vector in ${\mathbb{R}}^{3}$, arbitrarily moving in the direction of $\mathbf{\mathbf{d}}$ may, however, move the filament out of $\mathbf{\mathbf{r}}$. To define the translocation consistently, we define it $(u,v)$space by taking the translocation angle with respect to the $u$axis to also be $\theta $. Translocating along an angle $\theta $ with respect to the $u$axis in $(u,v)$space ensures that the filament remains on the surface, and the direction of translocation corresponds to that on the smooth patch defined locally on $\mathbf{\mathbf{r}}$.
As a discretetime flow in parametric space, and with suitable units of time so that the filament may reorient at every timestep, the 2D equation of filament motion is
where, in $(u,v)$space, ${X}_{n}$, ${\mathrm{\ell}}_{n}$, and $\theta}_{n$ are the position, step size, and translocation angle, respectively, of the filament at a timestep $n$. Here ${\theta}_{n}$ is the value of $\theta $ computed at the surface point corresponding to ${X}_{n}$ and assuming $\eta \sim \mathcal{\mathcal{N}}(0,{\sigma}^{2})$—that is, the angular noise is normallydistributed, with mean zero and variance ${\sigma}^{2}$. The strength of the noise may depend on factors such as the energetic difference of binding along the two principal directions, and we examine cases where $\sigma $ depends on the difference, $\mathrm{\Delta}c$, between principal curvatures below. ${\chi}_{n}$ is a random sign, which accounts for the possibility of both lefthanded and righthanded translocation, and may not substantially vary in $n$ if filament motion is processive (namely, the filament does not backtrack). We assume that ${\mathrm{\ell}}_{n}$ satisfies the integral equation
where $X(\tau )={X}_{n}+\tau (\mathrm{cos}{\theta}_{n},\mathrm{sin}{\theta}_{n})$ and $g$ denotes the metric tensor (with $u$ and $v$ corresponding to the indices $1$ and $2$, respectively), which relates it to a constant filament step size, $L$, on the surface. When $g$ is slowly varying, that is $g(X(\tau ))\approx g({X}_{n})$, as is expected in the limit of small ${\mathrm{\ell}}_{n}$, or in particular
where $\hat{Y}=(\mathrm{cos}{\theta}_{n},\mathrm{sin}{\theta}_{n})$, $I$ is the integrand of Equation (S17), and we discard higherorder terms in $\nabla I$ and ${\mathrm{\ell}}_{n}$, Equation (S17) simplifies to a linear equation for ${\mathrm{\ell}}_{n}$:
In general, we wish to determine the distribution of $X$ in $(u,v)$space, which would determine the distribution of the filament on the surface. Below, we introduce activation and deactivation dynamics, as discussed in the main text, and examine the statistics of $X$ on several surfaces.
2.1. Numerical solutions in the finite and continuum cases
Before discussing the implications of Equation (S16), we provide implementation details for the simulations and numerical calculations discussed below and in the main text. We numerically implemented the dynamics of Equation (S16) with Langevin simulations and verified our results with coarsegrained, continuum calculations corresponding to solutions of the corresponding FokkerPlanck equation. For the former, individual filaments were activated and deactivated as discussed in the main text, and their trajectories were simulated directly according to Equation (S16). The final positions of all filaments were pooled together to determine the corresponding filament concentrations. For the latter, we discretized $(u,v)$space uniformly into ${M}_{u}\times {M}_{v}$ rectangular elements, computed the transition matrix corresponding to Equation (S16) and a given set of parameter values, and multiplied a vector corresponding to the filament number, ${N}_{F}$, to model the dynamics. Deactivation effects were modeled by multiplying the vector of filament numbers after each timestep by ${e}^{\lambda}$, and filament concentrations were determined by dividing the vector of filament numbers by the corresponding surface area elements. Specific implementation details, such as geometric parameters, are summarized for each figure in this work in §3.
For a piecewise parametric surface, we allowed filaments to translocate between two subsurfaces while conserving the total step size, $L$, translocated at every timestep. For the piecewise geometries considered in this work, the filament coordinates were connected between subsurfaces in an obvious manner.
2.2. Cylinder
Consider a cylinder parameterized by $\mathbf{\mathbf{r}}=(a\mathrm{cos}u,a\mathrm{sin}u,av)$, where $a\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$, $u\in [0,2\pi ]$, $v\in [0,z]$, and $z$ is a large enough constant so that we do not consider filament translocation out of the cylinder. The nonvanishing components of the metric tensor are ${g}_{11}={g}_{22}={a}^{2}$. Hence, the step sizes in parameter space are identical at all timesteps $n$ and equal to ${\mathrm{\ell}}_{n}=\mathrm{\ell}=L/a$, and translocation occurs mainly along the $u$axis. Assuming processive filament motion, we recover a Pearsonlike random walk in a 2D plane with periodic boundary conditions in $u$, where the translocation angle satisfies ${\theta}_{n}\sim N(0,{\sigma}^{2})$, and ${\chi}_{n}=\pm 1$ for all $n$. Correlated Pearson random walks of a similar form have been studied by Kareiva and Shigesada in the context of insect movement (Kareiva and Shigesada, 1983). Denoting by ${U}_{N}$ and ${V}_{N}$ the displacement along the $u$ and $v$ coordinate, respectively, after $N$ steps, with $U}_{N}=\sum _{i=1}^{N}\ell \mathrm{cos}{\theta}_{i$ and $V}_{N}=\sum _{i=1}^{N}\ell \mathrm{sin}{\theta}_{i$, we find that
in agreement with simulations of Equation (S16) for parameter values relevant to B. subtilis (§3.1 and Figure 3B of the main text).
2.3. Torus and Helix
Let $\mathbf{\mathbf{r}}=((R+a\mathrm{sin}u)\mathrm{cos}(av/R),(R+a\mathrm{sin}u)\mathrm{sin}(av/R),a\mathrm{cos}u)$, where $u\in [0,2\pi ]$, $v\in [0,R\mathrm{\Phi}/a]$, $\mathrm{\Phi}$ is the axial subtended angle, and $0\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}a\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}R$. Here the nonvanishing elements of the metric tensor are ${g}_{11}={a}^{2}$ and ${g}_{22}={a}^{2}{(1+(a/R)\mathrm{sin}u)}^{2}$, and the step sizes depend on $u$. As the principal curvatures in the circumferential $u$direction, ${c}_{u}$, and in the axial $v$direction, ${c}_{v}$, satisfy ${c}_{u}=1/a\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}{c}_{v}=(R+a\mathrm{sin}u{)}^{1}\mathrm{sin}u$, translocation occurs mainly in the direction of $u$. Due to the dependence of the step sizes on $u$, $\u27e8{U}_{N}\u27e9$, $\u27e8{V}_{N}\u27e9$, $\mathrm{Var}({U}_{N})$, and $\mathrm{Var}({V}_{N})$ may differ from that of a cylinder and are analytically difficult to calculate. We focus instead on determining the filament concentration, $C}_{F$, which, unlike the case of a cylinder, will be nonuniform in $X$.
For intuition, we consider the case of no translocation noise. Consider the probability ${P}_{N}(X;Y)$ of observing the filament at any point $X$ in $(u,v)$space at a timestep $N$ given that it is initially at $Y$, and assume ${p}_{\mathrm{init}}(Y)$ to be proportional to $dA={a}^{2}(1+(a/R)\mathrm{sin}u)dudv$. Since $dA$ depends on $u$, the probability ${P}_{N}(X)$ of observing the filament at a timestep $N$, averaged over initial positions, is not uniform. Nevertheless, considering the dynamics of activating and deactivating filaments as above, we find that, in the limit of continuous time, the expected number of filaments at a coordinate $X=(u,v)$ and time $t$ is
where ${k}^{\prime}=k{p}_{\mathrm{init}}/dA$ and $\nu $ is a speed corresponding to ${U}_{N}$ in the continuum limit. As $\lambda \to \mathrm{\infty}$, corresponding to the limit in which filaments are instantaneously deactivated, Equation (S21) predicts a number enhancement on the outer edge of the form ${N}_{F}(X,t)=(2{a}^{2}{k}^{\prime}/\lambda )(1+(a/R)\mathrm{sin}u)$, while as $\lambda \to 0$, corresponding to the limit in which filaments persist indefinitely, Equation (S21) predicts an approximate uniform distribution $2{a}^{2}{k}^{\prime}(t+\frac{a}{R\nu}\mathrm{sin}(t\nu )\mathrm{sin}u)\approx 2{a}^{2}{k}^{\prime}t$ in the limit of large $t$. While the concentration of filaments, ${C}_{F}$, is uniform over the surface in the former case, the latter case corresponds to the ‘washingout’ of initial activation conditions and implies a larger value of ${C}_{F}$ on the inner edge. Namely, in the formal limit $\lambda \to 0$ followed by $t\to \mathrm{\infty}$, with $\lambda t\to 0$,
which is maximized at $u=\pi /2$, corresponding to the inner edge of the torus. The filament concentration at the inner edge becomes larger than that of the outer edge by a factor of $(1+a/R)/(1a/R)$, which depends only on the geometry of the torus. This effect may be interpreted as a ‘geometric focusing’ caused by both the filament number, ${N}_{F}$, becoming uniform in $(u,v)$coordinates and variations in the area element $dA$ in $u$ or $v$. Even for a range of noises, including the limit of large noise ($\sigma \to \mathrm{\infty}$), simulations show that this description is valid and agrees with Equation (S22) (Figure 4—figure supplement 1a).
While the geometric focusing effect is clearly applicable to other geometries, particularly those involving bent or sinusoidal surfaces (Takeuchi et al., 2005; Renner et al., 2013), we examine whether it applies to ellipsoids and other geometries considered in this work below. We also note that this effect is different from filaments translocating along a direction of curvature until it is ‘attracted,’ as is the case for a spherocylinder (§2.5) or a surface with small wavelength undulations (Figure 3F of the main text and §2.8). The main difference is that the dynamics corresponding to geometric focusing is recurrent, and hence, due to spatial variations in the area element $dA$, leads to ${C}_{F}$ depending on geometric features such as in Equation (S22). In contrast, the dynamics corresponding to an attractor is not. For instance, in the absence of noise and the limit of large processivity, ${C}_{F}$ is nonzero over a toroidal surface and the relative localization is determined by Equation (S22), but ${C}_{F}$ vanishes at the hemispherical caps on a spherocylinder, as discussed later (§2.5).
Note that characteristic values of MreB persistence, as summarized in Supplementary file 2, suggest that MreB filaments lie in the $\lambda \approx 0$ regime. The observation of enhanced MreB concentration at the inner edge of toroidal cells is quantitatively consistent with the geometric focusing effect discussed in this section, as explained in the main text and prior work (Wong et al., 2017).
Finally, we note that similar arguments as that above describe the case of a helix, which can be parameterized by
where $u\in [0,2\pi ]$, $v\in [0,R\mathrm{\Phi}/a]$, $\mathrm{\Phi}$ determines the extended length of the helix, $0\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}a\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}R$, and $\varphi $ is the helical pitch. Here the nonvanishing elements of the metric tensor are ${g}_{11}={a}^{2}$, ${g}_{12}={g}_{21}=\varphi {a}^{3}{(R\sqrt{{R}^{2}+{\varphi}^{2}})}^{1}$, and ${g}_{22}={a}^{2}(2({R}^{2}+{\varphi}^{2}{)}^{2}+{a}^{2}({R}^{2}+2{\varphi}^{2})$ $+aR(4({R}^{2}+{\varphi}^{2})\mathrm{sin}uaR\mathrm{cos}2u))(2{R}^{2}({R}^{2}+{\varphi}^{2}){)}^{1}$, and the step sizes again depend on $u$. As the principal curvatures in the circumferential direction, ${c}_{u}$, and in the axial direction, ${c}_{v}$, satisfy ${c}_{u}=1/a\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}{c}_{v}=({R}^{2}+{\varphi}^{2}+aR\mathrm{sin}u{)}^{1}R\mathrm{sin}u$, translocation occurs mainly in the direction of $u$. Like the case of a torus, we expect filaments to become uniformly distributed across circumferential hoops in the limit of infinite processivity. The predicted value of ${C}_{F}$ in this case, ${C}_{F}\propto 1/dA={({a}^{2}({R}^{2}+{\varphi}^{2}+aR\mathrm{sin}u)dudv)}^{1}R\sqrt{{R}^{2}+{\varphi}^{2}}$, is quantitatively consistent with numerical simulations and larger at the inner edge (Figure 4—figure supplement 1b–c).
2.4. Ellipsoid
Consider the surface $\mathbf{\mathbf{r}}=(a\mathrm{sin}u\mathrm{cos}v,a\mathrm{sin}u\mathrm{sin}v,b\mathrm{cos}u)$, where $a,b\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$, $u\in [0,\pi ]$, and $v\in [0,2\pi ]$. Here the nonvanishing components of the metric tensor are ${g}_{11}=({a}^{2}+{b}^{2}+({a}^{2}{b}^{2})\mathrm{cos}(2u))/2$ and ${g}_{22}={a}^{2}{\mathrm{sin}}^{2}u$, and $dA=\frac{1}{\sqrt{2}}a\mathrm{sin}u\sqrt{{a}^{2}+{b}^{2}+({a}^{2}{b}^{2})\mathrm{cos}(2u)}dudv$. When $b\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}a$, the principal curvatures in the circumferential $v$direction, ${c}_{v}$, and in the axial $u$direction, ${c}_{u}$, satisfy
with equality only at the poles ($u=0,\pi $). Hence, in the case of no noise, filaments always translocate in the direction of $v$, and similarly, the opposite is true when $b<a$.
As calculations involving a finite noise are analytically complex, we first consider the case of no noise. Here ${U}_{N}=0$ and ${V}_{N}^{u}=LN/(a\mathrm{sin}u)$, where ${U}_{N}$ is the displacement along the $u$coordinate after $N$ steps and ${V}_{N}^{u}$ is the displacement along the $v$coordinate after $N$ steps assuming that the filament remains on the curve of constant $u$. For an ensemble of filaments, considering similar activation and deactivation dynamics as above gives
where ${k}^{\prime}$ is the same as above and $\nu $ is a constant speed corresponding to ${V}_{N}={V}_{N}^{u}\mathrm{sin}u$ in the continuum limit. Hence, for any value of $\lambda $, ${C}_{F}$ is a constant, in agreement with numerics (Figure 4—figure supplement 2).
We next consider the case of a finite noise. We first note that the area element $dA$ vanishes at the tips, and here the resulting value of ${C}_{F}$ would diverge if ${N}_{F}$ were nonzero. To compare filament concentrations over regions with differently sized area elements, it is convenient to define the average concentration on the (sub)surface $S$ as
which is different from a direct averaging of ${C}_{F}$ over area elements as considered for a torus (the latter has the form of a harmonic sum). Numerical calculations accounting for an ensemble of filaments with similar activation and deactivation dynamics as above show that, intriguingly, ${\u27e8{C}_{F}\u27e9}_{S}$ is enhanced at the ellipsoidal poles ($S=\{X:u\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}\pi /4\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}\mathrm{o}\mathrm{r}\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}u\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}3\pi /4\}$) for a range of noises (Figure 4—figure supplement 2). Different from the case of zero noise, here the increased values of ${C}_{F}$ at the poles arise because filaments may randomly translocate to the poles and the area element is significantly smaller there than that away from the poles. Similar observations also hold in the limit of large noise, $\sigma \to \mathrm{\infty}$, and the limiting value of the ratio of ${\u27e8{C}_{F}\u27e9}_{S}$ at the tips versus the bulk is the ratio of the corresponding areas (Figure 4—figure supplement 2).
In the case of a finite noise, localization in a torus, a helix, and an ellipsoid arises because of ${N}_{F}$ becoming spatially homogeneous and spatial variations of the area elements. However, in the case of zero noise, we remark on a key difference between ${C}_{F}$ on a torus and an ellipsoid as follows. For a torus, spatially heterogeneous filament concentration arises due to the trajectories of individual filaments being closed orbits in the circumferential direction and the long persistence of filaments, which ‘washes out’ the initial distribution of filament position. For an ellipsoid, filaments do not move in the axial direction and the resultant filament concentration remains uniform on the surface, regardless of persistence. That a uniform distribution is maintained in the absence of noise on an ellipsoid is also different from the cases of a spherocylinder (§2.5) and a surface with small wavelength undulations (Figure 3F of the main text and §2.8), as discussed later. In general, ${C}_{F}$ would be uniform in the case of zero noise for any surface where the direction in which $dA$ varies does not coincide with the direction of filament translocation.
2.5. Spherocylinder
Consider the surface parameterized piecewise by ${\mathbf{\mathbf{r}}}_{1}=(a\mathrm{cos}u,a\mathrm{sin}u,av)$, where $a\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$, $u\in [0,2\pi ]$, $v\in [0,z]$ for some $z\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$, ${\mathbf{\mathbf{r}}}_{2}=(a\mathrm{sin}v\mathrm{cos}u,a\mathrm{sin}v\mathrm{sin}u,a\mathrm{cos}v+az)$, where $u\in [0,2\pi ]$ and $v\in [0,\pi /2]$, and ${\mathbf{\mathbf{r}}}_{3}={\mathbf{\mathbf{r}}}_{2}(0,0,az)$, where $u\in [0,2\pi ]$ and $v\in [\pi /2,\pi ]$. In the absence of noise, filaments translocate along the $u$ axis in the cylinder and a random angle in the hemispheres. This could allow for translocation out of the caps and into the cylindrical body. Indeed, in the absence of noise, it is clear that any filament in the hemispherical caps will eventually translocate into and remain at the cylindrical rims.
For the spherocylindrical surface described above, general calculations in the case of a finite noise are analytically complex due to the irregular geometry. In the presence of noise, however, the depletion of filaments at the hemispherical caps can be supported numerically. For parameter values relevant to MreB in E. coli and B. subtilis (Supplementary file 2), numerical calculations show that ${\u27e8{C}_{F}\u27e9}_{S}$ is larger by a multiplicative factor of ∼2.0 when $S$ is the cylindrical bulk than that when $S$ is a hemispherical endcap (Figure 3C of the main text).
2.6. Filament concentration is independent of Gaussian curvature
While filament movement depends on the principal curvatures, we wondered whether this implies that filament concentration always depend on the Gaussian curvature. Here we show that, for general surfaces, filament concentration is independent of Gaussian curvature under the dynamics considered in this work. Consider the parameterization $\mathbf{\mathbf{r}}=(\mathrm{sin}u,(1\mathrm{cos}u)\mathrm{cos}u,v)$, where $u\in [\pi ,\pi ]$, $v\in [0,z]$, and $z$ is a large enough constant so that we do not consider filament translocation out of the surface. Although the Gaussian curvature vanishes everywhere, any crosssection of the surface has regions of both negative and positive principal curvature, and the regions of positive principal curvature represent attractors to the filament dynamics: filaments translocating into such regions rarely translocate out.
While numerical simulations show that filaments are localized at regions where the largest principle direction coincides with the axial direction (Figure 4C in the main text and Figure 4—figure supplement 3a–b), analytical calculations are difficult to undertake due to the complicated geometry. Nevertheless, we may consider a similar 1D problem of a particle moving with velocity $\nu $ along a circular ring, which is parameterized without loss of generality by $u\in [0,2\pi ]$, and contains a single absorbing coordinate at $u=0$. Accounting for activation and deactivation dynamics as above, the filament number at $u=0$ at a time $t\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}{t}^{\ast}$ in the case of vanishing translocation noise ($\sigma =0$) and assuming $\nu \phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}2\pi \lambda$, so that the flux described below is nonzero, can be written as
where ${k}^{\prime}$ is defined above, the term in parenthesis is the filament flux into $u=0$ at time ${t}_{1}$, and ${t}^{*}=2\pi /\nu $ represents the maximal time needed for any filament to become absorbed. Further assuming a constant value of $dA(u,v)=dA$ for simplicity, direct evaluation of Equation (S29) yields
which quantitatively describes the dependence of the filament number at the absorbing coordinate as a function of the processivity, which is determined by $\lambda $, and other parameters. Note that, as $\lambda \to 0$, corresponding to the case of infinite processivity, ${N}_{F}(u=0)$ is predicted to diverge, while as $\lambda \to \mathrm{\infty}$, corresponding to the case of zero processivity, ${N}_{F}(u=0)\to 0$. It is straightforward to generalize Equation (S30) to the case of several absorbing points and different geometries, provided that similar simplifying assumptions can be employed as above. Importantly, these calculations can be compared to numerical calculations for the geometry considered in this section. While the value of ${N}_{F}$ at attracting regions may generally depend on the geometry of such regions, the dependence of ${N}_{F}$ on the processivity, $\lambda $, can be conveniently explored by defining the localization ratio, $\rho $, as the ratio between filament numbers at a certain value of $\lambda $ compared to that at $\lambda /2$:
Defining $\rho $ analogously for the geometry considered in this section, we find that numerical calculations of $\rho $ for the geometry considered here across a range of $\sigma $ are consistent with the theoretical prediction for $\rho $ based on the simplified model considered in this paragraph (Figure 4—figure supplement 3c). Hence, we conclude that Equation (S30) captures the dependence of ${N}_{F}$ on $\lambda $ for more general geometries.
2.7. Cylinder with a bulge
Building on work that has examined MreB motion in cells with membrane bulges (Hussain et al., 2018), here and below we consider bulges of positive Gaussian curvature on a cylindrical surface. We consider a class of cylinders with bulges parameterized piecewise by ${\mathbf{\mathbf{r}}}_{1}=(v,\mathrm{cos}u,\mathrm{sin}u)$, where $u\in [0,2\pi ]$ and $v\in [0,z]$ for some $z\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$ and $\sqrt{{(u\pi /2)}^{2}+{(vz/2)}^{2}}<b$ and ${\mathbf{\mathbf{r}}}_{2}=(b\mathrm{sin}v\mathrm{cos}u+z/2,b\mathrm{sin}v\mathrm{sin}u,sc\mathrm{cos}v+1)$ where $u\in [0,2\pi ]$ and $v\in [0,\pi /2]$ otherwise. We suppose $b,c\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$ and $b\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}z$, so that the intersection area constitutes a small fraction of the cylindrical body. Here $s=1$ if the bulge protrudes outward and $s=1$ if the bulge protrudes inward. For this class of parametric surfaces, three general cases of outward bulges ($s=1$) can be classified depending on whether $c=b$, $c\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}b$, or $c\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}b$. Upon computing the principal curvatures in each case, we find that the three cases correspond respectively to random translocation ($c=b$), in which case the translocation direction is random; polar translocation ($c\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}b$), in which case translocation occurs predominantly along the bulge $v$coordinate; and circumferential translocation ($c\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}b$), in which case translocation occurs predominantly along the bulge $u$coordinate (Figure 3—figure supplement 1a). We note that inward bulges of positive Gaussian curvature ($s=1$) exhibit the same behavior, with polar translocation when $c\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}b$ and circumferential translocation when $c\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}b$.
Figure 4B of the main text illustrates filament dynamics for both random and circumferential translocation where $b=c=0.5$ (spherical bulge) and $b=0.5,c=1$ (ellipsoidal bulge), in both the cases of zero and infinite processivity. We focus here on the latter geometry as the motion there is consistent with experiments (Hussain et al., 2018). When there is no noise, considering similar activation and deactivation dynamics as above shows that, similar to a spherocylinder (§2.5), the bulge attracts filaments in the case of small $\lambda $ or large processivity. As simulations demonstrate, this is also true in the case of a finite noise, where we find the bulge to contain larger numbers of filaments relative to the case of a uniform distribution on the surface (Figure 4B in the main text). For characteristic parameter values relevant to B. subtilis MreB, we find that ${\u27e8{C}_{F}\u27e9}_{S}$ is increased at the bulge (Figure 3E of the main text) and substantial localization may occur at the bulge neck, consistent with previous experiments (Figure 3—figure supplement 1b).
2.8. Cylinder with small wavelength undulations
Previous work has examined the correlation of MreB concentration with subcellularscale shape fluctuations in E. coli cells (Ursell et al., 2014). Based on the weak correlation observed between outer and inner contour curvatures in Ursell et al. (2014), the authors determined that the experimentally observed MreB enrichment was not caused by bending modes: in a circular torus, for instance, the correlation between outer and inner contour curvatures should be strictly negatively correlated. The authors concluded that short lengthscale, high magnitude fluctuations dominate experimental observations of cell shape.
To probe such a geometry, we consider the parameterization $\mathbf{\mathbf{r}}=((a+c\mathrm{sin}(Pv))\mathrm{sin}u,(a+c\mathrm{sin}(Pv))\mathrm{cos}u,v)$, for $u\in [0,2\pi ]$ and $v\in [0,2\pi ]$. Here $a$ denotes an average cylinder radius, we assume that $0\phantom{\rule{thinmathspace}{0ex}}<\phantom{\rule{thinmathspace}{0ex}}c\phantom{\rule{thinmathspace}{0ex}}\ll \phantom{\rule{thinmathspace}{0ex}}a$, and $P\phantom{\rule{thinmathspace}{0ex}}>\phantom{\rule{thinmathspace}{0ex}}0$ is a variable controlling the number of periods. For large wavelength undulations, $P$ is less than, or on the order of, unity. However, the Gaussian and mean curvatures are uncorrelated in this case, inconsistent with the positive correlation observed in Ursell et al., 2014 for cells growing in sinusoidallyshaped chambers, and the ranges of Gaussian and mean curvatures are significantly smaller than those measured for wildtype, unconfined, sinusoidallyconfined, thin, and wide cells in different experiments (Ursell et al., 2014; Shi et al., 2017; Bratton et al., 2018) (Figure 3—figure supplement 2a and Figure 3G of the main text). We therefore consider a geometry with short wavelength undulations, for which $P\phantom{\rule{thinmathspace}{0ex}}\ge \phantom{\rule{thinmathspace}{0ex}}1$ (Figure 3—figure supplement 2b). Numerical results for principal curvaturedependent translocation on this geometry, which is consistent with the predicted filament binding orientation (Figure 3—figure supplement 2c–d), are presented in Figure 3G in the main text.
2.9. Effects of principal curvaturedependent translocation noise and varying filament step size on model predictions
Prior experiments have shown that the variation in MreB trajectory direction is widthdependent, suggesting that the translocation noise may depend on the difference of principal curvatures, $\mathrm{\Delta}c$ (Hussain et al., 2018). As discussed in the main text, we may model this observation by letting $\sigma $ vary with $\mathrm{\Delta}c$: for simplicity, we set
(Figure 3—figure supplement 3a), but note that more complicated functional dependencies, such as a quadratic dependence of the form $\sigma =\beta {(\mathrm{\Delta}c)}^{2}$, do not significantly change our results (Figure 3—figure supplement 3b). Unless otherwise indicated (see §3), all simulations and calculations in this work pertaining to MreB assume Equation (S32), and we further note that, for the parameter values relevant to MreB translocation in E. coli and B. subtilis summarized in Supplementary file 2, taking a constant value of $\sigma =0.3$ also does not significantly change our results. Similarly, we verify that localization arises even for vary large step sizes ($L=2\mu \text{m}$), in which case MreB filaments realign infrequently (Figure 3—figure supplement 3b).
2.10. Effects of filament twist, flexural rigidity, and Gaussian curvaturedependent activation on model predictions
Recent work has shown that regions of negative Gaussian curvature can allow twisted filaments to bind with low elastic energy (Quint et al., 2016). In another study, Wang and Wingreen assumed that MreB assembles into bundles with significantly larger flexural rigidity than that of filaments considered in this work (Wang and Wingreen, 2013). In this section, we examine how our model predictions for MreB localization in E. coli change over ranges of three parameters: (1) the intrinsic filament twist, ${\omega}_{0}$, (2) the coupling, $\gamma $, of filament activation to Gaussian curvature, and (3) the filament flexural rigidity, $B$.
We first note that $\gamma $ may vary independently of ${\omega}_{0}$. Quantitatively, the rate of filament activation may not only depend on the parameters of the twist, but also on other cellular parameters (Wong et al., 2017). We assume that the filament activation rate per unit area, $k$, varies with $\gamma $ as
where ${k}_{0}$ and $\gamma $ are constants and $G(u,v)$ is the Gaussian curvature at the parametric coordinate $(u,v)$. When $\gamma =0$, we recover our original assumption that the activation rate per unit area is constant.
We now explore the effects of twist and Gaussian curvaturedependent activation on filament concentration. For simplicity, we consider a torus (Figure 4—figure supplement 1a), for which $G=\mathrm{sin}u{[a(R+a\mathrm{sin}u)]}^{1}$ (c.f. §2.3), but show later in Figure 3—figure supplement 3b that the results for the undulating geometry of Figure 3F of the main text also remain qualitatively similar. Due to twist, MreB filaments may move in a direction which deviates from the direction of largest curvature. In particular, while further work should verify the robustness of this equation for nonzero turgor pressures and intrinsic filament curvatures, Equation (12) of Quint et al. (2016),
predicts the angular deviation from the largest principal direction in a cylinder, ${\theta}_{0}$, due to twist in the limit that MreB filaments are fully bound to the membrane. Here $B$ is the flexural rigidity of a filament, $K$ is the elastic twist stiffness, $a$ is the cell radius, and ${\omega}_{0}$ is the intrinsic helical twist. Thus, upon knowing the values of ${k}_{0}$, ${\omega}_{0}$, $\gamma $, and $B$, we may compute $k$ and ${\theta}_{0}$ using Equation (S33) and (S34), from which we can determine ${C}_{F}$ via simulations similar to those above. Doing so for a large parameter range which includes characteristic values of $B$ (Supplementary file 1) and theoretically hypothesized values for MreB of $K=2000kT\cdot \text{nm}$ and ${\omega}_{0}a=1$ to $5$ (Quint et al., 2016) reveals the final ratio of filament concentrations to be quantitatively similar in all cases (Figure 4—figure supplement 4a–f). We note, in particular, that the effect of filament twist alone is small in all cases: this is because the biasing of translocation angles due to twist, as predicted by Equation (S34), is irrelevant for the toroidal geometry, where translocating along deviatory angles still traces out hoops (c.f. §2.3). In these simulations, to illustrate the range of translocation behavior we have assumed that binding is always energetically favorable; however, there will be an energetic cost of unwinding of the form
where $\omega $ is the bound filament twist (Quint et al., 2016). For the parameter values considered in this work (Supplementary file 1), binding becomes energetically unfavorable at large twists (${\omega}_{0}a\phantom{\rule{thinmathspace}{0ex}}\gtrsim \phantom{\rule{thinmathspace}{0ex}}100$), for which ${\theta}_{0}\approx \pi /4$.
Additionally, our previous measurements of MreB fluorescence in bent E. coli cells confined to toroidal microchambers (which may be modeled as torii with geometric parameters $a=1$ and $R=10$) shows MreB concentration to be enhanced at the inner edges by a factor of approximately $1.1$ (Wong et al., 2017). This suggests the strength of Gaussian curvature coupling for MreB, if indeed such coupling does exist, to be small, as shown in Figure 4—figure supplement 4g–h. As we have demonstrated previously (Wong et al., 2017), the small observed enhancement is consistent with processivity alone. Furthermore, as mentioned above, modeling characteristic parameter values of both filament twist and Gaussian curvaturedependent activation still results in localization for geometries other than a torus, such as that considered in Figure 4F of the main text (Figure 3—figure supplement 3b).
Finally, we were interested to determine if our model predictions were robust to variation in the filament flexural rigidity, which for the parameter values in Supplementary file 1 has a value of $B\approx 1.65\times {10}^{25}\text{J}\cdot \text{m}$. To explore this for a torus, we repeated the foregoing dynamical simulations across a range of flexural rigidities. We first considered a tenfold smaller value of $B=1.0\times {10}^{26}\text{J}\cdot \text{m}$, which, by Equation (S10), still predicts (1) the filament to bend to conform to the membrane and (2) the depth of the potential well corresponding to Figure 2B of the main text to be approximately $3kT$, a value that may be large enough to be robust to thermal noise and other sources of stochasticity (Figure 4—figure supplement 4k). We next considered a 100fold larger value of $B=1.5\times {10}^{23}\text{J}\cdot \text{m}$, approximately the flexural rigidity of a thick bundle with $r}_{f}=10\phantom{\rule{thinmathspace}{0ex}}\mathrm{n}\mathrm{m$ (Wang and Wingreen, 2013). For this value of $B$ and ${r}_{f}$, no twist, and the remaining parameters summarized in Supplementary file 1, Equation (S11) estimates the critical pressure to be $p}^{\ast}\approx 5\phantom{\rule{thinmathspace}{0ex}}\mathrm{a}\mathrm{t}\mathrm{m$, a value which is larger than characteristic estimates of turgor in E. coli. As the corresponding value of ${E}_{\mathrm{bend}}$ is less than ${E}_{\mathrm{int}}$, we anticipate that both the membrane and the bundle may bend (see Figure 2C of the main text). To show that it remains energetically favorable for thick bundles to bind to membranes and that the binding orientation along the largest principal direction is robust, we numerically solved the shape equation (Equation (S6)) and found that the free energy change due to binding, $\mathrm{\Delta}E$, decreases with the size of the domain $\mathrm{\Omega}$ and, for a given $\mathrm{\Omega}$, is minimal when the membrane bends to conform to the bundle (Figure 4—figure supplement 4k). In the limit $\mathrm{\Omega}\to U$, the Monge gauge assumption underlying Equation (S6) becomes invalid; nevertheless, $\mathrm{\Delta}E$ tends to the analytical expression of Equation (S10) under the condition that the minimizer, $R$, is close to the intrinsic radius of curvature of the bundle, ${R}_{s}$. In both this limit and the simulations of Figure 4—figure supplement 4k, binding remains energetically favorable and the binding orientation remains robust even for thick filaments, suggesting that translocating along directions of largest principal curvature remains relevant. Figure 4—figure supplement 4f,i–j shows simulation results for both values of $B$ compared to the value ($1.65\times {10}^{25}\text{J}\cdot \text{m}$) assumed in this work. We find that, in all cases, the model predictions remain quantitatively similar.
3. Summary of Figures
For convenience, here we summarize implementation details used to generate figures in this work.
3.1. Figure 3B of the main text
Here $a=1$ and $z=100$. Langevin simulations to generate trajectories were undertaken with ${10}^{5}$ activated filaments, $L=0.4$, $\sigma =0.3$, and $N$ determined by the number of steps needed to translocate one hoop. Filaments were activated at the center ($v=0$) so that none of them translocated beyond the range specified by $z$.
3.2. Figure 3C of the main text
Here $a=1$ and $z=4$. Langevin simulations to generate a representative trajectory were undertaken with a single activated filament and the parameters relevant to E. coli summarized in Supplementary file 2, except $N=500$. Note that we use the linear relation $\sigma =\alpha {(\mathrm{\Delta}c)}^{1}$, where $\mathrm{\Delta}c$ is the difference of principal curvatures and the value of $\alpha $ is provided in Supplementary file 2. Numerical calculations for ensemble dynamics were undertaken with the parameters relevant to E. coli summarized in Supplementary file 2. Each subsurface is discretized into $60\times 60$ bins.
3.3. Figure 3D of the main text
Here $a=1$, $R=10$, and $\mathrm{\Phi}=\pi $. Periodic boundary conditions in $v$ are assumed. Note that we use the linear relation $\sigma =\alpha {(\mathrm{\Delta}c)}^{1}$, where $\mathrm{\Delta}c$ is the difference of principal curvatures and the value of $\alpha $ is provided in Supplementary file 2. Numerical calculations for ensemble dynamics were undertaken with the parameters relevant to E. coli summarized in Supplementary file 2. The surface is discretized into $30\times 30$ bins.
3.4. Figure 3E of the main text
Here $z=4$, $b=0.5$, and $c=1$. Periodic boundary conditions in $v$ for the cylinder are assumed. Note that we use the linear relation $\sigma =\alpha {(\mathrm{\Delta}c)}^{1}$, where $\mathrm{\Delta}c$ is the difference of principal curvatures and the value of $\alpha $ is provided in Supplementary file 2. Numerical calculations for ensemble dynamics were undertaken with the parameters relevant to B. subtilis summarized in Supplementary file 2. The bulge is discretized into $20\times 20$ bins and the cylinder is discretized into $40\times 40$ bins.
3.5. Figure 3F of the main text
Here $c=0.1$, $P=4$, and $z=\pi $. Periodic boundary conditions in $v$ are assumed. Note that we use the linear relation $\sigma =\alpha {(\mathrm{\Delta}c)}^{1}$, where $\mathrm{\Delta}c$ is the difference of principal curvatures and the value of $\alpha $ is provided in Supplementary file 2. Langevin simulations to generate a representative trajectory were undertaken with a single activated filament and the parameters relevant to E. coli summarized in Supplementary file 2, except $N=100$. Langevin simulations and numerical calculations for ensemble dynamics were undertaken with the parameters relevant to E. coli summarized in Supplementary file 2. In the finite case, the surface is discretized into $25\times 25$ bins into which individual trajectories are aggregated. In the continuum case, the surface is discretized into $200\times 200$ bins.
3.6. Figure 4A of the main text
For the ellipsoid, $a=1$ and $b=2$. For the torus, $a=1$, $R=2$, and $\mathrm{\Phi}=2\pi $. For the helix, $a=1$, $R=2$, $\mathrm{\Phi}=2\pi $, $\varphi =1$, and periodic boundary conditions in $v$ are assumed. For all geometries, Langevin simulations to generate a representative trajectory were undertaken with a single activated filament and $L=0.4$, $\sigma =0.3$, and $N=300$. Numerical calculations for ensemble dynamics were undertaken with the parameters summarized in Supplementary file 2 but $\sigma =0.3$, $N$ large enough to correspond to a fixed point for the dynamics ($N={10}^{3}$), and either $\lambda =0$ (infinite processivity) or $\lambda =\mathrm{\infty}$ (zero processivity). The surfaces are discretized into $30\times 30$ bins.
3.7. Figure 4B of the main text
Here $z=4$, $b=0.5$, and $c=1$ or $c=0.5$. Periodic boundary conditions in $v$ for the cylinders are assumed. Langevin simulations to generate a representative trajectory were undertaken with a single activated filament and $L=0.4$, $\sigma =0.3$, and $N=30$. Numerical calculations for ensemble dynamics were undertaken with the parameters summarized in Supplementary file 2 but $\sigma =0.3$, $N$ large enough to correspond to a fixed point for the dynamics ($N={10}^{3}$), and either $\lambda =0$ (infinite processivity) or $\lambda =\mathrm{\infty}$ (zero processivity). The bulges are discretized into $20\times 20$ bins and the cylinders are discretized into $40\times 40$ bins.
3.8. Figure 4C of the main text
Here $z=4$ with periodic boundary conditions in $v$. Langevin simulations to generate a representative trajectory were undertaken with a single activated filament and $L=0.4$, $\sigma =0.3$, and $N=15$. Numerical calculations for ensemble dynamics were undertaken with the parameters summarized in Supplementary file 2 but $\sigma =0.3$, $N$ large enough to correspond to a fixed point for the dynamics ($N={10}^{3}$), and either $\lambda =0$ (infinite processivity) or $\lambda =\mathrm{\infty}$ (zero processivity). The surface is discretized into $60\times 60$ bins.
3.9. Figure 3—figure supplement 1
The numerical results shown in panel a are identical to those in Figure 4B of the main text, with an additional representative trajectory in the case $c=0.2$. For the simulation in the inset of panel a, $\sigma =0$. The numerical result shown in panel b is identical to that in Figure 3E of the main text.
3.10. Figure 3—figure supplement 3
The numerical results shown in panel b are identical to those in Figure 3G of the main text, with the exception of (1) a constant value of $\sigma =0.3$ and (2) a quadratic dependence of $\sigma $ on the difference of principal curvatures, $\sigma =\beta {(\mathrm{\Delta}c)}^{2}$, where the value of $\beta $ is provided in Supplementary file 2. A numerical result corresponding to Figure 3G of the main text, but with a step size of $L=2\mu \text{m}$, is also shown. Here the same parameters summarized in Supplementary file 2 apply, except the larger value of $L$ implies the following rescaling of parameters: $L=4$, $N=6$, and $\lambda =0.66$. Finally, a numerical result corresponding to Figure 3G of the main text, but for a nonzero filament twist of ${\omega}_{0}a=5$ and Gaussian curvaturedependent activation parameter of $\gamma /{k}_{0}=1$ (see also §2.10 and Figure 4—figure supplement 4h) is shown.
3.11. Figure 4—figure supplement 1
For the torus, $a=1$, $R=2$, and $\mathrm{\Phi}=2\pi $. For the helix, $a=1$, $R=2$, $\mathrm{\Phi}=2\pi $, and $\varphi =1$ unless otherwise stated, and periodic boundary conditions in $v$ are assumed. Numerical calculations for ensemble dynamics were undertaken with the parameters summarized in Supplementary file 2 but varying $\sigma $, $N$ large enough to correspond to a fixed point for the dynamics ($N={10}^{3}$), and either $\lambda =0$ (infinite processivity) or $\lambda =\mathrm{\infty}$ (zero processivity). In panel c, $\sigma $ is fixed at $\sigma =0$ while $\varphi $ varies. The surfaces are discretized into $30\times 30$ bins.
3.12. Figure 4—figure supplement 2
Here $a=1$ and $b=2$. Numerical calculations for ensemble dynamics were undertaken with the parameters summarized in Supplementary file 2 but varying $\sigma $, $N$ large enough to correspond to a fixed point for the dynamics ($N={10}^{3}$), and either $\lambda =0$ (infinite processivity) or $\lambda =\mathrm{\infty}$ (zero processivity). The ellipsoid is discretized into $30\times 30$ bins.
3.13. Figure 4—figure supplement 3
The numerical results shown in Figure 4—figure supplement 3a are identical to those in Figure 4C in the main text. The same details apply for Figure 4—figure supplement 3c, except that $\sigma $ and $\lambda $ are varied.
3.14. Figure 4—figure supplement 4
The numerical results shown in all panels use the parameter values relevant to E. coli as summarized in Supplementary file 2. Note that we use the linear relation $\sigma =\alpha {(\mathrm{\Delta}c)}^{1}$, where $\mathrm{\Delta}c$ is the difference of principal curvatures and the value of $\alpha $ is provided in Supplementary file 2. Generally, $a=1$, $R=2$, and $\mathrm{\Phi}=2\pi $ except for panels g and h, for which $R=10$. All simulated torii were discretized into $30\times 30$ bins.
References
 1

2
Getting into shape: how do rodlike bacteria control their geometry?Systems and Synthetic Biology 8:227–235.https://doi.org/10.1007/s1169301491439

3
Thermodynamics and mechanics of membrane curvature generation and sensing by proteins and lipidsAnnual Review of Physical Chemistry 62:483–506.https://doi.org/10.1146/annurev.physchem.012809.103450
 4
 5

6
FEM2D_POISSON finite element solution of Poisson’s equation on a triangulated regionAccessed December 27, 2015.
 7

8
FtsZ in bacterial cytokinesis: cytoskeleton and force generator all in oneMicrobiology and Molecular Biology Reviews 74:504–528.https://doi.org/10.1128/MMBR.0002110

9
Bacterial filament systems: toward understanding their emergent behavior and cellular functionsJournal of Biological Chemistry 290:17181–17189.https://doi.org/10.1074/jbc.R115.637876
 10
 11
 12

13
Elastic properties of lipid bilayers: theory and possible experimentsZeitschrift für Naturforschung C 28:693–703.https://doi.org/10.1515/znc1973111209
 14
 15
 16
 17
 18
 19

20
Patch curvatureAccessed December 27, 2015.
 21
 22
 23
 24

25
Making surface plots from scatter dataAccessed December 27, 2015.
 26

27
Curved FtsZ protofilaments generate bending forces on liposome membranesThe EMBO Journal 28:3476–3484.https://doi.org/10.1038/emboj.2009.277

28
MreB orientation correlates with cell diameter in Escherichia coliBiophysical Journal 111:1035–1043.https://doi.org/10.1016/j.bpj.2016.07.017
 29

30
Shape selection of surfacebound helical filaments: biopolymers on curved membranesBiophysical Journal 111:1575–1585.https://doi.org/10.1016/j.bpj.2016.08.017
 31
 32
 33

34
Motion of variablelength MreB filaments at the bacterial cell membrane influences cell morphologyMolecular Biology of the Cell 24:2340–2349.https://doi.org/10.1091/mbc.e12100728
 35
 36
 37
 38

39
Direct membrane binding by bacterial actin MreBMolecular Cell 43:478–487.https://doi.org/10.1016/j.molcel.2011.07.008
 40
 41
 42

43
Volume of a surface triangulationAccessed December 27, 2015.
 44

45
Controlling the shape of filamentous cells of Escherichia coliNano Letters 5:1819–1823.https://doi.org/10.1021/nl0507360
 46
 47
 48
 49
 50
 51

52
Cell shape can mediate the spatial organization of the bacterial cytoskeletonBiophysical Journal 104:541–552.https://doi.org/10.1016/j.bpj.2012.12.027

53
Determination of turgor pressure in Bacillus subtilis: a possible role for K+ in turgor regulationJournal of General Microbiology 136:2521–2526.https://doi.org/10.1099/00221287136122521
 54
 55
 56

57
Exo70 generates membrane curvature for morphogenesis and cell migrationDevelopmental Cell 26:266–278.https://doi.org/10.1016/j.devcel.2013.07.007

58
Instability and deformation of a spherical vesicle by pressurePhysical Review Letters 59:2486–2488.https://doi.org/10.1103/PhysRevLett.59.2486
 59

60
How proteins produce cellular membrane curvatureNature Reviews Molecular Cell Biology 7:9–19.https://doi.org/10.1038/nrm1784
Decision letter

Tâm MignotReviewing Editor; Aix Marseille UniversityCNRS UMR7283, France

Gisela StorzSenior Editor; National Institute of Child Health and Human Development, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Mechanics and dynamics of translocating MreB filaments on curved membranes" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Gisela Storz as the Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
The manuscript by Wong et al. seeks to link the molecular details of polymer binding and motion on curved surfaces to phenomena observed in the localization and motion of MreB in rodshaped bacteria. This is timely work and focuses on a multiscale process that is both fundamental to the way bacteria grow and firmly grounded in molecular biophysics. The authors present a very large number of analytical and numerical results and make a good effort to link these to experimental observations. Of particular note, Gaussian curvature localization (which seems to be somewhat controversial) is linked in this model not to an energetic preference of the proteins for these parts of the cell but rather to a consequence of the motion and the geometry via kinetics.
A series of important changes is requested before the manuscript could be considered for publication in eLife.
1) The authors need to do a better and more evenhanded job addressing the current theoretical literature that has attempted to explain MreB localization, and comparing their results to previous theoretical works. Specifically, the authors make very different assumptions/focus on different details of the polymers than several previous papers. For example, Quint, Gopinathan and Grason et al., 2016, show that ribbonlike polymers (of which MreB is likely to be one example due to its protofilament architecture) can have both intrinsic curvature and twist which can give rise to localization at specific mean and Gaussian curvatures. Wang and Wingreen, 2013, (which is barely discussed despite being likely the most relevant previous theoretical work on this subject) discuss curvature localization in detail and show that the specific pattern of polymer deformation (e.g. where most of the deformation is focused at the polymer tips) is highly important. Why did the authors of the manuscript under review choose the details of their model and how would it be different if they had added the details considered by others and shown to affect geometric localization? What if the localization or motion did depend on Gaussian curvature or if the filaments were slightly twisted ribbons? The authors specifically state in Appendix 1 "we have assumed the filament to be bent uniformly, but the case of a curvature which varies with position along the filament length can be considered similarly." But would it have made a difference (as suggested by Wang and Wingreen)? Would this fit the data better? Worse? The same but with more parameters? The manuscript is written as if there is no other theoretical work out there and that is just not true. The authors need to both motivate their decisions and compare their results quantitatively to those of previous models in the literature. Otherwise, the reader cannot put the work into context and it is unclear which model is more likely correct given available experimental data.
2) As shown in Figure 2C the model has broader implications in terms of filament proteins and membrane binding than MreB bound to cylindrical membranes and discussion of how it could apply to other biological systems should be included. In particular, the detailed mathematical model and all its generality must be explained in a context that fits the scope of eLife. The paper contains notation commonly used in mathematics (but not that common in biophysics). For example, after Equation 3, the angular noise is denoted by N(0, σ^2) where neither N or its arguments are defined nor given biophysical significance. This trend towards a formal mathematical approach is even more prevalent in the supplementary information. Another example, is the discussion in the main text of the processivity parameters λ and k in the subsection “Dynamics of translocation”. While they are used in some of the detailed calculations in the supplementary information, the main text with the application to MreB proteins shows results for either zero processivity for which the proteins cannot attain their minimal free energy configuration in regions of high membrane curvature or for infinite processivity, in which case the proteins can move to reside in the high curvature regions of the membrane – with some fluctuations about those positions. Thus, it was not clear why the dynamic calculation was performed in such detail since the main text only showed the results for the two limiting cases of kinetic restrictions and no kinetic restrictions.
3) The kinetic equation (Equation 3) seems to coarse grain over the proteinmembrane interactions and intrinsic noise to result in a Langevin equation with an effective noise that depends on the curvatures. It is not at all clear how this is obtained from the standard physical treatment of particle motion where the Langevin equation contains a term corresponding to the physical forces (in this case, due to the curvature energies of the membrane and protein) and another corresponding to white noise (which could also model cell activity if so desired). Had the standard treatment been used, the authors might have been able to predict the steadystate position and orientation as well as the fluctuations about that steadystate using a Boltzmann factor (with an effective noise) which would have provided more intuition than the kinetic simulations that eventually attain such a steadystate.
Specific points:
 The authors state that MreB motion is processive "in contrast to spontaneous or diffusive motion," but then state "we may model the trajectories of filaments as random walks." This seems contradictory.
 "translocates along the largest principal direction." Why was this chosen? Can the authors motivate this with a molecular mechanism via which enzymes/the polymers might know which direction has the largest principle curvature? This to me is the least well motivated part of the model in terms of linking realistic molecular biophysics to the global behavior.
 The supplement is already quite lengthy and detailed, but there remain statements in the text that either state something can be shown, that some parameter/choice is rather unimportant, or simply make an assertion and state "(not shown)". I find this extremely problematic. The authors should back up all of these statements (e.g. with a derivation, calculation, simulation, etc.). This reminded me of the professor who states, "it is trivial to show…". If it is obvious, then it should be easy to show. If it is hard to show, then the offhanded remark seems incorrect.
 "filaments are increasingly misaligned in wider B. subtilis cells". Define misaligned so the reader does not have to go to the referenced paper and figure out what this means.
 Figure 3: The '0' and '\infty' designations were confusing for B and D, particularly after reading the fifth paragraph of the subsection “Dynamics of translocation”. I think this is the most confusing part of the text and after reading the text, supplement, and figures/supplementary figures I admit to not fully grasping all of it. Perhaps it could be made more clear.
 Figure 4 inset: Why is this not compared to the experimental data (e.g. from Bratton et al., 2018) as the Mean curvature plot is? This would make this result (a key one for the paper) much stronger.
 Supplementary file 1, 2: Is it possible to put confidence intervals on these estimates from different papers? I understand the authors need to make choices for their modeling, but it would be useful to know how well the field agrees on these numbers.
 Figure 3—figure supplement 1: In the fluorescence image, it is not clear where the cell/bulge are, especially compared to the rest of the figure. Maybe add a schematic?
 Figure 3—figure supplement 2 and Appendix subsection “Cylinder with small wavelength undulations”: The authors state that since Ursell et al. showed a high degree of correlation between mean and Gaussian curvature, that the small wavelength undulating model is better. But in Ursell, this was done for the channelcurved sinusoidal cells and not normal rod cells. Could this make a difference such that the long wavelength model might be better? If it was, would this make a difference?
 Figure 3—figure supplement 3: Unless I missed it, this figure is not used in the text or supplement. Does it matter that the data clearly does not resemble the assumption?
https://doi.org/10.7554/eLife.40472.023Author response
A series of important changes is requested before the manuscript could be considered for publication in eLife.
1) The authors need to do a better and more evenhanded job addressing the current theoretical literature that has attempted to explain MreB localization, and comparing their results to previous theoretical works. Specifically, the authors make very different assumptions/focus on different details of the polymers than several previous papers. For example, Quint, Gopinathan and Grason, 2016, show that ribbonlike polymers (of which MreB is likely to be one example due to its protofilament architecture) can have both intrinsic curvature and twist which can give rise to localization at specific mean and Gaussian curvatures. Wang and Wingreen, 2013, (which is barely discussed despite being likely the most relevant previous theoretical work on this subject) discuss curvature localization in detail and show that the specific pattern of polymer deformation (e.g. where most of the deformation is focused at the polymer tips) is highly important.
We thank the reviewers for these comments. We have further explained these previous studies better, and also clarified how our study is different, as follows:
1) In the section on MreB binding, we highlight the work done by Wang and Wingreen in showing that thick, twisted MreB bundles can orient along membranes and have their lengths regulated by twist [Wang and Wingreen, 2013].
2) We also discuss the work by Quint et al. demonstrating that surface regions of negative Gaussian curvature can allow twisted filaments to bind with low elastic energy [Quint et al., 2016].
3) We then detail how our work is different from these studies, in that it includes the directional motion of MreB filaments along their lengths. This directional motion is shown in the new Figure 1—video 1. Specifically, we show how the observed cellular distributions of MreB filaments across different experiments [Hussain et al., 2018; Wong et al., 2017; Kawazura et al., 2017; Ursell et al., 2014; Shi et al., 2017] can be explained with translocation alone. In principle, we have shown that the simplest assumption of MreB motion can give rise to localization.
While we agree with the reviewers that it is interesting to model the combined effects of twist and motion for general filament systems, we note that there is no in vitro evidence that membranebound MreB filaments are twisted. Rather, in all experimental observations of membranebound filaments, they are flat and untwisted (discussed below). The idea that MreB filaments contain twist is based only on theory, first arising in Wang and Wingreen’s work [Wang and Wingreen, 2013], and later supported by molecular dynamics simulations [Colavin et al., Proc. Natl. Acad. Sci. USA 111, 3585 (2014)]. The only experimental data indicating that MreB filaments might have twist are from observations of in vitro polymerization of E. coli MreB in the absence of membranes, where it aggregates into nonphysiological bundles [Nurse and Marians, J. Biol. Chem. 288, 3469 (2013)]. These bundles, which appear to have some twist, are likely driven by aggregation of the hydrophobic domains [Esue et al., J. Biol. Chem. 280, 2628 (2005)]. However, when MreB is assembled in the presence of its native membrane substrate, three separate cryoEM studies have found MreB filaments are flat and untwisted [Salje et al., 2011; van den Ent et al., 2014; Hussain et al., 2018]. When MreB is assembled on flat, supported lipid bilayers, filaments have no twist and are short [Salje et al., 2011; van den Ent et al., 2014]. When MreB filaments are assembled inside liposomes, they become much longer and wrap around the inside of the liposome, but again show no twist.
Regardless, we agree it remains a possibility MreB filaments from different species might be twisted. Thus, we now include in the paper both a more thorough discussion of the previous work, as well as an exploration of how twist affects our model predictions. As shown below in our next response, our model predictions are largely robust to twist.
Why did the authors of the manuscript under review choose the details of their model and how would it be different if they had added the details considered by others and shown to affect geometric localization? What if the localization or motion did depend on Gaussian curvature or if the filaments were slightly twisted ribbons?
We thank the reviewers for this comment. As noted above, we initially chose the parameters of our model based on the untwisted filaments found in all experimental observations of membranebound MreB filaments to date. However, we also agree with the reviewers that modeling the different effects of twist and bending rigidity for general filament systems is intriguing. A general model of filamentmembrane binding might be applicable for other species of MreB.
Per the reviewers’ suggestion, we have added new work showing the possible effects of three parameters on model predictions: these parameters are the (1) intrinsic filament twist, (2) coupling of filament activation to Gaussian curvature, and (3) filament bending rigidity. A systematic exploration of these parameters has been introduced in a new subsection “Effects of filament twist, flexural rigidity, and Gaussian curvaturedependent activation on model predictions” of the Appendix, Figure 3—figure supplement 3, and Figure 4—figure supplement 4.
Briefly, we summarize our results as follows:
1) Previous work [Quint et al., 2016] has shown that the binding angles of filaments may become biased due to filament twist: Equation 12 of that work,
predicts the angular deviation from the largest principal direction in a cylinder, ${\theta}_{0}$, due to twist in the limit that MreB filaments are fully bound to the membrane. Here $C$ is the bending rigidity of a filament, $K$ is the elastic twist stiffness, $r$ is the cell radius, and ${\omega}_{0}$ is the intrinsic helical twist. While future work should verify the robustness of this equation for nonzero turgor pressures and intrinsic filament curvatures, we found that the corresponding deviation from motion along the largest principal direction does not significantly change our model predictions for all possible values of twist on a torus.
2) We then examined the influence of Gaussian curvaturecoupling on our model predictions. While Quint et al. show that twist may contribute to Gaussian curvaturedependent activation, the coupling of filament activation to Gaussian curvature can, in general, be influenced by other cellular processes and may be independent from twist. We therefore varied the coupling parameter independently. We found that, for a broad range of such couplings, the model predictions were quantitatively similar, with larger localization at the inner edge of a torus for larger Gaussian curvature couplings. Nevertheless, our prior experiments measuring MreB fluorescence intensity in cells confined to curved microchambers suggest the magnitude of Gaussian curvature coupling, if it exists, to be small [Wong et al., 2017], and our measurements in that work can be explained without Gaussian curvature coupling (Figure 4—figure supplement 4h). We have also shown that the results pertaining to a different geometry, that of a cylinder with small wavelength undulations, are largely robust to twist and Gaussian curvature coupling in Figure 3—figure supplement 3.
3) Finally, we showed that our results are largely robust to MreB filament rigidity by modeling both the binding and motion of MreB filaments with significantly smaller and larger bending rigidities. While the details of membrane binding differ among these cases, the ensuing motion and predictions of localization do not. The main difference we observed is that MreB bundles are not predicted to bend to conform to the membrane; rather, the membrane conforms to the bundle even for physiological values of turgor in E. coli. This is qualitatively consistent with Wang and Wingreen’s model, which predicts the largest deformation to occur at the tips [Wang and Wingreen, 2013].
To clarify these points in the main text, a discussion of the works mentioned above has been added to the first paragraph of the subsection “Mechanics of binding”, where we further detail our rationale for focusing on specific parameter values. We state that our results remain largely robust to twist, Gaussian curvature coupling, and filament rigidity in a new subsection “Dependence of localization on processivity and Gaussian curvature”, where we refer readers to the Appendix subsection “Effects of filament twist, flexural rigidity, and Gaussian curvaturedependent activation on model predictions” for a more detailed discussion.
The authors specifically state in Appendix 1 "we have assumed the filament to be bent uniformly, but the case of a curvature which varies with position along the filament length can be considered similarly." But would it have made a difference (as suggested by Wang and Wingreen)? Would this fit the data better? Worse? The same but with more parameters?
As noted above in our first response to point (1), Wang and Wingreen model MreB bundles that have a significantly larger bending rigidity. While it is interesting to explore the relation between different bundle shapes and their accompanying membrane deformations (see also our second response to point (1)), we note the following. First, these large MreB bundles are not physiological, as they are only observed in the absence of membranes (as discussed further in response to point (1)). Second, we note that all three in vitro cryoEM studies of MreB have observed membranebound MreB filaments to be bent uniformly along their lengths [Hussain et al., 2018; Salje et al., 2011; van den Ent et al., 2014]. This is also a natural result if MreB filaments deform against the membrane and are short, so that the membrane curvature can be assumed uniform. Third, we note that the goal of our work is different from that of Wang and Wingreen’s: we were less interested in exploring the different ways in which thick bundles can deform membranes and more focused on clarifying why MreB filaments (1) align along principal curvatures and (2) localize differently to different cellular regions.
Furthermore, for the case of MreB filaments – the case on which we focus in this work – a calculation based on Equation S2 in the Appendix shows that varying curvature with position along the filament will result in the same predictions of alignment along directions of largest curvatures, given that the filaments are, on average, more curved than the membrane. In particular, let ${\kappa}_{s}(\mathrm{\ell})$ denote the intrinsic filament curvature as a function of position along the filament length, $\mathrm{\ell}$, and $\kappa$ denote the deformed curvature. In cylindrical cells, $\kappa$ does not vary as a function of $\ell$ because it is most energetically favorable for the filament to bend completely to match the ambient membrane curvature, as we have shown in our work. The total bending energy of the filament is then
where $B$ is the bending rigidity (identical to the variable $C$ in Quint et al.’s work mentioned above) and $L}_{f$ is the filament length. When binding to an angle that deviates from the circumferential direction in a cylinder, the deformed curvature will be smaller: let $\kappa ={\kappa}_{0}{\kappa}^{\prime}$, where $\kappa}_{0$ is the curvature along the circumferential direction of a cylinder and ${\kappa}^{\mathrm{\prime}}\ge 0$ is a constant correction to $\kappa}_{0$ depending on the binding angle. Then, the difference in bending energies between binding in the direction of $\kappa$ as opposed the circumferential $\kappa}_{0$ direction is
which is larger than zero provided the filament is, on average, more curved than the cell (that is, the second term above is nonnegative). Hence, because the binding orientation is robust, our model predictions will remain the same.
As readers may be interested in calculations similar to the above, we have appended the foregoing discussion to the Appendix subsection “Preferred orientation of filament binding and binding phase diagram”. We have, however, decided to keep matters simple in the main text by focusing on the case of a uniformly bent MreB filament.
The manuscript is written as if there is no other theoretical work out there and that is just not true. The authors need to both motivate their decisions and compare their results quantitatively to those of previous models in the literature. Otherwise, the reader cannot put the work into context and it is unclear which model is more likely correct given available experimental data.
We appreciate the reviewers’ note and have made a serious attempt to address it, as detailed in response to point (1). Specifically, we have discussed the two theoretical studies of Quint et al. and Wang and Wingreen in the first paragraph of the subsection “Mechanics of binding” and “Dependence of localization on processivity and Gaussian curvature” in the main text, where we also detail the rationale behind the parameter choices in this work (see also the response to the seventh Specific point). As it remains possible that MreB from different species may be twisted and thick, we have shown that our results are largely robust to filament twist, rigidity, and Gaussian curvature coupling in a new subsection “Effects of filament twist, flexural rigidity, and Gaussian curvaturedependent activation on model predictions”, of the Appendix. We emphasize the main difference that indeed, we know of no prior work on similar curvaturebased, directed random walks used to model MreB motion in this work.
2) As shown in Figure 2C the model has broader implications in terms of filament proteins and membrane binding than MreB bound to cylindrical membranes and discussion of how it could apply to other biological systems should be included. In particular, the detailed mathematical model and all its generality must be explained in a context that fits the scope of eLife.
We appreciate the reviewers’ insight that our results can be cast in a broader context. To address this point, we have substantially revised the main text to make it more broadly applicable as follows:
1) We better highlight previous theoretical work and discuss how our approach is different from these studies.
2) We have reorganized the Results to introduce a new subsection “Dependence of localization on processivity and Gaussian curvature”, which explores our model in generality and discusses the applicability of our results to different filament systems. In particular, we have extended our model to explore the effects of twist, Gaussian curvature coupling, and different filament bending rigidities, as detailed in our second response to point (1).
3) Finally, we have substantially rewritten the Discussion as to discuss potential applications of our results to other contexts, including bacterial cytokinesis and membrane trafficking, growth, and movement in both prokaryotes and eukaryotes.
We hope that these changes better highlight the broad scope of our work.
The paper contains notation commonly used in mathematics (but not that common in biophysics). For example, after Equation 3, the angular noise is denoted by N(0, σ^2) where neither N or its arguments are defined nor given biophysical significance. This trend towards a formal mathematical approach is even more prevalent in the supplementary information.
We have clarified the meaning of the notation in the main text and, whenever possible, in the Appendix. This includes the following changes:
1) The offending phrase and its surrounding context has been edited to read “Here and below, we assume ${\xi}_{n}\sim N(0,{\sigma}^{2})$—that is, the angular noise is normallydistributed, with mean zero and a variance, $\sigma}^{2$, to be inferred from data—and note that the translocation noise, $\sigma$, may depend on quantities such as the principal curvatures, as discussed later”.
2) “where $(x,y)\in {\mathbb{R}}^{2}$” has been changed to “where $x$ and $y$ are real numbers”.
3) The meaning of the notation $\mathrm{\nabla}h$ is clarified as “the gradient of $h$”.
4) The phrase “$\mathbf{r}=r(u,v)\subset {\mathbb{R}}^{3}$” has been clarified. The phrase and its surrounding context now reads “We consider the membrane as a parametric surface, $r=r(u,v)$, embedded in threedimensional space ($\mathbb{R}}^{3$) with surface coordinates $u$ and $v$”.
5) The notation $O[{\left(\mathrm{\nabla}h\right)}^{2}]$ has been clarified as “where the bigO notation signifies$H\frac{1}{2}{\mathrm{\nabla}}^{2}h\le M(\mathrm{\nabla}h{)}^{2}$when $0\phantom{\rule{thinmathspace}{0ex}}<(\mathrm{\nabla}h{)}^{2}\phantom{\rule{thinmathspace}{0ex}}<\delta$ for some positive numbers $\delta$ and $M$”.
6) The notation $\mathrm{\partial}\mathrm{\Omega}$ has been clarified as “the boundary of $\mathrm{\Omega}$”.
7) The notation ${\xi}_{n}\sim N(0,{\sigma}^{2})$ is again clarified in the Appendix. The revised sentence reads “Here and below, ${\xi}_{n}\sim N(0,{\sigma}^{2})$—that is, the angular noise is normallydistributed, with mean 0 and variance $\sigma}^{2$”.
We thank the reviewers for pointing out this issue and hope that these clarifications will make our work more accessible to a nonmathematically inclined audience.
Another example, is the discussion in the main text of the processivity parameters λ and k in the subsection “Dynamics of translocation”. While they are used in some of the detailed calculations in the supplementary information, the main text with the application to MreB proteins shows results for either zero processivity for which the proteins cannot attain their minimal free energy configuration in regions of high membrane curvature or for infinite processivity, in which case the proteins can move to reside in the high curvature regions of the membrane – with some fluctuations about those positions. Thus, it was not clear why the dynamic calculation was performed in such detail since the main text only showed the results for the two limiting cases of kinetic restrictions and no kinetic restrictions.
We apologize for the confusion. The cases of zero and infinite processivity were explored to demonstrate limiting cases, as to first cast our results in a broader context and understand the underlying physics. To better motivate this exploration, we have significantly reorganized our Results as to first discuss the implications of our model to MreB (subsection “Implications to MreB localization”) and then explore our model in a broader context (subsection “Dependence of localization on processivity and Gaussian curvature”). The results for MreB, which involve a finite processivity (Supplementary file 2) and are neither in the limit of zero nor infinite processivity, are now shown in Figure 3. The subsequent exploration of the physics is now shown in Figure 4.
To further clarify the relevance of finite processivity to MreB, we have:
1) Appended a clarifying sentence to the legend of Figure 3: “In this figure, the localization of MreB filaments, which exhibit a finite processivity and are assumed to follow the parameter values summarized in Supplementary file 2, is shown,”
2) Revised a relevant sentence in the legend of Figure 4 to read “Note that cases of zero processivity correspond to uniform distributions and that we have considered the limiting cases of zero and infinite processivity here. Figure 3 shows numerical results for the case of a large, but finite, processivity relevant to MreB,” and changed relevant phrases in the main text to better emphasize the MreBspecific (Figure 3) and general (Figure 4) processivity values adopted, including “While we assume the values of λ and σ to be based on these measurements, we examine the effects of varying λ and σ in the following section” and “in the case of small λ corresponding to large processivity – a limiting case that is relevant to MreB”.
We hope this clarification also better motivates some of the calculations undertaken in the Appendix.
Finally, we note that processivity plays an important role in our model: in the limit of zero processivity, MreB filaments are assumed to be deactivated as soon as they are activated. If the activation of MreB filaments were simply proportional to area, as we have assumed in our work, then the binding orientation is inconsequential and no localization can occur, as shown in Figure 4. However, we would also like to emphasize that, at least for the model considered in our work, there is generally no explicit connection between processivity and (1) attaining a minimal free energy configuration or (2) moving to reside in highly curved membrane regions in our work. For (1), we hypothesize the minimal free energy configuration to be relevant for binding and orienting the filament only (subsection “Mechanics of binding”). The processivity is a variable that determines the lifetime of MreB filaments while moving and is not necessarily related to binding. As for (2), the movement of MreB filaments along the largest principal directions need not entail that they will end up residing in highly curved membrane regions. Whether or not such correlation manifests depends on the geometry at hand, as detailed more carefully in the new subsection “Dependence of localization on processivity and Gaussian curvature”. In other words, localization is not generically correlated with scalar measures of curvature. Figure 4C illustrates this point succinctly: in this geometry, motion along the largest principal directions results in accumulation of filaments, but the Gaussian curvature vanishes everywhere. In contrast, correlations between filament concentration and the Gaussian curvature exist for the geometry shown in Figure 3F.
3) The kinetic equation (Equation 3) seems to coarse grain over the proteinmembrane interactions and intrinsic noise to result in a Langevin equation with an effective noise that depends on the curvatures. It is not at all clear how this is obtained from the standard physical treatment of particle motion where the Langevin equation contains a term corresponding to the physical forces (in this case, due to the curvature energies of the membrane and protein) and another corresponding to white noise (which could also model cell activity if so desired). Had the standard treatment been used, the authors might have been able to predict the steadystate position and orientation as well as the fluctuations about that steadystate using a Boltzmann factor (with an effective noise) which would have provided more intuition than the kinetic simulations that eventually attain such a steadystate.
We thank the reviewers for this comment. The main point raised is the justification of using kinetic simulations that attain a steady state but do not relate directly to a Langevin equation describing an equilibrium process, whose steady state corresponds to the Boltzmann distribution. A crucial difference between our work and such an approach is that the dynamics here are inherently away from equilibrium, which implies that the steadystate solution will not generally be well described by the Boltzmann distribution. The nonequilibrium nature of the dynamics can be appreciated by observing that detailed balance is broken due to directional motion. To see this for the cylinder, we may consider the transitions in (X_{n}, χ_{n})space between two small area elements as below; as in the main text, X_{n} and χ_{n} respectively denote the position and direction of the filament at a timestep n.
Note that the steadystate probability of a filament in (X_{n}, χ_{n})space is uniform: it is uniform in position space due to symmetry and in direction space due to the hypothesis that χ_{n} is uniformly random (but does not vary in n). However, detailed balance does not hold in Author response image 1 because, for small deviatory angles between the two area elements, the transition rate corresponding to motion in the direction of χ_{n} is larger than the transition rate corresponding to motion in the direction of −χ_{n}. Hence, the dynamics is out of equilibrium. This argument applies across many of the simulations considered in our work. Intriguingly, the nonequilibrium dynamics implies the expenditure of energy, and prior studies have shown that MreB filaments require the energyintensive process of cell wall synthesis in order to move processively [Garner et al., 2011; van Teeffelen et al., 2011; DomínguezEscobar et al., 2011].
More broadly, while we have been able to obtain theoretical predictions for displacements and fluctuations in the simple case of an infinite cylinder (subsection “Cylinder” of the Appendix) and infer the filament concentrations for other limited geometries (subsections “Torus and Helix”, “Ellipsoid”, and “Filament concentration is independent of Gaussian curvature” of the Appendix), we have not been able to recast Equation 3 into a more accessible form that, even away from equilibrium, would permit a general solution. We therefore thought it best to describe the stochastic process of Equation 3 accessibly and resort to numerical simulation. We hope that this will provide a foundation for further theoretical work which builds on our model.
Specific points:
 The authors state that MreB motion is processive "in contrast to spontaneous or diffusive motion," but then state "we may model the trajectories of filaments as random walks." This seems contradictory.
We understand the reviewer’s confusion, as the motion of MreB filaments is processive and highly directional. We use the term “random walks” to designate any stochastic process which describes a path comprising a succession of random steps. To mitigate confusion, we now use the term “biased random walks” and clarify, in the following sentence, that “a ‘biased random walk’ refers to a succession of random steps which may be processive: while the meansquared displacement of a filament will be approximately quadratic, and not linear, in time, the processive motion we consider is random only because the translocation direction can deviate from directions of largest membrane curvature due to sources of stochasticity”.
 "translocates along the largest principal direction." Why was this chosen? Can the authors motivate this with a molecular mechanism via which enzymes/the polymers might know which direction has the largest principle curvature? This to me is the least well motivated part of the model in terms of linking realistic molecular biophysics to the global behavior.
The oriented motion of MreB filaments has been welldocumented in both E. coli and B. subtilis. In cells, MreB filaments move directionally, and this motion is oriented predominantly around the cell width and perpendicular to the long axis [Garner et al., 2011; van Teeffelen et al., 2011; DomínguezEscobar et al., 2011]. In our original work [Hussain et al., 2018], we worked to understand what caused MreB motion to be circumferentially oriented. We found that MreB orientation is most likely determined only by MreB filaments binding to membranes. The orientation arises as these highly inwardlycurved filaments bind to the membrane in a manner which maximizes their interaction energy – i.e., filaments orient along the direction of largest curvature. We demonstrated that MreB translocates along the direction of largest curvature in a variety of cell shapes, including rods, round cells, and bulged cells. Furthermore, we showed that MreB filaments also orient in directions determined by their binding orientations in both spheroplasted cells and liposomes (Author response image 2), in agreement with previous work [Olshausen et al., 2013]. In this work, we now explore the implications of filament orientation and directional motion on the cellular distribution of filaments.
As it appears this background was not evident in our initial draft, in the updated manuscript we have summarized the evidence for principal curvaturebased motion in a new Figure 1C and the section on Dynamics of translocation. We hope this improved discussion of previous work clarifies the rationale for our model.
If the reviewers are querying in regards to how the directional motion of MreB filaments arises, we note that the exact mechanism is not known, although this motion requires – and is perhaps driven by – cell wall synthesis [Garner et al., 2011; van Teeffelen et al., 2011; DomínguezEscobar et al., 2011]. We believe many possibilities exist, but have decided not to speculate upon them in the text.
 The supplement is already quite lengthy and detailed, but there remain statements in the text that either state something can be shown, that some parameter/choice is rather unimportant, or simply make an assertion and state "(not shown)". I find this extremely problematic. The authors should back up all of these statements (e.g. with a derivation, calculation, simulation, etc.). This reminded me of the professor who states, "it is trivial to show…". If it is obvious, then it should be easy to show. If it is hard to show, then the offhanded remark seems incorrect.
We thank the reviewers for this comment and have sought to address it by backing up all such statements as follows:
1) In the main text, we have replaced the only instance of the phrase “it can be shown” before “$C}_{F$ is uniform over the surface of an ellipsoid” with a reference to the Appendix, which indeed shows this via both calculations and simulations (subsection “Ellipsoid” of the Appendix and Figure 4—figure supplement 2).
2) In subsection “Model of a protein filament binding to a membrane” of the Appendix, we have removed the sentence “For the parameter values summarized in Supplementary file 2, it can be shown that l_{b}=L_{f} below,” as it is not necessary to claim this. We have instead mentioned in subsection “Finite element solutions of the shape equation” of the Appendix, that the numerical optimization runs over varying sizes of U, which corresponds to the bound length of the filament, and that they predict l_{b}=L_{f}. This result is consistent with the linear dependence on l_{b} of the free energy change due to binding in the case of physiological turgor pressures, as we have clarified following Equation S11 of the Appendix.
3) The claim that “our results do not significantly change for different b,” the fractional crosssectional binding region, has been supplemented with a reference to Figure 4 of [Hussain et al., 2018] which shows this in panels F and G.
4) In subsections “Torus and Helix” and “Ellipsoid” of the Appendix, claims that certain results hold in the limit of large noise (σ→∞) have been corroborated by showing the results of additional simulations in this limit in Figure 4—figure supplements 1 and 2.
5) In subsection “Cylinder with a bulge” of the Appendix, the claim that bulges attract filaments in the case of no noise and large processivity is represented by a simulation which shows this (inset of Figure 3—figure supplement 1A).
6) In subsection "Effects of principal curvaturedependent translocation noise and varying filament step size on model predictions”, claims that varying the dependence of σ on ∆c do not significantly change our results are corroborated with new simulation results in a revised Figure 3—figure supplement 3 (see also the response to the last Specific point).
We believe these revisions address all instances of unsupported assertions, both in the main text and in the Appendix.
 "filaments are increasingly misaligned in wider B. subtilis cells". Define misaligned so the reader does not have to go to the referenced paper and figure out what this means.
We thank the reviewers for drawing our attention to this point. We have revised the offending phrase and its surrounding context to read “As the distribution of MreB filament angles gradually becomes broader as B. subtilis cells become wider” and hope this clarifies our meaning.
 Figure 3: The '0' and '\infty' designations were confusing for B and D, particularly after reading the fifth paragraph of the subsection “Dynamics of translocation”. I think this is the most confusing part of the text and after reading the text, supplement, and figures/supplementary figures I admit to not fully grasping all of it. Perhaps it could be made more clear.
We thank the reviewers for this comment and hope that they find this point clarified in our third response to point (2). To clarify the definition of processivity, we have added the phrase “that is, the mean number of steps that a filament takes on the membrane surface before becoming inactive” upon first mention of processivity.
 Figure 4 inset: Why is this not compared to the experimental data (e.g. from Bratton et al., 2018) as the Mean curvature plot is? This would make this result (a key one for the paper) much stronger.
We thank the reviewers for pointing this out and have split the old Figure 4E into two panels, which are now shown in Figure 3G. Here, we individually compare our predictions with statistics from the works of Ursell et al. and Shi et al. [Ursell et al., 2014; Shi et al., 2017] for the mean curvature and Bratton et al. [Bratton et al., 2018] for the Gaussian curvature.
While the geometry we have assumed in Figure 3F exhibits a larger range of mean and Gaussian curvatures than is needed to explain both the mean and Gaussian curvature data, we note that Bratton et al.’s measurements of filament enrichment in wildtype E. coli cells are quantitatively consistent with our predictions and indicates a decreased concentration of MreB filaments at regions of positive Gaussian curvature. We therefore believe that the revised Figure 3G strengthens our comparison to previous work.
 Supplementary file 1, 2: Is it possible to put confidence intervals on these estimates from different papers? I understand the authors need to make choices for their modeling, but it would be useful to know how well the field agrees on these numbers.
We appreciate the reviewers’ point. To address it, we have performed a detailed literature search and appended different estimates of parameters in Supplementary files 1 and 2, along with relevant references. We hope that this additional information better contextualizes the values chosen.
 Figure 3—figure supplement 1: In the fluorescence image, it is not clear where the cell/bulge are, especially compared to the rest of the figure. Maybe add a schematic?
We thank the reviewers for pointing this out and have added a schematic to Figure 3—figure supplement 1.
 Figure 3—figure supplement 2 and Appendix 1 subsection “Cylinder with small wavelength undulations”: The authors state that since Ursell et al. showed a high degree of correlation between mean and Gaussian curvature, that the small wavelength undulating model is better. But in Ursell, this was done for the channelcurved sinusoidal cells and not normal rod cells. Could this make a difference such that the long wavelength model might be better? If it was, would this make a difference?
We thank the reviewers for this comment. The long wavelength model cannot make a difference for an additional reason: the ranges of Gaussian and mean curvatures predicted by such a model would be significantly smaller than those measured for wildtype, unconfined, sinusoidallyconfined, thin, and wide E. coli cells from the works of Ursell et al., Shi et al., and Bratton et al. [Ursell et al., 2014; Shi et al., 2017; Bratton et al., 2018]. For instance, the range of mean curvature predicted by the long wavelength model in Figure 3—figure supplement 2 is, in dimensionful units, [.994 μm^{1}, 1.012 μm^{1}]. The data for both confined and unconfined cells from Ursell et al.exhibit a much larger range of approximately [0 μm^{1}, 2 μm^{1}]. Thus, strong regions of curvature are necessary to recapitulate the experimental data.
To clarify this point further, we have introduced additional text to the subsection “Cylinder with small wavelength undulations” of the Appendix, and a relevant sentence now reads: “However, the Gaussian and mean curvatures are uncorrelated in this case, inconsistent with the positive correlation observed in [Ursell et al., 2014] for cells growing in sinusoidallyshaped chambers, and the ranges of Gaussian and mean curvatures are significantly smaller than those measured for wildtype, unconfined, sinusoidallyconfined, thin, and wide cells in different experiments (Figure 3—figure supplement 3A and Figure 3G of the main text).” We hope that this better motivates our consideration of a short wavelength geometry.
 Figure 3—figure supplement 3: Unless I missed it, this figure is not used in the text or supplement. Does it matter that the data clearly does not resemble the assumption?
Figure 3—figure supplement 3 is referenced in subsection “Effects of principal curvaturedependent translocation noise and varying filament step size on model predictions” of the Appendix, in addressing the possibility that the translocation noise, σ, depends on the difference of principal curvatures, ∆c (subsection “Implications to MreB localization”, first paragraph in the main text). The figure shows the empirically observed dependence of σon cell width, along with a linear fit. While the linear fit clearly does not capture all the variability in the data, the point is that similar dependencies (including the quadratic fit below) of σ on ∆c do not qualitatively change our results. To show this explicitly (and in the spirit of our response to the third Specific point), we have regenerated the theoretical predictions of Figure 3G, one of the main results of this work, for different fits of σ as a function of∆c. As discussed in response to points and (3), we have also done the same for different values of filament twist and Gaussian curvature coupling. For all cases, we have added new panels to Figure 3—figure supplement 3 showing the qualitative robustness: all assumptions still result in filament enhancement at regions of negative Gaussian curvature.
We thank the reviewers for giving us the opportunity to clarify this point.
https://doi.org/10.7554/eLife.40472.024Article and author information
Author details
Funding
National Science Foundation (DGE1144152)
 Felix Wong
Quantitative Biology Initiative at Harvard
 Felix Wong
National Institutes of Health (DP2AI11792301)
 Ethan C Garner
Smith Family Award
 Ethan C Garner
Searle Scholar Fellowship
 Ethan C Garner
Volkswagen Foundation
 Ethan C Garner
 Ariel Amir
Materials Research and Engineering Center at Harvard
 Ariel Amir
Kavli Institute for Bionano Science and Technology at Harvard
 Ariel Amir
Alfred P. Sloan Foundation
 Ariel Amir
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
FW was supported by the National Science Foundation Graduate Research Fellowship under grant no. DGE1144152 and the Quantitative Biology Initiative at Harvard. ECG was supported by the National Institutes of Health under grant no. DP2AI11792301, the Smith Family Award, and the Searle Scholar Fellowship. AA was supported by the Materials Research and Engineering Center at Harvard, the Kavli Institute for Bionano Science and Technology at Harvard, and the Alfred P Sloan Foundation. ECG and AA were supported by the Volkswagen Foundation. We thank Carl Wivagg, Saman Hussain, Ned Wingreen, and Siyuan (Steven) Wang for discussions and Sven van Teeffelen, Jie Lin, and PoYi Ho for comments on the manuscript.
Senior Editor
 Gisela Storz, National Institute of Child Health and Human Development, United States
Reviewing Editor
 Tâm Mignot, Aix Marseille UniversityCNRS UMR7283, France
Publication history
 Received: August 2, 2018
 Accepted: February 15, 2019
 Accepted Manuscript published: February 18, 2019 (version 1)
 Version of Record published: May 7, 2019 (version 2)
Copyright
© 2019, Wong et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 2,092
 Page views

 459
 Downloads

 5
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.