Alignment and actuation of liquid crystals via 3D confinement and two-photon laser printing.

Liquid crystalline (LC) materials are especially suited for the preparation of active three-dimensional (3D) and 4D microstructures using two-photon laser printing. To achieve the desired actuation, the alignment of the LCs has to be controlled during the printing process. In most cases studied before, the alignment relied on surface modifications and complex alignment patterns and concomitant actuation were not possible. Here, we introduce a strategy for spatially aligning LC domains in three-dimensional space by using 3D-printed polydimethylsiloxane-based microscaffolds as confinement barriers, which induce the desired director field. The director field resulting from the boundary conditions is calculated with Landau de Gennes theory and validated by comparing experimentally measured and theoretically predicted birefringence patterns. We demonstrate our procedures for structures of varying complexity and then employed them to fabricate 4D microstructures that show the desired actuation. Overall, we obtain excellent agreement between theory and experiment. This opens the door for rational design of functional materials for 4D (micro)printing in the future.

DNA microbeads for spatio-temporally controlled morphogen release within organoids.

Organoids are transformative in vitro model systems that mimic features of the corresponding tissue in vivo. However, across tissue types and species, organoids still often fail to reach full maturity and function because biochemical cues cannot be provided from within the organoid to guide their development. Here we introduce nanoengineered DNA microbeads with tissue mimetic tunable stiffness for implementing spatio-temporally controlled morphogen gradients inside of organoids at any point in their development. Using medaka retinal organoids and early embryos, we show that DNA microbeads can be integrated into embryos and organoids by microinjection and erased in a non-invasive manner with light. Coupling a recombinant surrogate Wnt to the DNA microbeads, we demonstrate the spatio-temporally controlled morphogen release from the microinjection site, which leads to morphogen gradients resulting in the formation of retinal pigmented epithelium while maintaining neuroretinal cell types. Thus, we bioengineered retinal organoids to more closely mirror the cell type diversity of in vivo retinae. Owing to the facile, one-pot fabrication process, the DNA microbead technology can be adapted to other organoid systems for improved tissue mimicry.

A geometrical theory of gliding motility based on cell shape and surface flow.

Gliding motility proceeds with little changes in cell shape and often results from actively driven surface flows of adhesins binding to the extracellular environment. It allows for fast movement over surfaces or through tissue, especially for the eukaryotic parasites from the phylum apicomplexa, which includes the causative agents of the widespread diseases malaria and toxoplasmosis. We have developed a fully three-dimensional active particle theory which connects the self-organized, actively driven surface flow over a fixed cell shape to the resulting global motility patterns. Our analytical solutions and numerical simulations show that straight motion without rotation is unstable for simple shapes and that straight cell shapes tend to lead to pure rotations. This suggests that the curved shapes of Plasmodium sporozoites and Toxoplasma tachyzoites are evolutionary adaptations to avoid rotations without translation. Gliding motility is also used by certain myxo- or flavobacteria, which predominantly move on flat external surfaces and with higher control of cell surface flow through internal tracks. We extend our theory for these cases. We again find a competition between rotation and translation and predict the effect of internal track geometry on overall forward speed. While specific mechanisms might vary across species, in general, our geometrical theory predicts and explains the rotational, circular, and helical trajectories which are commonly observed for microgliders. Our theory could also be used to design synthetic microgliders.

The role of the nucleus for cell mechanics: an elastic phase field approach.

The nucleus of eukaryotic cells typically makes up around 30% of the cell volume and has significantly different mechanics, which can make it effectively up to ten times stiffer than the surrounding cytoplasm. Therefore it is an important element for cell mechanics, but a quantitative understanding of its mechanical role during whole cell dynamics is largely missing. Here we demonstrate that elastic phase fields can be used to describe dynamical cell processes in adhesive or confining environments in which the nucleus acts as a stiff inclusion in an elastic cytoplasm. We first introduce and verify our computational method and then study several prevalent cell-mechanical measurement methods. For cells on adhesive patterns, we find that nuclear stress is shielded by the adhesive pattern. For cell compression between two parallel plates, we obtain force-compression curves that allow us to extract an effective modulus for the cell-nucleus composite. For micropipette aspiration, the effect of the nucleus on the effective modulus is found to be much weaker, highlighting the complicated interplay between extracellular geometry and cell mechanics that is captured by our approach. We also show that our phase field approach can be used to investigate the effects of Kelvin-Voigt-type viscoelasticity and cortical tension.

