Showing papers in "arXiv: Biological Physics in 2019"

Physical Models of Collective Cell Migration

Abstract: Collective cell migration is a key driver of embryonic development, wound healing, and some types of cancer invasion. Here we provide a physical perspective of the mechanisms underlying collective cell migration. We begin with a catalogue of the cell-cell and cell-substrate interactions that govern cell migration, which we classify into positional and orientational interactions. We then review the physical models that have been developed to explain how these interactions give rise to collective cellular movement. These models span the sub-cellular to the supracellular scales, and they include lattice models, phase fields models, active network models, particle models, and continuum models. For each type of model, we discuss its formulation, its limitations, and the main emergent phenomena that it has successfully explained. These phenomena include flocking and fluid-solid transitions, as well as wetting, fingering, and mechanical waves in spreading epithelial monolayers. We close by outlining remaining challenges and future directions in the physics of collective cell migration.

Active wetting of epithelial tissues

Abstract: Development, regeneration and cancer involve drastic transitions in tissue morphology. In analogy with the behavior of inert fluids, some of these transitions have been interpreted as wetting transitions. The validity and scope of this analogy are unclear, however, because the active cellular forces that drive tissue wetting have been neither measured nor theoretically accounted for. Here we show that the transition between 2D epithelial monolayers and 3D spheroidal aggregates can be understood as an active wetting transition whose physics differs fundamentally from that of passive wetting phenomena. By combining an active polar fluid model with measurements of physical forces as a function of tissue size, contractility, cell-cell and cell-substrate adhesion, and substrate stiffness, we show that the wetting transition results from the competition between traction forces and contractile intercellular stresses. This competition defines a new intrinsic lengthscale that gives rise to a critical size for the wetting transition in tissues, a striking feature that has no counterpart in classical wetting. Finally, we show that active shape fluctuations are dynamically amplified during tissue dewetting. Overall, we conclude that tissue spreading constitutes a prominent example of active wetting --- a novel physical scenario that may explain morphological transitions during tissue morphogenesis and tumor progression.

Lotka-Volterra representation of general nonlinear systems

Abstract: In this paper we elaborate on the structure of the Generalized Lotka-Volterra (GLV) form for nonlinear differential equations. We discuss here the algebraic properties of the GLV family, such as the invariance under quasimonomial transformations and the underlying structure of classes of equivalence. Each class possesses a unique representative under the classical quadratic Lotka-Volterra form. We show how other standard modelling forms of biological interest, such as S-systems or mass-action systems are naturally embedded into the GLV form, which thus provides a formal framework for their comparison, and for the establishment of transformation rules. We also focus on the issue of recasting of general nonlinear systems into the GLV format. We present a procedure for doing so, and point at possible sources of ambiguity which could make the resulting Lotka-Volterra system dependent on the path followed. We then provide some general theorems that define the operational and algorithmic framework in which this is not the case.

Three-dimensional view of ultrafast dynamics in photoexcited bacteriorhodopsin

Abstract: Bacteriorhodopsin (bR) is a light-driven proton pump. We use time-resolved crystallography at an X-ray free-electron laser to follow the structural changes in multiphoton-excited bR from 250 femtoseconds to 10 picoseconds. Quantum chemistry and ultrafast spectroscopy allow identifying a sequential two-photon absorption process, leading to excitation of a tryptophan residue flanking the retinal chromophore, as a first manifestation of multi-photon effects. We resolve distinct stages in the structural dynamics of the all-trans retinal in photoexcited bR to a highly twisted 13-cis conformation. Other active site sub-picosecond rearrangements include correlated vibrational motions of the electronically excited retinal chromophore, the surrounding amino acids and water molecules as well as their hydrogen bonding network. These results show that this extended photo-active network forms an electronically and vibrationally coupled system in bR, and most likely in all retinal proteins.

Terradynamically streamlined shapes in animals and robots enhances traversability through densely cluttered terrain.

