An endogenous rate of time preference, the Penrose effect, and dynamic optimality of environmental quality.

In the present paper, the endogenous theory of time preference is extended to analyze those processes of capital accumulation and changes in environmental quality that are dynamically optimum with respect to the intertemporal preference ordering of the representative individual of the society in question. The analysis is carried out within the conceptual framework of the dynamic analysis of environmental quality, as has been developed by a number of economists for specific cases of the fisheries and forestry commons. The duality principles on intertemporal preference ordering and capital accumulation are extended to the situation where processes of capital accumulation are subject to the Penrose effect, which exhibit the marginal decrease in the effect of investment in private and social overhead capital upon the rate at which capital is accumulated. The dynamically optimum time-path of economic activities is characterized by the proportionality of two systems of imputed, or efficient, prices, one associated with the given intertemporal ordering and another associated with processes of accumulation of private and social overhead capital. It is particularly shown that the dynamically optimality of the processes of capital accumulation involving both private and social overhead capital is characterized by the conditions that are identical with those involving private capital, with the role of social overhead capital only indirectly exhibited.

Dynamic mapping at the laminar level of odor-elicited responses in rat olfactory bulb by functional MRI

We have applied functional MRI (fMRI) based on blood oxygenation level-dependent (BOLD) image-contrast to map odor-elicited olfactory responses at the laminar level in the rat olfactory bulb (OB) elicited by iso-amyl acetate (10−2 dilution of saturated vapor) with spatial and temporal resolutions of 220×220×1,000 μm and 36 s. The laminar structure of the OB was clearly depicted by high-resolution in vivo anatomical MRI with spatial resolution of 110×110×1,000 μm. In repeated BOLD fMRI measurements, highly significant (P < 0.001) foci were located in the outer layers of both OBs. The occurrence of focal OB activity within a domain at the level of individual glomeruli or groups of glomeruli was corroborated on an intra- and inter-animal basis under anesthetized conditions with this noninvasive method. The dynamic studies demonstrated that the odor-elicited BOLD activations were highly reproducible on a time scale of minutes, whereas over tens of minutes the activations sometimes varied slowly. We found large BOLD signal (ΔS/S = 10–30%) arising from the olfactory nerve layer, which is devoid of synapses and composed of unmyelinated fibers and glial cells. Our results support previous studies with other methods showing that odors elicit activity within glomerular layer domains in the mammalian OB, and extend the analysis to shorter time periods at the level of individual glomeruli or groups of glomeruli. With further improvement, BOLD fMRI should be ideal for systematic analysis of the functional significance of individual glomeruli in olfactory information encoding and of spatiotemporal processing within the olfactory system.

γ-Aminobutyric acid type B receptor-dependent burst-firing in thalamic neurons: A dynamic clamp study

Synchronized network responses in thalamus depend on phasic inhibition originating in the thalamic reticular nucleus (nRt) and are mediated by the neurotransmitter γ-aminobutyric acid (GABA). A suggested role for intra-nRt connectivity in inhibitory phasing remains controversial. Recently, functional GABA type B (GABAB) receptors were demonstrated on nRt cells, and the slow time course of the GABAB synaptic response seems ideally suited to deinactivate low-threshold calcium channels. This promotes burst firing, a characteristic feature of synchronized responses. Here we investigate GABAB-mediated rebound burst firing in thalamic cells. Whole-cell current-clamp recordings were obtained from nRt cells and somatosensory thalamocortical relay cells in rat brain slices. Synthetic GABAB inhibitory postsynaptic potentials, generated by a hybrid computer–neuron synapse (dynamic clamp), triggered rebound low-threshold calcium spikes in both cell types when peak inhibitory postsynaptic potential hyperpolarization was greater than −92 mV. The threshold inhibitory postsynaptic potential conductance for rebound burst generation was comparable in nRt (7 nS) and thalamocortical (5 nS) cells. However, burst onset in nRt (1 s) was considerably delayed compared with thalamocortical (0.6 s) cells. Thus, GABAB inhibitory postsynaptic potentials can elicit low-threshold calcium spikes in both relay and nRt neurons, but the resultant oscillation frequency would be faster for thalamocortical–nRt networks (3 Hz) than for nRt–nRt networks (1–2 Hz). We conclude, therefore, that fast (>2 Hz) GABAB-dependent thalamic oscillations are maintained primarily by reciprocal connections between excitatory and inhibitory cells. These findings further indicate that when oscillatory neural networks contain both recurrent and reciprocal inhibition, then distinct population frequencies may result when one or the other type of inhibition is favored.

The Drosophila melanogaster Homologue of the Xeroderma Pigmentosum D Gene Product Is Located in Euchromatic Regions and Has a Dynamic Response to UV Light-induced Lesions in Polytene Chromosomes

