Structural polymorphism of the cytoskeleton: a model of linker-assisted filament aggregation

Charge Regulation Triggers Condensation of Short Oligopeptides to Polyelectrolytes

Electrostatic interactions between charged macromolecules are ubiquitous in biological systems, and they are important also in materials design. Attraction between oppositely charged molecules is often interpreted as if the molecules had a fixed charge, which is not affected by their interaction. Less commonly, charge regulation is invoked to interpret such interactions, i.e., a change of the charge state in response to a change of the local environment. Although some theoretical and simulation studies suggest that charge regulation plays an important role in intermolecular interactions, experimental evidence supporting such a view is very scarce. In the current study, we used a model system, composed of a long polyanion interacting with cationic oligolysines, containing up to 8 lysine residues. We showed using both simulations and experiments that while these lysines are only weakly charged in the absence of the polyanion, they charge up and condense on the polycations if the pH is close to the pKa of the lysine side chains. We show that the lysines coexist in two distinct populations within the same solution: (1) practically nonionized and free in solution; (2) highly ionized and condensed on the polyanion. Using this model system, we demonstrate under what conditions charge regulation plays a significant role in the interactions of oppositely charged macromolecules and generalize our findings beyond the specific system used here.

DNA-protamine condensates under low salt conditions: molecular dynamics simulation with a simple coarse-grained model focusing on electrostatic interactions

Protamine, a small, strongly positively-charged protein, plays a key role in achieving chromatin condensation inside sperm cells and is also involved in the formulation of nanoparticles for gene therapy and packaging of mRNA-based vaccines against viral infection and cancer. The detailed mechanisms of such condensations are still poorly understood especially under low salt conditions where electrostatic interaction predominates. Our previous study, with a refined coarse-grained model in full consideration of the long-range electrostatic interactions, has demonstrated the crucial role of electrostatic interaction in protamine-controlled reversible DNA condensation. Therefore, we herein pay our attention only to the electrostatic interaction and devise a coarser-grained bead-spring model representing the right linear charge density on protamine and DNA chains but treating other short-range interactions as simply as possible, which would be suitable for real-scale simulations. Effective pair potential calculations and large-scale molecular dynamics simulations using this extremely simple model reproduce the phase behaviour of DNA in a wide range of protamine concentrations under low salt conditions, again revealing the importance of the electrostatic interaction in this process and providing a detailed nanoscale picture of bundle formation mediated by a charge disproportionation mechanism. Our simulations also show that protamine length alters DNA overcharging and in turn redissolution thresholds of DNA condensates, revealing the important role played by entropies and correlated fluctuations of condensing agents and thus offering an additional opportunity to design tailored nanoparticles for gene therapy. The control mechanism of DNA-protamine condensates will also provide a better microscopic picture of biomolecular condensates, i.e., membraneless organelles arising from liquid-liquid phase separation, that are emerging as key principles of intracellular organization. Such condensates controlled by post-translational modification of protamine, in particular phosphorylation, or by variations in protamine length from species to species may also be responsible for the chromatin-nucleoplasm patterning observed during spermatogenesis in several vertebrate and invertebrate species.

Organization and dynamics of cross-linked actin filaments in confined environments

The organization of the actin cytoskeleton is impacted by the interplay between physical confinement, features of cross-linking proteins, and deformations of semiflexible actin filaments. Some cross-linking proteins preferentially bind filaments in parallel, although others bind more indiscriminately. However, a quantitative understanding of how the mode of binding influences the assembly of actin networks in confined environments is lacking. Here we employ coarse-grained computer simulations to study the dynamics and organization of semiflexible actin filaments in confined regions upon the addition of cross-linkers. We characterize how the emergent behavior is influenced by the system shape, the number and type of cross-linking proteins, and the length of filaments. Structures include isolated clusters of filaments, highly connected filament bundles, and networks of interconnected bundles and loops. Elongation of one dimension of the system promotes the formation of long bundles that align with the elongated axis. Dynamics are governed by rapid cross-linking into aggregates, followed by a slower change in their shape and connectivity. Cross-linking decreases the average bending energy of short or sparsely connected filaments by suppressing shape fluctuations. However, it increases the average bending energy in highly connected networks because filament bundles become deformed, and small numbers of filaments exhibit long-lived, highly unfavorable configurations. Indiscriminate cross-linking promotes the formation of high-energy configurations due to the increased likelihood of unfavorable, difficult-to-relax configurations at early times. Taken together, this work demonstrates physical mechanisms by which cross-linker binding and physical confinement impact the emergent behavior of actin networks, which is relevant both in cells and in synthetic environments.

Kinetochore- and chromosome-driven transition of microtubules into bundles promotes spindle assembly

Mitotic spindle assembly is crucial for chromosome segregation and relies on bundles of microtubules that extend from the poles and overlap in the middle. However, how these structures form remains poorly understood. Here we show that overlap bundles arise through a network-to-bundles transition driven by kinetochores and chromosomes. STED super-resolution microscopy reveals that PRC1-crosslinked microtubules initially form loose arrays, which become rearranged into bundles. Kinetochores promote microtubule bundling by lateral binding via CENP-E/kinesin-7 in an Aurora B-regulated manner. Steric interactions between the bundle-associated chromosomes at the spindle midplane drive bundle separation and spindle widening. In agreement with experiments, theoretical modeling suggests that bundles arise through competing attractive and repulsive mechanisms. Finally, perturbation of overlap bundles leads to inefficient correction of erroneous kinetochore-microtubule attachments. Thus, kinetochores and chromosomes drive coarsening of a uniform microtubule array into overlap bundles, which promote not only spindle formation but also chromosome segregation fidelity.

