Aggregative cycles evolve as a solution to conflicts in social investment

Multicellular organization is particularly vulnerable to conflicts between different cell types when the body forms from initially isolated cells, as in aggregative multicellular microbes. Like other functions of the multicellular phase, coordinated collective movement can be undermined by conflicts between cells that spend energy in fuelling motion and ‘cheaters’ that get carried along. The evolutionary stability of collective behaviours against such conflicts is typically addressed in populations that undergo extrinsically imposed phases of aggregation and dispersal. Here, via a shift in perspective, we propose that aggregative multicellular cycles may have emerged as a way to temporally compartmentalize social conflicts. Through an eco-evolutionary mathematical model that accounts for individual and collective strategies of resource acquisition, we address regimes where different motility types coexist. Particularly interesting is the oscillatory regime that, similarly to life cycles of aggregative multicellular organisms, alternates on the timescale of several cell generations phases of prevalent solitary living and starvation-triggered aggregation. Crucially, such self-organized oscillations emerge as a result of evolution of cell traits associated to conflict escalation within multicellular aggregates.     

Fluctuations and correlations in crystalline protein dynamics: a simulation analysis of staphylococcal nuclease

Understanding collective motions in protein crystals is likely to furnish insight into functional protein dynamics and will improve models for refinement against diffraction data. Here, four 10 ns molecular dynamics simulations of crystalline Staphylococcal nuclease are reported and analyzed in terms of fluctuations and correlations in atomic motion. The simulation-derived fluctuations strongly correlate with, but are slightly higher than, the values derived from the experimental B-factors. Approximately 70% of the atomic fluctuations are due to internal protein motion. For 65% of the protein atoms the internal fluctuations converge on the nanosecond timescale. Convergence is much slower for the elements of the interatomic displacement correlation matrix--of these, >80% converge within 1 ns for interatomic distances less, approximately <6 A, but only 10% for separations approximately =12 A. Those collective motions that converged on the nanosecond timescale involve mostly correlations within the beta-barrel or between alpha-helices of the protein. The R-factor with the experimental x-ray diffuse scattering for the crystal, which is determined by the displacement variance-covariance matrix, decreases to 8% after 10 ns simulation. Both the number of converged correlation matrix elements and the R-factor depend logarithmically on time, consistent with a model in which the number of energy minima sampled depends exponentially on the maximum energy barrier crossed. The logarithmic dependence is also extrapolated to predict a convergence time for the whole variance-covariance matrix of approximately 1 micros.     

Maintenance of Motility Bias during Cyanobacterial Phototaxis

Signal transduction in bacteria is complex, ranging across scales from molecular signal detectors and effectors to cellular and community responses to stimuli. The unicellular, photosynthetic cyanobacterium Synechocystis sp. PCC6803 transduces a light stimulus into directional movement known as phototaxis. This response occurs via a biased random walk toward or away from a directional light source, which is sensed by intracellular photoreceptors and mediated by Type IV pili. It is unknown how quickly cells can respond to changes in the presence or directionality of light, or how photoreceptors affect single-cell motility behavior. In this study, we use time-lapse microscopy coupled with quantitative single-cell tracking to investigate the timescale of the cellular response to various light conditions and to characterize the contribution of the photoreceptor TaxD1 (PixJ1) to phototaxis. We first demonstrate that a community of cells exhibits both spatial and population heterogeneity in its phototactic response. We then show that individual cells respond within minutes to changes in light conditions, and that movement directionality is conferred only by the current light directionality, rather than by a long-term memory of previous conditions. Our measurements indicate that motility bias likely results from the polarization of pilus activity, yielding variable levels of movement in different directions. Experiments with a photoreceptor (taxD1) mutant suggest a supplementary role of TaxD1 in enhancing movement directionality, in addition to its previously identified role in promoting positive phototaxis. Motivated by the behavior of the taxD1 mutant, we demonstrate using a reaction-diffusion model that diffusion anisotropy is sufficient to produce the observed changes in the pattern of collective motility. Taken together, our results establish that single-cell tracking can be used to determine the factors that affect motility bias, which can then be coupled with biophysical simulations to connect changes in motility behaviors at the cellular scale with group dynamics.     

THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Predicting the conformational changes in proteins that are relevant for substrate binding is an ongoing challenge in the aim of elucidating the functional states of proteins. The motions that are induced by protein-ligand interactions are governed by the protein global modes. Our measurements indicate that the detected changes in the global backbone motion of the enzyme upon binding reflect a shift from the large-scale collective dominant mode in the unbound state towards a functional twisting deformation that assists in closing the binding cleft. Correlated motion in lysozyme has been implicated in enzyme function in previous studies, but detailed characterization of the internal fluctuations that enable the protein to explore the ensemble of conformations that ultimately foster large-scale conformational change is yet unknown. For this reason, we use THz spectroscopy to investigate the picosecond time scale binding modes and collective structural rearrangements that take place in hen egg white lysozyme (HEWL) when bound by the inhibitor (NAG) ₃. These protein thermal motions correspond to fluctuations that have a role in both selecting and sampling from the available protein intrinsic conformations that communicate function. Hence, investigation of these fast, collective modes may provide knowledge about the mechanism leading to the preferred binding process in HEWL-(NAG) ₃. Specifically, in this work we find that the picosecond time scale hydrogen-bonding rearrangements taking place in the protein hydration shell with binding modify the packing density within the hydrophobic core on a local level. These localized, intramolecular contact variations within the protein core appear to facilitate the large cooperative movements within the interfacial region separating the α- and β- domain that mediate binding. The THz time-scale fluctuations identified in the protein-ligand system may also reveal a molecular mechanism for substrate recognition.     

A study of low-frequency microturbulence by CO2-laser collective scattering in the FT-2 tokamak

Results are presented from studies of small-scale plasma density fluctuations in the FT-2 tokamak by the method of far-forward CO₂-laser collective scattering. The frequency and wavenumber spectra of fluctuations are measured using parallel k analysis at various positions of the scattering volume in the plane of the minor cross section of the torus. The data obtained are interpreted using numerical simulations. In phenomenological models, plasma fluctuations are substituted by a superposition of two-dimensional noninteracting cells with Gaussian profiles. A comparison of the calculated and experimental spectra shows that plasma fluctuations should be described based on the concept of strong microturbulence. The poloidal rotation velocity and the characteristic scale length of the scattering fluctuations, as well as the radial position of the region where they are located, are determined. The diffusion coefficient of the cells introduced in the model turns out to be close to the thermal diffusivity determined from the electron energy balance in the ohmic phase of the discharge.     

Collective oscillations in developing cells: insights from simple systems

From hormonal secretion to gene expression, multicellular dynamics are rich in oscillatory regulation. When organized in space and time, periodic cell-cell signaling can give rise to long-range coordination of gene expression and cell movement in tissues. Lack of synchrony of the oscillations on the other hand can serve as a source of initial divergence of cell fate in stem cells. How properties of individual cells can account for collective rhythmic behaviors at the tissue level remains elusive in most cases. Recently, studies in chemical reactions, synthetic gene circuits, yeast and social amoeba Dictyostelium have greatly enhanced our view of collective oscillations in cell populations. From these relatively simple systems, a unified view of how excitable and oscillatory regulations could be tuned and coupled to give rise to tissue-level oscillations is emerging. The review focuses on recent progress in cyclic adenosine monophosphate oscillations in Dictyostelium and highlights similarities and differences with other systems. We will see that the autonomy of single-cell level oscillations and different ways in which cells are coupled influence how group-level information can be encoded in collective oscillations.     

Minimal Network Topologies for Signal Processing during Collective Cell Chemotaxis

