Researcher profile

Ulrich S. Schwarz

Ulrich S. Schwarz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

First-passage times in complex energy landscapes: a case study with nonmuscle myosin II assembly

Complex energy landscapes often arise in biological systems, e.g. for protein folding, biochemical reactions or intracellular transport processes. Their physical effects are often reflected in the first-passage times arising from these energy landscapes. However, their calculation is notoriously challenging and it is often difficult to identify the most relevant features of a given energy landscape. Here we show how this can be achieved by coarse-graining the Fokker-Planck equation to a master equation and decomposing its first-passage times in an iterative process. We apply this method to the electrostatic interaction between two rods of nonmuscle myosin II (NM2), which is the main molecular motor for force generation in nonmuscle cells. Energy landscapes are computed directly from the amino acid sequences of the three different isoforms. Our approach allows us to identify the most relevant energy barriers for their self-assembly into nonmuscle myosin II minifilaments and how they change under force. In particular, we find that antiparallel configurations are more stable than parallel ones, but also show more changes under mechanical loading. Our work demonstrates the rich dynamics that can be expected for NM2-assemblies under mechanical load and in general shows how one can identify the most relevant energy barriers in complex energy landscapes.

preprint2020arXiv

Rheology of mixed motor ensembles

The rheology of biological cells is not only determined by their cytoskeletal networks, but also by the molecular motors that crosslink and contract them. Recently it has been found that the assemblies of myosin II molecular motors in non-muscle cells are mixtures of fast and slow motor variants. Using computer simulations and an analytical mean field theory of a crossbridge model for myosin II motors, here we show that such motor ensembles effectively behave as active Maxwell elements. We calculate storage and loss moduli as a function of the model parameters and show that the rheological properties cross over from viscous to elastic as one increases the ratio of slow to fast motors. This suggests that cells tune their mechanical properties by regulating the composition of their myosin assemblies.

preprint2019arXiv

Mechanosensitive Self-Assembly of Myosin II Minifilaments

Self-assembly and force generation are two central processes in biological systems that usually are considered in separation. However, the signals that activate non-muscle myosin II molecular motors simultaneously lead to self-assembly into myosin II minifilaments as well as progression of the motor heads through the crossbridge cycle. Here we investigate theoretically the possible effects of coupling these two processes. Our assembly model, which builds upon a consensus architecture of the minifilament, predicts a critical aggregation concentration at which the assembly kinetics slows down dramatically. The combined model predicts that increasing actin filament concentration and force both lead to a decrease in the critical aggregation concentration. We suggest that due to these effects, myosin II minifilaments in a filamentous context might be in a critical state that reacts faster to varying conditions than in solution. We finally compare our model to experiments by simulating fluorescence recovery after photobleaching.

preprint2016arXiv

Measuring cellular traction forces on non-planar substrates

Animal cells use traction forces to sense the mechanics and geometry of their environment. Measuring these traction forces requires a workflow combining cell experiments, image processing and force reconstruction based on elasticity theory. Such procedures have been established before mainly for planar substrates, in which case one can use the Green's function formalism. Here we introduce a worksflow to measure traction forces of cardiac myofibroblasts on non-planar elastic substrates. Soft elastic substrates with a wave-like topology were micromolded from polydimethylsiloxane (PDMS) and fluorescent marker beads were distributed homogeneously in the substrate. Using feature vector based tracking of these marker beads, we first constructed a hexahedral mesh for the substrate. We then solved the direct elastic boundary volume problem on this mesh using the finite element method (FEM). Using data simulations, we show that the traction forces can be reconstructed from the substrate deformations by solving the corresponding inverse problem with a L1-norm for the residue and a L2-norm for 0th order Tikhonov regularization. Applying this procedure to the experimental data, we find that cardiac myofibroblast cells tend to align both their shapes and their forces with the long axis of the deformable wavy substrate.

preprint2015arXiv

Role of dynamic capsomere supply for viral capsid self-assembly