The XPD/ERCC2/Rad3 gene is required for excision repair of UV-damaged DNA and is an important component of nucleotide excision repair. Mutations in the XPD gene generate the cancer-prone syndrome, xeroderma pigmentosum, Cockayne’s syndrome, and trichothiodystrophy. XPD has a 5′- to 3′-helicase activity and is a component of the TFIIH transcription factor, which is essential for RNA polymerase II elongation. We present here the characterization of the Drosophila melanogaster XPD gene (DmXPD). DmXPD encodes a product that is highly related to its human homologue. The DmXPD protein is ubiquitous during development. In embryos at the syncytial blastoderm stage, DmXPD is cytoplasmic. At the onset of transcription in somatic cells and during gastrulation in germ cells, DmXPD moves to the nuclei. Distribution analysis in polytene chromosomes shows that DmXPD is highly concentrated in the interbands, especially in the highly transcribed regions known as puffs. UV-light irradiation of third-instar larvae induces an increase in the signal intensity and in the number of sites where the DmXPD protein is located in polytene chromosomes, indicating that the DmXPD protein is recruited intensively in the chromosomes as a response to DNA damage. This is the first time that the response to DNA damage by UV-light irradiation can be visualized directly on the chromosomes using one of the TFIIH components.

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.

“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Solvent dependence of dynamic transitions in protein solutions

A transition as a function of increasing temperature from harmonic to anharmonic dynamics has been observed in globular proteins by using spectroscopic, scattering, and computer simulation techniques. We present here results of a dynamic neutron scattering analysis of the solvent dependence of the picosecond-time scale dynamic transition behavior of solutions of a simple single-subunit enzyme, xylanase. The protein is examined in powder form, in D2O, and in four two-component perdeuterated single-phase cryosolvents in which it is active and stable. The scattering profiles of the mixed solvent systems in the absence of protein are also determined. The general features of the dynamic transition behavior of the protein solutions follow those of the solvents. The dynamic transition in all of the mixed cryosolvent–protein systems is much more gradual than in pure D2O, consistent with a distribution of energy barriers. The differences between the dynamic behaviors of the various cryosolvent protein solutions themselves are remarkably small. The results are consistent with a picture in which the picosecond-time scale atomic dynamics respond strongly to melting of pure water solvent but are relatively invariant in cryosolvents of differing compositions and melting points.

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

Signals and signs in the nervous system: The dynamic anatomy of electrical activity is probably information-rich

The dichotomy between two groups of workers on neuroelectrical activity is retarding progress. To study the interrelations between neuronal unit spike activity and compound field potentials of cell populations is both unfashionable and technically challenging. Neither of the mutual disparagements is justified: that spikes are to higher functions as the alphabet is to Shakespeare and that slow field potentials are irrelevant epiphenomena. Spikes are not the basis of the neural code but of multiple codes that coexist with nonspike codes. Field potentials are mainly information-rich signs of underlying processes, but sometimes they are also signals for neighboring cells, that is, they exert influence. This paper concerns opportunities for new research with many channels of wide-band (spike and slow wave) recording. A wealth of structure in time and three-dimensional space is different at each scale—micro-, meso-, and macroactivity. The depth of our ignorance is emphasized to underline the opportunities for uncovering new principles. We cannot currently estimate the relative importance of spikes and synaptic communication vs. extrasynaptic graded signals. In spite of a preponderance of literature on the former, we must consider the latter as probably important. We are in a primitive stage of looking at the time series of wide-band voltages in the compound, local field, potentials and of choosing descriptors that discriminate appropriately among brain loci, states (functions), stages (ontogeny, senescence), and taxa (evolution). This is not surprising, since the brains in higher species are surely the most complex systems known. They must be the greatest reservoir of new discoveries in nature. The complexity should not deter us, but a dose of humility can stimulate the flow of imaginative juices.

Modeling and optimization of populations subject to time-dependent mutation.

It has become clear that many organisms possess the ability to regulate their mutation rate in response to environmental conditions. So the question of finding an optimal mutation rate must be replaced by that of finding an optimal mutation schedule. We show that this task cannot be accomplished with standard population-dynamic models. We then develop a "hybrid" model for populations experiencing time-dependent mutation that treats population growth as deterministic but the time of first appearance of new variants as stochastic. We show that the hybrid model agrees well with a Monte Carlo simulation. From this model, we derive a deterministic approximation, a "threshold" model, that is similar to standard population dynamic models but differs in the initial rate of generation of new mutants. We use these techniques to model antibody affinity maturation by somatic hypermutation. We had previously shown that the optimal mutation schedule for the deterministic threshold model is phasic, with periods of mutation between intervals of mutation-free growth. To establish the validity of this schedule, we now show that the phasic schedule that optimizes the deterministic threshold model significantly improves upon the best constant-rate schedule for the hybrid and Monte Carlo models.