Nonlinear control of heart rate variability in human infants.

Nonlinear analyses of infant heart rhythms reveal a marked rise in the complexity of the electrocardiogram with maturation. We find that normal mature infants (gestation greater than or equal to 35 weeks) have complex and distinctly nonlinear heart rhythms (consistent with recent reports for healthy adults) but that such nonlinearity is lacking in preterm infants (gestation > or = to 27 weeks) where parasympathetic-sympathetic interaction and function are presumed to be less well developed. Our study further shows that infants with clinical brain death and those treated with atropine exhibit a similar lack of nonlinear feedback control. These three lines of evidence support the hypothesis championed by Goldberger et al. [Goldberger, A.L., Rigney, D.R. & West, B.J. (1990) Sci. Am. 262, 43-49] that autonomic nervous system control underlies the nonlinearity and possible chaos of normal heart rhythms. This report demonstrates the acquisition of nonlinear heart rate dynamics and possible chaos in developing human infants and its loss in brain death and with the administration of atropine. It parallels earlier work documenting changes in the variability of heart rhythms in each of these cases and suggests that nonlinearity may provide additional power in characterizing physiological states.

Variation in per capita interaction strength: Thresholds due to nonlinear dynamics and nonequilibrium conditions

I measured the strength of interaction between a marine herbivore and its growing resource over a realistic range of absolute and relative abundances. The herbivores (hermit crabs: Pagurus spp.) have slow and/or weak functional and numerical responses to epiphytic diatoms (Isthmia nervosa), which show logistic growth in the absence of consumers. By isolating this interaction in containers in the field, I mimicked many of the physical and biological variables characteristic of the intertidal while controlling the densities of focal species. The per capita effects of consumers on the population dynamics of their resource (i.e., interaction strength) were defined by using the relationship between hermit crab density and proportional change in the resource. When this relationship is fit by a Weibull function, a single parameter distinguishes constant interaction strength from one that varies as a function of density. Constant interaction strength causes the proportion of diatoms to fall linearly or proportionally as hermit crab density increases whereas per capita effects that increase with density cause an accelerating decline. Although many mathematical models of species interactions assume linear dynamics and invariant parameters, at least near equilibrium, the per capita effects of hermit crabs on diatoms varied substantially, apparently crossing a threshold from weak to strong when consumption exceeded resource production. This threshold separates a domain of coexistence from one of local extinction of the resource. Such thresholds may help explain trophic cascades, resource compensation, and context-dependent interaction strengths, while indicating a way to predict trophic effects, despite nonlinearities, as a function of vital rates.

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.

Nonlinear-dynamical arrhythmia control in humans

Nonlinear-dynamical control techniques, also known as chaos control, have been used with great success to control a wide range of physical systems. Such techniques have been used to control the behavior of in vitro excitable biological tissue, suggesting their potential for clinical utility. However, the feasibility of using such techniques to control physiological processes has not been demonstrated in humans. Here we show that nonlinear-dynamical control can modulate human cardiac electrophysiological dynamics by rapidly stabilizing an unstable target rhythm. Specifically, in 52/54 control attempts in five patients, we successfully terminated pacing-induced period-2 atrioventricular-nodal conduction alternans by stabilizing the underlying unstable steady-state conduction. This proof-of-concept demonstration shows that nonlinear-dynamical control techniques are clinically feasible and provides a foundation for developing such techniques for more complex forms of clinical arrhythmia.

Nonlinear dynamics in ventricular fibrillation.

Electrogram recordings of ventricular fibrillation appear complex and possibly chaotic. However, sequences of beat-to-beat intervals obtained from these recordings are generally short, making it difficult to explicitly demonstrate nonlinear dynamics. Motivated by the work of Sugihara on atmospheric dynamics and the Durbin-Watson test for nonlinearity, we introduce a new statistical test that recovers significant dynamical patterns from smoothed lag plots. This test is used to show highly significant nonlinear dynamics in a stable canine model of ventricular fibrillation.