Abstract: Many animals, modern aircraft, and underwater vehicles use streamlined body shapes that reduce fluid dynamic drag to achieve fast and effective locomotion in air and water. Similarly, numerous small terrestrial animals move through cluttered terrain where 3-D, multi-component obstacles like grass, shrubs, vines, and leaf litter resist motion, but it is unknown whether their body shape plays a major role in traversal. Few ground vehicles or terrestrial robots have used body shape to effectively traverse cluttered terrain. Here, we challenged forest-floor-dwelling discoid cockroaches possessing a thin, rounded body to traverse tall, narrowly spaced, vertical, grass-like compliant beams. Animals displayed high traversal performance (79 +/- 12% probability and 3.4 +/- 0.7 s time). Although we observed diverse traversal strategies, cockroaches primarily (48 +/- 9 % probability) used a novel roll maneuver, allowing them to rapidly traverse obstacle gaps narrower than half body width (2.0 +/- 0.5 s traversal time). Reduction of body roundness by addition of artificial shells nearly inhibited roll maneuvers and decreased traversal performance. Inspired by this discovery, we added a thin, rounded exoskeletal shell to a legged robot with a nearly cuboidal body, common to many existing terrestrial robots. Without adding sensory feedback or changing the open-loop control, the rounded shell enabled the robot to traverse beam obstacles with gaps narrower than shell width via body roll. Terradynamically streamlined shapes can reduce terrain resistance and enhance traversability by assisting effective body reorientation via distributed mechanical feedback. Our findings highlight the need to consider body shape to improve robot mobility in real-world terrain often filled with clutter, and to develop better locomotor-ground contact models to understand interaction with complex terrain.

Machine learning for protein folding and dynamics

Abstract: Many aspects of the study of protein folding and dynamics have been affected by the recent advances in machine learning. Methods for the prediction of protein structures from their sequences are now heavily based on machine learning tools. The way simulations are performed to explore the energy landscape of protein systems is also changing as force-fields are started to be designed by means of machine learning methods. These methods are also used to extract the essential information from large simulation datasets and to enhance the sampling of rare events such as folding/unfolding transitions. While significant challenges still need to be tackled, we expect these methods to play an important role on the study of protein folding and dynamics in the near future. We discuss here the recent advances on all these fronts and the questions that need to be addressed for machine learning approaches to become mainstream in protein simulation.

OnionNet: a multiple-layer inter-molecular contact based convolutional neural network for protein-ligand binding affinity prediction

Abstract: Computational drug discovery provides an efficient tool helping large scale lead molecules screening. One of the major tasks of lead discovery is identifying molecules with promising binding affinities towards a target, a protein in general. The accuracies of current scoring functions which are used to predict the binding affinity are not satisfactory enough. Thus, machine learning (ML) or deep learning (DL) based methods have been developed recently to improve the scoring functions. In this study, a deep convolutional neural network (CNN) model (called OnionNet) is introduced and the features are based on rotation-free element-pair specific contacts between ligands and protein atoms, and the contacts were further grouped in different distance ranges to cover both the local and non-local interaction information between the ligand and the protein. The prediction power of the model is evaluated and compared with other scoring functions using the comparative assessment of scoring functions (CASF-2013) benchmark and the v2016 core set of PDBbind database. When compared to a previous CNN-based scoring function, our model shows improvements of 0.08 and 0.16 in the correlations (R) and standard deviations (SD) of regression, respectively, between the predicted binding affinities and the experimental measured binding affinities. The robustness of the model is further explored by predicting the binding affinities of the complexes generated from docking simulations instead of experimentally determined PDB structures.

Bacteria display optimal transport near surfaces -- bacteria as intermittent active chiral particles: trapped by hydrodynamics, escaping by adhesion

