The fusion of cell aggregates widely exists during biological processes such as development, tissue regeneration, and tumor invasion. Cellular spheroids (spherical cell aggregates) are commonly used to study this phenomenon. In previous studies, with approximated assumptions and measurements, researchers found that the fusion of two spheroids with some cell type is similar to the coalescence of two liquid droplets. However, with more accurate measurements focusing on the overall shape evolution in this process, we find that even in the previously-regarded liquid-like regime, the fusion process of spheroids can be very different from regular liquid coalescence. We conduct numerical simulations using both standard particulate models and vertex models with both Molecular Dynamics and Brownian Dynamics. The simulation results show that the difference between spheroids and regular liquid droplets is caused by the microscopic overdamped dynamics of each cell rather than the topological cell-cell interactions in the vertex model. Our research reveals the necessity of a new continuum theory for “liquid” with microscopically overdamped components, such as cellular and colloidal systems. Detailed analysis of our simulation results of different system sizes provides the basis for developing the new theory.

We present new kinematic bending measures and quadratic energies for isotropic elastic plates and shells, with certain desirable features not present in commonly employed models in mechanics and soft matter. These are justified both by simple physical arguments related to the through-thickness variation in strain, and through a detailed reduction from a three-dimensional energy quadratic in stretch. The measure of plate bending is a dilation-invariant surface tensor that couples stretch and curvature in a natural extension of primitive generalized bending strains for straight rods. The extension to naturally-curved rods and shells, for which the pure stretching of a curved rest configuration is not a dilation, contrasts with previous ad hoc postulated forms. Our results provide a clean basis for simple models of low-dimensional elastic systems, and should enable more accurate probing of the structure of singularities in soft sheets and membranes.

Self-organization is a universal concept spanning numerous disciplines including mathematics, physics and biology. Chromosomes are self-organizing polymers that fold into orderly, hierarchical and yet dynamic structures. In the past decade, advances in experimental biology have provided a means to reveal information about chromosome connectivity, allowing us to directly use this information from experiments to generate 3D models of individual genes, chromosomes and even genomes. In this talk I will present a novel data-driven modeling approach and discuss a number of possibilities that this method holds. I will discuss a detailed study of the time-evolution of X chromosome inactivation, highlighting both global and local properties of chromosomes that result in topology-driven dynamical arrest and present and characterize a novel type of motion we discovered in knots that may have applications to nanoscale materials and machines.

Macromolecular protein complexes perform essential biological functions across life forms. A fundamental, though yet unsolved question in biology is how the function of such complexes is regulated by intracellular or extracellular signals. For instance, we have little understanding of how forces affect multi-protein machines whose function is often mechanical in nature. We address this question by studying the bacterial flagellar motor, a large complex that powers swimming motility in many bacteria. This rotary motor autonomously adapts to changes in mechanical load by adding or removing force-generating ‘stator’ units that power rotation. In the bacterium Escherichia coli, up to 11 units drive the motor at high load while all the units are released at low load. We manipulate motor load using electrorotation, a technique in which a rapidly rotating electric field applies an external torque on the motor. This allows us to change motor load at will and measure the resulting stator dynamics at single-unit resolution. We found that the force generated by the stator units controls their unbinding, forming a feedback loop that leads to autoregulation of the assembly. We complemented our experiments with theoretical models that provide insight into the underlying molecular interactions. Torque-dependent remodeling takes place within seconds, making it a highly responsive control mechanism, one that is mediated by the mechano-chemical tuning of protein interactions.

There is a recent surge of interest in the behavior of active particles that can at the same time align their direction of movement and synchronize their oscillations, known as swarmalators. While analytical and numerical models of such systems are now abundant, no real-life examples have been shown to date. I will present an experimental investigation of the collective motion of the nematode Turbatrix aceti, which self-propel by body undulation. I will show that under favorable conditions these nematodes can synchronize their body oscillations, forming striking traveling metachronal waves which, similar to the case of beating cilia, produce strong fluid flows. I will demonstrate that the location and strength of this collective state can be controlled through the shape of the confining structure; in our case the contact angle of a droplet. This opens a way for producing controlled work such as on-demand flows or displacement of objects. I will illustrate this by a practical example: showing that the force generated by the collectively moving nematodes is sufficient to change the mode of evaporation of fluid droplets, by counteracting the surface-tension force, which allow us to estimate its strength.

