Source author record

Herbert Levine

Herbert Levine appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

27 published item(s)

preprint2021arXiv

Collective Motility, Mechanical Waves, and Durotaxis in Cell Clusters

Epithelial cell clusters often move collectively on a substrate. Mechanical signals play a major role in organizing this behavior. There are a number of experimental observations in these systems which await a comprehensive explanation. These include: the internal strains are tensile even for clusters that expand by proliferation; the tractions on the substrate are often confined to the edges of the cluster; there can exist density waves within the cluster; clusters can exhibit collective durotaxis when individual cells show no effect; and for cells in an annulus there is a transition between expanding clusters with proliferation and the case where cells fill the annulus and rotate around it. We formulate a mechanical model to examine these effects. We use a molecular clutch picture which allows "stalling" -- inhibition of cell contraction by external forces. Stalled cells are passive from a physical point of view and the un-stalled cells are active. By attaching cells to the substrate and to each other, and taking into account contact inhibition of locomotion, we get a simple picture for many of these findings as well as predictions that could be tested. SI text/figures included, SI movies at https://rice.box.com/s/xiy3mwsfj3203lfu7pk0udfklexcgsew

preprint2021arXiv

Contact map dependence of a T cell receptor binding repertoire

The T cell arm of the adaptive immune system provides the host protection against unknown pathogens by discriminating between host and foreign material. This discriminatory capability is achieved by the creation of a repertoire of cells each carrying a T cell receptor (TCR) specific to non-self antigens displayed as peptides bound to the major histocompatibility complex (pMHC). The understanding of the dynamics of the adaptive immune system at a repertoire level is complex, due to both the nuanced interaction of a TCR-pMHC pair and to the number of different possible TCR-pMHC pairings, making computationally exact solutions currently unfeasible. To gain some insight into this problem, we study an affinity-based model for TCR-pMHC binding in which a crystal structure is used to generate a distance-based contact map that weights the pairwise amino acid interactions. We find that the TCR-pMHC binding energy distribution strongly depends both on the number of contacts and the repeat structure allowed by the topology of the contact map of choice; this in turn influences T cell recognition probability during negative selection, with higher variances leading to higher survival probabilities. In addition, we quantify the degree to which neoantigens with mutations in sites with higher contacts are recognized at a higher rate.

preprint2021arXiv

Understanding cytoskeletal avalanches using mechanical stability analysis

Eukaryotic cells are mechanically supported by a polymer network called the cytoskeleton, which consumes chemical energy to dynamically remodel its structure. Recent experiments in vivo have revealed that this remodeling occasionally happens through anomalously large displacements, reminiscent of earthquakes or avalanches. These cytoskeletal avalanches might indicate that the cytoskeleton's structural response to a changing cellular environment is highly sensitive, and they are therefore of significant biological interest. However, the physics underlying "cytoquakes" is poorly understood. Here, we use agent-based simulations of cytoskeletal self-organization to study fluctuations in the network's mechanical energy. We robustly observe non-Gaussian statistics and asymmetrically large rates of energy release compared to accumulation in a minimal cytoskeletal model. The large events of energy release are found to correlate with large, collective displacements of the cytoskeletal filaments. We also find that the changes in the localization of tension and the projections of the network motion onto the vibrational normal modes are asymmetrically distributed for energy release and accumulation. These results imply an avalanche-like process of slow energy storage punctuated by fast, large events of energy release involving a collective network rearrangement. We further show that mechanical instability precedes cytoquake occurrence through a machine learning model that dynamically forecasts cytoquakes using the vibrational spectrum as input. Our results provide the first connection between the cytoquake phenomenon and the network's mechanical energy and can help guide future investigations of the cytoskeleton's structural susceptibility.

preprint2020arXiv

Compression stiffening of fibrous networks with stiff inclusions

Tissues commonly consist of cells embedded within a fibrous biopolymer network. Whereas cell-free reconstituted biopolymer networks typically soften under applied uniaxial compression, various tissues, including liver, brain, and fat, have been observed to instead stiffen when compressed. The mechanism for this compression stiffening effect is not yet clear. Here, we demonstrate that when a material composed of stiff inclusions embedded in a fibrous network is compressed, heterogeneous rearrangement of the inclusions can induce tension within the interstitial network, leading to a macroscopic crossover from an initial bending-dominated softening regime to a stretching-dominated stiffening regime, which occurs before and independently of jamming of the inclusions. Using a coarse-grained particle-network model, we first establish a phase diagram for compression-driven, stretching-dominated stress propagation and jamming in uniaxially compressed 2- and 3-dimensional systems. Then, we demonstrate that a more detailed computational model of stiff inclusions in a subisostatic semiflexible fiber network exhibits quantitative agreement with the predictions of our coarse-grained model as well as qualitative agreement with experiments.