Abstract: The near-surface swimming patterns of bacteria are strongly determined by the hydrodynamic interactions between bacteria and the surface, which trap bacteria in smooth circular trajectories that lead to inefficient surface exploration. Here, we show by combining experiments and a data-driven mathematical model that surface exploration of enterohemorrhagic Escherichia coli (EHEC) -- a pathogenic strain of E. coli causing serious illnesses such as bloody diarrhea -- results from a complex interplay between motility and transient surface adhesion events. These events allow EHEC to break the smooth circular trajectories and regulate their transport properties by the use stop-adhesion events that lead to a characteristic intermittent motion on surfaces. We find that the experimentally measured frequency of stop-adhesion events in EHEC is located at the value predicted by the developed mathematical model that maximizes bacterial surface diffusivity. We indicate that these results and the developed model apply to other bacterial strains on different surfaces, which suggests that swimming bacteria use transient adhesion to regulate surface motion.

Brillouin Light Scattering: Applications in Biomedical Sciences

Abstract: Brillouin spectroscopy and imaging are emerging techniques in analytical science, biophotonics and biomedicine. They are based on Brillouin light scattering from acoustic waves or phonons in the GHz range, providing a nondestructive contactless probe of the mechanics on a micro-scale. Novel approaches and applications of these techniques to the field of biomedical sciences are discussed highlighting the theoretical foundations and experimental methods that have been developed to date, acknowledging that this is a fast moving field and a comprehensive account of the relevant literature is critically assessed here.

Multiscale modeling meets machine learning: What can we learn?

Abstract: Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learning to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

Magnetic resonance imaging of mean cell size in human breast tumors

Abstract: Purpose: Cell size is a fundamental characteristic of all tissues, and changes in cell size in cancer reflect tumor status and response to treatments, such as apoptosis and cell cycle arrest. Unfortunately, cell size can only be obtained by pathologic evaluation of the tumor in the current standard of care. Previous imaging approaches can be implemented on only animal MRI scanners or require relatively long acquisition times that are undesirable for clinical imaging. There is a need to develop cell size imaging for clinics. Experimental Design: We propose a new method, IMPULSED (Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion) that can characterize mean cell sizes in solid tumors. We report the use of combined sequences with different gradient waveforms on human MRI and analytical equations that link DWI signals of real gradient waveforms and specific microstructural parameters such as cell size. We also describe comprehensive validations using computer simulations, cell experiments in vitro, and animal experiments in vivo and finally demonstrate applications in pre-operative breast cancer patients. Results: With fast acquisitions (~ 7 mins), IMPULSED can provide high-resolution (1.3 mm in-plane) mapping of mean cell size of human tumors in vivo on currently-available 3T MRI scanners. All validations suggest IMPULSED provide accurate and reliable measurements of mean cell size. Conclusion: The proposed IMPULSED method can assess cell size variations in the tumor of breast cancer patients, which may have the potential to assess early response to neoadjuvant therapy.

Timing of fungal spore release dictates survival during atmospheric transport

Abstract: The fungi disperse spores to move across landscapes and spore liberation takes different patterns. While many species release spores intermittently, others release spores at specific times of day or night according to intrinsic rhythms. Despite intriguing evidence of diurnal rhythms, why the timing of spore liberation would matter to a fungus remains an open question. Here we use state-of-the-art numerical simulations of atmospheric transport with meteorological data to follow the trajectory of many spores released in the open atmosphere at different times of day, during different seasons and at different locations across North America. While individual spores follow unpredictable trajectories due to turbulence, in the aggregate patterns emerge: statistically, spores released during the day fly for several days, while spores released at night return to the ground within a few hours. Differences are caused by intense turbulence during the day and weak turbulence at night. The pattern is widespread but its reliability varies, for example, day/night patterns are stronger in southern regions, where temperatures are warmer. Results provide a set of testable hypotheses explaining intermittent and regular patterns of spore release as strategies to maximize spore survival in the air. Species with short lived spores reproducing where there is strong and regular turbulence during the day, for example in Mexico, will maximize survival by routinely releasing spores at night. Where cycles are weak, for example in Canada during spring, there will be no benefit to releasing spores at the same time every day. We also challenge the perception of atmospheric dispersal as risky, wasteful, and beyond control of a sporocarp; our data suggest the timing of spore liberation may be finely tuned by a fungus to maximize fitness during atmospheric transport.

