Mechanical behavior in living cells consistent with the tensegrity model

Alternative models of cell mechanics depict the living cell as a simple mechanical continuum, porous filament gel, tensed cortical membrane, or tensegrity network that maintains a stabilizing prestress through incorporation of discrete structural elements that bear compression. Real-time microscopic analysis of cells containing GFP-labeled microtubules and associated mitochondria revealed that living cells behave like discrete structures composed of an interconnected network of actin microfilaments and microtubules when mechanical stresses are applied to cell surface integrin receptors. Quantitation of cell tractional forces and cellular prestress by using traction force microscopy confirmed that microtubules bear compression and are responsible for a significant portion of the cytoskeletal prestress that determines cell shape stability under conditions in which myosin light chain phosphorylation and intracellular calcium remained unchanged. Quantitative measurements of both static and dynamic mechanical behaviors in cells also were consistent with specific a priori predictions of the tensegrity model. These findings suggest that tensegrity represents a unified model of cell mechanics that may help to explain how mechanical behaviors emerge through collective interactions among different cytoskeletal filaments and extracellular adhesions in living cells.

A general two-process model describes the hydrogen exchange behavior of RNase A in unfolding conditions.

When NMR hydrogen exchange was used previously to monitor the kinetics of RNase A unfolding, some peptide NH protons were found to show EX2 exchange (detected by base catalysis) in addition to the expected EX1 exchange, whose rate is limited by the kinetic unfolding process. In earlier work, two groups showed independently that a restricted two-process model successfully fits published hydrogen exchange rates of native RNase A in the range 0-0.7 M guanidinium chloride. We find that this model predicts properties that are very different from the observed properties of the EX2 exchange reactions of RNase A in conditions where guanidine-induced unfolding takes place. The model predicts that EX2 exchange should be too fast to measure by the technique used, whereas it is readily measurable. Possible explanations for the contradiction are considered here, and we show that removing the restriction from the earlier two-process model is sufficient to resolve the contradiction; instead of specifying that exchange caused by global unfolding occurs by the EX2 mechanism, we allow it to occur by the general mechanism, which includes both the EX1 and EX2 cases. It is logical to remove this restriction because global unfolding of RNase A is known to give rise to EX1 exchange in these unfolding conditions. Resolving the contradiction makes it possible to determine whether populated unfolding intermediates contribute to the EX2 exchange, and this question is considered elsewhere. The results and simulations indicate that moderate or high denaturant concentrations readily give rise to EX1 exchange in native proteins. Earlier studies showed that hydrogen exchange in native proteins typically occurs by the EX2 mechanism but that high temperatures or pH values above 7 may give rise to EX1 exchange. High denaturant concentrations should be added to the list of variables likely to cause EX1 exchange.

Self-organized defensive behavior in honeybees

We investigated the defensive behavior of honeybees under controlled experimental conditions. During an attack on two identical targets, the spatial distribution of stings varied as a function of the total number of stings, evincing the classic “pitchfork bifurcation” phenomenon of nonlinear dynamics. The experimental results support a model of defensive behavior based on a self-organizing mechanism. The model helps to explain several of the characteristic features of the honeybee defensive response: (i) the ability of the colony to localize and focus its attack, (ii) the strong variability between different hives in the intensity of attack, as well as (iii) the variability observed within the same hive, and (iv) the ability of the colony to amplify small differences between the targets.

A model of excitation and adaptation in bacterial chemotaxis

Bacterial chemotaxis is widely studied because of its accessibility and because it incorporates processes that are important in a number of sensory systems: signal transduction, excitation, adaptation, and a change in behavior, all in response to stimuli. Quantitative data on the change in behavior are available for this system, and the major biochemical steps in the signal transduction/processing pathway have been identified. We have incorporated recent biochemical data into a mathematical model that can reproduce many of the major features of the intracellular response, including the change in the level of chemotactic proteins to step and ramp stimuli such as those used in experimental protocols. The interaction of the chemotactic proteins with the motor is not modeled, but we can estimate the degree of cooperativity needed to produce the observed gain under the assumption that the chemotactic proteins interact directly with the motor proteins.