preprint2020arXiv

Decoding the mechanisms underlying cell-fate decision-making during stem cell differentiation by Random Circuit Perturbation

Stem cells can precisely and robustly undergo cellular differentiation and lineage commitment, referred to as stemness. However, how the gene network underlying stemness regulation reliably specifies cell fates is not well understood. To address this question, we applied a recently developed computational method, Random Circuit Perturbation (RACIPE), to a nine-component gene regulatory network (GRN) governing stemness, from which we identified fifteen robust gene states. Among them, four out of the five most probable gene states exhibit gene expression patterns observed in single mouse embryonic cells at 32-cell and 64-cell stages. These gene states can be robustly predicted by the stemness GRN but not by randomized versions of the stemness GRN. Strikingly, we found a hierarchical structure of the GRN with the Oct4/Cdx2 motif functioning as the first decision-making module followed by Gata6/Nanog. We propose that stem cell populations, instead of being viewed as all having a specific cellular state, can be regarded as a heterogeneous mixture including cells in various states. Upon perturbations by external signals, stem cells lose the capacity to access certain cellular states, thereby becoming differentiated. The findings demonstrate that the functions of the stemness GRN is mainly determined by its well-evolved network topology rather than by detailed kinetic parameters.

preprint2020arXiv

Insights from Graph Theory on the Morphologies of Actomyosin Networks with Multilinkers

Quantifying the influence of microscopic details on the dynamics of development of the overall structure of a filamentous network is important in a number of biologically relevant contexts, but it is not obvious what order parameters can be used to adequately describe this complex process. In this paper, we investigated the role of multivalent actin-binding proteins (ABPs) in reorganizing actin filaments into higher-order complex networks via a computer model of semiflexible filaments. We characterize the importance of local connectivity among actin filaments as well as the global features of actomyosin networks. We first map the networks into local graph representations and then, using principles from network-theory order parameters, combine properties from these representations to gain insight on the heterogeneous morphologies of actomyosin networks at a global level. We find that ABPs with a valency greater than two promote filament bundles and large filament clusters to a much greater extent than bivalent multilinkers. We also show that active myosin-like motor proteins promote the formation of dendritic branches from a stalk of actin bundles. Our work motivates future studies to embrace network theory as a tool to characterize complex morphologies of actomyosin detected by experiments, leading to a quantitative understanding of the role of ABPs in manipulating the self-assembly of actin filaments into unique architectures that underlie the structural scaffold of a cell relating to its mobility and shape.

preprint2020arXiv

Leader-cell-driven epithelial sheet fingering

Collective cell migration is crucial in many biological processes such as wound healing, tissue morphogenesis, and tumor progression. The leading front of a collective migrating epithelial cell layer often destabilizes into multicellular finger-like protrusions, each of which is guided by a leader cell at the fingertip. Here, we develop a subcellular-element-based model of this fingering instability, which incorporates leader cells and other related properties of a monolayer of epithelial cells. Our model recovers multiple aspects of the dynamics, especially the traction force patterns and velocity fields, observed in experiments on MDCK cells. Our model predicts the necessity of the leader cell and its minimal functions for the formation and maintenance of a stable finger pattern. Meanwhile, our model allows for an analysis of the role of supra-cellular actin cable on the leading front, predicting that while this observed structure helps maintain the shape of the finger, it is not required in order to form a finger. In addition, we also study the driving instability in the context of continuum active fluid model, which justifies some of our assumptions in the computational approach. In particular, we show that in our model no finger protrusions would emerge in a phenotypically homogenous active fluid and hence the role of the leader cell and its followers are often critical.

preprint2020arXiv

Ordered Hexagonal Patterns via Notch-Delta Signaling

