Seminars

Through recent advances in microscopy, we now have an unprecedented view of the brain and body of the fruit fly Drosophila melanogaster. We now know the connectivity at single neuron resolution across the whole brain. How do we translate these new measurements into a deeper understanding of how the brain processes sensory information and produces behavior? I will describe two computational efforts to model the brain and the body of the fruit fly. First, I will describe a new modeling method which makes highly accurate predictions of neural activity in the fly visual system as measured in the living brain, using only measurements of its connectivity from a dead brain [1], joint work with Jakob Macke. Second, I will describe a whole body physics simulation of the fruit fly which can accurately reproduce its locomotion behaviors, both flight and walking [2], joint work with Google DeepMind.

Myopia is predicted to affect 50% of all people worldwide by 2050, and is a risk factor for significant, potentially blinding ocular pathologies, such as retinal detachment and glaucoma. Thus, there is significant motivation to better understand the process of myopigenesis and to develop effective anti-myopigenic treatments. In nearly all cases of human myopia, scleral remodeling is an obligate step in the axial elongation that characterizes the condition. Here I will describe the development of a biomechanical assay based on transient unconfined compression of scleral samples. By treating the scleral as a poroelastic material, one can determine scleral biomechanical properties from extremely small samples, such as obtained from the mouse eye. These properties provide proxy measures of scleral remodeling, and have allowed us to identify all-trans retinoic acid (atRA) as a myopigenic stimulus in mice. I will also describe nascent collaborative work on modeling the transport of atRA in the eye.

To effectively move in complex and changing environments, animals must control locomotor speed and gait, while precisely coordinating and adapting limb movements to the terrain. The underlying neuronal control is facilitated by circuits in the spinal cord, which integrate supraspinal commands and afferent feedback signals to produce coordinated rhythmic muscle activations necessary for stable locomotion. I will present a series of computational models investigating dynamics of central neuronal interactions as well as a neuromechanical model that integrates neuronal circuits with a model of the musculoskeletal system. These models closely reproduce speed-dependent gait expression and experimentally observed changes following manipulation of multiple classes of genetically-identified neuronal populations. I will discuss the utility of these models in providing experimentally testable predictions for future studies.

The eye is a multi-component biological system, where mechanics, optics, transport phenomena and chemical reactions are strictly interlaced, characterized by the typical bio-variability in sizes and material properties. The eye’s response to external action is patient-specific and it can be predicted only by a customized approach, that accounts for the multiple physics and for the intrinsic microstructure of the tissues, developed with the aid of forefront means of computational biomechanics. Our activity in the last years has been devoted to the development of a comprehensive model of the cornea that aims at being entirely patient-specific. While the geometrical aspects are fully under control, given the sophisticated diagnostic machinery able to provide a fully three-dimensional images of the eye, the major difficulties are related to the characterization of the tissues, which require the setup of in-vivo tests to complement the well documented results of in-vitro tests. The interpretation of in-vivo tests is very complex, since the entire structure of the eye is involved and the characterization of the single tissue is not trivial. The availability of micromechanical models constructed from detailed images of the eye represents an important support for the characterization of the corneal tissues, especially in the case of pathologic conditions. In this presentation I will provide an overview of the research developed in our group in terms of computational models and experimental approaches developed for the human cornea.

The cytoskeleton has the remarkable ability to self-organize into active materials which underlie diverse cellular processes ranging from motility to cell division. Actomyosin is a canonical example of an active material, which generates cellularscale contractility in part through the forces exerted by myosin motors on actin filaments. While the molecular players underlying actomyosin contractility have been well characterized, how cellular-scale deformation in disordered actomyosin networks emerges from filament-scale interactions is not well understood. In this talk, I’ll present work done in collaboration with Sebastian Fürthauer and Nikta Fakhri addressing this question in vivo using the meiotic surface contraction wave seen in oocytes of the bat star Patiria miniata as a model system. By perturbing actin polymerization, we find that the cellular deformation rate is a nonmonotonic function of cortical actin density peaked near the wild type density. To understand this, we develop an active fluid model coarse-grained from filament-scale interactions and find quantitative agreement with the measured data. The model makes further predictions, including the surprising prediction that deformation rate decreases with increasing motor concentration. We test these predictions through protein overexpression and find quantitative agreement. Taken together, this work is an important step for bridging the molecular and cellular length scales for cytoskeletal networks in vivo.

