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 two- and three-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.

Normal stress anisotropy and marginal stability in athermal elastic networks.

Hydrogels of semiflexible biopolymers such as collagen have been shown to contract axially under shear strain, in contrast to the axial dilation observed for most elastic materials. Recent work has shown that this behavior can be understood in terms of the porous, two-component nature and consequent time-dependent compressibility of hydrogels. The apparent normal stress measured by a torsional rheometer reflects only the tensile contribution of the axial component σzz on long (compressible) timescales, crossing over to the first normal stress difference, N1 = σxx - σzz at short (incompressible) times. While the behavior of N1 is well understood for isotropic viscoelastic materials undergoing affine shear deformation, biopolymer networks are often anisotropic and deform nonaffinely. Here, we numerically study the normal stresses that arise under shear in subisostatic, athermal semiflexible polymer networks. We show that such systems exhibit strong deviations from affine behavior and that these anomalies are controlled by a rigidity transition as a function of strain.

Cell motility, contact guidance, and durotaxis.

Mechanical properties of the substrate play a vital role in cell motility. In particular, cells have been shown to migrate along aligned fibers in the substrate (contact guidance) and up stiffness gradients (durotaxis). Here we present a simple mechanical model for cell migration coupled to substrate properties, by placing a simulated cell on a lattice mimicking biopolymer gels or hydrogels. In our model cells attach to the substrate via focal adhesions (FAs). As the cells contract, forces are generated at the FAs, determining their maturation and detachment. At the same time, the cell was also allowed to move and rotate to maintain force and torque balance. Our model, in which the cells only have access to information regarding forces acting at the FAs, without a prior knowledge of the substrate stiffness or geometry, is able to reproduce both contact guidance and durotaxis.

Stress relaxation in F-actin solutions by severing.

Networks of filamentous actin (F-actin) are important for the mechanics of most animal cells. These cytoskeletal networks are highly dynamic, with a variety of actin-associated proteins that control cross-linking, polymerization and force generation in the cytoskeleton. Inspired by recent rheological experiments on reconstituted solutions of dynamic actin filaments, we report a theoretical model that describes stress relaxation behavior of these solutions in the presence of severing proteins. We show that depending on the kinetic rates of assembly, disassembly, and severing, one can observe both length-dependent and length-independent relaxation behavior.

Hindrances to precise recovery of cellular forces in fibrous biopolymer networks.

How cells move through the three-dimensional extracellular matrix (ECM) is of increasing interest in attempts to understand important biological processes such as cancer metastasis. Just as in motion on flat surfaces, it is expected that experimental measurements of cell-generated forces will provide valuable information for uncovering the mechanisms of cell migration. However, the recovery of forces in fibrous biopolymer networks may suffer from large errors. Here, within the framework of lattice-based models, we explore possible issues in force recovery by solving the inverse problem: how can one determine the forces cells exert to their surroundings from the deformation of the ECM? Our results indicate that irregular cell traction patterns, the uncertainty of local fiber stiffness, the non-affine nature of ECM deformations and inadequate knowledge of network topology will all prevent the precise force determination. At the end, we discuss possible ways of overcoming these difficulties.

Micromechanics of cellularized biopolymer networks.

Collagen gels are widely used in experiments on cell mechanics because they mimic the extracellular matrix in physiological conditions. Collagen gels are often characterized by their bulk rheology; however, variations in the collagen fiber microstructure and cell adhesion forces cause the mechanical properties to be inhomogeneous at the cellular scale. We study the mechanics of type I collagen on the scale of tens to hundreds of microns by using holographic optical tweezers to apply pN forces to microparticles embedded in the collagen fiber network. We find that in response to optical forces, particle displacements are inhomogeneous, anisotropic, and asymmetric. Gels prepared at 21 °C and 37 °C show qualitative difference in their micromechanical characteristics. We also demonstrate that contracting cells remodel the micromechanics of their surrounding extracellular matrix in a strain- and distance-dependent manner. To further understand the micromechanics of cellularized extracellular matrix, we have constructed a computational model which reproduces the main experiment findings.

Growth feedback as a basis for persister bistability.

A small fraction of cells in many bacterial populations, called persisters, are much less sensitive to antibiotic treatment than the majority. Persisters are in a dormant metabolic state, even while remaining genetically identical to the actively growing cells. Toxin and antitoxin modules in bacteria are believed to be one possible cause of persistence. A two-gene operon, HipBA, is one of many chromosomally encoded toxin and antitoxin modules in Escherichia coli and the HipA7 allelic variant was the first validated high-persistence mutant. Here, we present a stochastic model that can generate bistability of the HipBA system, via the reciprocal coupling of free HipA to the cellular growth rate. The actively growing state and the dormant state each correspond to a stable state of this model. Fluctuations enable transitions from one to the other. This model is fully in agreement with experimental data obtained with synthetic promoter constructs.

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. Nat. Phys., 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.

RACIPE: a computational tool for modeling gene regulatory circuits using randomization.