Spatially Heterogeneous Dynamics of Cells in a Growing Tumor Spheroid: Comparison Between Theory and Experiments

Abstract: Collective cell movement, characterized by multiple cells that are in contact for substantial periods of time and undergo correlated motion, plays a central role in cancer and embryogenesis. Recent imaging experiments have provided time-dependent traces of individual cells, thus providing an unprecedented picture of tumor spheroid growth. By using simulations of a minimal cell model, we analyze the experimental data that map the movement of cells in fibrosarcoma tumor spheroid embedded in a collagen matrix. Both simulations and experiments show that cells in the core of the spheroid exhibit subdiffusive glassy dynamics (mean square displacement, $\Delta(t) \approx t^{\alpha}$ with $\alpha 1)$. The motion of most of the cells near the periphery undergo highly persistent and correlated directional motion, thus explaining the observed superdiffusive behavior. The $\alpha$ values for cells in the core and periphery, extracted from simulations and experiments are in near quantitative agreement with each other, which is surprising given that no parameter in the model was used to fit the measurements. The qualitatively different dynamics of cells in the core and periphery is captured by the fourth order susceptibility, introduced to characterize metastable states in glass forming systems. Analyses of the velocity autocorrelation of individual cells show remarkable spatial heterogeneity with no two cells exhibiting similar behavior. The prediction that $\alpha$ should depend on the location of the cells in the tumor is amenable to experimental test. The highly heterogeneous dynamics of cells in the tumor spheroid provides a plausible mechanism for the origin of intratumor heterogeneity.

\textit{C. elegans} collectively forms dynamical networks

Abstract: Understanding physical rules underlying collective motions requires perturbation of controllable parameters in self-propelled particles. However, controlling parameters in animals is generally not easy, which makes collective behaviours of animals elusive. Here, we report an experimental system in which a conventional model animal, \textit {Caenorhabditis elegans}, collectively forms dynamical networks of bundle-shaped aggregates. We investigate the dependence of our experimental system on various extrinsic parameters (material of substrate, ambient humidity and density of worms). Taking advantage of well-established \textit {C.~elegans} genetics, we also control the intrinsic parameters (genetically determined motility) by mutations and by forced neural activation via optogenetics. Furthermore, we develop a minimal agent-based model that reproduces the dynamical network formation and its dependence on the parameters, suggesting that the key factors are alignment of worms after collision and smooth turning. Our findings imply that the concepts of active matter physics may help us to understand biological functions of animal groups.

Nanopore fabrication via tip-controlled local breakdown using an atomic force microscope

Abstract: The dielectric breakdown approach for forming nanopores has greatly accelerated the pace of research in solid-state nanopore sensing, enabling inexpensive formation of nanopores via a bench top setup. Here we demonstrate the potential of tip controlled dielectric breakdown (TCLB) to fabricate pores 100$\times$ faster, with high scalability and nanometre positioning precision. A conductive atomic force microscope (AFM) tip is brought into contact with a nitride membrane positioned above an electrolyte reservoir. Application of a voltage pulse at the tip leads to the formation of a single nanoscale pore. Pores are formed precisely at the tip position with a complete suppression of multiple pore formation. In addition, our approach greatly accelerates the electric breakdown process, leading to an average pore fabrication time on the order of 10 ms, at least 2 orders of magnitude shorter than achieved by classic dielectric breakdown approaches. With this fast pore writing speed we can fabricate over 300 pores in half an hour on the same membrane.

High-resolution Markov state models for the dynamics of Trp-cage miniprotein constructed over slow folding modes identified by state-free reversible VAMPnets.