Many developmental processes in biology utilize Notch-Delta signaling to construct an ordered pattern of cellular differentiation. This signaling modality is based on nearest-neighbor contact, as opposed to the more familiar mechanism driven by the release of diffusible ligands. Here, exploiting this "juxtacrine" property, we present an exact treatment of the pattern formation problem via a system of nine coupled ordinary differential equations. The possible patterns that are realized for realistic parameters can be analyzed by considering a co-dimension 2 pitchfork bifurcation of this system. This analysis explains the observed prevalence of hexagonal patterns with high Delta at their center, as opposed to those with central high Notch levels. We show that outside this range of parameters, in particular for low cis-coupling, a novel kind of pattern is produced, where high Delta cells have high Notch as well. It also suggests that the biological system is only weakly first order, so that an additional mechanism is required to generate the observed defect-free patterns. We construct a simple strategy for producing such defect-free patterns.

preprint2019arXiv

Biological Regulatory Networks are Minimally Frustrated

Characterization of the differences between biological and random networks can reveal the design principles that enable the robust realization of crucial biological functions including the establishment of different cell types. Previous studies, focusing on identifying topological features that are present in biological networks but not in random networks, have, however, provided few functional insights. We use a Boolean modeling framework and ideas from spin glass literature to identify functional differences between five real biological networks and random networks with similar topological features. We show that minimal frustration is a fundamental property that allows biological networks to robustly establish cell types and regulate cell fate choice, and this property can emerge in complex networks via Darwinian evolution. The study also provides clues regarding how the regulation of cell fate choice can go awry in a disease like cancer and lead to the emergence of aberrant cell types.

preprint2016arXiv

Modeling delayed processes in biological systems

Delayed processes are ubiquitous in biological systems and are often characterized by delay differential equations (DDEs) and their extension to include stochastic effects. DDEs do not explicitly incorporate intermediate states associated with a delayed process but instead use an estimated average delay time. In an effort to examine the validity of this approach, we study systems with significant delays by explicitly incorporating intermediate steps. We show by that such explicit models often yield significantly different equilibrium distributions and transition times as compared to DDEs with deterministic delay values. Additionally, different explicit models with qualitatively different dynamics can give rise to the same DDEs revealing important ambiguities. We also show that DDE-based predictions of oscillatory behavior may fail for the corresponding explicit model.

preprint2016arXiv

Nonlinear self-adapting wave patterns

We propose a new type of traveling wave pattern, one that can adapt to the size of physical system in which it is embedded. Such a system arises when the initial state has an instability that extends down to zero wavevector, connecting at that point to two symmetry modes of the underlying dynamical system. The Min system of proteins in E. coli is such as system with the symmetry emerging from the global conservation of two proteins, MinD and MinE. For this and related systems, traveling waves can adiabatically deform as the system is increased in size without the increase in node number that would be expected for an oscillatory version of a Turing instability containing an allowed wavenumber band with a finite minimum.

preprint2015arXiv

Collective signal processing in cluster chemotaxis: roles of adaptation, amplification, and co-attraction in collective guidance

Single eukaryotic cells commonly sense and follow chemical gradients, performing chemotaxis. Recent experiments and theories, however, show that even when single cells do not chemotax, clusters of cells may, if their interactions are regulated by the chemoattractant. We study this general mechanism of "collective guidance" computationally with models that integrate stochastic dynamics for individual cells with biochemical reactions within the cells, and diffusion of chemical signals between the cells. We show that if clusters of cells use the well-known local excitation, global inhibition (LEGI) mechanism to sense chemoattractant gradients, the speed of the cell cluster becomes non-monotonic in the cluster's size - clusters either larger or smaller than an optimal size will have lower speed. We argue that the cell cluster speed is a crucial readout of how the cluster processes chemotactic signal; both amplification and adaptation will alter the behavior of cluster speed as a function of size. We also show that, contrary to the assumptions of earlier theories, collective guidance does not require persistent cell-cell contacts and strong short range adhesion to function. If cell-cell adhesion is absent, and the cluster cohesion is instead provided by a co-attraction mechanism, e.g. chemotaxis toward a secreted molecule, collective guidance may still function. However, new behaviors, such as cluster rotation, may also appear in this case. Together, the combination of co-attraction and adaptation allows for collective guidance that is robust to varying chemoattractant concentrations while not requiring strong cell-cell adhesion.

preprint2015arXiv

Generic Criticality in Ecological and Neuronal Networks

We investigate the dynamics of two models of biological networks with purely suppressive interactions between the units; species interacting via niche competition and neurons via inhibitory synaptic coupling. In both of these cases, power-law scaling of the density of states with probability arises without any fine-tuning of the model parameters. These results argue against the increasingly popular notion that non-equilibrium living systems operate at special critical points, driven by there by evolution so as to enable adaptive processing of input data.

preprint2015arXiv

Implications of the hybrid epithelial/mesenchymal phenotype in metastasis