Fluorescence and photobleaching dynamics of single light-harvesting complexes

Single light-harvesting complexes LH-2 from Rhodopseudomonas acidophila were immobilized on various charged surfaces under physiological conditions. Polarized light experiments showed that the complexes were situated on the surface as nearly upright cylinders. Their fluorescence lifetimes and photobleaching properties were obtained by using a confocal fluorescence microscope with picosecond time resolution. Initially all molecules fluoresced with a lifetime of 1 ± 0.2 ns, similar to the bulk value. The photobleaching of one bacteriochlorophyll molecule from the 18-member assembly caused the fluorescence to switch off completely, because of trapping of the mobile excitations by energy transfer. This process was linear in light intensity. On continued irradiation the fluorescence often reappeared, but all molecules did not show the same behavior. Some LH-2 complexes displayed a variation of their quantum yields that was attributed to photoinduced confinement of the excited states and thereby a diminution of the superradiance. Others showed much shorter lifetimes caused by excitation energy traps that are only ≈3% efficient. On repeated excitation some molecules entered a noisy state where the fluorescence switched on and off with a correlation time of ≈0.1 s. About 490 molecules were examined.

Experimental evidence for a behavior-mediated trophic cascade in a terrestrial food chain

Predators of herbivorous animals can affect plant populations by altering herbivore density, behavior, or both. To test whether the indirect effect of predators on plants arises from density or behavioral responses in a herbivore population, we experimentally examined the dynamics of terrestrial food chains comprised of old field plants, leaf-chewing grasshoppers, and spider predators in Northeast Connecticut. To separate the effects of predators on herbivore density from the effects on herbivore behavior, we created two classes of spiders: (i) risk spiders that had their feeding mouth parts glued to render them incapable of killing prey and (ii) predator spiders that remained unmanipulated. We found that the effect of predators on plants resulted from predator-induced changes in herbivore behavior (shifts in activity time and diet selection) rather than from predator-induced changes in grasshopper density. Neither predator nor risk spiders had a significant effect on grasshopper density relative to a control. This demonstrates that the behavioral response of prey to predators can have a strong impact on the dynamics of terrestrial food chains. The results make a compelling case to examine behavioral as well as density effects in theoretical and empirical research on food chain dynamics.

Birefringence Imaging Directly Reveals Architectural Dynamics of Filamentous Actin in Living Growth ConesV⃞

We have investigated the dynamic behavior of cytoskeletal fine structure in the lamellipodium of nerve growth cones using a new type of polarized light microscope (the Pol-Scope). Pol-Scope images display with exquisite resolution and definition birefringent fine structures, such as filaments and membranes, without having to treat the cell with exogenous dyes or fluorescent labels. Furthermore, the measured birefringence of protein fibers in the thin lamellipodial region can be interpreted in terms of the number of filaments in the bundles. We confirmed that birefringent fibers are actin-based using conventional fluorescence-labeling methods. By recording movies of time-lapsed Pol-Scope images, we analyzed the creation and dynamic composition of radial fibers, filopodia, and intrapodia in advancing growth cones. The strictly quantitative information available in time-lapsed Pol-Scope images confirms previously deduced behavior and provides new insight into the architectural dynamics of filamentous actin.

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.

Modeling the dynamics of human hair cycles by a follicular automaton

The hair follicle cycle successively goes through the anagen, catagen, telogen, and latency phases, which correspond, respectively, to hair growth, arrest, shedding, and absence before a new anagen phase is initiated. Experimental observations collected over a period of 14 years in a group of 10 male volunteers, alopecic and nonalopecic, allowed us to determine the characteristics of scalp hair follicle cycles. On the basis of these observations, we propose a follicular automaton model to simulate the dynamics of human hair cycles. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states anagen, telogen, and latency, and the subsequent return to anagen. The transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a variance. The follicular automaton model accounts both for the dynamical transitions observed in a single follicle and for the behavior of an ensemble of independently cycling follicles. Thus, the model successfully reproduces the evolution of the fractions of follicle populations in each of the three phases, which fluctuate around steady-state or slowly drifting values. We apply the follicular automaton model to the study of spatial patterns of follicular growth that result from a spatially heterogeneous distribution of parameters such as the mean duration of anagen phase. When considering that follicles die or miniaturize after going through a critical number of successive cycles, the model can reproduce the evolution to hair patterns similar to well known types of diffuse or androgenetic alopecia.