Abstract: State-free reversible VAMPnets (SRVs) are a neural network-based framework capable of learning the leading eigenfunctions of the transfer operator of a dynamical system from trajectory data. In molecular dynamics simulations, these data-driven collective variables (CVs) capture the slowest modes of the dynamics and are useful for enhanced sampling and free energy estimation. In this work, we employ SRV coordinates as a feature set for Markov state model (MSM) construction. Compared to the current state of the art, MSMs constructed from SRV coordinates are more robust to the choice of input features, exhibit faster implied timescale convergence, and permit the use of shorter lagtimes to construct higher kinetic resolution models. We apply this methodology to study the folding kinetics and conformational landscape of the Trp-cage miniprotein. Folding and unfolding mean first passage times are in good agreement with prior literature, and a nine macrostate model is presented. The unfolded ensemble comprises a central kinetic hub with interconversions to several metastable unfolded conformations and which serves as the gateway to the folded ensemble. The folded ensemble comprises the native state, a partially unfolded intermediate "loop" state, and a previously unreported short-lived intermediate that we were able to resolve due to the high time-resolution of the SRV-MSM. We propose SRVs as an excellent candidate for integration into modern MSM construction pipelines.

Interfacing Graphene-Based Materials With Neural Cells

Abstract: The scientific community has witnessed an exponential increase in the applications of graphene and graphene-based materials in a wide range of fields. For what concerns neuroscience, the interest raised by these materials is two-fold. On one side, nanosheets made of graphene or graphene derivatives (graphene oxide, or its reduced form) can be used as carriers for drug delivery. Here, an important aspect is to evaluate their toxicity, which strongly depends on flake composition, chemical functionalization and dimensions. On the other side, graphene can be exploited as a substrate for tissue engineering. In this case, conductivity is probably the most relevant amongst the various properties of the different graphene materials, as it may allow to instruct and interrogate neural networks, as well as to drive neural growth and differentiation. In this review, we try to give a comprehensive view of the accomplishments and new challenges of the field, as well as which in our view are the most exciting directions to take in the immediate future. These include the need to engineer multifunctional nanoparticles able to cross the blood-brain-barrier to reach neural cells, and to achieve on-demand delivery of specific drugs. We describe the state-of-the-art in the use of graphene materials to engineer three-dimensional scaffolds to drive neuronal growth and regeneration in vivo, and the possibility of using graphene as a component of hybrid composites/multi-layer organic electronics devices. Last but not least, we address the need of an accurate theoretical modeling of the interface between graphene and biological material, by modeling the interaction of graphene with proteins and cell membranes at the nanoscale, and describing the physical mechanism(s) of charge transfer by which the various graphene materials can influence the excitability and physiology of neural cells.

MnO2-gated Nanoplatforms with Targeted Controlled Drug Release and Contrast-Enhanced MRI Properties: from 2D Cell Culture to 3D Biomimetic Hydrogels.

Abstract: Multifunctional nanomaterials combining diagnosis and therapeutic properties have attracted a considerable attention in cancer research. Yet some important challenges are still to be faced, including an optimal coupling between these two types of properties that would be effective within complex biological tissues. To address these points, we have prepared novel nanoplatforms associating controlled drug delivery of doxorubicin and Magnetic Resonance Imaging (MRI) contrast-enhancement that exhibit high specificity towards cancer cells compared to normal cells and evaluated them both in 2D cultures and within 3D tissue-like biomimetic matrices. Methods: Nanoplatforms were prepared from hollow silica nanoparticles coated with MnO2 nanosheets and conjugated with the AS1411 aptamer as a targeting agent. They were fully characterized from a chemical and structural point of view as well as for drug release and MRI signal enhancement. Standard two-dimensional monolayer cultures were performed using HeLa and Normal Human Dermal Fibroblasts (NHDF) cells to testify targeting and cytotoxicity. Cellularized type I collagen-based hydrogels were also used to study nanoparticles behavior in 3D biomimetic environments. Results: The as-established nanoplatforms can enter HeLa cells, leading to the dissociation of the MnO2 nanosheets into Mn 2+ that enhanced T1 magnetic resonance signals and concomitantly release doxorubicin, both effects being markedly more significant than in the presence of NHDFs. Moreover, particles functionality and specificity were preserved when the cells were immobilized within type I collagen-based fibrillar hydrogels. Conclusion: The use of MnO2 nanosheets as glutathione-sensitive coatings of drug loaded nanoparticles together with surface conjugation with a targeting aptamer offers an effective strategy to obtain efficient and specific nanotheranostic systems for cancer research, both in 2D and 3D. The here-described tissue-like models should be easy to implement and could constitute an interesting intermediate validation step for newly-developed theranostic nanoparticles before in vivo evaluation.