The rising and setting of the sun is accompanied by changes in both the irradiance and the spectral distribution of the sky. Since the discovery of intrinsically photosensitive retinal ganglion cells (ipRGCs) 20 years ago, considerable progress has been made in understanding melanopsin's contributions to encoding irradiance. Much less is known about the cone inputs to ipRGCs and how they could encode changes in the color of the sky. I will summarize our recent connectomic investigation into the cone-opponent inputs to primate ipRGCs and the implications of this work on our understanding of circadian photoentrainment and the evolution of color vision.

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.

Networks of stiff biopolymers are an omnipresent structural motif in cells and tissues. A prominent modeling framework for describing biopolymer network mechanics is rigidity percolation theory. This theory describes model networks as nodes joined by randomly placed, springlike bonds. Increasing the amount of bonds in a network results in an abrupt, dramatic increase in elastic moduli above a certain threshold – an example of a mechanical phase transition. While homogeneous networks are well studied, many tissues are made of disparate components and exhibit spatial fluctuations in the concentrations of their constituents. In this talk, I will first discuss recent work in which we explained the structural basis of the shear mechanics of healthy and chemically degraded cartilage by coupling a rigidity percolation framework with a background gel. Our model takes into account collagen concentration, as well as the concentration of peptidoglycans in the surrounding polyelectrolyte gel, to produce a structureproperty relationship that describes the shear mechanics of both sound and diseased cartilage. I will next discuss the introduction of structural correlation in constructing networks, such that sparse and dense patches emerge. I find moderate correlation allows a network to become rigid with fewer bonds, while this benefit is partly erased by excessive correlation. We explain this phenomenon through analysis of the spatial fluctuations in strained networks’ displacement fields. Finally, I will address our work’s implications for non-invasive diagnosis of pathology, as well as rational design of prostheses and novel soft materials.

Tissue folding is a ubiquitous shape change event during development whereby a cell sheet bends into a curved 3D structure. This mechanical process is remarkably robust, and the correct final form is almost always achieved despite internal fluctuations and external perturbations inherent in living systems. While many genetic and molecular strategies that lead to robust development have been established, much less is known about how mechanical patterns and movements are ensured at the population level. I will describe how quantitative imaging, physical modeling and concepts from network science can uncover collective interactions that govern tissue patterning and shape change. Actin and myosin are two important cytoskeletal proteins involved in the force generation and movement of cells. Both parts of this talk will be about the spontaneous organization of actomyosin networks and their role in collective tissue dynamics. First, I will present how out-of-plane curvature can trigger the global alignment of actin fibers and a novel transition from collective to individual cell migration in culture. I will then describe how tissue-scale cytoskeletal patterns can guide tissue folding in the early fruit fly embryo. I will show that actin and myosin organize into a network that spans a domain of the embryo that will fold. Redundancy in this supracellular network encodes the tissue’s intrinsic robustness to mechanical and molecular perturbations during folding.

Hydra is an extraordinary creature. Continuously replacing itself, it can live indefinitely, performing a stable repertoire of reasonably sophisticated behaviors. This remarkable stability under plasticity may be due to the uniform nature of its nervous system, which consists of two apparently noncommunicating nerve net layers. We use modeling to understand the role of active muscles and biomechanics interact with neural activity to shape Hydra behaviour. We will discuss our findings and thoughts on how this simple nervous system may self-organize to produce purposeful behavior.