Many viruses rely on the self-assembly of their capsids to protect and transport their genomic material. For many viral systems, in particular for human viruses like hepatitis B, adeno or human immunodeficiency virus, that lead to persistent infections, capsomeres are continuously produced in the cytoplasm of the host cell while completed capsids exit the cell for a new round of infection. Here we use coarse-grained Brownian dynamics simulations of a generic patchy particle model to elucidate the role of the dynamic supply of capsomeres for the reversible self-assembly of empty T1 icosahedral virus capsids. We find that for high rates of capsomere influx only a narrow range of bond strengths exists for which a steady state of continuous capsid production is possible. For bond strengths smaller and larger than this optimal value, the reaction volume becomes crowded by small and large intermediates, respectively. For lower rates of capsomere influx a broader range of bond strengths exists for which a steady state of continuous capsid production is established, although now the production rate of capsids is smaller. Thus our simulations suggest that the importance of an optimal bond strength for viral capsid assembly typical for in vitro conditions can be reduced by the dynamic influx of capsomeres in a cellular environment.

preprint2015arXiv

Stochastic dynamics of virus capsid formation: direct versus hierarchical self-assembly

In order to replicate within their cellular host, many viruses have developed self-assembly strategies for their capsids which are sufficiently robust as to be reconstituted in vitro. Mathematical models for virus self-assembly usually assume that the bonds leading to cluster formation have constant reactivity over the time course of assembly (direct assembly). In some cases, however, binding sites between the capsomers have been reported to be activated during the self-assembly process (hierarchical assembly). In order to study possible advantages of such hierarchical schemes for icosahedral virus capsid assembly, we use Brownian dynamics simulations of a patchy particle model that allows us to switch binding sites on and off during assembly. For T1 viruses, we implement a hierarchical assembly scheme where inter-capsomer bonds become active only if a complete pentamer has been assembled. We find direct assembly to be favorable for reversible bonds allowing for repeated structural reorganizations, while hierarchical assembly is favorable for strong bonds with small dissociation rate, as this situation is less prone to kinetic trapping. However, at the same time it is more vulnerable to monomer starvation during the final phase. Increasing the number of initial monomers does have only a weak effect on these general features. The differences between the two assembly schemes become more pronounced for more complex virus geometries, as shown here for T3 viruses, which assemble through homogeneous pentamers and heterogeneous hexamers in the hierarchical scheme. In order to complement the simulations for this more complicated case, we introduce a master equation approach that agrees well with the simulation results.

preprint2015arXiv

Traction force microscopy on soft elastic substrates: a guide to recent computational advances

The measurement of cellular traction forces on soft elastic substrates has become a standard tool for many labs working on mechanobiology. Here we review the basic principles and different variants of this approach. In general, the extraction of the substrate displacement field from image data and the reconstruction procedure for the forces are closely linked to each other and limited by the presence of experimental noise. We discuss different strategies to reconstruct cellular forces as they follow from the foundations of elasticity theory, including two- versus three-dimensional, inverse versus direct and linear versus non-linear approaches. We also discuss how biophysical models can improve force reconstruction and comment on practical issues like substrate preparation, image processing and the availability of software for traction force microscopy.

preprint2014arXiv

Dynamics of Cell Shape and Forces on Micropatterned Substrates Predicted by a Cellular Potts Model

Micropatterned substrates are often used to standardize cell experiments and to quantitatively study the relation between cell shape and function. Moreover, they are increasingly used in combination with traction force microscopy on soft elastic substrates. To predict the dynamics and steady states of cell shape and forces without any a priori knowledge of how the cell will spread on a given micropattern, here we extend earlier formulations of the two-dimensional cellular Potts model. The third dimension is treated as an area reservoir for spreading. To account for local contour reinforcement by peripheral bundles, we augment the cellular Potts model by elements of the tension-elasticity model. We first parameterize our model and show that it accounts for momentum conservation. We then demonstrate that it is in good agreement with experimental data for shape, spreading dynamics, and traction force patterns of cells on micropatterned substrates. We finally predict shapes and forces for micropatterns that have not yet been experimentally studied.