Label-free biochemical quantitative phase imaging with mid-infrared photothermal effect

Abstract: Label-free optical imaging is valuable in biology and medicine with its non-destructive property and reduced optical and chemical damages. Quantitative phase (QPI) and molecular vibrational imaging (MVI) are the two most successful label-free methods, providing morphology and biochemistry, respectively, that have pioneered numerous applications along their independent technological maturity over the past few decades. However, the distinct label-free contrasts are inherently complementary and difficult to integrate due to the use of different light-matter interactions. Here, we present a unified imaging scheme that realizes simultaneous and in-situ acquisition of MV-fingerprint contrasts of single cells in the framework of QPI utilizing the mid-infrared photothermal effect. The fully label-free and robust integration of subcellular morphology and biochemistry would have important implications, especially for studying complex and fragile biological phenomena such as drug delivery, cellular diseases and stem cell development, where long-time observation of unperturbed cells are needed under low phototoxicity.

Dynamics of a large population of red blood cells under shear flow

Abstract: An exhaustive description of the dynamics under shear flow of a large number of red blood cells in dilute regime is proposed, which highlights and takes into account the dispersion in cell properties within a given blood sample. Physiological suspending fluid viscosity is considered, a configuration surprisingly seldom considered in experimental studies, as well as a more viscous fluid that is a reference in the literature. Stable and unstable flipping motions well described by Jeffery orbits or modified Jeffery orbits are identified, as well as transitions to and from tank-treading motion in the more viscous suspending fluid case. Hysteresis loops upon shear rate increase or decrease are highlighted for the transitions between unstable and stable orbits as well as for the transition between flipping and tank-treading. We identify which of the characteristic parameters of motion and of the transition thresholds depend on flow stress only or also on suspending fluid viscosity.

Geometric phase and dimensionality reduction in locomoting living systems

Abstract: The apparent ease with which animals move requires the coordination of their many degrees of freedom to manage and properly utilize environmental interactions. Identifying effective strategies for locomotion has proven challenging, often requiring detailed models that generalize poorly across modes of locomotion, body morphologies, and environments. We present the first biological application of a gauge-theory-based geometric framework for movement, originally proposed by Wilczek and Shapere nearly $40$ years ago, to describe self-deformation-driven movements through dissipative environments. Using granular resistive force theory to model environmental forces and principal components analysis to identify a low-dimensional space of animal postures and dynamics, we show that our approach captures key features of how a variety of animals, from undulatory swimmers and slitherers to sidewinding rattlesnakes, coordinate body movements and leverage environmental interactions to generate locomotion. Our results demonstrate that this geometric approach is a powerful and general framework that enables the discovery of effective control strategies, which could be further augmented by physiologically-relevant parameters and constraints to provide a deeper understanding of locomotion in a wide variety of biological systems and environments.

A mechanical model reveals that non-axisymmetric buckling lowers the energy barrier associated with membrane neck constriction

Abstract: Membrane neck formation is essential for scission, which, as recent experiments on tubules have demonstrated, can be location dependent. The diversity of biological machinery that can constrict a neck such as dynamin, actin, ESCRTs and BAR proteins, and the range of forces and deflection over which they operate, suggest that the constriction process is functionally mechanical and robust to changes in biological environment. In this study, we used a mechanical model of the lipid bilayer to systematically investigate the influence of location, symmetry constraints, and helical forces on membrane neck constriction. Simulations from our model demonstrated that the energy barriers associated with constriction of a membrane neck are location-dependent. Importantly, if symmetry restrictions are relaxed, then the energy barrier for constriction is dramatically lowered and the membrane buckles at lower values of forcing parameters. Our simulations also show that constriction due to helical proteins further reduces the energy barrier for neck formation compared to cylindrical proteins. These studies establish that despite different molecular mechanisms of neck formation in cells, the mechanics of constriction naturally leads to a loss of symmetry that can lower the energy barrier to constriction.