Theory of Weakly Polydisperse Cytoskeleton Filaments

Cytoskeleton filaments have the extraordinary ability to change conformations dynamically in response to alterations of the number density of actins/tubulin, the number density and type of binding agents, and the electrolyte concentration. This property is crucial for eukaryotic cells to achieve specific biological functions in different cellular compartments. Conventional approaches to biopolymers’ solution break down for cytoskeleton filaments because they entail several approximations to treat their polyelectrolyte and mechanical properties. In this article, we introduce a novel density functional theory for polydisperse, semiflexible cytoskeleton filaments. The approach accounts for the equilibrium polymerization kinetics, length and orientation filament distributions, as well as the electrostatic interaction between filaments and the electrolyte. This is essential for cytoskeleton polymerization in different cell compartments generating filaments of different lengths, sometimes long enough to become semiflexible. We characterized the thermodynamics properties of actin filaments in electrolyte aqueous solutions. We calculated the free energy, pressure, chemical potential, and second virial coefficient for each filament conformation. We also calculated the phase diagram of actin filaments’ solution and compared with the corresponding results in in vitro experiments.

Interplay between Brownian motion and cross-linking kinetics controls bundling dynamics in actin networks

Morphology changes in cross-linked actin networks are important in cell motility, division, and cargo transport. Here, we study the transition from a weakly cross-linked network of actin filaments to a heavily cross-linked network of actin bundles through microscopic Brownian dynamics simulations. We show that this transition occurs in two stages: first, a composite bundle network of small and highly aligned bundles evolves from cross-linking of individual filaments and, second, small bundles coalesce into the clustered bundle state. We demonstrate that Brownian motion speeds up the first stage of this process at a faster rate than the second. We quantify the time to reach the composite bundle state and show that it strongly increases as the mesh size increases only when the concentration of cross-links is small and that it remains roughly constant if we decrease the relative ratio of cross-linkers as we increase the actin concentration. Finally, we examine the dependence of the bundling timescale on filament length, finding that shorter filaments bundle faster because they diffuse faster.

Stratification relieves constraints from steric hindrance in the generation of compact actomyosin asters at the membrane cortex

Recent in vivo studies reveal that several membrane proteins are driven to form nanoclusters by active contractile flows arising from localized dynamic patterning of F-actin and myosin at the cortex. Since myosin-II assemble as minifilaments with tens of myosin heads, one might worry that steric considerations would obstruct the emergence of nanoclustering. Using coarse-grained, agent-based simulations that account for steric constraints, we find that the patterns exhibited by actomyosin in two dimensions, do not resemble the steady-state patterns in our in vitro reconstitution of actomyosin on a supported bilayer. We perform simulations in a thin rectangular slab, separating the layer of actin filaments from myosin-II minifilaments. This recapitulates the observed features of in vitro patterning. Using super resolution microscopy, we find evidence for such stratification in our in vitro system. Our study suggests that molecular stratification may be an important organizing feature of the cortical cytoskeleton in vivo.

Short and Soft: Multidomain Organization, Tunable Dynamics, and Jamming in Suspensions of Grafted Colloidal Cylinders with a Small Aspect Ratio

The yet virtually unexplored class of soft colloidal rods with a small aspect ratio is investigated and shown to exhibit a very rich phase and dynamic behavior, spanning from liquid to nearly melt state. Instead of the nematic order, these short and soft nanocylinders alter their organization with increasing concentration from isotropic liquid with random orientation to small domains with preferred local orientation and eventually a multidomain arrangement with a local orientational order. The latter gives rise to a kinetically suppressed state akin to structural glass with detectable terminal relaxation, which, on further increasing concentration, reveals features of hexagonally packed order as in ordered block copolymers. The respective dynamic response comprises four regimes, all above the overlapping concentration of 0.02 g/mL:(I) from 0.03 to 0.1 g/mol, the system undergoes a liquid-to-solidlike transition with a structural relaxation time that grows by 4 orders of magnitude. (II) From 0.1 to 0.2 g/mL, a dramatic slowing-down is observed and is accompanied by an evolution from isotropic to a multidomain structure. (III) Between 0.2 and 0.6 g/mol, the suspensions exhibit signatures of shell interpenetration and jamming, with the colloidal plateau modulus depending linearly on concentration. (IV) At 0.74 g/mL, in the densely jammed state, the viscoelastic signature of hexagonally packed cylinders from microphase-separated block copolymers is detected. These properties set short and soft nanocylinders apart from long colloidal rods (with a large aspect ratio) and provide insights for fundamentally understanding the physics in this intermediate soft colloidal regime and for tailoring the flow properties of nonspherical soft colloids.

Distinct Binding Modes of Vinculin Isoforms Underlie Their Functional Differences

Vinculin and its splice isoform metavinculin play key roles in regulating cellular morphology, motility, and force transduction. Vinculin is distinct from metavinculin in its ability to bundle filamentous actin (F-actin). To elucidate the molecular basis for these differences, we employed computational and experimental approaches. Results from these analyses suggest that the C terminus of both vinculin and metavinculin form stable interactions with the F-actin surface. However, the metavinculin tail (MVt) domain contains a 68 amino acid insert, with helix 1 (H1) sequestered into a globular subdomain, which protrudes from the F-actin surface and prevents actin bundling by sterically occluding actin filaments. Consistent with our model, deletion and selective point mutations within the MVt H1 disrupt this protruding structure, and facilitate actin bundling similar to vinculin tail (Vt) domain.