BACKGROUND: One of the major challenges in traditional mathematical modeling of gene regulatory circuits is the insufficient knowledge of kinetic parameters. These parameters are often inferred from existing experimental data and/or educated guesses, which can be time-consuming and error-prone, especially for large networks. RESULTS: We present a user-friendly computational tool for the community to use our newly developed method named random circuit perturbation (RACIPE), to explore the robust dynamical features of gene regulatory circuits without the requirement of detailed kinetic parameters. Taking the network topology as the only input, RACIPE generates an ensemble of circuit models with distinct randomized parameters and uniquely identifies robust dynamical properties by statistical analysis. Here, we discuss the implementation of the software and the statistical analysis methods of RACIPE-generated data to identify robust gene expression patterns and the functions of genes and regulatory links. Finally, we apply the tool on coupled toggle-switch circuits and a published circuit of B-lymphopoiesis. CONCLUSIONS: We expect our new computational tool to contribute to a more comprehensive and unbiased understanding of mechanisms underlying gene regulatory networks. RACIPE is a free open source software distributed under (Apache 2.0) license and can be downloaded from GitHub ( https://github.com/simonhb1990/RACIPE-1.0 ).

Phenomenological modeling of durotaxis.

Cells exhibit qualitatively different behaviors on substrates with different rigidities. The fact that cells are more polarized on the stiffer substrate motivates us to construct a two-dimensional cell with the distribution of focal adhesions dependent on substrate rigidities. This distribution affects the forces exerted by the cell and thereby determines its motion. Our model reproduces the experimental observation that the persistence time is higher on the stiffer substrate. This stiffness-dependent persistence will lead to durotaxis, the preference in moving towards stiffer substrates. This propensity is characterized by the durotaxis index first defined in experiments. We derive and validate a two-dimensional corresponding Fokker-Planck equation associated with our model. Our approach highlights the possible role of the focal adhesion arrangement in durotaxis.

Stress-induced plasticity of dynamic collagen networks.

The structure and mechanics of tissues is constantly perturbed by endogenous forces originated from cells, and at the same time regulate many important cellular functions such as migration, differentiation, and growth. Here we show that 3D collagen gels, major components of connective tissues and extracellular matrix (ECM), are significantly and irreversibly remodeled by cellular traction forces, as well as by macroscopic strains. To understand this ECM plasticity, we develop a computational model that takes into account the sliding and merging of ECM fibers. We have confirmed the model predictions with experiment. Our results suggest the profound impacts of cellular traction forces on their host ECM during development and cancer progression, and suggest indirect mechanical channels of cell-cell communications in 3D fibrous matrices.The structure and mechanics of tissues is constantly perturbed by endogenous forces originated from cells. Here the authors show that 3D collagen gels, major components of connective tissues and extracellular matrix, are significantly and irreversibly remodelled by cellular traction forces and by macroscopic strains.

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

Alignment and nonlinear elasticity in biopolymer gels.

We present a Landau-type theory for the nonlinear elasticity of biopolymer gels with a part of the order parameter describing induced nematic order of fibers in the gel. We attribute the nonlinear 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.

A formalized design process for bacterial consortia that perform logic computing.

The concept of microbial consortia is of great attractiveness in synthetic biology. Despite of all its benefits, however, there are still problems remaining for large-scaled multicellular gene circuits, for example, how to reliably design and distribute the circuits in microbial consortia with limited number of well-behaved genetic modules and wiring quorum-sensing molecules. To manage such problem, here we propose a formalized design process: (i) determine the basic logic units (AND, OR and NOT gates) based on mathematical and biological considerations; (ii) establish rules to search and distribute simplest logic design; (iii) assemble assigned basic logic units in each logic operating cell; and (iv) fine-tune the circuiting interface between logic operators. We in silico analyzed gene circuits with inputs ranging from two to four, comparing our method with the pre-existing ones. Results showed that this formalized design process is more feasible concerning numbers of cells required. Furthermore, as a proof of principle, an Escherichia coli consortium that performs XOR function, a typical complex computing operation, was designed. The construction and characterization of logic operators is independent of "wiring" and provides predictive information for fine-tuning. This formalized design process provides guidance for the design of microbial consortia that perform distributed biological computation.

Automated design of genetic toggle switches with predetermined bistability.

Synthetic biology aims to rationally construct biological devices with required functionalities. Methods that automate the design of genetic devices without post-hoc adjustment are therefore highly desired. Here we provide a method to predictably design genetic toggle switches with predetermined bistability. To accomplish this task, a biophysical model that links ribosome binding site (RBS) DNA sequence to toggle switch bistability was first developed by integrating a stochastic model with RBS design method. Then, to parametrize the model, a library of genetic toggle switch mutants was experimentally built, followed by establishing the equivalence between RBS DNA sequences and switch bistability. To test this equivalence, RBS nucleotide sequences for different specified bistabilities were in silico designed and experimentally verified. Results show that the deciphered equivalence is highly predictive for the toggle switch design with predetermined bistability. This method can be generalized to quantitative design of other probabilistic genetic devices in synthetic biology.