A hierarchical neuronal network for planning behavior

Planning a goal-directed sequence of behavior is a higher function of the human brain that relies on the integrity of prefrontal cortical areas. In the Tower of London test, a puzzle in which beads sliding on pegs must be moved to match a designated goal configuration, patients with lesioned prefrontal cortex show deficits in planning a goal-directed sequence of moves. We propose a neuronal network model of sequence planning that passes this test and, when lesioned, fails in a way that mimics prefrontal patients’ behavior. Our model comprises a descending planning system with hierarchically organized plan, operation, and gesture levels, and an ascending evaluative system that analyzes the problem and computes internal reward signals that index the correct/erroneous status of the plan. Multiple parallel pathways connecting the evaluative and planning systems amend the plan and adapt it to the current problem. The model illustrates how specialized hierarchically organized neuronal assemblies may collectively emulate central executive or supervisory functions of the human brain.

Mouse model of Sanfilippo syndrome type B produced by targeted disruption of the gene encoding α-N-acetylglucosaminidase

The Sanfilippo syndrome type B is an autosomal recessive disorder caused by mutation in the gene (NAGLU) encoding α-N-acetylglucosaminidase, a lysosomal enzyme required for the stepwise degradation of heparan sulfate. The most serious manifestations are profound mental retardation, intractable behavior problems, and death in the second decade. To generate a model for studies of pathophysiology and of potential therapy, we disrupted exon 6 of Naglu, the homologous mouse gene. Naglu−/− mice were healthy and fertile while young and could survive for 8–12 mo. They were totally deficient in α-N-acetylglucosaminidase and had massive accumulation of heparan sulfate in liver and kidney as well as secondary changes in activity of several other lysosomal enzymes in liver and brain and elevation of gangliosides GM2 and GM3 in brain. Vacuolation was seen in many cells, including macrophages, epithelial cells, and neurons, and became more prominent with age. Although most vacuoles contained finely granular material characteristic of glycosaminoglycan accumulation, large pleiomorphic inclusions were seen in some neurons and pericytes in the brain. Abnormal hypoactive behavior was manifested by 4.5-mo-old Naglu−/− mice in an open field test; the hyperactivity that is characteristic of affected children was not observed even in younger mice. In a Pavlovian fear conditioning test, the 4.5-mo-old mutant mice showed normal response to context, indicating intact hippocampal-dependent learning, but reduced response to a conditioning tone, perhaps attributable to hearing impairment. The phenotype of the α-N-acetylglucosaminidase-deficient mice is sufficiently similar to that of patients with the Sanfilippo syndrome type B to make these mice a good model for study of pathophysiology and for development of therapy.

Social transmission of reproductive behavior increases frequency of inherited disorders in a young-expanding population

The observation of high frequencies of certain inherited disorders in the population of Saguenay–Lac Saint Jean can be explained in terms of the variance and the correlation of effective family size (EFS) from one generation to the next. We have shown this effect by using the branching process approach with real demographic data. When variance of EFS is included in the model, despite its profound effect on mutant allele frequency, any mutant introduced in the population never reaches the known carrier frequencies (between 0.035 and 0.05). It is only when the EFS correlation between generations is introduced into the model that we can explain the rise of the mutant alleles. This correlation is described by a c parameter that reflects the dependency of children’s EFS on their parents’ EFS. The c parameter can be considered to reflect social transmission of demographic behavior. We show that such social transmission dramatically reduces the effective population size. This could explain particular distributions in allele frequencies and unusually high frequency of certain inherited disorders in some human populations.

The Dynamic Behavior of Individual Microtubules Associated with Chromosomes In Vitro