Synthetic Transient Crosslinks Program the Mechanics of Soft, Biopolymer-Based Materials

Actin networks are adaptive materials enabling dynamic and static functions of living cells. A central element for tuning their underlying structural and mechanical properties is the ability to reversibly connect, i.e., transiently crosslink, filaments within the networks. Natural crosslinkers, however, vary across many parameters. Therefore, systematically studying the impact of their fundamental properties like size and binding strength is unfeasible since their structural parameters cannot be independently tuned. Herein, this problem is circumvented by employing a modular strategy to construct purely synthetic actin crosslinkers from DNA and peptides. These crosslinkers mimic both intuitive and noncanonical mechanical properties of their natural counterparts. By isolating binding affinity as the primary control parameter, effects on structural and dynamic behaviors of actin networks are characterized. A concentration-dependent triphasic behavior arises from both strong and weak crosslinkers due to emergent structural polymorphism. Beyond a certain threshold, strong binding leads to a nonmonotonic elastic pulse, which is a consequence of self-destruction of the mechanical structure of the underlying network. The modular design also facilitates an orthogonal regulatory mechanism based on enzymatic cleaving. This approach can be used to guide the rational design of further biomimetic components for programmable modulation of the properties of biomaterials and cells.

The dynamics of filament assembly define cytoskeletal network morphology

The actin cytoskeleton is a key component in the machinery of eukaryotic cells, and it self-assembles out of equilibrium into a wide variety of biologically crucial structures. Although the molecular mechanisms involved are well characterized, the physical principles governing the spatial arrangement of actin filaments are not understood. Here we propose that the dynamics of actin network assembly from growing filaments results from a competition between diffusion, bundling and steric hindrance, and is responsible for the range of observed morphologies. Our model and simulations thus predict an abrupt dynamical transition between homogeneous and strongly bundled networks as a function of the actin polymerization rate. This suggests that cells may effect dramatic changes to their internal architecture through minute modifications of their nonequilibrium dynamics. Our results are consistent with available experimental data.

The dynamics of actin cytoskeleton is essential to the function of living cells. Here, Foffano et al. describe a nonequilibrium filament model to mimic the formation of cytoskeleton and pinpoint the key role played by the actin entanglement during the transition from homogeneous to bundled networks.

Semiflexible Biopolymers in Bundled Arrangements

Bundles and networks of semiflexible biopolymers are key elements in cells, lending them mechanical integrity while also enabling dynamic functions. Networks have been the subject of many studies, revealing a variety of fundamental characteristics often determined via bulk measurements. Although bundles are equally important in biological systems, they have garnered much less scientific attention since they have to be probed on the mesoscopic scale. Here, we review theoretical as well as experimental approaches, which mainly employ the naturally occurring biopolymer actin, to highlight the principles behind these structures on the single bundle level.

Electrostatics and depletion determine competition between 2D nematic and 3D bundled phases of rod-like DNA nanotubes

Rod-like particles form solutions of technological and biological importance. In particular, biofilaments such as actin and microtubules are known to form a variety of phases, both in vivo and in vitro, whose appearance can be controlled by depletion, confinement, and electrostatic interactions. Here, we utilize DNA nanotubes to undertake a comprehensive study of the effects of those interactions on two particular rod-like phases: a 2D nematic phase consisting of aligned rods pressed against a glass surface, and a 3D bundled network phase. We experimentally measure the stability of these two phases over a range of depletant concentrations and ionic strengths, finding that the 2D phase is slightly more stable than the 3D phase. We formulate a quantitative model of phase stability based on consideration of pairwise rod-rod and rod-surface interactions; notably, we include a careful accounting of solution electrostatics interactions using an effective-charge strategy. The model is relatively simple and contains no free parameters, yet predicts phase boundaries in good agreement with the experiment. Our results indicate that electrostatic interactions, rather than depletion, are largely responsible for the enhanced stability of the 2D phase. This work provides insight into the polymorphism of rod-like solutions, indicating why certain phases appear, and providing a means (and a predictive model) for controlling those phases.

Bundling in semiflexible polymers: A theoretical overview

Supramolecular assemblies of polymers are key modules to sustain the structure of cells and their function. The main elements of these assemblies are charged semiflexible polymers (polyelectrolytes) generally interacting via a long(er)-range repulsion and a short(er)-range attraction. The most common supramolecular structure formed by these polymers is the bundle. In the present paper, we critically review some recent theoretical and computational advances on the problem of bundle formation, and point a few promising directions for future work.

Computational analysis of morphologies and phase transitions of cross-linked, semi-flexible polymer networks

The eukaryotic cytoskeleton is a protein fibre network mainly consisting of the semi-flexible biopolymer F-actin, microtubules and intermediate filaments. It is well known to exhibit a pronounced structural polymorphism, which enables intracellular processes such as cell adhesion, cell motility and cell division. We present a computational study on cross-linked networks of semi-flexible polymers, which offers a detailed analysis of the network structure and phase transitions from one morphology to another. We elaborate the morphological differences, their mechanical implications and the order of the observed phase transitions. Finally, we present a perspective on how the information gained in our simulations can be exploited in order to build both flexible and accurate, microstructurally informed, homogenized constitutive models of the cytoskeleton.