After more than 4 million deaths worldwide, the ongoing vaccination to conquer the COVID-19 disease is now competing with the emergence of increasingly contagious mutations, repeatedly supplanting earlier strains. Following the near-absence of historical examples of the long-time evolution of infectious diseases under similar circumstances, models are crucial to exemplify possible scenarios. Accordingly, in the present work we systematically generalize the popular susceptible-infected-recovered model to account for mutations leading to repeatedly occurring new strains, which we coarse grain based on tools from statistical mechanics to derive a model predicting the most likely outcomes. The model predicts that mutations can induce a super exponential growth of infection numbers at early times, which self-amplify to giant infection waves which are caused by a positive feedback loop between infection numbers and mutations and lead to a simultaneous infection of the majority of the population. At later stages -- if vaccination progresses too slowly -- mutations can interrupt an ongoing decrease of infection numbers and can cause infection revivals which occur as single waves or even as whole wave trains featuring alternative periods of decreasing and increasing infection numbers. Our results might be useful for discussions regarding the importance of a release of vaccine-patents to reduce the risk of mutation-induced infection revivals but also to coordinate the release of measures following a downwards trend of infection numbers.

Dictyostelium Discoideum (DD) are a fascinating single-cellular organism. When nutrients are plentiful, the DD cells act as autonomous individuals foraging their local vicinity. At the onset of starvation, a few (<0.1%) cells begin communicating with others by emitting a spike in the chemoattractant protein cyclic-AMP. Nearby cells sense the chemical gradient and respond by moving toward it and emitting a cyclic-AMP spike of their own. Cyclic-AMP activity increases over time, and eventually a spiral wave emerges, attracting hundreds of thousands of cells to an aggregation center. How DD cells go from autonomous individuals to a collective entity remains an open question for more than 60 years--a question whose answer would shed light on the emergence of multi-cellular life. Recently, trans-scale imaging has allowed the ability to sense the cyclic-AMP activity at both cell and colony levels. Using both the images as well as toy simulation models, this research aims to clarify whether the activity at the colony level is in fact initiated by a few cells, which may be deemed "leader" or "pacemaker" cells. In this talk, I will demonstrate the use of information-theoretic techniques to classify leaders and followers based on trajectory data, as well as to infer the domain of interaction of leader cells. We validate the techniques on toy models where leaders and followers are known, and then try to answer the question in real data--do leader cells drive collective behavior in DD colonies?

Many microorganisms and cells function in complex (non-Newtonian) fluids, which are mixtures of different materials and exhibit both viscous and elastic stresses. For example, mammalian sperm swim through cervical mucus on their journey through the female reproductive tract, and they must penetrate the viscoelastic gel outside the ovum to fertilize. In micro-scale swimming the dynamics emerge from the coupled interactions between the complex rheology of the surrounding media and the passive and active body dynamics of the swimmer. We use computational models of swimmers in viscoelastic fluids to investigate and provide mechanistic explanations for emergent swimming behaviors. I will discuss how flexible filaments (such as flagella) can store energy from a viscoelastic fluid to gain stroke boosts due to fluid elasticity. I will also describe 3D simulations of model organisms such as C. Reinhardtii and mammalian sperm, where we use experimentally measured stroke data to separate naturally coupled stroke and fluid effects. We explore why strokes that are adapted to Newtonian fluid environments might not do well in viscoelastic environments.

The settling speeds of small biological particles influence the distribution of organisms such as plants, corals, and phytoplankton, but these speeds are difficult to quantify without magnification. In this talk, I highlight my novel method,  using dynamic scaling principles and 3D printed models to solve this problem. Dynamic scaling involves creating models with differ in size to the original system and match the physical forces acting upon the model to the original system. I discuss the methodology behind the technique and show how it differs to previous works using dynamically scaled models. I show the flexibility of the technique and suggest how it can be applied to other free-falling particles (e.g. seeds and spores).

At small length-scales, capillary effects are significant, and thus the mechanics of soft material interfaces may be dominated by solid surface stresses or liquid surface tensions. The balance between surface and bulk properties is described by an elasto-capillary length-scale, in which equilibrium interfacial energies are constant. However, at small length-scales in biological materials, including living cells and tissues, interfacial energies are not constant but are actively regulated and driven far from equilibrium. Thus, the balance between surface and bulk properties depends upon the distance from equilibrium. Here, we model the spreading (wetting) of living cell aggregates as ‘active droplets’, with a non-equilibrium surface tension that depends upon internal stress generated by the actomyosin cytoskeleton. Depending upon the extent of activity, droplet surface properties adapt to the mechanics of their surroundings. The impact of this adaptation challenges contemporary models of interfacial mechanics, including extensively used models of contact mechanics and wetting.