Understanding cell-fate decisions during tumorigenesis and metastasis is a major challenge in modern cancer biology. One canonical cell-fate decision that cancer cells undergo is Epithelial-to-Mesenchymal Transition (EMT) and its reverse Mesenchymal-to-Epithelial Transition (MET). While transitioning between these two phenotypes - epithelial and mesenchymal - cells can also attain a hybrid epithelial/mesenchymal (i.e. partial or intermediate EMT) phenotype. Cells in this phenotype have mixed epithelial (e.g. adhesion) and mesenchymal (e.g. migration) properties, thereby allowing them to move collectively as clusters of Circulating Tumor Cells (CTCs). If these clusters enter the circulation, they can be more apoptosis-resistant and more capable of initiating metastatic lesions than cancer cells moving individually with wholly mesenchymal phenotypes, having undergo a complete EMT. Here, we review the operating principles of the core regulatory network for EMT/MET that acts as a three-way switch giving rise to three distinct phenotypes - epithelial, mesenchymal and hybrid epithelial/mesenchymal. We further characterize this hybrid E/M phenotype in terms of its capabilities in terms of collective cell migration, tumor-initiation, cell-cell communication, and drug resistance. We elucidate how the highly interconnected coupling between these modules coordinates cell-fate decisions among a population of cancer cells in the dynamic tumor, hence facilitating tumor-stoma interactions, formation of CTC clusters, and consequently cancer metastasis. Finally, we discuss the multiple advantages that the hybrid epithelial/mesenchymal phenotype have as compared to a complete EMT phenotype and argue that these collectively migrating cells are the primary 'bad actors' of metastasis.

preprint2015arXiv

Nonlinear elasticity of disordered fiber networks

Disordered biopolymer gels have striking mechanical properties including strong nonlinearities. In the case of athermal gels (such as collagen-I) the nonlinearity has long been associated with a crossover from a bending dominated to a stretching dominated regime of elasticity. The physics of this crossover is related to the existence of a central-force isostatic point and to the fact that for most gels the bending modulus is small. This crossover induces scaling behavior for the elastic moduli. In particular, for linear elasticity such a scaling law has been demonstrated [Broedersz et al. Nature Physics, 2011 7, 983]. In this work we generalize the scaling to the nonlinear regime with a two-parameter scaling law involving three critical exponents. We test the scaling law numerically for two disordered lattice models, and find a good scaling collapse for the shear modulus in both the linear and nonlinear regimes. We compute all the critical exponents for the two lattice models and discuss the applicability of our results to real systems.

preprint2014arXiv

Alignment and Nonlinear Elasticity in Biopolymer Gels

We present a Landau type theory for the non-linear elasticity of biopolymer gels with a part of the order parameter describing induced nematic order of fibers in the gel. We attribute the non-linear elastic behavior of these materials to fiber alignment induced by strain. We suggest an application to contact guidance of cell motility in tissue. We compare our theory to simulation of a disordered lattice model for biopolymers. We treat homogeneous deformations such as simple shear, hydrostatic expansion, and simple extension, and obtain good agreement between theory and simulation. We also consider a localized perturbation which is a simple model for a contracting cell in a medium.

preprint2014arXiv

Properties of Cooperatively Induced Phases in Sensing Models

A large number of eukaryotic cells are able to directly detect external chemical gradients with great accuracy and the ultimate limit to their sensitivity has been a topic of debate for many years. Previous work has been done to understand many aspects of this process but little attention has been paid to the possibility of emergent sensing states. Here we examine how cooperation between sensors existing in a two dimensional network, as they do on the cell's surface, can both enhance and fundamentally alter the response of the cell to a spatially varying signal. We show that weakly interacting sensors linearly amplify the sensors response to an external gradient while a network of strongly interacting sensors form a collective non-linear response with two separate domains of active and inactive sensors forming what have called a "1/2-state" . In our analysis we examine the cell's ability to sense the direction of a signal and pay special attention to the substantially different behavior realized in the strongly interacting regime.

preprint2014arXiv

Scaling solution in the large population limit of the general asymmetric stochastic Luria-Delbrück evolution process