Designing highly tunable semiflexible filament networks

Semiflexible polymers can generate a range of filamentous networks significantly different in structure from those seen in conventional polymer solutions. Our coarse-grained simulations with an implicit cross-linker potential show that networks of branching bundles, knotted morphologies, and structural chirality can be generated by a generalized approach independent of specific cross-linkers. Network structure depends primarily on filament flexibility and separation, with significant connectivity increase after percolation. Results should guide the design of engineered semiflexible polymers.

Direct evidence for structural transition promoting shear thinning in cylindrical colloid assemblies

In this paper, we describe stimuli-responsive hydrogels prepared from a rigid rod-like polyelectrolyte 'imogolite' and a dicarboxylic acid. The hydrogel exhibited thixotropy in response to mechanical shock within the order of seconds or sub-seconds. Here, using the latest structural/rheological characterisation techniques, the relationship between the structural transition processes and the shear thinning was estimated. The evidence obtained by the experiments revealed for the first time the direct relationship between the microscopic structural change and the macroscopic thixotropic behavior that have been extensively discussed. The thixotropic hydrogel has the hierarchical architecture in the combination of imogolite and dicarboxylic acid, i.e., sheathed nanotubes/hydroclusters of cross-bridged nanotubes/frameworks. The formation and disintegration of the network structure upon resting and agitating, respectively, were the origin of gel/sol transition (thixotropy), although the hydroclusters of cross-bridged nanotubes were maintained throughout the transition.

Alignment of nematic and bundled semiflexible polymers in cell-sized confinement

The finite size of cells poses severe spatial constraints on the network of semiflexible filaments called the cytoskeleton, a main determinant of cell shape. At the same time, the high packing density of cytoskeletal filaments poses mutual packing constraints. Here we investigate the competition between excluded volume interactions in the bulk and surface packing constraints on the orientational ordering of confined actin filaments as a function of filament density and the presence of crosslinks. We grow fluorescently labeled actin filaments in shallow (thickness dz 3 μm), rectangular microchambers with a systematically varied length (dy between 5 and 100 μm) and in-plane aspect ratio (dx/dy between 1 and 10). We determine the nematic director field by image analysis of fluorescence confocal images. We find that high-density (nematic) solutions respond sensitively to changes in the size and aspect ratio of the chambers. In small chambers (dy ≤ 20 μm), filaments align parallel to the long walls as soon as the aspect ratio is ≥1.5, indicating that surface-induced ordering dominates. In larger chambers, the filaments instead align along the chamber diagonal, indicating that bulk packing constraints dominate. The nematic order parameter is maximal in small and highly anisometric chambers. In contrast to the nematic solutions, low-density (isotropic) solutions are rather insensitive to confinement. Bundled actin solutions behave similarly to nematic solutions, but are less well-ordered. Our observations imply that the orientational order of actin filaments in flat confining geometries is primarily determined by a balance between bulk and surface packing constraints with a minimal effect of the enthalpic cost of filament bending. Our assay provides an interesting platform for the future reconstitution of more complex, active cytoskeletal systems with actively treadmilling filaments or molecular motors.

Properties of intermediate filament networks assembled from keratin 8 and 18 in the presence of Mg²+

The mechanical properties of epithelial cells are modulated by structural changes in keratin intermediate filament networks. To investigate the relationship between network architecture and viscoelasticity, we assembled keratin filaments from recombinant keratin proteins 8 (K8) and 18 (K18) in the presence of divalent ions (Mg(2+)). We probed the viscoelastic modulus of the network by tracking the movement of microspheres embedded in the network during assembly, and studied the network architecture using scanning electron microscopy. Addition of Mg(2+) at physiological concentrations (<1 mM) resulted in networks whose structure was similar to that of keratin networks in epithelial cells. Moreover, the elastic moduli of networks assembled in vitro were found to be within the same magnitude as those measured in keratin networks of detergent-extracted epithelial cells. These findings suggest that Mg(2+)-induced filament cross-linking represents a valid model for studying the cytoskeletal mechanics of keratin networks.

Assembly kinetics determine the architecture of α-actinin crosslinked F-actin networks

The actin cytoskeleton is organized into diverse meshworks and bundles that support many aspects of cell physiology. Understanding the self-assembly of these actin-based structures is essential for developing predictive models of cytoskeletal organization. Here we show that the competing kinetics of bundle formation with the onset of dynamic arrest arising from filament entanglements and cross-linking determine the architecture of reconstituted actin networks formed with α-actinin cross-links. Cross-link mediated bundle formation only occurs in dilute solutions of highly mobile actin filaments. As actin polymerization proceeds, filament mobility and bundle formation are arrested concomitantly. By controlling the onset of dynamic arrest, perturbations to actin assembly kinetics dramatically alter the architecture of biochemically identical samples. Thus, the morphology of reconstituted F-actin networks is a kinetically determined structure similar to those formed by physical gels and glasses. These results establish mechanisms controlling the structure and mechanics in diverse semi-flexible biopolymer networks.

Supramolecular cellular filament systems: how and why do they form?