Cell-cell communication plays an important role in collective cell migration. However, it remains unclear how cells in a group cooperatively process external signals to determine the group’s direction of motion. Although the topology of signaling pathways is vitally important in single-cell chemotaxis, the signaling topology for collective chemotaxis has not been systematically studied. Here, we combine mathematical analysis and simulations to find minimal network topologies for multicellular signal processing in collective chemotaxis. We focus on border cell cluster chemotaxis in the Drosophila egg chamber, in which responses to several experimental perturbations of the signaling network are known. Our minimal signaling network includes only four elements: a chemoattractant, the protein Rac (indicating cell activation), cell protrusion, and a hypothesized global factor responsible for cell-cell interaction. Experimental data on cell protrusion statistics allows us to systematically narrow the number of possible topologies from more than 40,000,000 to only six minimal topologies with six interactions between the four elements. This analysis does not require a specific functional form of the interactions, and only qualitative features are needed; it is thus robust to many modeling choices. Simulations of a stochastic biochemical model of border cell chemotaxis show that the qualitative selection procedure accurately determines which topologies are consistent with the experiment. We fit our model for all six proposed topologies; each produces results that are consistent with all experimentally available data. Finally, we suggest experiments to further discriminate possible pathway topologies.     

Nanosecond Motions in Proteins Impose Bounds on the Timescale Distributions of Local Dynamics

We elucidate the physics of protein dynamical transition via 10-100-ns molecular dynamics simulations at temperatures spanning 160-300 K. By tracking the energy fluctuations, we show that the protein dynamical transition is marked by a crossover from nonstationary to stationary processes that underlie the dynamics of protein motions. A two-timescale function captures the nonexponential character of backbone structural relaxations. One timescale is attributed to the collective segmental motions and the other to local relaxations. The former is well defined by a single-exponential, nanosecond decay, operative at all temperatures. The latter is described by a set of processes that display a distribution of timescales. Although their average remains on the picosecond timescale, the distribution is markedly contracted at the onset of the transition. It is shown that the collective motions impose bounds on timescales spanned by local dynamical processes. The nonstationary character below the transition implicates the presence of a collection of substates whose interactions are restricted. At these temperatures, a wide distribution of local-motion timescales, extending beyond that of nanoseconds, is observed. At physiological temperatures, local motions are confined to timescales faster than nanoseconds. This relatively narrow window makes possible the appearance of multiple channels for the backbone dynamics to operate.     

Synchronizing a multicellular system by external input: an artificial control strategy

MOTIVATION: Although there are significant advances on elucidating the collective behaviors on biological organisms in recent years, the essential mechanisms by which the collective rhythms arise remain to be fully understood, and further how to synchronize multicellular networks by artificial control strategy has not yet been well explored. RESULTS: A control strategy is developed to synchronize gene regulatory networks in a multicellular system when spontaneous synchronization cannot be achieved. We first construct an impulsive control system to model the process of periodically injecting coupling substances with constant or random impulsive control amounts into the common extracellular medium, and further study its effects on the dynamics of individual cells. We derive the threshold of synchronization induced by the periodic substance input. Therefore, we can synchronize the multicellular network to a specific collective behavior by changing the frequency and amplitude of the periodic stimuli. Moreover, a two-stage scheme is proposed to facilitate the synchronization in this paper. We show that the presence of the external input may also initiate different dynamics. The multicellular network of coupled repressilators is used to show the effectiveness of the proposed method. The results not only provide a perspective to understand the interactions between external stimuli and intrinsic physiological rhythms, but also may lead to development of realistic artificial control strategy and medical therapy. AVAILABILITY: CONTACT: aihara@sat.t.u-tokyo.ac.jp.     

Estimates of the local parameters of plasma density fluctuations by reflectometry measurements

The question is considered of how to estimate the parameters of local plasma density fluctuations from reflectometry measurements made by probing the plasma with an extraordinary electromagnetic wave. In the geometrical-optics approximation, a formula is derived that relates the fluctuation amplitude of the phase of the reflected signal to the amplitude of local plasma density fluctuations and the range of its applicability is considered. The spectral sensitivity of reflectometry measurements in a reflection region of finite dimensions to poloidal perturbations with wavenumbers k _⊥ ? k ₀ is estimated by the phase-screen method, and the expressions obtained are compared with the results of numerical simulations. Based on the relationships derived, an algorithm is proposed for recovering the amplitude of the local plasma density fluctuations from the fluctuations in the reflected reflectometer signal. The results obtained are compared with the results of the full-wave simulations of the reflection of microwaves from a turbulent plasma. Finally, an example is given of how to recover the data on the amplitude of the local plasma density fluctuations in the T-10 tokamak plasma.     

