Control of time-dependent biological processes by temporally patterned input

Temporal patterning of biological variables, in the form of oscillations and rhythms on many time scales, is ubiquitous. Altering the temporal pattern of an input variable greatly affects the output of many biological processes. We develop here a conceptual framework for a quantitative understanding of such pattern dependence, focusing particularly on nonlinear, saturable, time-dependent processes that abound in biophysics, biochemistry, and physiology. We show theoretically that pattern dependence is governed by the nonlinearity of the input–output transformation as well as its time constant. As a result, only patterns on certain time scales permit the expression of pattern dependence, and processes with different time constants can respond preferentially to different patterns. This has implications for temporal coding and decoding, and allows differential control of processes through pattern. We show how pattern dependence can be quantitatively predicted using only information from steady, unpatterned input. To apply our ideas, we analyze, in an experimental example, how muscle contraction depends on the pattern of motorneuron firing.

Time-dependent vascular regression and permeability changes in established human tumor xenografts induced by an anti-vascular endothelial growth factor/vascular permeability factor antibody

The hyperpermeability of tumor vessels to macromolecules, compared with normal vessels, is presumably due to vascular endothelial growth factor/vascular permeability factor (VEGF/VPF) released by neoplastic and/or host cells. In addition, VEGF/VPF is a potent angiogenic factor. Removal of this growth factor may reduce the permeability and inhibit tumor angiogenesis. To test these hypotheses, we transplanted a human glioblastoma (U87), a human colon adenocarcinoma (LS174T), and a human melanoma (P-MEL) into two locations in immunodeficient mice: the cranial window and the dorsal skinfold chamber. The mice bearing vascularized tumors were treated with a bolus (0.2 ml) of either a neutralizing antibody (A4.6.1) (492 μg/ml) against VEGF/VPF or PBS (control). We found that tumor vascular permeability to albumin in antibody-treated groups was lower than in the matched controls and that the effect of the antibody was time-dependent and influenced by the mode of injection. Tumor vascular permeability did not respond to i.p. injection of the antibody until 4 days posttreatment. However, the permeability was reduced within 6 h after i.v. injection of the same amount of antibody. In addition to the reduction in vascular permeability, the tumor vessels became smaller in diameter and less tortuous after antibody injections and eventually disappeared from the surface after four consecutive treatments in U87 tumors. These results demonstrate that tumor vascular permeability can be reduced by neutralization of endogenous VEGF/VPF and suggest that angiogenesis and the maintenance of integrity of tumor vessels require the presence of VEGF/VPF in the tissue microenvironment. The latter finding reveals a new mechanism of tumor vessel regression—i.e., blocking the interactions between VEGF/VPF and endothelial cells or inhibiting VEGF/VPF synthesis in solid tumors causes dramatic reduction in vessel diameter, which may block the passage of blood elements and thus lead to vascular regression.

Temperature, ordering, and equilibrium with time-dependent confining forces

The concepts of temperature and equilibrium are not well defined in systems of particles with time-varying external forces. An example is a radio frequency ion trap, with the ions laser cooled into an ordered solid, characteristic of sub-mK temperatures, whereas the kinetic energies associated with the fast coherent motion in the trap are up to 7 orders of magnitude higher. Simulations with 1,000 ions reach equilibrium between the degrees of freedom when only aperiodic displacements (secular motion) are considered. The coupling of the periodic driven motion associated with the confinement to the nonperiodic random motion of the ions is very small at low temperatures and increases quadratically with temperature.

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.

Detection of Ca2+-Dependent Transglutaminase Activity in Root and Leaf Tissue of Monocotyledonous and Dicotyledonous Plants

Protein extracted from root and leaf tissue of the dicotyledonous plants pea (Pisum sativum) and broad bean (Vicia faba) and the monocotyledonous plants wheat (Triticum aestivum) and barley (Hordeum vulgare) were shown to catalyze the incorporation of biotin-labeled cadaverine into microtiter-plate-bound N′,N′-dimethylcasein and the cross-linking of biotin-labeled casein to microtiter-plate-bound casein in a Ca2+-dependent manner. The cross-linking of biotinylated casein and the incorporation of biotin-labeled cadaverine into N′,N′-dimethylcasein were time-dependent reactions with a pH optimum of 7.9. Transglutaminase activity was shown to increase over a 2-week growth period in both the roots and leaves of pea. The product of transglutaminase's protein-cross-linking activity, ε-(γ-glutamyl)-lysine isodipeptide, was detected in root and shoot protein from pea, broad bean, wheat, and barley by cation-exchange chromatography. The presence of the isodipeptide was confirmed by reversed-phase chromatography. Hydrolysis of the isodipeptide after cation-exchange chromatography confirmed the presence of glutamate and lysine.

A joint hazard and time scaling model to compare survival curves.

To provide a more general method for comparing survival experience, we propose a model that independently scales both hazard and time dimensions. To test the curve shape similarity of two time-dependent hazards, h1(t) and h2(t), we apply the proposed hazard relationship, h12(tKt)/ h1(t) = Kh, to h1. This relationship doubly scales h1 by the constant hazard and time scale factors, Kh and Kt, producing a transformed hazard, h12, with the same underlying curve shape as h1. We optimize the match of h12 to h2 by adjusting Kh and Kt. The corresponding survival relationship S12(tKt) = [S1(t)]KtKh transforms S1 into a new curve S12 of the same underlying shape that can be matched to the original S2. We apply this model to the curves for regional and local breast cancer contained in the National Cancer Institute's End Results Registry (1950-1973). Scaling the original regional curves, h1 and S1 with Kt = 1.769 and Kh = 0.263 produces transformed curves h12 and S12 that display congruence with the respective local curves, h2 and S2. This similarity of curve shapes suggests the application of the more complete curve shapes for regional disease as templates to predict the long-term survival pattern for local disease. By extension, this similarity raises the possibility of scaling early data for clinical trial curves according to templates of registry or previous trial curves, projecting long-term outcomes and reducing costs. The proposed model includes as special cases the widely used proportional hazards (Kt = 1) and accelerated life (KtKh = 1) models.

Implications of fault constitutive properties for earthquake prediction.

The rate- and state-dependent constitutive formulation for fault slip characterizes an exceptional variety of materials over a wide range of sliding conditions. This formulation provides a unified representation of diverse sliding phenomena including slip weakening over a characteristic sliding distance Dc, apparent fracture energy at a rupture front, time-dependent healing after rapid slip, and various other transient and slip rate effects. Laboratory observations and theoretical models both indicate that earthquake nucleation is accompanied by long intervals of accelerating slip. Strains from the nucleation process on buried faults generally could not be detected if laboratory values of Dc apply to faults in nature. However, scaling of Dc is presently an open question and the possibility exists that measurable premonitory creep may precede some earthquakes. Earthquake activity is modeled as a sequence of earthquake nucleation events. In this model, earthquake clustering arises from sensitivity of nucleation times to the stress changes induced by prior earthquakes. The model gives the characteristic Omori aftershock decay law and assigns physical interpretation to aftershock parameters. The seismicity formulation predicts large changes of earthquake probabilities result from stress changes. Two mechanisms for foreshocks are proposed that describe observed frequency of occurrence of foreshock-mainshock pairs by time and magnitude. With the first mechanism, foreshocks represent a manifestation of earthquake clustering in which the stress change at the time of the foreshock increases the probability of earthquakes at all magnitudes including the eventual mainshock. With the second model, accelerating fault slip on the mainshock nucleation zone triggers foreshocks.