Modelling cell shape in 3D structured environments: A quantitative comparison with experiments.

Cell shape plays a fundamental role in many biological processes, including adhesion, migration, division and development, but it is not clear which shape model best predicts three-dimensional cell shape in structured environments. Here, we compare different modelling approaches with experimental data. The shapes of single mesenchymal cells cultured in custom-made 3D scaffolds were compared by a Fourier method with surfaces that minimize area under the given adhesion and volume constraints. For the minimized surface model, we found marked differences to the experimentally observed cell shapes, which necessitated the use of more advanced shape models. We used different variants of the cellular Potts model, which effectively includes both surface and bulk contributions. The simulations revealed that the Hamiltonian with linear area energy outperformed the elastic area constraint in accurately modelling the 3D shapes of cells in structured environments. Explicit modelling the nucleus did not improve the accuracy of the simulated cell shapes. Overall, our work identifies effective methods for accurately modelling cellular shapes in complex environments.

The role of flow in the self-assembly of dragline spider silk proteins.

Hydrodynamic flow in the spider duct induces conformational changes in dragline spider silk proteins (spidroins) and drives their assembly, but the underlying physical mechanisms are still elusive. Here we address this challenging multiscale problem with a complementary strategy of atomistic and coarse-grained molecular dynamics simulations with uniform flow. The conformational changes at the molecular level were analyzed for single-tethered spider silk peptides. Uniform flow leads to coiled-to-stretch transitions and pushes alanine residues into ß sheet and poly-proline II conformations. Coarse-grained simulations of the assembly process of multiple semi-flexible block copolymers using multi-particle collision dynamics reveal that the spidroins aggregate faster but into low-order assemblies when they are less extended. At medium-to-large peptide extensions (50%-80%), assembly slows down and becomes reversible with frequent association and dissociation events, whereas spidroin alignment increases and alanine repeats form ordered regions. Our work highlights the role of flow in guiding silk self-assembly into tough fibers by enhancing alignment and kinetic reversibility, a mechanism likely relevant also for other proteins whose function depends on hydrodynamic flow.

Membrane Tension Inhibits Lipid Mixing by Increasing the Hemifusion Stalk Energy.

Fusion of biological membranes is fundamental in various physiological events. The fusion process involves several intermediate stages with energy barriers that are tightly dependent on the mechanical and physical properties of the system, one of which is membrane tension. As previously established, the late stages of fusion, including hemifusion diaphragm and pore expansions, are favored by membrane tension. However, a current understanding of how the energy barrier of earlier fusion stages is affected by membrane tension is lacking. Here, we apply a newly developed experimental approach combining micropipette-aspirated giant unilamellar vesicles and optically trapped membrane-coated beads, revealing that membrane tension inhibits lipid mixing. We show that lipid mixing is 6 times slower under a tension of 0.12 mN/m compared with tension-free membranes. Furthermore, using continuum elastic theory, we calculate the dependence of the hemifusion stalk formation energy on membrane tension and intermembrane distance and find the increase in the corresponding energy barrier to be 1.6 kBT in our setting, which can explain the increase in lipid mixing time delay. Finally, we show that tension can be a significant factor in the stalk energy if the pre-fusion intermembrane distance is on the order of several nanometers, while for membranes that are tightly docked, tension has a negligible effect.

Force propagation between epithelial cells depends on active coupling and mechano-structural polarization.