Filaments (slender, microscopic elastic bodies) are prevalent in biological and industrial settings. In the biological case, the filaments are often active, in that they are driven internally by motor proteins, with the prime examples being cilia and flagella. For cilia in particular, which can appear in dense arrays, their resulting motions are coupled through the surrounding fluid, as well as through surfaces to which they are attached. In this talk, I present numerical simulations exploring the coordinated motion of active filaments and how it depends on the driving force, density of filaments, as well as the attached surface. In particular, we find that when the surface is spherical, its topology introduces local defects in coordinated motion which can then feedback and alter the global state. This is particularly true when the surface is not held fixed and is free to move in the surrounding fluid. These simulations take advantage of a computational framework we developed for fully 3D filament motion that combines unit quaternions, implicit geometric time integration, quasi-Newton methods, and fast, matrix-free methods for hydrodynamic interactions and it will also be presented.

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 journey of development begins with sperm swimming through the female reproductive tract en-route to the egg. In order to successfully complete this journey sperm must beat a single flagellum, propelling themselves through a wide range of fluids, from liquified semen to viscous cervical mucus. It is well-known that the beating tail is driven by an array of 9 microtubule doublets surrounding a central pair, with interconnecting dynein motors generating shear forces and driving elastic wave propagation. Despite this knowledge, the exact mechanism by which coordination of these motors drives oscillating waves along the flagellum remains unknown; hypothesised mechanisms include curvature control, sliding control, and geometric clutch. In this talk we will discuss the mechanisms of flagellar bending, and present a simple model of active curvature that is able to produce many of the various sperm waveforms that are seen experimentally, including those in low and high viscosity fluids and after a cell has ‘hyperactivated’ (a chemical process thought to be key for fertilization). We will show comparisons between these simulated waveforms and sperm that have been experimentally tracked, and discuss methods for fitting simulated mechanistic parameters to these real cells.

Embryo morphogenesis is impacted by dynamic changes in tissue material properties, which have been proposed to occur via processes akin phase transitions (PTs). Here, we show that rigidity percolation provides a simple and robust theoretical framework to predict material/structural PTs of embryonic tissues from local cell connectivity. By using percolation theory, combined with directly monitoring dynamic changes in tissue rheology and cell contact mechanics, we demonstrate that the zebrafish blastoderm undergoes a genuine rigidity PT, brought about by a small reduction in adhesion-dependent cell connectivity below a critical value. We quantitatively predict and experimentally verify hallmarks of PTs, including power-law exponents and associated discontinuities of macroscopic observables at criticality. Finally, we show that this uniform PT depends on blastoderm cells undergoing meta-synchronous divisions causing random and, consequently, uniform changes in cell connectivity. Collectively, our theoretical and experimental findings reveal the structural basis of material PTs in an organismal context.

'Shape' of microorganisms are diverse. However, we sometimes approximate them as a sphere or a spheroid when we mathematically model the hydrodynamics of motile and non-motile cells. Such a geometrical simplification can be theoretically validated for motions in a linear background flow, since the dynamics, known as the Jeffery orbit, only contain a single geometric parameter, called the Bretherton constant. In this talk, we generalise the Jeffery equations for a chiral axisymmetric object using the low-Reynolds-number hydrokinetic symmetry and then demonstrate that the dynamics of a certain type of chiral object in a fluid flow are characterised by a new chiral parameter in addition to the Bretherton constant. We also discuss how the generalised Jeffery orbits are applied to biased locomotion of bacteria in a bulk shear flow and we will share the idea of hydrodynamic `shape' of microorganisms to simplify the description of their dynamics.

Microorganisms can swim in a variety of environments, interacting with chemicals and other proteins in the fluid. In this talk, we will highlight recent computational methods and results for swimming efficiency and hydrodynamic interactions of swimmers in different fluid environments. Sperm are modeled via a centerline representation where forces are solved for using elastic rod theory. The method of regularized Stokeslets is used to solve the fluid-structure interaction where emergent swimming speeds can be compared to asymptotic analysis. In the case of fluids with extra proteins or cells that may act as friction, swimming speeds may be enhanced, and attraction may not occur. We will also highlight how parameter estimation techniques can be utilized to infer fluid and/or swimmer properties.

A key challenge in modern soft matter is to identify the principles that govern the organisation and functionality in non-equilibrium systems. Current research efforts largely focus on non-equilibrium processes that occur either at the single-molecule scale (e.g. protein and DNA conformations under driving forces), or at the scale of whole tissues, organisms, and active colloidal and microscopic objects. However, the range of the scales in-between — from molecules to large-scaled molecular assemblies that consume energy and perform work — remains under-explored. This is, nevertheless, the scale that is crucial for the function of a living cell, where molecular self-assembly driven far from equilibrium produces mechanical work needed for cell reshaping, transport, motility, division, and healing. Today I will discuss physical modelling of active elastic filaments, called ESCRT-III filaments, that dynamically assemble and disassemble on cell membranes. This dynamic assembly changes the filaments’ shape and mechanical properties and leads to the remodelling and cutting of cells. I will present a range of experimental comparisons of our simulation results: from ESCRT-III-driven trafficking in eukaryotes to division of evolutionary simple archaeal cells.

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 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.