Bridging scales in a multiscale pattern-forming system

Microscopic interactions control a structural transition in active mixtures of microtubules and molecular motors

Microtubules and molecular motors are essential components of the cellular cytoskeleton, driving fundamental processes in vivo, including chromosome segregation and cargo transport. When reconstituted in vitro, these cytoskeletal proteins serve as energy-consuming building blocks to study the self-organization of active matter. Cytoskeletal active gels display rich emergent dynamics, including extensile flows, locally contractile asters, and bulk contraction. However, it is unclear how the protein-protein interaction kinetics set their contractile or extensile nature. Here, we explore the origin of the transition from extensile bundles to contractile asters in a minimal reconstituted system composed of stabilized microtubules, depletant, adenosine 5'-triphosphate (ATP), and clusters of kinesin-1 motors. We show that the microtubule-binding and unbinding kinetics of highly processive motor clusters set their ability to end-accumulate, which can drive polarity sorting of the microtubules and aster formation. We further demonstrate that the microscopic time scale of end-accumulation sets the emergent time scale of aster formation. Finally, we show that biochemical regulation is insufficient to fully explain the transition as generic aligning interactions through depletion, cross-linking, or excluded volume interactions can drive bundle formation despite end-accumulating motors. The extensile-to-contractile transition is well captured by a simple self-assembly model where nematic and polar aligning interactions compete to form either bundles or asters. Starting from a five-dimensional organization phase space, we identify a single control parameter given by the ratio of the different component concentrations that dictates the material-scale organization. Overall, this work shows that the interplay of biochemical and mechanical tuning at the microscopic level controls the robust self-organization of active cytoskeletal materials.

Forceful patterning: theoretical principles of mechanochemical pattern formation

Biological pattern formation is essential for generating and maintaining spatial structures from the scale of a single cell to tissues and even collections of organisms. Besides biochemical interactions, there is an important role for mechanical and geometrical features in the generation of patterns. We review the theoretical principles underlying different types of mechanochemical pattern formation across spatial scales and levels of biological organization.

Coarsening and wavelength selection far from equilibrium: A unifying framework based on singular perturbation theory

Intracellular protein patterns are described by (nearly) mass-conserving reaction-diffusion systems. While these patterns initially form out of a homogeneous steady state due to the well-understood Turing instability, no general theory exists for the dynamics of fully nonlinear patterns. We develop a unifying theory for nonlinear wavelength-selection dynamics in (nearly) mass-conserving two-component reaction-diffusion systems independent of the specific mathematical model chosen. Previous work has shown that these systems support an extremely broad band of stable wavelengths, but the mechanism by which a specific wavelength is selected has remained unclear. We show that an interrupted coarsening process selects the wavelength at the threshold to stability. Based on the physical intuition that coarsening is driven by competition for mass and interrupted by weak source terms that break strict mass conservation, we develop a singular perturbation theory for the stability of stationary patterns. The resulting closed-form analytical expressions enable us to quantitatively predict the coarsening dynamics and the final pattern wavelength. We find excellent agreement with numerical results throughout the diffusion- and reaction-limited regimes of the dynamics, including the crossover region. Further, we show how, in these limits, the two-component reaction-diffusion systems map to generalized Cahn-Hilliard and conserved Allen-Cahn dynamics, therefore providing a link to these two fundamental scalar field theories. The systematic understanding of the length-scale dynamics of fully nonlinear patterns in two-component systems provided here builds the basis to reveal the mechanisms underlying wavelength selection in multicomponent systems with potentially several conservation laws.

Protein-Based Patterning to Spatially Functionalize Biomimetic Membranes