Cell Cycle-Dependent Changes in Microtubule Dynamics in Living Cells Expressing Green Fluorescent Protein-α TubulinV⃞

LLCPK-1 cells were transfected with a green fluorescent protein (GFP)-α tubulin construct and a cell line permanently expressing GFP-α tubulin was established (LLCPK-1α). The mitotic index and doubling time for LLCPK-1α were not significantly different from parental cells. Quantitative immunoblotting showed that 17% of the tubulin in LLCPK-1α cells was GFP-tubulin; the level of unlabeled tubulin was reduced to 82% of that in parental cells. The parameters of microtubule dynamic instability were compared for interphase LLCPK-1α and parental cells injected with rhodamine-labeled tubulin. Dynamic instability was very similar in the two cases, demonstrating that LLCPK-1α cells are a useful tool for analysis of microtubule dynamics throughout the cell cycle. Comparison of astral microtubule behavior in mitosis with microtubule behavior in interphase demonstrated that the frequency of catastrophe increased twofold and that the frequency of rescue decreased nearly fourfold in mitotic compared with interphase cells. The percentage of time that microtubules spent in an attenuated state, or pause, was also dramatically reduced, from 73.5% in interphase to 11.4% in mitosis. The rates of microtubule elongation and rapid shortening were not changed; overall dynamicity increased 3.6-fold in mitosis. Microtubule release from the centrosome and a subset of differentially stable astral microtubules were also observed. The results provide the first quantitative measurements of mitotic microtubule dynamics in mammalian cells.

Molecular dynamics of a grid-mounted molecular dipolar rotor in a rotating electric field

Classical molecular dynamics is applied to the rotation of a dipolar molecular rotor mounted on a square grid and driven by rotating electric field E(ν) at T ≃ 150 K. The rotor is a complex of Re with two substituted o-phenanthrolines, one positively and one negatively charged, attached to an axial position of Rh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{2}^{4+}}}\end{equation*}\end{document} in a [2]staffanedicarboxylate grid through 2-(3-cyanobicyclo[1.1.1]pent-1-yl)malonic dialdehyde. Four regimes are characterized by a, the average lag per turn: (i) synchronous (a < 1/e) at E(ν) = |E(ν)| > Ec(ν) [Ec(ν) is the critical field strength], (ii) asynchronous (1/e < a < 1) at Ec(ν) > E(ν) > Ebo(ν) > kT/μ, [Ebo(ν) is the break-off field strength], (iii) random driven (a ≃ 1) at Ebo(ν) > E(ν) > kT/μ, and (iv) random thermal (a ≃ 1) at kT/μ > E(ν). A fifth regime, (v) strongly hindered, W > kT, Eμ, (W is the rotational barrier), has not been examined. We find Ebo(ν)/kVcm−1 ≃ (kT/μ)/kVcm−1 + 0.13(ν/GHz)1.9 and Ec(ν)/kVcm−1 ≃ (2.3kT/μ)/kVcm−1 + 0.87(ν/GHz)1.6. For ν > 40 GHz, the rotor behaves as a macroscopic body with a friction constant proportional to frequency, η/eVps ≃ 1.14 ν/THz, and for ν < 20 GHz, it exhibits a uniquely molecular behavior.

Linear and nonlinear pathways of spectral information transmission in the cochlear nucleus