All cells, from simple bacteria to complex human tissues, rely on extensive networks of protein fibers to help maintain their proper form and function. These filament systems usually do not operate as single filaments, but form complex suprastructures, which are essential for specific cellular functions. Here, we describe the progress in determining the architectures of molecular filamentous suprastructures, the principles leading to their formation, and the mechanisms by which they may facilitate function. The complex eukaryotic cytoskeleton is tightly regulated by a large number of actin- or microtubule-associated proteins. In contrast, recently discovered bacterial actins and tubulins have few associated regulatory proteins. Hence, the quest to find basic principles that govern the formation of filamentous suprastructures is simplified in bacteria. Three common principles, which have been probed extensively during evolution, can be identified that lead to suprastructures formation: cationic counterion fluctuations; self-association into liquid crystals; and molecular crowding. The underlying physics of these processes will be discussed with respect to physiological circumstance.

Cooperative crosslink (un)binding in slowly driven bundles of semiflexible filaments

Combining simulations and theory I study the interplay between bundle elastic degrees of freedom and crosslink binding propensity. By slowly driving bundles into a deformed configuration, and depending on the mechanical stiffness of the crosslinking agent, the binding affinity is shown to display a sudden and discontinuous drop. This indicates a cooperative unbinding process that involves the crossing of a free-energy barrier. Choosing the proper crosslinker therefore not only allows us to change the composite elastic properties of the bundle but also the relevant time scales which can be tuned from the single crosslink binding rate to the much longer escape time over the free-energy barrier.

Polymorphism of highly cross-linked F-actin networks: probing multiple length scales

The assembly properties of F-actin filaments in the presence of different biological cross-linker concentrations and types have been investigated using a combined approach of fluorescence confocal microscopy and coarse-grained molecular dynamics simulation. In particular for highly cross-linked regimes, new network morphologies are observed. Complex network formation and the details of the resulting structure are strongly dependent on the ratio of cross-linkers to actin monomers and cross-linker shape but only weakly dependent on overall actin concentration and filament length. The work presented here may help to provide some fundamental understanding of how excessive cross-linkers interact with the actin filament solution, creating different structures in the cell under high cross-linker concentrations. F-actin is not only of biological importance but also, as an example of a semiflexible polymer, has attracted significant interest in its physical behavior. In combination with different cross-linkers semiflexible filaments may provide new routes to bio-materials development and act as the inspiration for new hierarchical network-based materials.

Self-organization of muscle cell structure and function

The organization of muscle is the product of functional adaptation over several length scales spanning from the sarcomere to the muscle bundle. One possible strategy for solving this multiscale coupling problem is to physically constrain the muscle cells in microenvironments that potentiate the organization of their intracellular space. We hypothesized that boundary conditions in the extracellular space potentiate the organization of cytoskeletal scaffolds for directed sarcomeregenesis. We developed a quantitative model of how the cytoskeleton of neonatal rat ventricular myocytes organizes with respect to geometric cues in the extracellular matrix. Numerical results and in vitro assays to control myocyte shape indicated that distinct cytoskeletal architectures arise from two temporally-ordered, organizational processes: the interaction between actin fibers, premyofibrils and focal adhesions, as well as cooperative alignment and parallel bundling of nascent myofibrils. Our results suggest that a hierarchy of mechanisms regulate the self-organization of the contractile cytoskeleton and that a positive feedback loop is responsible for initiating the break in symmetry, potentiated by extracellular boundary conditions, is required to polarize the contractile cytoskeleton.

How muscle is organized impacts its function. However, understanding how muscle organizes is challenging, as the process occurs over several length scales. We approach this multiscale coupling problem by constraining the overall shapes of muscle cells to indirectly control the organization of their intracellular space. We hypothesized the cellular boundary conditions direct the organization of cytoskeletal scaffolds. We developed a model of how the cytoskeleton of cardiomyocytes organizes with respect to boundary cues. Our computational and experimental results to control myocyte shape indicated that distinct muscle architectures arise from two main organizational mechanisms: the interaction between actin fibers, premyofibrils and focal adhesions, as well as cooperative alignment and parallel bundling of more mature myofibrils. We show that a hierarchy of processes regulate the self-organization of cardiomyocytes. Our results suggest that a symmetry break, due to the boundary conditions imposed on the cell, is responsible for polarization of the contractile cytoskeletal organization.

Orientational order and glassy states in networks of semiflexible polymers

Motivated by the structure of networks of cross-linked cytoskeletal biopolymers, we study orientationally ordered phases in two-dimensional networks of randomly cross-linked semiflexible polymers. We consider permanent cross-links which prescribe a finite angle and treat them as quenched disorder in a semimicroscopic replica field theory. Starting from a fluid of un-cross-linked polymers and small polymer clusters (sol) and increasing the cross-link density, a continuous gelation transition occurs. In the resulting gel, the semiflexible chains either display long-range orientational order or are frozen in random directions depending on the value of the crossing angle, the cross-link concentration, and the stiffness of the polymers. A crossing angle θ~2π/M leads to long-range M-fold orientational order, for example, "hexatic" or "tetratic" for θ=60° or 90°, respectively. The transition to the orientationally ordered state is discontinuous and the critical cross-link density, which is higher than that of the gelation transition, depends on the bending stiffness of the polymers and the cross-link angle: The higher the stiffness and the lower the M, the lower is the critical number of cross-links. In between the sol and the long-range ordered state, we always expect a gel which is a statistically isotropic amorphous solid with random positional and random orientational localization of the participating polymers.

Structural reorganization of parallel actin bundles by crosslinking proteins: incommensurate states of twist