preprint2014arXiv

Geometrical model for malaria parasite migration in structured environments

Malaria is transmitted to vertebrates via a mosquito bite, during which rod-like and crescent-shaped parasites, called sporozoites, are injected into the skin of the host. Searching for a blood capillary to penetrate, sporozoites move quickly in locally helical trajectories, that are frequently perturbed by interactions with the extracellular environment. Here we present a theoretical analysis of the active motility of sporozoites in a structured environment. The sporozoite is modelled as a self-propelled rod with spontaneous curvature and bending rigidity. It interacts with hard obstacles through collision rules inferred from experimental observation of two-dimensional sporozoite movement in pillar arrays. Our model shows that complex motion patterns arise from the geometrical shape of the parasite and that its mechanical flexibility is crucial for stable migration patterns. Extending the model to three dimensions reveals that a bent and twisted rod can associate to cylindrical obstacles in a manner reminiscent of the association of sporozoites to blood capillaries, supporting the notion of a prominent role of cell shape during malaria transmission.

preprint2014arXiv

Stochastic dynamics and mechanosensitivity of myosin II minifilaments

Tissue cells are in a state of permanent mechanical tension that is maintained mainly by myosin II minifilaments, which are bipolar assemblies of tens of myosin II molecular motors contracting actin networks and bundles. Here we introduce a stochastic model for myosin II minifilaments as two small myosin II motor ensembles engaging in a stochastic tug-of-war. Each of the two ensembles is described by the parallel cluster model that allows us to use exact stochastic simulations and at the same time to keep important molecular details of the myosin II cross-bridge cycle. Our simulation and analytical results reveal a strong dependence of myosin II minifilament dynamics on environmental stiffness that is reminiscent of the cellular response to substrate stiffness. For small stiffness, minifilaments form transient crosslinks exerting short spikes of force with negligible mean. For large stiffness, minifilaments form near permanent crosslinks exerting a mean force which hardly depends on environmental elasticity. This functional switch arises because dissociation after the power stroke is suppressed by force (catch bonding) and because ensembles can no longer perform the power stroke at large forces. Symmetric myosin II minifilaments perform a random walk with an effective diffusion constant which decreases with increasing ensemble size, as demonstrated for rigid substrates with an analytical treatment.

preprint2014arXiv

Studying protein assembly with reversible Brownian dynamics of patchy particles

Assembly of protein complexes like virus shells, the centriole, the nuclear pore complex or the actin cytoskeleton is strongly determined by their spatial structure. Moreover it is becoming increasingly clear that the reversible nature of protein assembly is also an essential element for their biological function. Here we introduce a computational approach for the Brownian dynamics of patchy particles with anisotropic assemblies and fully reversible reactions. Different particles stochastically associate and dissociate with microscopic reaction rates depending on their relative spatial positions. The translational and rotational diffusive properties of all protein complexes are evaluated on-the-fly. Because we focus on reversible assembly, we introduce a scheme which ensures detailed balance for patchy particles. We then show how the macroscopic rates follow from the microscopic ones. As an instructive example, we study the assembly of a pentameric ring structure, for which we find excellent agreement between simulation results and a macroscopic kinetic description without any adjustable parameters. This demonstrates that our approach correctly accounts for both the diffusive and reactive processes involved in protein assembly.

preprint2013arXiv

Mesoscopic model for filament orientation in growing actin networks: the role of obstacle geometry

