computation
Latest
Towards model-based control of active matter: active nematics and oscillator networks
The richness of active matter's spatiotemporal patterns continues to capture our imagination. Shaping these emergent dynamics into pre-determined forms of our choosing is a grand challenge in the field. To complicate matters, multiple dynamical attractors can coexist in such systems, leading to initial condition-dependent dynamics. Consequently, non-trivial spatiotemporal inputs are generally needed to access these states. Optimal control theory provides a general framework for identifying such inputs and represents a promising computational tool for guiding experiments and interacting with various systems in soft active matter and biology. As an exemplar, I first consider an extensile active nematic fluid confined to a disk. In the absence of control, the system produces two topological defects that perpetually circulate. Optimal control identifies a time-varying active stress field that restructures the director field, flipping the system to its other attractor that rotates in the opposite direction. As a second, analogous case, I examine a small network of coupled Belousov-Zhabotinsky chemical oscillators that possesses two dominant attractors, two wave states of opposing chirality. Optimal control similarly achieves the task of attractor switching. I conclude with a few forward-looking remarks on how the same model-based control approach might come to bear on problems in biology.
Reconstruct cellular dynamics from single cell data
Recent advances of single cell techniques catalyzed quantitative studies on the dynamics of cell phenotypic transitions (CPT) emerging as a new field. However, fixed cell-based approaches have fundamental limits on revealing temporal information, and fluorescence-based live cell imaging approaches are technically challenging for multiplex long-term imaging. To tackle the challenges, we developed an integrated experimental/computational platform for reconstructing single cell phenotypic transition dynamics. Experimentally, we developed a live-cell imaging platform to record the phenotypic transition path of A549 VIM-RFP reporter cell line and unveil parallel paths of epithelial-to-mesenchymal transition (EMT). Computationally, we modified a finite temperature string method to reconstruct the reaction coordinate from the paths, and reconstruct a corresponding quasi-potential, which reveals that the EMT process resembles a barrier-less relaxation process. Our work demonstrates the necessity of extracting dynamical information of phenotypic transitions and the existence of a unified theoretical framework describing transition and relaxation dynamics in systems with and without detailed balance.
Coordinated motion of active filaments on spherical surfaces
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.
Microorganism locomotion in viscoelastic fluids
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.
Sperm Navigation: from hydrodynamic interactions to parameter estimation
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.
Simons-Emory Workshop on Neural Dynamics: What could neural dynamics have to say about neural computation, and do we know how to listen?
Speakers will deliver focused 10-minute talks, with periods reserved for broader discussion on topics at the intersection of neural dynamics and computation. Organizer & Moderator: Chethan Pandarinath - Emory University and Georgia Tech Speakers & Discussants: Adrienne Fairhall - U Washington Mehrdad Jazayeri - MIT John Krakauer - John Hopkins Francesca Mastrogiuseppe - Gatsby / UCL Abigail Person - U Colorado Abigail Russo - Princeton Krishna Shenoy - Stanford Saurabh Vyas - Columbia
Endless forms most beautiful: how to program materials using geometry, topology and singularities
The dream of programmable matter is to create materials whose physical properties (shape, moduli, response to perturbations, etc.) can be changed on the fly. For many years, my group has been thinking about how to program flat sheets that fold up into three dimensional shapes, most recently by exploiting the principles of origami design. Unfortunately, a combinatorial explosion of folding pathways makes robust folding particularly challenging. In this talk, I will discuss how this pluripotency arises from the topology of the configuration space. This suggests a broader understanding of a larger class of materials spanning from folding forms to spring networks to mechanical structures that perform computational logic.
Is there universality in biology?
It is sometimes said that there are two reasons why physics is so successful as a science. One is that it deals with very simple problems. The other is that it attempts to account only for universal aspects of systems at a desired level of description, with lower level phenomena subsumed into a small number of adjustable parameters. It is a widespread belief that this approach seems unlikely to be useful in biology, which is intimidatingly complex, where “everything has an exception”, and where there are a huge number of undetermined parameters. I will try to argue, nonetheless, that there are important, experimentally-testable aspects of biology that exhibit universality, and should be amenable to being tackled from a physics perspective. My suggestion is that this can lead to useful new insights into the existence and universal characteristics of living systems. I will try to justify this point of view by contrasting the goals and practices of the field of condensed matter physics with materials science, and then by extension, the goals and practices of the newly emerging field of “Physics of Living Systems” with biology. Specific biological examples that I will discuss include the following: Universal patterns of gene expression in cell biology Universal scaling laws in ecosystems, including the species-area law, Kleiber’s law, Paradox of the Plankton Universality of the genetic code Universality of thermodynamic utilization in microbial communities Universal scaling laws in the tree of life The question of what can be learned from studying universal phenomena in biology will also be discussed. Universal phenomena, by their very nature, shed little light on detailed microscopic levels of description. Yet there is no point in seeking idiosyncratic mechanistic explanations for phenomena whose explanation is found in rather general principles, such as the central limit theorem, that every microscopic mechanism is constrained to obey. Thus, physical perspectives may be better suited to answering certain questions such as universality than traditional biological perspectives. Concomitantly, it must be recognized that the identification and understanding of universal phenomena may not be a good answer to questions that have traditionally occupied biological scientists. Lastly, I plan to talk about what is perhaps the central question of universality in biology: why does the phenomenon of life occur at all? Is it an inevitable consequence of the laws of physics or some special geochemical accident? What methodology could even begin to answer this question? I will try to explain why traditional approaches to biology do not aim to answer this question, by comparing with our understanding of superconductivity as a physical phenomenon, and with the theory of universal computation. References Nigel Goldenfeld, Tommaso Biancalani, Farshid Jafarpour. Universal biology and the statistical mechanics of early life. Phil. Trans. R. Soc. A 375, 20160341 (14 pages) (2017). Nigel Goldenfeld and Carl R. Woese. Life is Physics: evolution as a collective phenomenon far from equilibrium. Ann. Rev. Cond. Matt. Phys. 2, 375-399 (2011).
“Biophysics of Structural Plasticity in Postsynaptic Spines”
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.
Building a synthetic cell: Understanding the clock design and function
Clock networks containing the same central architectures may vary drastically in their potential to oscillate, raising the question of what controls robustness, one of the essential functions of an oscillator. We computationally generate an atlas of oscillators and found that, while core topologies are critical for oscillations, local structures substantially modulate the degree of robustness. Strikingly, two local structures, incoherent and coherent inputs, can modify a core topology to promote and attenuate its robustness, additively. The findings underscore the importance of local modifications to the performance of the whole network. It may explain why auxiliary structures not required for oscillations are evolutionary conserved. We also extend this computational framework to search hidden network motifs for other clock functions, such as tunability that relates to the capabilities of a clock to adjust timing to external cues. Experimentally, we developed an artificial cell system in water-in-oil microemulsions, within which we reconstitute mitotic cell cycles that can perform self-sustained oscillations for 30 to 40 cycles over multiple days. The oscillation profiles, such as period, amplitude, and shape, can be quantitatively varied with the concentrations of clock regulators, energy levels, droplet sizes, and circuit design. Such innate flexibility makes it crucial to studying clock functions of tunability and stochasticity at the single-cell level. Combined with a pressure-driven multi-channel tuning setup and long-term time-lapse fluorescence microscopy, this system enables a high-throughput exploration in multi-dimension continuous parameter space and single-cell analysis of the clock dynamics and functions. We integrate this experimental platform with mathematical modeling to elucidate the topology-function relation of biological clocks. With FRET and optogenetics, we also investigate spatiotemporal cell-cycle dynamics in both homogeneous and heterogeneous microenvironments by reconstructing subcellular compartments.
“Models for Liquid-liquid Phase Separation of Intrinsically Disordered Proteins”
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.
Continuum modelling of active fluids beyond the generalised Taylor dispersion
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.
Flow, fluctuate and freeze: Epithelial cell sheets as soft active matter
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.
Dynamics of microbiota communities during physical perturbation
The consortium of microbes living in and on our bodies is intimately connected with human biology and deeply influenced by physical forces. Despite incredible gains in describing this community, and emerging knowledge of the mechanisms linking it to human health, understanding the basic physical properties and responses of this ecosystem has been comparatively neglected. Most diseases have significant physical effects on the gut; diarrhea alters osmolality, fever and cancer increase temperature, and bowel diseases affect pH. Furthermore, the gut itself is comprised of localized niches that differ significantly in their physical environment, and are inhabited by different commensal microbes. Understanding the impact of common physical factors is necessary for engineering robust microbiota members and communities; however, our knowledge of how they affect the gut ecosystem is poor. We are investigating how changes in osmolality affect the host and the microbial community and lead to mechanical shifts in the cellular environment. Osmotic perturbation is extremely prevalent in humans, caused by the use of laxatives, lactose intolerance, or celiac disease. In our studies we monitored osmotic shock to the microbiota using a comprehensive and novel approach, which combined in vivo experiments to imaging, physical measurements, computational analysis and highly controlled microfluidic experiments. By bridging several disciplines, we developed a mechanistic understanding of the processes involved in osmotic diarrhea, linking single-cell biophysical changes to large-scale community dynamics. Our results indicate that physical perturbations can profoundly and permanently change the competitive and ecological landscape of the gut, and affect the cell wall of bacteria differentially, depending on their mechanical characteristics.
Inaugural Simons-Emory Symposium On Motor Control: "What tools are we missing to understand motor control? What could we learn if we had them
This is the inaugural symposium of the Simons-Emory International Consortium on Motor Control, and speakers will deliver 10 minute talks (each followed by 10 minutes of discussion) addressing this question: "What tools are we missing to understand motor control, and what could we learn if we had them?”
computation coverage
15 items