We construct a coarse-grained model of parallel actin bundles crosslinked by compact globular bundling proteins, such as fascin and espin, necessary components of filopodial and mechanosensory bundles. Consistent with structural observations of bundles, we find that the optimal geometry for crosslinking is overtwisted, requiring a coherent structural change of the helical geometry of the filaments. We study the linker-dependent thermodynamic transition of bundled actin filaments from their native state to the overtwisted state and map out the "twist-state" phase diagram in terms of the availability as well as the flexibility of crosslinker proteins. We predict that the transition from the uncrosslinked to fully crosslinked state is highly sensitive to linker flexibility: flexible crosslinking smoothly distorts the twist state of bundled filaments, while rigidly crosslinked bundles undergo a phase transition, rapidly overtwisting filaments over a narrow range of free crosslinker concentrations. Additionally, we predict a rich spectrum of intermediate structures, composed of alternating domains of sparsely bound (untwisted) and strongly bound (overtwisted) filaments. This model reveals that subtle differences in crosslinking agents themselves modify not only the detailed structure of parallel actin bundles, but also the thermodynamic pathway by which they form.

Electrostatics of strongly charged biological polymers: ion-mediated interactions and self-organization in nucleic acids and proteins

Charges on biological polymers in physiologically relevant solution conditions are strongly screened by water and salt solutions containing counter-ions. However, the entropy of these counterions can result in surprisingly strong interactions between charged objects in water despite short screening lengths, via coupling between osmotic and electrostatic interactions. Widespread work in theory, experiment, and computation has been carried out to gain a fundamental understanding of the rich, yet sometimes counterintuitive, behavior of these polyelectrolyte systems. Examples of polyelectrolyte association in biology include DNA packaging and RNA folding, as well as aggregation and self-organization phenomena in different disease states.

Suprastructures and dynamic properties of Mycobacterium tuberculosis FtsZ

Tuberculosis causes the most death in humans by any bacterium. Drug targeting of bacterial cytoskeletal proteins requires detailed knowledge of the various filamentous suprastructures and dynamic properties. Here, we have investigated by high resolution electron microscopy the assembly of cell division protein and microtubule homolog FtsZ from Mycobacterium tuberculosis (MtbFtsZ) in vitro in the presence of various monovalent salts, crowding agents and polycations. Supramolecular structures, including two-dimensional rings, three-dimensional toroids, and multistranded helices formed in the presence of molecular crowding, were similar to those observed by fluorescence microscopy in bacteria in vivo. Dynamic properties of MtbFtsZ filaments were visualized by light scattering and real time total internal reflection fluorescence microscopy. Interestingly, MtbFtsZ revealed a form of dynamic instability at steady state. Cation-induced condensation phenomena of bacterial cytomotive polymers have not been investigated in any detail, although it is known that many bacteria can contain high amounts of polycations, which may modulate the prokaryotic cytoskeleton. We find that above a threshold concentration of polycations which varied with the valence of the cation, ionic strength, and pH, MtbFtsZ mainly formed sheets. The general features of these cation-induced condensation phenomena could be explained in the framework of the Manning condensation theory. Chirality and packing defects limited the dimensions of sheets and toroids at steady state as predicted by theoretical models. In first approximation simple physical principles seem to govern the formation of MtbFtsZ suprastructures.

From branched networks of actin filaments to bundles

Cross-linking proteins can mediate the emergence of rigid bundles from a dense branched network of actin filaments. To enable their binding, the filaments must first bend towards each other. We derive an explicit criterion for the onset of bundling, in terms of the initial length of filaments L, their spacing b, and cross-linker concentration f, reflecting the balance between bending and binding energies. Our model system contains actin, the branching complex Arp2/3 and the bundling protein fascin. In the first distinct stage, during which only actin and Arp2/3 are active, an entangled aster-like mesh of actin filaments is formed. Tens of seconds later, when filaments at the aster periphery are long and barely branched, a sharp transition takes place into a star-like structure, marking the onset of bundling. Now fascin and actin govern bundle growth; Arp2/3 plays no role. Using kinetic Monte Carlo simulations we calculate the temporal evolution of b and L, and predict the onset of bundling as a function of f. Our predictions are in good qualitative agreement with several new experiments that are reported herein and demonstrate how f controls the aster-star transition and bundle length. We also present two models for aster growth corresponding to different experimental realizations. The first treats filament and bundle association as an irreversible sequence of elongation-association steps. The second, applicable for low f, treats bundling as a reversible self-assembly process, where the optimal bundle size is dictated by the balance between surface and bending energies. Finally, we discuss the relevance of our conclusions for the lamellipodium to filopodia transition in living cells, noting that bundles are more likely nucleated by "tip complex" cross-linkers (e.g. mDia2 or Ena/VASP), whereas fascin is mainly involved in bundle maintenance.

Molecular mechanism of bundle formation by the bacterial actin ParM

The actin homolog ParM plays a microtubule-like role in segregating DNA prior to bacterial cell division. Fluorescence and cryo-electron microscopy have shown that ParM forms filament bundles between separating DNA plasmids in vivo. Given the lack of ParM bundling proteins it remains unknown how ParM bundles form at the molecular level. Here we show using time-lapse TIRF microscopy, under in vitro molecular crowding conditions, that ParM-bundle formation consists of two distinct phases. At the onset of polymerization bundle thickness and shape are determined in the form of nuclei of short helically disordered filaments arranged in a liquid-like lattice. These nuclei then undergo an elongation phase whereby they rapidly increase in length. At steady state, ParM bundles fuse into one single large aggregate. This behavior had been predicted by theory but has not been observed for any other cytomotive biopolymer, including F-actin. We employed electron micrographs of ParM rafts, which are 2-D analogs of 3-D bundles, to identify the main molecular interfilament contacts within these suprastructures. The interface between filaments is similar for both parallel and anti-parallel orientations and the distribution of filament polarity is random within a bundle. We suggest that the interfilament interactions are not due to the interactions of specific residues but rather to long-range, counter ion mediated, electrostatic attractive forces. A randomly oriented bundle ensures that the assembly is rigid and that DNA may be captured with equal efficiency at both ends of the bundle via the ParR binding protein.