Propulsion by growing actin networks is a universal mechanism used in many different biological systems. Although the core molecular machinery for actin network growth is well preserved in most cases, the geometry of the propelled obstacle can vary considerably. In recent years, filament orientation distribution has emerged as an important observable characterizing the structure and dynamical state of the growing network. Here we derive several continuum equations for the orientation distribution of filaments growing behind stiff obstacles of various shapes and validate the predicted steady state orientation patterns by stochastic computer simulations based on discrete filaments. We use an ordinary differential equation approach to demonstrate that for flat obstacles of finite size, two fundamentally different orientation patterns peaked at either +35/-35 or +70/0/-70 degrees exhibit mutually exclusive stability, in agreement with earlier results for flat obstacles of very large lateral extension. We calculate and validate phase diagrams as a function of model parameters and show how this approach can be extended to obstacles with piecewise straight contours. For curved obstacles, we arrive at a partial differential equation in the continuum limit, which again is in good agreement with the computer simulations. In all cases, we can identify the same two fundamentally different orientation patterns, but only within an appropriate reference frame, which is adjusted to the local orientation of the obstacle contour. Our results suggest that two fundamentally different network architectures compete with each other in growing actin networks, irrespective of obstacle geometry, and clarify how simulated and electron tomography data have to be analyzed for non-flat obstacle geometries.

preprint2013arXiv

Physics of adherent cells

One of the most unique physical features of cell adhesion to external surfaces is the active generation of mechanical force at the cell-material interface. This includes pulling forces generated by contractile polymer bundles and networks, and pushing forces generated by the polymerization of polymer networks. These forces are transmitted to the substrate mainly by focal adhesions, which are large, yet highly dynamic adhesion clusters. Tissue cells use these forces to sense the physical properties of their environment and to communicate with each other. The effect of forces is intricately linked to the material properties of cells and their physical environment. Here a review is given of recent progress in our understanding of the role of forces in cell adhesion from the viewpoint of theoretical soft matter physics and in close relation to the relevant experiments.

preprint2013arXiv

Stochastic dynamics of small ensembles of non-processive molecular motors: the parallel cluster model

Non-processive molecular motors have to work together in ensembles in order to generate appreciable levels of force or movement. In skeletal muscle, for example, hundreds of myosin II molecules cooperate in thick filaments. In non-muscle cells, by contrast, small groups with few tens of non-muscle myosin II motors contribute to essential cellular processes such as transport, shape changes or mechanosensing. Here we introduce a detailed and analytically tractable model for this important situation. Using a three-state crossbridge model for the myosin II motor cycle and exploiting the assumptions of fast power stroke kinetics and equal load sharing between motors in equivalent states, we reduce the stochastic reaction network to a one-step master equation for the binding and unbinding dynamics (parallel cluster model) and derive the rules for ensemble movement. We find that for constant external load, ensemble dynamics is strongly shaped by the catch bond character of myosin II, which leads to an increase of the fraction of bound motors under load and thus to firm attachment even for small ensembles. This adaptation to load results in a concave force-velocity relation described by a Hill relation. For external load provided by a linear spring, myosin II ensembles dynamically adjust themselves towards an isometric state with constant average position and load. The dynamics of the ensembles is now determined mainly by the distribution of motors over the different kinds of bound states. For increasing stiffness of the external spring, there is a sharp transition beyond which myosin II can no longer perform the power stroke. Slow unbinding from the pre-power-stroke state protects the ensembles against detachment.

preprint2013arXiv

Unifying autocatalytic and zeroth order branching models for growing actin networks

The directed polymerization of actin networks is an essential element of many biological processes, including cell migration. Different theoretical models considering the interplay between the underlying processes of polymerization, capping and branching have resulted in conflicting predictions. One of the main reasons for this discrepancy is the assumption of a branching reaction that is either first order (autocatalytic) or zeroth order in the number of existing filaments. Here we introduce a unifying framework from which the two established scenarios emerge as limiting cases for low and high filament number. A smooth transition between the two cases is found at intermediate conditions. We also derive a threshold for the capping rate, above which autocatalytic growth is predicted at sufficiently low filament number. Below the threshold, zeroth order characteristics are predicted to dominate the dynamics of the network for all accessible filament numbers. Together, this allows cells to grow stable actin networks over a large range of different conditions.