The bottom-up reconstitution of proteins for their modular engineering into synthetic cellular systems can reveal hidden protein functions in vitro. This is particularly evident for the bacterial Min proteins, a paradigm for self-organizing reaction-diffusion systems that displays an unexpected functionality of potential interest for bioengineering: the directional active transport of any diffusible cargo molecule on membranes. Here, the MinDE protein system is reported as a versatile surface patterning tool for the rational design of synthetically assembled 3D systems. Employing two-photon lithography, microswimmer-like structures coated with tailored lipid bilayers are fabricated and demonstrate that Min proteins can uniformly pattern bioactive molecules on their surface. Moreover, it is shown that the MinDE system can form stationary patterns inside lipid vesicles, which allow the targeting and distinctive clustering of higher-order protein structures on their inner leaflet. Given their facile use and robust function, Min proteins thus constitute a valuable molecular toolkit for spatially patterned functionalization of artificial biosystems like cell mimics and microcarriers.

Pattern selection by material aging: Modeling chemical gardens in two and three dimensions

Chemical gardens are complex, often macroscopic, structures formed by precipitation reactions. Their thin walls compartmentalize the system and adjust in size and shape if the volume of the interior reactant solution is increased by osmosis or active injection. Spatial confinement to a thin layer is known to result in various patterns including self-extending filaments and flower-like patterns organized around a continuous, expanding front. Here, we describe a cellular automaton model for this type of self-organization, in which each lattice site is occupied by one of the two reactants or the precipitate. Reactant injection causes the random replacement of precipitate and generates an expanding near-circular precipitate front. If this process includes an age bias favoring the replacement of fresh precipitate, thin-walled filaments arise and grow-like in the experiments-at the leading tip. In addition, the inclusion of a buoyancy effect allows the model to capture various branched and unbranched chemical garden shapes in two and three dimensions. Our results provide a model of chemical garden structures and highlight the importance of temporal changes in the self-healing membrane material.

Geometry-induced patterns through mechanochemical coupling

Intracellular protein patterns regulate a variety of vital cellular processes such as cell division and motility, which often involve dynamic cell-shape changes. These changes in cell shape may in turn affect the dynamics of pattern-forming proteins, hence leading to an intricate feedback loop between cell shape and chemical dynamics. While several computational studies have examined the rich resulting dynamics, the underlying mechanisms are not yet fully understood. To elucidate some of these mechanisms, we explore a conceptual model for cell polarity on a dynamic one-dimensional manifold. Using concepts from differential geometry, we derive the equations governing mass-conserving reaction-diffusion systems on time-evolving manifolds. Analyzing these equations mathematically, we show that dynamic shape changes of the membrane can induce pattern-forming instabilities in parts of the membrane, which we refer to as regional instabilities. Deformations of the local membrane geometry can also (regionally) suppress pattern formation and spatially shift already existing patterns. We explain our findings by applying and generalizing the local equilibria theory of mass-conserving reaction-diffusion systems. This allows us to determine a simple onset criterion for geometry-induced pattern-forming instabilities, which is linked to the phase-space structure of the reaction-diffusion system. The feedback loop between membrane shape deformations and reaction-diffusion dynamics then leads to a surprisingly rich phenomenology of patterns, including oscillations, traveling waves, and standing waves, even if these patterns do not occur in systems with a fixed membrane shape. Our paper reveals that the local conformation of the membrane geometry acts as an important dynamical control parameter for pattern formation in mass-conserving reaction-diffusion systems.

Modelling network formation in folded protein hydrogels by cluster aggregation kinetics

Globular folded protein-based hydrogels are becoming increasingly attractive due to their specific biological functionality, as well as their responsiveness to stimuli. By modelling folded proteins as colloids, there are rich opportunities to explore network formation mechanisms in protein hydrogels that negate the need for computationally expensive simulations which capture the full complexity of proteins. Here we present a kinetic lattice-based model which simulates the formation of irreversibly chemically crosslinked, folded protein-based hydrogels. We identify the critical point of gel percolation, explore the range of network regimes covering diffusion-limited to reaction-limited cluster aggregation (DLCA and RLCA, respectively) network formation mechanisms and predict the final network structure, fractal dimensions and final gel porosity. We reveal a crossover between DLCA and RLCA mechanisms as a function of protein volume fraction and show how the final network structure is governed by the structure at the percolation point, regardless of the broad variation of non-percolating cluster masses observed across all systems. An analysis of the pore size distribution in the final network structures reveals that, approaching RLCA, gels have larger maximal pores than the DLCA counterparts for both volume fractions studied. This general kinetic model and the analysis tools generate predictions of network structure and concurrent porosity over a broad range of experimentally controllable parameters that are consistent with current expectations and understanding of experimental results.