Modeling mechanical inhomogeneities in small populations of proliferating monolayers and spheroids

Understanding the mechanical behavior of multicellular monolayers and spheroids is fundamental to tissue culture, organism development, and the early stages of tumor growth. Proliferating cells in monolayers and spheroids experience mechanical forces as they grow and divide and local inhomogeneities in the mechanical microenvironment can cause individual cells within the multicellular system to grow and divide at different rates. This differential growth, combined with cell division and reorganization, leads to residual stress. Multiple different modeling approaches have been taken to understand and predict the residual stresses that arise in growing multicellular systems, particularly tumor spheroids. Here, we show that by using a mechanically robust agent-based model constructed with the peridynamic framework, we gain a better understanding of residual stresses in multicellular systems as they grow from a single cell. In particular, we focus on small populations of cells (1–100 s) where population behavior is highly stochastic and prior investigation has been limited. We compare the average strain energy density of cells in monolayers and spheroids using different growth and division rules and find that, on average, cells in spheroids have a higher strain energy density than cells in monolayers. We also find that cells in the interior of a growing spheroid are, on average, in compression. Finally, we demonstrate the importance of accounting for stochastic fluctuations in the mechanical environment, particularly when the cellular response to mechanical cues is nonlinear. The results presented here serve as a starting point for both further investigation with agent-based models, and for the incorporation of major findings from agent-based models into continuum scale models when explicit representation of individual cells is not computationally feasible.     

Backbone dynamics of alamethicin bound to lipid membranes: spin-echo electron paramagnetic resonance of TOAC-spin labels

Alamethicin F50/5 is a hydrophobic peptide that is devoid of charged residues and that induces voltage-dependent ion channels in lipid membranes. The peptide backbone is likely to be involved in the ion conduction pathway. Electron spin-echo spectroscopy of alamethicin F50/5 analogs in which a selected Aib residue (at position n = 1, 8, or 16) is replaced by the TOAC amino-acid spin label was used to study torsional dynamics of the peptide backbone in association with phosphatidylcholine bilayer membranes. Rapid librational motions of limited angular amplitude were observed at each of the three TOAC sites by recording echo-detected spectra as a function of echo delay time, 2τ. Simulation of the time-resolved spectra, combined with conventional EPR measurements of the librational amplitude, shows that torsional fluctuations of the peptide backbone take place on the subnanosecond to nanosecond timescale, with little temperature dependence. Associated fluctuations in polar fields from the peptide could facilitate ion permeation.     

单细胞技术进展及其在植物研究中的应用

多细胞生物体的生存依赖于不同类型细胞特异性的功能分工,不同类型的细胞尽管基因组相同,但有其独特的发育过程和应对环境变化的能力。生物学的一大挑战就是揭示基因如何在正确的位置、正确的时间表达到正确的水平,最近出现了很多通过细胞类型特异性方法研究单细胞组学的工具,这些新技术使我们能通过空前分辨率,理解多细胞生物体内不同类型的单个细胞基因表达特点及其适应环境变化的机制。单细胞样品的获取一直是单细胞研究的一大技术瓶颈,因此本文将以如何获得起始材料为重点,探讨单细胞研究的样品标记、单细胞分离及获取、组学数据分析和结果验证等技术方法及其在植物研究中的应用。     

Enhanced persistence and collective migration in cooperatively aligning cell clusters