Cell-generated forces play a major role in coordinating the large-scale behavior of cell assemblies, in particular during development, wound healing, and cancer. Mechanical signals propagate faster than biochemical signals, but can have similar effects, especially in epithelial tissues with strong cell-cell adhesion. However, a quantitative description of the transmission chain from force generation in a sender cell, force propagation across cell-cell boundaries, and the concomitant response of receiver cells is missing. For a quantitative analysis of this important situation, here we propose a minimal model system of two epithelial cells on an H-pattern ('cell doublet'). After optogenetically activating RhoA, a major regulator of cell contractility, in the sender cell, we measure the mechanical response of the receiver cell by traction force and monolayer stress microscopies. In general, we find that the receiver cells show an active response so that the cell doublet forms a coherent unit. However, force propagation and response of the receiver cell also strongly depend on the mechano-structural polarization in the cell assembly, which is controlled by cell-matrix adhesion to the adhesive micropattern. We find that the response of the receiver cell is stronger when the mechano-structural polarization axis is oriented perpendicular to the direction of force propagation, reminiscent of the Poisson effect in passive materials. We finally show that the same effects are at work in small tissues. Our work demonstrates that cellular organization and active mechanical response of a tissue are key to maintain signal strength and lead to the emergence of elasticity, which means that signals are not dissipated like in a viscous system, but can propagate over large distances.

Force generation in human blood platelets by filamentous actomyosin structures.

Blood platelets are central elements of the blood clotting response after wounding. Upon vessel damage, they bind to the surrounding matrix and contract the forming thrombus, thus helping to restore normal blood circulation. The hemostatic function of platelets is directly connected to their mechanics and cytoskeletal organization. The reorganization of the platelet cytoskeleton during spreading occurs within minutes and leads to the formation of contractile actomyosin bundles, but it is not known if there is a direct correlation between the emerging actin structures and the force field that is exerted to the environment. In this study, we combine fluorescence imaging of the actin structures with simultaneous traction force measurements in a time-resolved manner. In addition, we image the final states with superresolution microscopy. We find that both the force fields and the cell shapes have clear geometrical patterns defined by stress fibers. Force generation is localized in a few hotspots, which appear early during spreading, and, in the mature state, anchor stress fibers in focal adhesions. Moreover, we show that, for a gel stiffness in the physiological range, force generation is a very robust mechanism and we observe no systematic dependence on the amount of added thrombin in solution or fibrinogen coverage on the substrate, suggesting that force generation after platelet activation is a threshold phenomenon that ensures reliable thrombus contraction in diverse environments.

The Ebola virus VP40 matrix layer undergoes endosomal disassembly essential for membrane fusion.

Ebola viruses (EBOVs) assemble into filamentous virions, whose shape and stability are determined by the matrix viral protein 40 (VP40). Virus entry into host cells occurs via membrane fusion in late endosomes; however, the mechanism of how the remarkably long virions undergo uncoating, including virion disassembly and nucleocapsid release into the cytosol, remains unknown. Here, we investigate the structural architecture of EBOVs entering host cells and discover that the VP40 matrix disassembles prior to membrane fusion. We reveal that VP40 disassembly is caused by the weakening of VP40-lipid interactions driven by low endosomal pH that equilibrates passively across the viral envelope without a dedicated ion channel. We further show that viral membrane fusion depends on VP40 matrix integrity, and its disassembly reduces the energy barrier for fusion stalk formation. Thus, pH-driven structural remodeling of the VP40 matrix acts as a molecular switch coupling viral matrix uncoating to membrane fusion during EBOV entry.

High curvature promotes fusion of lipid membranes: Predictions from continuum elastic theory.