Mitotic movements of chromosomes are usually coupled to the elongation and shortening of the microtubules to which they are bound. The lengths of kinetochore-associated microtubules change by incorporation or loss of tubulin subunits, principally at their chromosome-bound ends. We have reproduced aspects of this phenomenon in vitro, using a real-time assay that displays directly the movements of individual chromosome-associated microtubules as they elongate and shorten. Chromosomes isolated from cultured Chinese hamster ovary cells were adhered to coverslips and then allowed to bind labeled microtubules. In the presence of tubulin and GTP, these microtubules could grow at their chromosome-bound ends, causing the labeled segments to move away from the chromosomes, even in the absence of ATP. Sometimes a microtubule would switch to shortening, causing the direction of movement to change abruptly. The link between a microtubule and a chromosome was mechanically strong; 15 pN of tension was generally insufficient to detach a microtubule, even though it could add subunits at the kinetochore–microtubule junction. The behavior of the microtubules in vitro was regulated by the chromosomes to which they were bound; the frequency of transitions from polymerization to depolymerization was decreased, and the speed of depolymerization-coupled movement toward chromosomes was only one-fifth the rate of shortening for microtubules free in solution. Our results are consistent with a model in which each microtubule interacts with an increasing number of chromosome-associated binding sites as it approaches the kinetochore.

Fire–diffuse–fire model of dynamics of intracellular calcium waves

When Ca2+ is released from internal stores in living cells, the resulting wave of increased concentration can travel without deformation (continuous propagation) or with burst-like behavior (saltatory propagation). We analyze the “fire–diffuse–fire” model in order to illuminate the differences between these two modes of propagation. We show that the Ca2+ release wave in immature Xenopus oocytes and cardiac myocytes is saltatory, whereas the fertilization wave in the mature oocyte is continuous.

A model for individual and collective cell movement in Dictyostelium discoideum

The cellular slime mold Dictyostelium discoideum is a widely used model system for studying a variety of basic processes in development, including cell–cell signaling, signal transduction, pattern formation, cell motility, and the movement of tissue-like aggregates of cells. Many aspects of cell motion are poorly understood, including how individual cell behavior produces the collective motion of cells observed within the mound and slug. Herein, we describe a biologically realistic model for motile D. discoideum cells that can generate active forces, that interact via surface molecules, and that can detect and respond to chemotactic signals. We model the cells as deformable viscoelastic ellipsoids and incorporate signal transduction and cell–cell signaling by using a previously developed model. The shape constraint restricts the admissible deformations but makes the simulation of a large number of interacting cells feasible. Because the model is based on known processes, the parameters can be estimated or measured experimentally. We show that this model can reproduce the observations on the chemotactic behavior of single cells, streaming during aggregation, and the collective motion of an aggregate of cells driven by a small group of pacemakers. The model predicts that the motion of two-dimensional slugs [Bonner, J. T. (1998) Proc. Natl. Acad. Sci. USA 95, 9355–9359] results from the same behaviors that are exhibited by individual cells; it is not necessary to invoke different mechanisms or behaviors. Our computational experiments also suggest previously uncharacterized phenomena that may be experimentally observable.

A model for amplification of hair-bundle motion by cyclical binding of Ca2+ to mechanoelectrical-transduction channels

Amplification of auditory stimuli by hair cells augments the sensitivity of the vertebrate inner ear. Cell-body contractions of outer hair cells are thought to mediate amplification in the mammalian cochlea. In vertebrates that lack these cells, and perhaps in mammals as well, active movements of hair bundles may underlie amplification. We have evaluated a mathematical model in which amplification stems from the activity of mechanoelectrical-transduction channels. The intracellular binding of Ca2+ to channels is posited to promote their closure, which increases the tension in gating springs and exerts a negative force on the hair bundle. By enhancing bundle motion, this force partially compensates for viscous damping by cochlear fluids. Linear stability analysis of a six-state kinetic model reveals Hopf bifurcations for parameter values in the physiological range. These bifurcations signal conditions under which the system’s behavior changes from a damped oscillatory response to spontaneous limit-cycle oscillation. By varying the number of stereocilia in a bundle and the rate constant for Ca2+ binding, we calculate bifurcation frequencies spanning the observed range of auditory sensitivity for a representative receptor organ, the chicken’s cochlea. Simulations using prebifurcation parameter values demonstrate frequency-selective amplification with a striking compressive nonlinearity. Because transduction channels occur universally in hair cells, this active-channel model describes a mechanism of auditory amplification potentially applicable across species and hair-cell types.