Most cells possess the capacity to locomote. Alone or collectively, this allows them to adapt, to rearrange, and to explore their surroundings. The biophysical characterization of such motile processes, in health and in disease, has so far focused mostly on two limiting cases: single-cell motility on the one hand and the dynamics of confluent tissues such as the epithelium on the other. The in-between regime of clusters, composed of relatively few cells moving as a coherent unit, has received less attention. Such small clusters are, however, deeply relevant in development but also in cancer metastasis. In this work, we use cellular Potts models and analytical active matter theory to understand how the motility of small cell clusters changes with N, the number of cells in the cluster. Modeling and theory reveal our two main findings: cluster persistence time increases with N, whereas the intrinsic diffusivity decreases with N. We discuss a number of settings in which the motile properties of more complex clusters can be analytically understood, revealing that the focusing effects of small-scale cooperation and cell-cell alignment can overcome the increased bulkiness and internal disorder of multicellular clusters to enhance overall migrational efficacy. We demonstrate this enhancement for small-cluster collective durotaxis, which is shown to proceed more effectively than for single cells. Our results may provide some novel, to our knowledge, insights into the connection between single-cell and large-scale collective motion and may point the way to the biophysical origins of the enhanced metastatic potential of small tumor cell clusters.     

Backbone dynamics of the natively unfolded pro-peptide of subtilisin by heteronuclear NMR relaxation studies

The dynamics of the natively unfolded form of the pro-peptide of subtilisin (PPS) have been characterized at two different pHs (6.0 and 3.0) by ¹⁵N relaxation experiments. ¹⁵N relaxation data is obtained at multiple field strengths and a detailed comparison of spectral density mapping, the model free approach and the recently proposed Cole–Cole model free (CC-MF) analysis is presented. The CC-MF analysis provides a better fit to the observed magnetic field dependence of ¹⁵N relaxation data of unfolded PPS than conventional model free approaches and shows that fluctuations in R₂ may be accounted for by a distribution of correlation times on the nanosecond timescale. A new parameter derives from the analysis and represents the width of the distribution function and the heterogeneity of the dynamics on the nanosecond timescale at a particular site. Particularly interesting is the observation that is sensitive to pH changes and that PPS samples a wider distribution of nanosecond time scale motions at less acidic pHs than at more acidic pHs. These results suggest that PPS experiences a higher degree of correlated motion at pH 6.0 and that electrostatic interactions may be important for inducing correlated motions on the nanosecond timescale in unfolded PPS.     

Localization Microscopy Reveals Expression-Dependent Parameters of Chromatin Nanostructure

A combined approach of 2D high-resolution localization light microscopy and statistical methods is presented to infer structural features and density fluctuations at the nuclear nanoscale. Hallmarks of nuclear nanostructure are found on the scale below 100 nm for both human fibroblast and HeLa cells. Mechanical measures were extracted as a quantitative tool from the histone density fluctuations inside the cell to obtain structural fluctuations on the scale of several micrometers. Results show that different mechanisms of expression of the same nuclear protein type lead to significantly different patterns on the nanoscale and to pronounced differences in the detected compressibility of chromatin. The observed fluctuations, including the experimental evidence for dynamic looping, are consistent with a recently proposed chromatin model.     

Cell jamming: Collective invasion of mesenchymal tumor cells imposed by tissue confinement

Background

Cancer invasion is a multi-step process which coordinates interactions between tumor cells with mechanotransduction towards the surrounding matrix, resulting in distinct cancer invasion strategies. Defined by context, mesenchymal tumors, including melanoma and fibrosarcoma, develop either single-cell or collective invasion modes, however, the mechanical and molecular programs underlying such plasticity of mesenchymal invasion programs remain unclear.

Methods

To test how tissue anatomy determines invasion mode, spheroids of MV3 melanoma and HT1080 fibrosarcoma cells were embedded into 3D collagen matrices of varying density and stiffness and analyzed for migration type and efficacy with matrix metalloproteinase (MMP)-dependent collagen degradation enabled or pharmacologically inhibited.

Results

With increasing collagen density and dependent on proteolytic collagen breakdown and track clearance, but independent of matrix stiffness, cells switched from single-cell to collective invasion modes. Conversion to collective invasion included gain of cell-to-cell junctions, supracellular polarization and joint guidance along migration tracks.

Conclusions