Salt-induced aggregation of stiff polyelectrolytes

Molecular dynamics simulation techniques are used to study the process of aggregation of highly charged stiff polyelectrolytes due to the presence of multivalent salt. The dominant kinetic mode of aggregation is found to be the case of one end of one polyelectrolyte meeting others at right angles, and the kinetic pathway to bundle formation is found to be similar to that of flocculation dynamics of colloids as described by Smoluchowski. The aggregation process is found to favor the formation of finite bundles of 10-11 filaments at long times. Comparing the distribution of the cluster sizes with the Smoluchowski formula suggests that the energy barrier for the aggregation process is negligible. Also, the formation of long-lived metastable structures with similarities to the raft-like structures of actin filaments is observed within a range of salt concentration.

Controllable side-by-side and end-to-end assembly of Au nanorods by lyotropic chromonic materials

We present a simple and universal technique for assembling gold nanorods (NRs) using self-assembled stacks of lyotropic chromonic materials, without covalent bonding between NRs and the linking agent. The anisotropic electrostatic interaction between the chromonic stacks and NRs allows one to achieve either side-by-side or end-to-end assembly, depending on the surface charge of NRs. The assembled superstructures are stable within an extended temperature range; the degree of NR aggregation can be controlled by a number of factors influencing the self-assembly of chromonic materials, such as the concentration and pH of the solution.

Energetics and dynamics of global integrals modeling interaction between stiff filaments

The attractive and spacing interaction between pairs of filaments via cross-linkers, e.g. myosin oligomers connecting actin filaments, is modeled by global integral kernels for negative binding energies between two intersecting stiff and long rods in a (projected) two-dimensional situation, for simplicity. Whereas maxima of the global energy functional represent intersection angles of 'minimal contact' between the filaments, minima are approached for energy values tending to -infinity, representing the two degenerate states of parallel and anti-parallel filament alignment. Standard differential equations of negative gradient flow for such energy functionals show convergence of solutions to one of these degenerate equilibria in finite time, thus called 'super-stable' states. By considering energy variations under virtual rotation or translation of one filament with respect to the other, integral kernels for the resulting local forces parallel and orthogonal to the filament are obtained. For the special modeling situation that these variations only activate 'spring forces' in direction of the cross-links, explicit formulas for total torque and translational forces are given and calculated for typical examples. Again, the two degenerate alignment states are locally 'super-stable' equilibria of the assumed over-damped dynamics, but also other stable states of orthogonal arrangement and different asymptotic behavior can occur. These phenomena become apparent if stochastic perturbations of the local force kernels are implemented as additive Gaussian noise induced by the cross-link binding process with appropriate scaling. Then global filament dynamics is described by a certain type of degenerate stochastic differential equations yielding asymptotic stationary processes around the alignment states, which have generalized, namely bimodal Gaussian distributions. Moreover, stochastic simulations reveal characteristic sliding behavior as it is observed for myosin-mediated interaction between actin filaments. Finally, the forgoing explicit and asymptotic analysis as well as numerical simulations are extended to the more realistic modeling situation with filaments of finite length.

Packing defects and the width of biopolymer bundles

The formation of bundles composed of actin filaments and cross-linking proteins is an essential process in the maintenance of the cells' cytoskeleton. It has also been recreated by in-vitro experiments, where actin networks are routinely produced to mimic and study the cellular structures. It has been observed that these bundles seem to have a well-defined width distribution, which has not been adequately described theoretically. We propose here that packing defects of the filaments, quenched and random, contribute an effective repulsion that counters the cross-linking adhesion energy and leads to a well-defined bundle width. This is a two-dimensional strain-field version of the classic Rayleigh instability of charged droplets.

Thickness distribution of actin bundles in vitro

Bundles of filamentous actin form the primary building blocks of a broad range of cytoskeletal structures, including filopodia, stereocilia and microvilli. In each case, the cell uses specific associated proteins to tailor the dynamics, dimensions and mechanical properties of the bundles to suit a specific cellular function. While the length distribution of actin bundles was extensively studied, almost nothing is known about the thickness distribution. Here, we use high-resolution cryo-TEM to measure the thickness distribution of actin/fascin bundles, in vitro. We find that the thickness distribution has a prominent peak, with an exponential tail, supporting a scenario of an initial fast formation of a disc-like nucleus of short actin filaments, which only later elongates. The bundle thicknesses at steady state are found to follow the distribution of the initial nuclei indicating that no lateral coalescence occurs. Our results show that the distribution of bundles thicknesses can be controlled by monitoring the initial nucleation process. In vivo, this is done by using specific regulatory proteins complexes.

Anisotropic random networks of semiflexible polymers