The fusion of lipid membranes progresses through a series of hemifusion intermediates with two significant energy barriers related to the formation of stalk and fusion pore, respectively. These energy barriers determine the speed and success rate of many critical biological processes, including the fusion of highly curved membranes, for example synaptic vesicles and enveloped viruses. Here we use continuum elastic theory of lipid monolayers to determine the relationship between membrane shape and energy barriers to fusion. We find that the stalk formation energy decreases with curvature by up to 31 kBT in a 20-nm-radius vesicle compared with planar membranes and by up to 8 kBT in the fusion of highly curved, long, tubular membranes. In contrast, the fusion pore formation energy barrier shows a more complicated behavior. Immediately after stalk expansion to the hemifusion diaphragm, the fusion pore formation energy barrier is low (15-25 kBT) due to lipid stretching in the distal monolayers and increased tension in highly curved vesicles. Therefore, the opening of the fusion pore is faster. However, these stresses relax over time due to lipid flip-flop from the proximal monolayer, resulting in a larger hemifusion diaphragm and a higher fusion pore formation energy barrier, up to 35 kBT. Therefore, if the fusion pore fails to open before significant lipid flip-flop takes place, the reaction proceeds to an extended hemifusion diaphragm state, which is a dead-end configuration in the fusion process and can be used to prevent viral infections. In contrast, in the fusion of long tubular compartments, the surface tension does not accumulate due to the formation of the diaphragm, and the energy barrier for pore expansion increases with curvature by up to 11 kBT. This suggests that inhibition of polymorphic virus infection could particularly target this feature of the second barrier.

IFITM3 blocks influenza virus entry by sorting lipids and stabilizing hemifusion.

Interferon-induced transmembrane protein 3 (IFITM3) inhibits the entry of numerous viruses through undefined molecular mechanisms. IFITM3 localizes in the endosomal-lysosomal system and specifically affects virus fusion with target cell membranes. We found that IFITM3 induces local lipid sorting, resulting in an increased concentration of lipids disfavoring viral fusion at the hemifusion site. This increases the energy barrier for fusion pore formation and the hemifusion dwell time, promoting viral degradation in lysosomes. In situ cryo-electron tomography captured IFITM3-mediated arrest of influenza A virus membrane fusion. Observation of hemifusion diaphragms between viral particles and late endosomal membranes confirmed hemifusion stabilization as a molecular mechanism of IFITM3. The presence of the influenza fusion protein hemagglutinin in post-fusion conformation close to hemifusion sites further indicated that IFITM3 does not interfere with the viral fusion machinery. Collectively, these findings show that IFITM3 induces lipid sorting to stabilize hemifusion and prevent virus entry into target cells.

Grand canonical Brownian dynamics simulations of adsorption and self-assembly of SAS-6 rings on a surface.

The Spindle Assembly Abnormal Protein 6 (SAS-6) forms dimers, which then self-assemble into rings that are critical for the nine-fold symmetry of the centriole organelle. It has recently been shown experimentally that the self-assembly of SAS-6 rings is strongly facilitated on a surface, shifting the reaction equilibrium by four orders of magnitude compared to the bulk. Moreover, a fraction of non-canonical symmetries (i.e., different from nine) was observed. In order to understand which aspects of the system are relevant to ensure efficient self-assembly and selection of the nine-fold symmetry, we have performed Brownian dynamics computer simulation with patchy particles and then compared our results with the experimental ones. Adsorption onto the surface was simulated by a grand canonical Monte Carlo procedure and random sequential adsorption kinetics. Furthermore, self-assembly was described by Langevin equations with hydrodynamic mobility matrices. We find that as long as the interaction energies are weak, the assembly kinetics can be described well by coagulation-fragmentation equations in the reaction-limited approximation. By contrast, larger interaction energies lead to kinetic trapping and diffusion-limited assembly. We find that the selection of nine-fold symmetry requires a small value for the angular interaction range. These predictions are confirmed by the experimentally observed reaction constant and angle fluctuations. Overall, our simulations suggest that the SAS-6 system works at the crossover between a relatively weak binding energy that avoids kinetic trapping and a small angular range that favors the nine-fold symmetry.

Clathrin coats partially preassemble and subsequently bend during endocytosis.

Eukaryotic cells use clathrin-mediated endocytosis to take up a large range of extracellular cargo. During endocytosis, a clathrin coat forms on the plasma membrane, but it remains controversial when and how it is remodeled into a spherical vesicle. Here, we use 3D superresolution microscopy to determine the precise geometry of the clathrin coat at large numbers of endocytic sites. Through pseudo-temporal sorting, we determine the average trajectory of clathrin remodeling during endocytosis. We find that clathrin coats assemble first on flat membranes to 50% of the coat area before they become rapidly and continuously bent, and this mechanism is confirmed in three cell lines. We introduce the cooperative curvature model, which is based on positive feedback for curvature generation. It accurately describes the measured shapes and dynamics of the clathrin coat and could represent a general mechanism for clathrin coat remodeling on the plasma membrane.