The SARS-CoV-2 pandemic is awash in data, including daily, spatially-resolved COVID case data, virus sequence data, patients `omics data, and mobility data. Journals are now also awash in studies that make use of quantitative modeling approaches to gain insight into the geographic spread of SARS-CoV-2 and its temporal dynamics, as well as studies that predict the impact of control strategies on SARS-CoV-2 circulation. Some, but by no means all, of these studies are informed by the massive amounts of available data. Some, but by no means all, of these studies have been useful — in that their predictions revealed something beyond simple back of the envelope calculations. To summarize some of these findings, in this symposium, we will address questions such as: What do we want from models of disease spread? What can and should be predicted? Which data are the most useful for predictions? When do we need mechanistic models? What have we learned about how to model disease spread from unmet and/or conflicting predictions? The workshop speakers will explore these questions from different perspectives on what data need to be considered and how models can be evaluated. As at other TMLS workshops, each speaker will deliver a 10-minute talk with ample time set aside for moderated questions/discussion. We expect the talks to be provocative and bold, while respecting different perspectives.

The coordinated migration of groups of cells underlies many biological processes, including embryo development, wound healing and cancer metastasis. In many of these situations, tissues are able to tune themselves between liquid-like states, where cells flow collectively as in a liquid, and solid-like states that can support shear stresses. In this talk I will describe mesoscopic models of cell assemblies inspired by active matter physics to examine the roles of cell motility, cell crowding and the interplay of contractility and adhesion in controlling the rheological state of biological tissue.

Understanding the transport of microorganisms and self-propelled particles in porous media has important consequences in human health as well as for microbial ecology. In this work, we explore models for the dispersion of active particles in both periodic and random porous media. In a first problem, we analyze the long-time transport properties in a dilute system of active Brownian particles swimming in a periodic lattice in the presence of an external flow. Using generalized Taylor dispersion theory, we calculate the mean transport velocity and dispersion dyadic and explain their dependence on flow strength, swimming activity and geometry. In a second approach, we address the case of run-and-tumble particles swimming through unstructured porous media composed of randomly distributed circular pillars. There, we show that the long-time dispersion is described by a universal hindrance function that depends on the medium porosity and ratio of the swimmer run length to the pillar size. An asymptotic expression for the hindrance function is derived in dilute media, and its extension to semi-dilute and dense media is obtained using stochastic simulations. We conclude by discussing the role of hydrodynamic interactions and swimmer concentration effects.

The ability of the brain to encode and store information depends on the plastic nature of the individual synapses. The increase and decrease in synaptic strength, mediated through the structural plasticity of the spine, are important for learning, memory, and cognitive function. Dendritic spines are small structures that contain the synapse. They come in a variety of shapes (stubby, thin, or mushroom-shaped) and a wide range of sizes that protrude from the dendrite. These spines are the regions where the postsynaptic biochemical machinery responds to the neurotransmitters. Spines are dynamic structures, changing in size, shape, and number during development and aging. While spines and synapses have inspired neuromorphic engineering, the biophysical events underlying synaptic and structural plasticity of single spines remain poorly understood. Our current focus is on understanding the biophysical events underlying structural plasticity. I will discuss recent efforts from my group — first, a systems biology approach to construct a mathematical model of biochemical signaling and actin-mediated transient spine expansion in response to calcium influx caused by NMDA receptor activation and a series of spatial models to study the role of spine geometry and organelle location within the spine for calcium and cyclic AMP signaling. Second, I will discuss how mechanics of membrane-cytoskeleton interactions can give insight into spine shape region. And I will conclude with some new efforts in using reconstructions from electron microscopy to inform computational domains. I will conclude with how geometry and mechanics plays an important role in our understanding of fundamental biological phenomena and some general ideas on bio-inspired engineering.

Intrinsically disordered proteins (IDPs), lack of a well-defined folded structure, have been recently shown to be critical to forming membrane-less organelles via liquid-liquid phase separation (LLPS). Due to the flexible conformations of IDPs, it could be challenging to investigate IDPs with solely experimental techniques. Computational models can therefore provide complementary views at several aspects, including the fundamental physics underlying LLPS and the sequence determinants contributing to LLPS. In this presentation, I will start with our coarse-grained computational framework that can help generate sequence dependent phase diagrams. The coarse-grained model further led to the development of a polymer model with empirical parameters to quickly predict LLPS of IDPs. At last, I will show our preliminary efforts on addressing molecular interactions within LLPS of IDPs using all-atom explicit-solvent simulations.

Many of us came to biology from physics. There we have been trained on such classic examples as muon g-2, where experimental data and theoretical predictions agree to many significant digits. Now, working in biology, we routinely hear that it is messy, most details matter, and that the best hope for theory in biology is to be semi-qualitative, predict general trends, and to forgo the hope of ever making quantitative predictions with the precision that we are used to in physics. Colloquially, we should be satisfied even if data and models differ so much that plotting them on the same plot makes little sense. However, some of us won’t be satisfied by this. So can we take theory in biology seriously and predict experimental outcomes within (small) error bars? Certainly, we won’t be able to predict everything, but this is never required, even in traditional physics. But we should be able to choose some features of data that are nontrivial and interesting, and focus on them. We also should be able to find different classes of models --- maybe even null models --- that match biology better, and thus allow for a better agreement. It is even possible that large-dimensional datasets of modern high-throughput experiments, and the ensuing “more is different” statistical physics style models will make quantitative, precise theory easier. To explore the role of quantitative theory in biology, in this workshop, eight speakers will address some of the following general questions based on their specific work in different corners of biology: Which features of biological data are predictable? Which types of models are best suited to making quantitative predictions in different fields? Should theorists interested in quantitative predictions focus on different questions, not typically asked by biologists? Do large, multidimensional datasets make theories (and which theories?) more or less likely to succeed? This will be an unapologetically theoretical physics workshop — we won’t focus on a specific subfield of biology, but will explore these questions across the fields, hoping that the underlying theoretical frameworks will help us find the missing connections.

The Smoluchowski equation has often been used as the starting point of many continuum models of active suspensions. However, its six-dimensional nature depending on time, space and orientation requires a huge computational cost, fundamentally limiting its use for large-scale problems, such as mixing and transport of active fluids in turbulent flows. Despite the singular nature in strain-dominant flows, the generalised Taylor dispersion (GTD) theory (Frankel & Brenner 1991, J. Fluid Mech. 230:147-181) has been understood to be one of the most promising ways to reduce the Smoluchowski equation into an advection-diffusion equation, the mean drift and diffusion tensor of which rely on ‘local’ flow information only. In this talk, we will introduce an exact transformation of the Smoluchowski equation into such an advection-diffusion equation requiring only local flow information. Based on this transformation, a new advection-diffusion equation will subsequently be proposed by taking an asymptotic analysis in the limit of small particle velocity. With several examples, it will be demonstrated that the new advection-diffusion model, non-singular in strain-dominant flows, outperforms the GTD theory.

Epithelial cell sheets form a fundamental role in the developing embryo, and also in adult tissues including the gut and the cornea of the eye. Soft and active matter provides a theoretical and computational framework to understand the mechanics and dynamics of these tissues.I will start by introducing the simplest useful class of models, active brownian particles (ABPs), which incorporate uncoordinated active crawling over a substrate and mechanical interactions. Using this model, I will show how the extended ’swirly’ velocity fluctuations seen in sheets on a substrate can be understood using a simple model that couples linear elasticity with disordered activity. We are able to quantitatively match experiments using in-vitro corneal epithelial cells.Adding a different source of activity, cell division and apoptosis, to such a model leads to a novel 'self-melting' dense fluid state. Finally, I will discuss a direct application of this simple particle-based model to the steady-state spiral flow pattern on the mouse cornea.

Recently a powerful example of a replicating nano-machine entered our society. In principle, it’s just a normal disease, that one attempts to model with 3 or 4 simple coupled equations with 2 important parameters: a timescale, and a replication factor (the famous R0). Then one tries to guess how changes in society change R0 and perhaps adopt some more or less strong lock-down measures. However, this virus has more “personality” than that. It behaves differently in different persons, and persons behave differently. Presumably, only a few of us infect a lot, while most do not infect so much. This assumption is supported by the observation that couples living together only infect each other with about 15 percent probability, indicating that most infected people are not really infectious. I will discuss this and other aspects of Covid-19 in the perspective of models that describe heterogeneous individuals in a society. In particular, we suggest that limiting superspreading opportunities is a cost-effective strategy to mitigate Covid-19.

We have developed a robotic platform allowing us to monitor cytokines dynamics (including IL-2, IFN-g, TNF, IL-6) of immune cells in vitro, with unprecedented resolution. To understand the complex emerging dynamics, we use interpretable machine learning techniques to build a generative model of cytokine response. We discover that, surprisingly, immune activity is encoded into one global parameter, encoding ligand antigenic properties and to a less extent ligand quantity. Based on this we build a simple interpretable model which can fully explain the broad variability of cytokines dynamics. We validate our approach using different lines of cells and different ligands. Two processes are identified, connected to timing and intensity of cytokine response, which we successfully modulate using drugs or by changing conditions such as initial T cell numbers. Our work reveals a simple "cytokine code", which can be used to better understand immune response in different contexts including immunotherapy. More generally, it reveals how robotic platforms and machine learning can be leveraged to build and validate systems biology models.