Universal dynamics of biological pattern formation in spatio-temporal morphogen variations

In biological systems, chemical signals termed morphogens self-organize into patterns that are vital for many physiological processes. As observed by Turing in 1952, these patterns are in a state of continual development, and are usually transitioning from one pattern into another. How do cells robustly decode these spatio-temporal patterns into signals in the presence of confounding effects caused by unpredictable or heterogeneous environments? Here, we answer this question by developing a general theory of pattern formation in spatio-temporal variations of ‘pre-pattern’ morphogens, which determine gene-regulatory network parameters. Through mathematical analysis, we identify universal dynamical regimes that apply to wide classes of biological systems. We apply our theory to two paradigmatic pattern-forming systems, and predict that they are robust with respect to non-physiological morphogen variations. More broadly, our theoretical framework provides a general approach to classify the emergent dynamics of pattern-forming systems based on how the bifurcations in their governing equations are traversed.

Directing Min protein patterns with advective bulk flow

The Min proteins constitute the best-studied model system for pattern formation in cell biology. We theoretically predict and experimentally show that the propagation direction of in vitro Min protein patterns can be controlled by a hydrodynamic flow of the bulk solution. We find downstream propagation of Min wave patterns for low MinE:MinD concentration ratios, upstream propagation for large ratios, but multistability of both propagation directions in between. Whereas downstream propagation can be described by a minimal model that disregards MinE conformational switching, upstream propagation can be reproduced by a reduced switch model, where increased MinD bulk concentrations on the upstream side promote protein attachment. Our study demonstrates that a differential flow, where bulk flow advects protein concentrations in the bulk, but not on the surface, can control surface-pattern propagation. This suggests that flow can be used to probe molecular features and to constrain mathematical models for pattern-forming systems.

Can One Series of Self-Organized Nanoripples Guide Another Series of Self-Organized Nanoripples during Ion Bombardment: From the Perspective of Power Spectral Density Entropy?

Ion bombardment (IB) is a promising nanofabrication tool for self-organized nanostructures. When ions bombard a nominally flat solid surface, self-organized nanoripples can be induced on the irradiated target surface, which are called intrinsic nanoripples of the target material. The degree of ordering of nanoripples is an outstanding issue to be overcome, similar to other self-organization methods. In this study, the IB-induced nanoripples on bilayer systems with enhanced quality are revisited from the perspective of guided self-organization. First, power spectral density (PSD) entropy is introduced to evaluate the degree of ordering of the irradiated nanoripples, which is calculated based on the PSD curve of an atomic force microscopy image (i.e., the Fourier transform of the surface height. The PSD entropy can characterize the degree of ordering of nanoripples). The lower the PSD entropy of the nanoripples is, the higher the degree of ordering of the nanoripples. Second, to deepen the understanding of the enhanced quality of nanoripples on bilayer systems, the temporal evolution of the nanoripples on the photoresist (PR)/antireflection coating (ARC) and Au/ARC bilayer systems are compared with those of single PR and ARC layers. Finally, we demonstrate that a series of intrinsic IB-induced nanoripples on the top layer may act as a kind of self-organized template to guide the development of another series of latent IB-induced nanoripples on the underlying layer, aiming at improving the ripple ordering. The template with a self-organized nanostructure may alleviate the critical requirement for periodic templates with a small period of ~100 nm. The work may also provide inspiration for guided self-organization in other fields.