One of the most popular models for quantitatively understanding the emergence of drug resistance both in bacterial colonies and in malignant tumors was introduced long ago by Luria and Delbrück. Here, individual resistant mutants emerge randomly during the birth events of an exponentially growing sensitive population. A most interesting limit of this process occurs when the population size $N$ is large and mutation rates are low, but not necessarily small compared to $1/N$. Here we provide a scaling solution valid in this limit, making contact with the theory of Levy $α$-stable distributions, in particular one discussed long ago by Landau. One consequence of this association is that moments of the distribution are highly misleading as far as characterizing typical behavior. A key insight that enables our solution is that working in the fixed population size ensemble is not the same as working in a fixed time ensemble. Some of our results have been presented previously in shortened form.

preprint2013arXiv

Periodic migration in a physical model of cells on micropatterns

We extend a model for the morphology and dynamics of a crawling eukaryotic cell to describe cells on micropatterned substrates. This model couples cell morphology, adhesion, and cytoskeletal flow in response to active stresses induced by actin and myosin. We propose that protrusive stresses are only generated where the cell adheres, leading to the cell's effective confinement to the pattern. Consistent with experimental results, simulated cells exhibit a broad range of behaviors, including steady motion, turning, bipedal motion, and periodic migration, in which the cell crawls persistently in one direction before reversing periodically. We show that periodic motion emerges naturally from the coupling of cell polarization to cell shape by reducing the model to a simplified one-dimensional form that can be understood analytically.

preprint2012arXiv

How input fluctuations reshape the dynamics of a biological switching system

An important task in quantitative biology is to understand the role of stochasticity in biochemical regulation. Here, as an extension of our recent work [Phys. Rev. Lett. 107, 148101 (2011)], we study how input fluctuations affect the stochastic dynamics of a simple biological switch. In our model, the on transition rate of the switch is directly regulated by a noisy input signal, which is described as a nonnegative mean-reverting diffusion process. This continuous process can be a good approximation of the discrete birth-death process and is much more analytically tractable. Within this new setup, we apply the Feynman-Kac theorem to investigate the statistical features of the output switching dynamics. Consistent with our previous findings, the input noise is found to effectively suppress the input-dependent transitions. We show analytically that this effect becomes significant when the input signal fluctuates greatly in amplitude and reverts slowly to its mean.

preprint2011arXiv

Noise Effects in Nonlinear Biochemical Signaling

It has been generally recognized that stochasticity can play an important role in the information processing accomplished by reaction networks in biological cells. Most treatments of that stochasticity employ Gaussian noise even though it is a priori obvious that this approximation can violate physical constraints, such as the positivity of chemical concentrations. Here, we show that even when such nonphysical fluctuations are rare, an exact solution of the Gaussian model shows that the model can yield unphysical results. This is done in the context of a simple incoherent-feedforward model which exhibits perfect adaptation in the deterministic limit. We show how one can use the natural separation of time scales in this model to yield an approximate model, that is analytically solvable, including its dynamical response to an environmental change. Alternatively, one can employ a cutoff procedure to regularize the Gaussian result.

preprint2011arXiv

Propagating mode-I fracture in amorphous materials using the continuous random network (CRN) model

We study propagating mode-I fracture in two dimensional amorphous materials using atomistic simulations. We used the continuous random network (CRN) model of an amorphous material, creating samples using a two dimensional analogue of the WWW (Wooten, Winer & Weaire) Monte-Carlo algorithm. For modeling fracture, molecular-dynamics simulations were run on the resulting samples. The results of our simulations reproduce the main experimental features. In addition to achieving a steady-state crack under a constant driving displacement (which had not yet been achieved by other atomistic models for amorphous materials), the runs show micro-branching, which increases with driving, transitioning to macro-branching for the largest drivings. Beside the qualitative visual similarity of the simulated cracks to experiment, the simulation also succeeds in explaining the experimentally observed oscillations of the crack velocity.

preprint2011arXiv

Self-assisted Amoeboid Navigation in Complex Environments

Background: Living cells of many types need to move in response to external stimuli in order to accomplish their functional tasks; these tasks range from wound healing to immune response to fertilization. While the directional motion is typically dictated by an external signal, the actual motility is also restricted by physical constraints, such as the presence of other cells and the extracellular matrix. The ability to successfully navigate in the presence of obstacles is not only essential for organisms, but might prove relevant in the study of autonomous robotic motion. Methodology/principal findings: We study a computational model of amoeboid chemotactic navigation under differing conditions, from motion in an obstacle-free environment to navigation between obstacles and finally to moving in a maze. We use the maze as a simple stand-in for a motion task with severe constraints, as might be expected in dense extracellular matrix. Whereas agents using simple chemotaxis can successfully navigate around small obstacles, the presence of large barriers can often lead to agent trapping. We further show that employing a simple memory mechanism, namely secretion of a repulsive chemical by the agent, helps the agent escape from such trapping. Conclusions/significance: Our main conclusion is that cells employing simple chemotactic strategies will often be unable to navigate through maze-like geometries, but a simple chemical marker mechanism (which we refer to as "self-assistance") significantly improves success rates. This realization provides important insights into mechanisms that might be employed by real cells migrating in complex environments as well as clues for the design of robotic navigation strategies. The results can be extended to more complicated multi-cellular systems and can be used in the study of mammalian cell migration and cancer metastasis.