Simulating 3D Cell Shape with the Cellular Potts Model.

Computer simulations have become a widely used method for the field of mechanobiology. An important question is whether one can predict the shape and forces of cells as a function of the extracellular environment. Different types of models have been described before to simulate cell and tissue shapes in structured environments. In this chapter, we give a brief overview of commonly used models and then describe the Cellular Potts Model, a lattice-based modelling framework, in more detail. We provide a hands-on guide on how to build a model that simulates the shape of a single cell on a micropattern in three dimensions in different open source software packages using the Cellular Potts framework. A simulation is set up with an initial configuration of generalized cells that change shape and position due to an energy function that incorporates cellular volume and surface area constraints as well as interaction energies between the generalized cells.

An Improved Method for Assessing Antigen Presentation on the Surface of Plasmodium falciparum-Infected Erythrocytes by Immuno-Electron Microscopy.

Immuno-electron microscopy can detect and localize antigens in cells or tissues at a resolution of several nanometers. In the case of P. falciparum-infected erythrocytes, immuno-EM studies are frequently hampered by the electron-dense nature of the hemoglobin and access of antibodies to antigenic sites, particularly if the targeted protein is presented on the host cell surface or lies in proximity to the host cell cytoskeleton. Here, we describe an improved immuno-EM protocol that overcomes these problems. The improved signal to noise ratio and the enhanced access to antigenic sites now allows one to obtain information regarding target density and distribution and, hence, additional insights into the architecture and function of parasite-induced, or -affected, structures.

A particle-based computational model to analyse remodelling of the red blood cell cytoskeleton during malaria infections.

Red blood cells can withstand the harsh mechanical conditions in the vasculature only because the bending rigidity of their plasma membrane is complemented by the shear elasticity of the underlying spectrin-actin network. During an infection by the malaria parasite Plasmodium falciparum, the parasite mines host actin from the junctional complexes and establishes a system of adhesive knobs, whose main structural component is the knob-associated histidine rich protein (KAHRP) secreted by the parasite. Here we aim at a mechanistic understanding of this dramatic transformation process. We have developed a particle-based computational model for the cytoskeleton of red blood cells and simulated it with Brownian dynamics to predict the mechanical changes resulting from actin mining and KAHRP-clustering. Our simulations include the three-dimensional conformations of the semi-flexible spectrin chains, the capping of the actin protofilaments and several established binding sites for KAHRP. For the healthy red blood cell, we find that incorporation of actin protofilaments leads to two regimes in the shear response. Actin mining decreases the shear modulus, but knob formation increases it. We show that dynamical changes in KAHRP binding affinities can explain the experimentally observed relocalization of KAHRP from ankyrin to actin complexes and demonstrate good qualitative agreement with experiments by measuring pair cross-correlations both in the computer simulations and in super-resolution imaging experiments.

Asynchronous nuclear cycles in multinucleated Plasmodium falciparum facilitate rapid proliferation.

Malaria-causing parasites proliferate within erythrocytes through schizogony, forming multinucleated stages before cellularization. Nuclear multiplication does not follow a strict geometric 2n progression, and each proliferative cycle produces a variable number of progeny. Here, by tracking nuclei and DNA replication, we show that individual nuclei replicate their DNA at different times, despite residing in a shared cytoplasm. Extrapolating from experimental data using mathematical modeling, we provide strong indication that a limiting factor exists, which slows down the nuclear multiplication rate. Consistent with this prediction, our data show that temporally overlapping DNA replication events were significantly slower than partially overlapping or nonoverlapping events. Our findings suggest the existence of evolutionary pressure that selects for asynchronous DNA replication, balancing available resources with rapid pathogen proliferation.