Flagellar number governs bacterial spreading and transport efficiency

Abstract: Peritrichous bacteria synchronize and bundle their flagella to actively swim while disruption of the bundle leads to tumbling. It is still not known whether the number of flagella represents an evolutionary adaptation towards optimizing bacterial navigation. Here, we study the swimming dynamics of differentially flagellated Bacillus subtilis strains in a quasi-two-dimensional system. We find that decreasing the number of flagella N reduces the average turning angle between two successive run phases and enhances the duration and directional persistence of the run phase. As a result, having less flagella is beneficial for long-distance transport and fast spreading, while having a lot of flagella is advantageous for the processes which require localization and slow dynamics, such as biofilm formation. We develop a two-state random walk model that incorporates spontaneous switchings between the states and yields exact analytical expressions for transport properties, in remarkable agreement with experiments. The results of numerical simulations based on our two-state model suggest that the efficiency of searching and exploring the environment is optimized at intermediate values of N. The optimal choice of N, for which the search time is minimized, decreases with increasing the size of the environment in which the bacteria swim.

Pairing, waltzing and scattering of chemotactic active colloids.

Abstract: An interacting pair of chemotactic (anti-chemotactic) active colloids, that can rotate their axes of self-propulsion to align {parallel (anti-parallel)} to a chemical gradient, shows dynamical behaviour that varies from bound states to scattering. The underlying two-body interactions are purely dynamical, non-central, non-reciprocal, and controlled by changing the catalytic activity and phoretic mobility. Mutually chemotactic colloids trap each other in a final state of fixed separation; the resulting `active dimer' translates. A second type of bound state is observed where the polar axes undergo periodic cycles leading to phase-synchronised circular motion around a common point. These bound states are formed depending on initial conditions and can unbind on increasing the speed of self propulsion. Mutually anti-chemotactic swimmers always scatter apart. We also classify the fixed points underlying the bound states, and the bifurcations leading to transitions from one type of bound state to another, for the case of a single swimmer in the presence of a localised source of solute.

Collective Cell Mechanics of Small-Organoid Morphologies

Abstract: The study of organoids, artificially grown cell aggregates with the functionality and small-scale anatomy of real organs, is one of the most active areas of research in biology and biophysics, yet the basic physical origins of their different morphologies remain poorly understood. Here we propose a mechanistic theory of small-organoid morphologies. Using a 3D surface-tension-based vertex model, we reproduce the characteristic shapes, ranging from branched and budded structures to invaginated shapes. We find that the formation of branched morphologies relies strongly on junctional activity, enabling temporary aggregations of topological defects in cell packing. To elucidate our numerical results, we develop an effective elasticity theory, which allows one to estimate the apico-basal polarity from the organoid-scale modulation of cell height. Our work provides a generic interpretation of the observed small-organoid morphologies, highlighting the role of physical factors such as the differential surface tension, cell rearrangements, and tissue growth.

High Sensitivity and Specificity Biomechanical Imaging by Stimulated Brillouin Scattering Microscopy.

Abstract: Noncontact label-free biomechanical imaging is a crucial tool for unraveling the mechanical properties of biological systems, which play critical roles in the fields of engineering, physics, biology and medicine; yet, it represents a significant challenge in microscopy. Spontaneous Brillouin microscopy meets this challenge, but often requires long acquisition times or lacks high specificity for detecting biomechanical constituents with highly overlapping Brillouin bands. We developed stimulated Brillouin scattering (SBS) microscopy that provides intrinsic noncontact biomechanical contrast and generates mechanical cross-sectional images inside large specimens, with high mechanical specificity and pixel dwell times that are >10-fold improved over those of spontaneous Brillouin microscopy. We used SBS microscopy in different biological applications, including the quantification of the high-frequency complex longitudinal modulus of the pharyngeal region of live wild-type Caenorhabditis elegans nematodes, imaging of the variations in the high-frequency viscoelastic response to osmotic stress in the head of living worms, and in vivo mechanical contrast mesoscopy of developing nematodes

Cells exploit a phase transition to mechanically remodel the fibrous extracellular matrix

Abstract: Through mechanical forces, biological cells remodel the surrounding collagen network, generating striking deformation patterns. Tethers-tracts of high densification and fiber alignment-form between cells, thinner bands emanate from cell clusters. While tethers facilitate cell migration and communication, how they form is unclear. Combining modeling, simulation and experiment, we show that tether formation is a densification phase transition of the extracellular matrix, caused by buckling instability of network fibers under cell-induced compression, featuring unexpected similarities with martensitic microstructures. Multiscale averaging yields a two-phase, bistable continuum energy landscape for fibrous collagen, with a densified/aligned second phase. Simulations predict strain discontinuities between the undensified and densified phase, which localizes within tethers as experimentally observed. In our experiments, active particles induce similar localized patterns as cells. This shows how cells exploit an instability to mechanically remodel the extracellular matrix simply by contracting, thereby facilitating mechanosensing, invasion and metastasis.

Brillouin-derived viscoelastic parameters of hydrogel tissue models

Abstract: Many problems in mechanobiology urgently require characterisation of the micromechanical properties of the fibrous proteins of cells and tissues. Brillouin light scattering has been proposed as a new optical elastography technique to meet this need, but the information contained in the Brillouin spectrum is still a matter of debate. Here we investigate this question using gelatin as a model system in which the macroscopic physical properties can be manipulated to mimic all the relevant biological states of matter, ranging from the liquid to the gel and the glassy phase. We demonstrate that Brillouin spectroscopy is a powerful tool to reveal both the elastic and viscous properties of biomaterials that are central to their biological functions.

Magic numbers in polymer phase separation -- the importance of being rigid

Abstract: Cells possess non-membrane-bound bodies, many of which are now understood as phase-separated condensates. One class of such condensates is composed of two polymer species, where each consists of repeated binding sites that interact in a one-to-one fashion with the binding sites of the other polymer. Previous biologically-motivated modeling of such a two-component system surprisingly revealed that phase separation is suppressed for certain combinations of numbers of binding sites. This phenomenon, dubbed the "magic-number effect", occurs if the two polymers can form fully-bonded small oligomers by virtue of the number of binding sites in one polymer being an integer multiple of the number of binding sites of the other. Here we use lattice-model simulations and analytical calculations to show that this magic-number effect can be greatly enhanced if one of the polymer species has a rigid shape that allows for multiple distinct bonding conformations. Moreover, if one species is rigid, the effect is robust over a much greater range of relative concentrations of the two species. Our findings advance our understanding of the fundamental physics of two-component polymer-based phase-separation and suggest implications for biological and synthetic systems.

Centering and symmetry breaking in confined contracting actomyosin networks

Abstract: Centering and decentering of cellular components is essential for internal organization of cells and their ability to perform basic cellular functions such as division and motility. How cells achieve proper localization of their components is still not well-understood, especially in large cells such as oocytes. Here, we study actin-based positioning mechanisms in artificial cells with persistently contracting actomyosin networks, generated by encapsulating cytoplasmic Xenopus egg extracts into cell-sized water-in-oil droplets. We observe size-dependent localization of the contraction center, with a symmetric configuration in larger cells and a polar one in smaller cells. In the symmetric state, the contraction center is actively centered, via a hydrodynamic mechanism based on Darcy friction between the contracting network and the surrounding cytoplasm. During symmetry breaking, transient attachments to the cell boundary drive the contraction center to a polar location near the droplet boundary. Our findings demonstrate a robust, yet tunable, mechanism for subcellular localization.