preprint2010arXiv

Spatial and Temporal Correlates of Vesicular Release at Hippocampal Synapses

We develop a spatially explicit biophysical model of the hippocampal CA3-CA1 presynaptic bouton to study local calcium dynamics leading to vesicle fusion. A kinetic model with two calcium sensors is formulated specifically for the CA3-CA1 synapse. The model includes a sensor for fast synchronous release that lasts a few tens of milliseconds and a sensor for slow asynchronous release that lasts a few hundred milliseconds. We show that a variety of extant data on CA3-CA1 synapse can be accounted for consistently only when a refractory period of the order of few milliseconds between releases is introduced. Including a second sensor for asynchronous release that has a slow unbinding site and therefore an embedded long memory, is shown to play a role in short-term plasticity by facilitating release. For synchronous release mediated by Synaptotagmin II a third time scale is revealed in addition to the fast and slow release. This third time scale corresponds to "stimulus-correlated super-fast" neurotransmitter release. Our detailed spatial simulation indicates that all three-time scales of neurotransmitter release are an emergent property of the calcium sensor and independent of synaptic ultrastructure. Furthermore, it allows us to identify features of synaptic transmission that are universal and those that are modulated by structure.

preprint2009arXiv

Signal processing in local neuronal circuits based on activity-dependent noise and competition

We study the characteristics of weak signal detection by a recurrent neuronal network with plastic synaptic coupling. It is shown that in the presence of an asynchronous component in synaptic transmission, the network acquires selectivity with respect to the frequency of weak periodic stimuli. For non-periodic frequency-modulated stimuli, the response is quantified by the mutual information between input (signal) and output (network's activity), and is optimized by synaptic depression. Introducing correlations in signal structure resulted in the decrease of input-output mutual information. Our results suggest that in neural systems with plastic connectivity, information is not merely carried passively by the signal; rather, the information content of the signal itself might determine the mode of its processing by a local neuronal circuit.

preprint2007arXiv

microRNAs may sharpen spatial expression patterns

The precise layout of gene expression patterns is a crucial step in development. Formation of a sharp boundary between high and low expression domains requires a genetic mechanism which is both sensitive and robust to fluctuations, a demand that may not be easily achieved by morphogens alone. Recently it has been demonstrated that small RNAs (and, in particular, microRNAs) play many roles in embryonic development. While some RNAs are essential for embryogenesis, others are limited to fine-tuning a predetermined gene expression pattern. Here we explore the possibility that small RNAs participate in sharpening a gene expression profile that was crudely established by a morphogen. To this end we study a model where small RNAs interact with a target gene and diffusively move from cell to cell. Though diffusion generally smears spatial expression patterns, we find that intercellular mobility of small RNAs is actually critical in sharpening the interface between target expression domains in a robust manner. We discuss the applicability of our results, as examples, to the case of leaf polarity establishment in maize and Hox patterning in the early {\it Drosophila} embryo. Our findings point out the functional significance of some mechanistic properties, such as mobility of small RNAs and the irreversibility of their interactions. These properties are yet to be established directly for most classes of small RNAs. An indirect yet simple experimental test of the proposed mechanism is suggested in some detail.

preprint1995arXiv

Aggregation Patterns in Stressed Bacteria

We study the formation of spot patterns seen in a variety of bacterial species when the bacteria are subjected to oxidative stress due to hazardous byproducts of respiration. Our approach consists of coupling the cell density field to a chemoattractant concentration as well as to nutrient and waste fields. The latter serves as a triggering field for emission of chemoattractant. Important elements in the proposed model include the propagation of a front of motile bacteria radially outward form an initial site, a Turing instability of the uniformly dense state and a reduction of motility for cells sufficiently far behind the front. The wide variety of patterns seen in the experiments is explained as being due the variation of the details of the initiation of the chemoattractant emission as well as the transition to a non-motile phase.