At the level of the cochlear nucleus (CN), the auditory pathway divides into several parallel circuits, each of which provides a different representation of the acoustic signal. Here, the representation of the power spectrum of an acoustic signal is analyzed for two CN principal cells—chopper neurons of the ventral CN and type IV neurons of the dorsal CN. The analysis is based on a weighting function model that relates the discharge rate of a neuron to first- and second-order transformations of the power spectrum. In chopper neurons, the transformation of spectral level into rate is a linear (i.e., first-order) or nearly linear function. This transformation is a predominantly excitatory process involving multiple frequency components, centered in a narrow frequency range about best frequency, that usually are processed independently of each other. In contrast, type IV neurons encode spectral information linearly only near threshold. At higher stimulus levels, these neurons are strongly inhibited by spectral notches, a behavior that cannot be explained by level transformations of first- or second-order. Type IV weighting functions reveal complex excitatory and inhibitory interactions that involve frequency components spanning a wider range than that seen in choppers. These findings suggest that chopper and type IV neurons form parallel pathways of spectral information transmission that are governed by two different mechanisms. Although choppers use a predominantly linear mechanism to transmit tonotopic representations of spectra, type IV neurons use highly nonlinear processes to signal the presence of wide-band spectral features.

Ultrafast diffraction and structural dynamics: The nature of complex molecules far from equilibrium

Studies of molecular structures at or near their equilibrium configurations have long provided information on their geometry in terms of bond distances and angles. Far-from-equilibrium structures are relatively unknown—especially for complex systems—and generally, neither their dynamics nor their average geometries can be extrapolated from equilibrium values. For such nonequilibrium structures, vibrational amplitudes and bond distances play a central role in phenomena such as energy redistribution and chemical reactivity. Ultrafast electron diffraction, which was developed to study transient molecular structures, provides a direct method for probing the nature of complex molecules far from equilibrium. Here we present our ultrafast electron diffraction observations of transient structures for two cyclic hydrocarbons. At high internal energies of ≈4 eV, these molecules display markedly different behavior. For 1,3,5-cycloheptatriene, excitation results in the formation of hot ground-state structures with bond distances similar to those of the initial structure, but with nearly three times the average vibrational amplitude. Energy is redistributed within 5 ps, but with a negative temperature characterizing the nonequilibrium population. In contrast, the ring-opening reaction of 1,3-cyclohexadiene is shown to result in hot structures with a C—C bond distance of over 1.7 Å, which is 0.2 Å away from any expected equilibrium value. Even up to 400 ps, energy remains trapped in large-amplitude motions comprised of torsion and asymmetric stretching. These studies promise a new direction for studying structural dynamics in nonequilibrium complex systems.

Spindle Dynamics and the Role of γ-Tubulin in Early Caenorhabditis elegans EmbryosV⃞

γ-Tubulin is a ubiquitous and highly conserved component of centrosomes in eukaryotic cells. Genetic and biochemical studies have demonstrated that γ-tubulin functions as part of a complex to nucleate microtubule polymerization from centrosomes. We show that, as in other organisms, Caenorhabditis elegans γ-tubulin is concentrated in centrosomes. To study centrosome dynamics in embryos, we generated transgenic worms that express GFP::γ-tubulin or GFP::β-tubulin in the maternal germ line and early embryos. Multiphoton microscopy of embryos produced by these worms revealed the time course of daughter centrosome appearance and growth and the differential behavior of centrosomes destined for germ line and somatic blastomeres. To study the role of γ-tubulin in nucleation and organization of spindle microtubules, we used RNA interference (RNAi) to deplete C. elegans embryos of γ-tubulin. γ-Tubulin (RNAi) embryos failed in chromosome segregation, but surprisingly, they contained extensive microtubule arrays. Moderately affected embryos contained bipolar spindles with dense and long astral microtubule arrays but with poorly organized kinetochore and interpolar microtubules. Severely affected embryos contained collapsed spindles with numerous long astral microtubules. Our results suggest that γ-tubulin is not absolutely required for microtubule nucleation in C. elegans but is required for the normal organization and function of kinetochore and interpolar microtubules.