The density of the extracellulair matrix (ECM) determines the invasion mode of mesenchymal tumor cells. Whereas fibrillar, high porosity ECM enables single-cell dissemination, dense matrix induces cell–cell interaction, leader–follower cell behavior and collective migration as an obligate protease-dependent process.

General significance

These findings establish plasticity of cancer invasion programs in response to ECM porosity and confinement, thereby recapitulating invasion patterns of mesenchymal tumors in vivo. The conversion to collective invasion with increasing ECM confinement supports the concept of cell jamming as a guiding principle for melanoma and fibrosarcoma cells into dense tissue.This article is part of a Special Issue entitled Matrix-mediated cell behaviour and properties.     

Persistent cell-autonomous circadian oscillations in fibroblasts revealed by six-week single-cell imaging of PER2::LUC bioluminescence

Biological oscillators naturally exhibit stochastic fluctuations in period and amplitude due to the random nature of molecular reactions. Accurately measuring the precision of noisy oscillators and the heterogeneity in period and strength of rhythmicity across a population of cells requires single-cell recordings of sufficient length to fully represent the variability of oscillations. We found persistent, independent circadian oscillations of clock gene expression in 6-week-long bioluminescence recordings of 80 primary fibroblast cells dissociated from PER2::LUC mice and kept in vitro for 6 months. Due to the stochastic nature of rhythmicity, the proportion of cells appearing rhythmic increases with the length of interval examined, with 100% of cells found to be rhythmic when using 3-week windows. Mean period and amplitude are remarkably stable throughout the 6-week recordings, with precision improving over time. For individual cells, precision of period and amplitude are correlated with cell size and rhythm amplitude, but not with period, and period exhibits much less cycle-to-cycle variability (CV 7.3%) than does amplitude (CV 37%). The time series are long enough to distinguish stochastic fluctuations within each cell from differences among cells, and we conclude that the cells do exhibit significant heterogeneity in period and strength of rhythmicity, which we measure using a novel statistical metric. Furthermore, stochastic modeling suggests that these single-cell clocks operate near a Hopf bifurcation, such that intrinsic noise enhances the oscillations by minimizing period variability and sustaining amplitude.     

Nonequilibrium Brownian Motion beyond the Effective Temperature

The condition of thermal equilibrium simplifies the theoretical treatment of fluctuations as found in the celebrated Einstein’s relation between mobility and diffusivity for Brownian motion. Several recent theories relax the hypothesis of thermal equilibrium resulting in at least two main scenarios. With well separated timescales, as in aging glassy systems, equilibrium Fluctuation-Dissipation Theorem applies at each scale with its own “effective” temperature. With mixed timescales, as for example in active or granular fluids or in turbulence, temperature is no more well-defined, the dynamical nature of fluctuations fully emerges and a Generalized Fluctuation-Dissipation Theorem (GFDT) applies. Here, we study experimentally the mixed timescale regime by studying fluctuations and linear response in the Brownian motion of a rotating intruder immersed in a vibro-fluidized granular medium. Increasing the packing fraction, the system is moved from a dilute single-timescale regime toward a denser multiple-timescale stage. Einstein’s relation holds in the former and is violated in the latter. The violation cannot be explained in terms of effective temperatures, while the GFDT is able to impute it to the emergence of a strong coupling between the intruder and the surrounding fluid. Direct experimental measurements confirm the development of spatial correlations in the system when the density is increased.     

Coherent Neutron Scattering and Collective Dynamics in the Protein,GFP

Collective dynamics are considered to be one of the major properties of soft materials, including biological macromolecules. We present coherent neutron scattering studies of the low-frequency vibrations, the so-called boson peak, in fully deuterated green fluorescent protein (GFP). Our analysis revealed unexpectedly low coherence of the atomic motions in GFP. This result implies a low amount of in-phase collective motion of the secondary structural units contributing to the boson peak vibrations and fast conformational fluctuations on the picosecond timescale. These observations are in contrast to earlier studies of polymers and glass-forming systems, and suggest that random or out-of-phase motions of the β-strands contribute greater than two-thirds of the intensity to the low-frequency vibrational spectra of GFP.