Visual constraints in foraging bumblebees: Flower size and color affect search time and flight behavior

In optimal foraging theory, search time is a key variable defining the value of a prey type. But the sensory-perceptual processes that constrain the search for food have rarely been considered. Here we evaluate the flight behavior of bumblebees (Bombus terrestris) searching for artificial flowers of various sizes and colors. When flowers were large, search times correlated well with the color contrast of the targets with their green foliage-type background, as predicted by a model of color opponent coding using inputs from the bees' UV, blue, and green receptors. Targets that made poor color contrast with their backdrop, such as white, UV-reflecting ones, or red flowers, took longest to detect, even though brightness contrast with the background was pronounced. When searching for small targets, bees changed their strategy in several ways. They flew significantly slower and closer to the ground, so increasing the minimum detectable area subtended by an object on the ground. In addition, they used a different neuronal channel for flower detection. Instead of color contrast, they used only the green receptor signal for detection. We relate these findings to temporal and spatial limitations of different neuronal channels involved in stimulus detection and recognition. Thus, foraging speed may not be limited only by factors such as prey density, flight energetics, and scramble competition. Our results show that understanding the behavioral ecology of foraging can substantially gain from knowledge about mechanisms of visual information processing.

Forming classes by stimulus frequency: Behavior and theory

Visual classification is the way we relate to different images in our environment as if they were the same, while relating differently to other collections of stimuli (e.g., human vs. animal faces). It is still not clear, however, how the brain forms such classes, especially when introduced with new or changing environments. To isolate a perception-based mechanism underlying class representation, we studied unsupervised classification of an incoming stream of simple images. Classification patterns were clearly affected by stimulus frequency distribution, although subjects were unaware of this distribution. There was a common bias to locate class centers near the most frequent stimuli and their boundaries near the least frequent stimuli. Responses were also faster for more frequent stimuli. Using a minimal, biologically based neural-network model, we demonstrate that a simple, self-organizing representation mechanism based on overlapping tuning curves and slow Hebbian learning suffices to ensure classification. Combined behavioral and theoretical results predict large tuning overlap, implicating posterior infero-temporal cortex as a possible site of classification.

Neutral behavior of shared polymorphism

Several cases have been described in the literature where genetic polymorphism appears to be shared between a pair of species. Here we examine the distribution of times to random loss of shared polymorphism in the context of the neutral Wright–Fisher model. Order statistics are used to obtain the distribution of times to loss of a shared polymorphism based on Kimura’s solution to the diffusion approximation of the Wright–Fisher model. In a single species, the expected absorption time for a neutral allele having an initial allele frequency of ½ is 2.77 N generations. If two species initially share a polymorphism, that shared polymorphism is lost as soon as either of two species undergoes fixation. The loss of a shared polymorphism thus occurs sooner than loss of polymorphism in a single species and has an expected time of 1.7 N generations. Molecular sequences of genes with shared polymorphism may be characterized by the count of the number of sites that segregate in both species for the same nucleotides (or amino acids). The distribution of the expected numbers of these shared polymorphic sites also is obtained. Shared polymorphism appears to be more likely at genetic loci that have an unusually large number of segregating alleles, and the neutral coalescent proves to be very useful in determining the probability of shared allelic lineages expected by chance. These results are related to examples of shared polymorphism in the literature.