preprint2012arXiv

Force localization in contracting cell layers

Epithelial cell layers on soft elastic substrates or pillar arrays are commonly used as model systems for investigating the role of force in tissue growth, maintenance and repair. Here we show analytically that the experimentally observed localization of traction forces to the periphery of the cell layers does not necessarily imply increased local cell activity, but follows naturally from the elastic problem of a finite-sized contractile layer coupled to an elastic foundation. For homogeneous contractility, the force localization is determined by one dimensionless parameter interpolating between linear and exponential force profiles for the extreme cases of very soft and very stiff substrates, respectively. If contractility is sufficiently increased at the periphery, outward directed displacements can occur at intermediate positions, although the edge itself still retracts. We also show that anisotropic extracellular stiffness leads to force localization in the stiffer direction, as observed experimentally.

preprint2012arXiv

Stochastic force generation by small ensembles of myosin II motors

Forces in the actin cytoskeleton are generated by small groups of non-processive myosin II motors for which stochastic effects are highly relevant. Using a crossbridge model with the assumptions of fast powerstroke kinetics and equal load sharing between equivalent states, we derive a one-step master equation for the activity of a finite-sized ensemble of mechanically coupled myosin II motors. For constant external load, this approach yields analytical results for duty ratio and force-velocity relation as a function of ensemble size. We find that stochastic effects cannot be neglected for ensemble sizes below 15. The one-step master equation can be used also for efficient computer simulations with linear elastic external load and reveals the sequence of build-up of force and ensemble rupture that is characteristic for reconstituted actomyosin contractility.

preprint2011arXiv

Viscoelastic response of contractile filament bundles

The actin cytoskeleton of adherent tissue cells often condenses into filament bundles contracted by myosin motors, so-called stress fibers, which play a crucial role in the mechanical interaction of cells with their environment. Stress fibers are usually attached to their environment at the endpoints, but possibly also along their whole length. We introduce a theoretical model for such contractile filament bundles which combines passive viscoelasticity with active contractility. The model equations are solved analytically for two different types of boundary conditions. A free boundary corresponds to stress fiber contraction dynamics after laser surgery and results in good agreement with experimental data. Imposing cyclic varying boundary forces allows us to calculate the complex modulus of a single stress fiber.

preprint2010arXiv

Role of anisotropy for protein-protein encounter

Protein-protein interactions comprise both transport and reaction steps. During the transport step, anisotropy of proteins and their complexes is important both for hydrodynamic diffusion and accessibility of the binding site. Using a Brownian dynamics approach and extensive computer simulations, we quantify the effect of anisotropy on the encounter rate of ellipsoidal particles covered with spherical encounter patches. We show that the encounter rate $k$ depends on the aspect ratios $ξ$ mainly through steric effects, while anisotropic diffusion has only a little effect. Calculating analytically the crossover times from anisotropic to isotropic diffusion in three dimensions, we find that they are much smaller than typical protein encounter times, in agreement with our numerical results.

preprint2009arXiv

Stochastic simulations of cargo transport by processive molecular motors

We use stochastic computer simulations to study the transport of a spherical cargo particle along a microtubule-like track on a planar substrate by several kinesin-like processive motors. Our newly developed adhesive motor dynamics algorithm combines the numerical integration of a Langevin equation for the motion of a sphere with kinetic rules for the molecular motors. The Langevin part includes diffusive motion, the action of the pulling motors, and hydrodynamic interactions between sphere and wall. The kinetic rules for the motors include binding to and unbinding from the filament as well as active motor steps. We find that the simulated mean transport length increases exponentially with the number of bound motors, in good agreement with earlier results. The number of motors in binding range to the motor track fluctuates in time with a Poissonian distribution, both for springs and cables being used as models for the linker mechanics. Cooperativity in the sense of equal load sharing only occurs for high values for viscosity and attachment time.