Motivated by the organization of cross-linked cytoskeletal biopolymers, we present a semimicroscopic replica field theory for the formation of anisotropic random networks of semiflexible polymers. The networks are formed by introducing random permanent cross-links which fix the orientations of the corresponding polymer segments to align with one another. Upon increasing the cross-link density, we obtain a continuous gelation transition from a fluid phase to a gel where a finite fraction of the system gets localized at random positions. For sufficiently stiff polymers, this positional localization is accompanied by a continuous isotropic-to-nematic (IN) transition occurring at the same cross-link density. As the polymer stiffness decreases, the IN transition becomes first order, shifts to a higher cross-link density, and is preceded by an amorphous solid where the average polymer orientations freeze in random directions.

Aggregation kinetics of stiff polyelectrolytes in the presence of multivalent salt

Using molecular dynamics simulations, the kinetics of bundle formation for stiff polyelectrolytes such as actin is studied in the solution of multivalent salt. The dominant kinetic mode of aggregation is found to be the case of one end of one rod meeting others at a right angle due to electrostatic interactions. The kinetic pathway to bundle formation involves a hierarchical structure of small clusters forming initially and then feeding into larger clusters, which is reminiscent of the flocculation dynamics of colloids. For the first few cluster sizes, the Smoluchowski formula for the time evolution of the cluster size gives a reasonable account of the results of our simulation without a single fitting parameter. The description using the Smoluchowski formula provides evidence for the aggregation time scale to be controlled by diffusion, with no appreciable energy barrier to overcome.

Chirality and equilibrium biopolymer bundles

We use continuum theory to show that chirality is a key thermodynamic control parameter for the aggregation of biopolymers: chirality produces a stable disperse phase of hexagonal bundles under moderately poor solvent conditions, as has been observed in in vitro studies of F actin [O. Pelletier et al., Phys. Rev. Lett. 91, 148102 (2003)]. The large characteristic radius of these chiral bundles is not determined by a mysterious long-range molecular interaction but by in-plane shear elastic stresses generated by the interplay between a chiral torque and an unusual, but universal, nonlinear gauge term in the strain tensor of ordered chains that is imposed by rotational invariance.

Polyamine-induced bundling of F-actin

To better understand the mechanism of actin filament (F-actin) bundling by polyamines, we have measured the onset of bundling as a function of polyamine concentration. Samples were centrifuged at low speeds to separate bundles from unbundled actin, and the relative amounts of actin in the pellet and supernatant were determined via gel electrophoresis, yielding a description of the bundling transition as a function of actin and polyamine concentrations. These experiments were carried out for two different polyamines, spermine (tetravalent) and spermidine (trivalent). We found that the threshold concentration of polyamine needed to bundle actin is independent of both actin concentration and Mg2+ concentration over a wide range in Mg2+ concentration. We also find that spermine in F-actin bundles is essentially invisible in solution-phase proton NMR, suggesting that it is bound so tightly to F-actin that it is immobilized.

Density-functional theory for complex fluids

Density-functional theory (DFT) and its variations provide a fruitful approach to the computational modeling of the microscopic structures and phase behavior of soft-condensed matter. The methodology takes deep root in quantum mechanics but shares a mathematical similarity with a number of classical approaches in statistical mechanics. This review discusses different strategies commonly used to formulate the free-energy functional of complex fluids for either phenomena-oriented applications or as a generic description of the thermodynamic nonideality owing to various components of intermolecular forces. We emphasize the connections among different schemes of DFT approximations, their underlying assumptions, and inherent limitations. We also address extensions of equilibrium DFT to phenomenological theories for the dynamic properties of complex fluids and for the kinetics of phase transitions. In addition, we highlight the recent literature concerning applications of DFT to diverse static and time-dependent phenomena in complex fluids.

Simulation and theory of self-assembly and network formation in reversibly cross-linked equilibrium polymers

A simulation model of hard spheres capable of reversible assembly into chains, which then may reversibly cross-link into networks, has been studied through grand canonical Monte Carlo simulation. Effects of varying intra- and interchain bond strengths, chain flexibilities, and restrictions on cross-linking angle were investigated. Observations including chain-length distributions and phase separation could be captured in most cases using a simple model theory. The coupling of chain growth to cross-linking was shown to be highly sensitive to the treatment of cross-linking by chain ends. In some systems, ladderlike domains of several cross-links joining two chains were common, resulting from cooperativity in the cross-linking. Extended to account for this phenomenon, the model theory predicts that such cooperativity will suppress phase separation in weakly polymerizing chains and at high cross-link concentration. In the present model, cross-linking stabilizes the isotropic phase with respect to the nematic phase, causing a shift in the isotropic-nematic transition to higher monomer concentration than in simple equilibrium polymers.

Lattice model for self-assembly with application to the formation of cytoskeletal-like structures

We introduce a stochastic approach for self-assembly in systems far from equilibrium. The building blocks are represented by a lattice of discrete variables (Potts-like spins), and physically meaningful mechanisms are obtained by restricting transitions through spatially local rules based on experimental data. We use the method to study nucleation of filopodia-like bundles in a system consisting of purified actin, fascin, actin-related protein 2/3 , and beads coated with Wiskott-Aldrich syndrome protein. Consistent with previous speculation based on static experimental images, we find that bundles derive from Lambda-precursor-like patterns of spins on the lattice. The ratcheting of the actin network relative to the surface that represents beads plays an important role in determining the number and orientation of bundles due to the fact that branching is the primary means for generating barbed ends pointed in directions that allow rapid filament growth. By enabling the de novo formation of coexisting morphologies without the computational cost of explicit representation of proteins, the approach introduced complements earlier models of cytoskeletal behavior in vitro and in vivo.