Application of matched filtering to identify behavioral modulation of brain oscillations

Brain oscillations modulated by motor behaviors are coupled to steady-state and other potentially unrelated to movement oscillations, with energy in the same frequency bands as the signals of interest. We applied matched filtering, a quasi-optimum signal detection technique, to decouple and extract movement-related signals from local field potentials (LFPs) recorded in monkey motor cortical areas during the execution of a visually instructed reach-out task. Using a matched-filterbank, we examined coupling and interference of pre-movement and initial steady-state oscillations with movement-induced signals. Once these signal contributions were eliminated, we were able to identify significant correlations of the residual signals with behavioral parameters, which appeared attenuated by pre-movement signal interference in the raw LFPs. Specifically, the maximum and minimum amplitudes of filtered LFPs were directly modulated by peak movement velocity and micro-movements, respectively, identified in recorded hand velocity profiles. In addition, we identified phase correlations between signals during the delay (when the instructional cue was presented) and movement intervals, as well as modulation of LFP phase by movement direction. For pairs of orthogonal movement directions, corresponding LFP signals were consistently out of phase. Finally, β-band energy which is typically reduced during movement execution, possibly partly due to destructive interference between the modulated by behavior signal and unrelated oscillations, appeared to be recovered in the filtered signals.     

LaDIVA: A neurocomputational model providing laryngeal motor control for speech acquisition and production

Many voice disorders are the result of intricate neural and/or biomechanical impairments that are poorly understood. The limited knowledge of their etiological and pathophysiological mechanisms hampers effective clinical management. Behavioral studies have been used concurrently with computational models to better understand typical and pathological laryngeal motor control. Thus far, however, a unified computational framework that quantitatively integrates physiologically relevant models of phonation with the neural control of speech has not been developed. Here, we introduce LaDIVA, a novel neurocomputational model with physiologically based laryngeal motor control. We combined the DIVA model (an established neural network model of speech motor control) with the extended body-cover model (a physics-based vocal fold model). The resulting integrated model, LaDIVA, was validated by comparing its model simulations with behavioral responses to perturbations of auditory vocal fundamental frequency (f_o) feedback in adults with typical speech. LaDIVA demonstrated capability to simulate different modes of laryngeal motor control, ranging from short-term (i.e., reflexive) and long-term (i.e., adaptive) auditory feedback paradigms, to generating prosodic contours in speech. Simulations showed that LaDIVA’s laryngeal motor control displays properties of motor equivalence, i.e., LaDIVA could robustly generate compensatory responses to reflexive vocal f_o perturbations with varying initial laryngeal muscle activation levels leading to the same output. The model can also generate prosodic contours for studying laryngeal motor control in running speech. LaDIVA can expand the understanding of the physiology of human phonation to enable, for the first time, the investigation of causal effects of neural motor control in the fine structure of the vocal signal.     

Aging and immortality in a cell proliferation model

We investigate a model of cell division in which the length of telomeres within a cell regulates its proliferative potential. At each division, telomeres undergo a systematic length decrease as well as a superimposed fluctuation due to exchange of telomere DNA between the two daughter cells. A cell becomes senescent when one or more of its telomeres become shorter than a critical length. We map this telomere dynamics onto a biased branching-diffusion process with an absorbing boundary condition whenever any telomere reaches the critical length. Using first-passage ideas, we find a phase transition between finite lifetime and immortality (infinite proliferation) of the cell population as a function of the influence of telomere shortening, fluctuations, and cell division.     

A pharmacological investigation of the biphasic nature of the contractile response of rabbit and rat vas deferens to field stimulation

It has been demonstrated previously with the vas deferens of the guinea-pig that the first and second phases of the contractile response to motor nerve stimulation are preferentially antagonized by the P₂-purinoceptor antagonist arylazido aminopropionyl ATP (ANAPP₃), and the α₁-adrenoceptor antagonist prazosin, respectively. We have now investigated the effect of the two antagonists on the biphasic contraction in the vas deferens of two other species; rabbit and rat. ANAPP₃, in a concentration which antagonized responses to exogenously applied ATP but not those to exogenous norepinephrine, preferentially reduced the initial phasic response of the rabbit vas deferens to motor nerve stimulation without significantly reducing the secondary, tonic phase of the response. Prazosin had the opposite effect; antagonizing the response to norepinephrine but not to ATP and reducing the tonic response to motor nerve stimulation without significantly reducing the initial phasic response. Results obtained with the rat vas deferens were similar. The present results combined with previous findings suggest that ATP and norepinephrine act as cotransmitters in the vas deferens of several species.     

Non-Gaussian noise-optimized intracellular cytosolic calcium oscillations

We have numerically studied the effect of a particular kind of non-Gaussian colored noise (NGN), characterized by the deviation q from Gaussian noise (q = 1), on intracellular cytosolic calcium (Ca²⁺) oscillations. It is found that, as q is increased, the Ca²⁺ oscillation regularity increases and reaches a best performance at an optimal q, and then decreases with further increasing q, which represents the occurrence of coherence resonance, i.e., the most regular Ca²⁺ oscillations. Similar phenomena occur for different values of noise intensity and correlation time of the NGN. This phenomenon of deviation-optimized Ca²⁺ oscillations show that, external non-Gaussian noises of different types can enhance and even optimize the intracellular Ca²⁺ oscillations. This result provides new insights into the constructive roles and potential applications of non-Gaussian noises in intracellular cytosolic Ca²⁺ oscillations.     

Allosteric kinetics of pyruvate kinase of Saccharomyces carlsbergensis

The allosteric model of Monod et al. (1965) has been used to analyse the steadystate kinetics of pyruvate kinase from Saccharomyces carlsbergensis. The dissociation constants for the substrate phosphoenolpyruvate, the inhibitor ATP as well as the activator fructose-1, 6-diphosphate from the R and T state were calculated using a series of computer programs. On the basis of a crucial relation (derived in the Appendix), which correlates the Hill coefficient and the half-saturating concentration of substrate saturation curves with the parameters of the model of Monod et al., it is possible to differentiate between exclusive and non-exclusive ligand binding. On the other hand, this relation makes it possible to fit the experimental data to an extended model assuming only partially concerted transitions in each enzyme molecule.The physical data of yeast pyruvate kinase point to a tetrameric structure, whereas the steady-state kinetics favour a trimeric one. This discrepancy in the number of protomers can be overcome by the use of an extended model, which permits the occurrence of hybrid states R_tT_n?t. The introduction of one symmetrical hybrid state R₂T₂ into the model explains the kinetic data of yeast pyruvate kinase on the basis of four, probably identical, protomers. The equilibrium constants between the states are given.In the Appendix the derivation of the equation describing the occurrence of hybrid states is reported.     

The cerebellum and initiation of movement: the stretch reflex

Studies of the stretch reflex in decerebrate cats indicate a phase advance of peak sinusoidal tension in steady-state cycles between 0.1 and 10 Hz. This phase advance is reduced in acute and chronic cerebellectomy, as shown in previous investigations. Also, the augmentation of muscle peak tension in initial sinusoidal stretch cycles at 0.5-5 Hz has been found to be reduced during the time of reflex and motor instability in the several months following cerebellar ablation. This report shows the increased amplitude and phase lead of integrated electromyographic activity in initiating sinusoidal stretch cycles in the decerebrate cat. These reflex aspects are demonstrated in relation to the discharge of neurons in the dorsal spinocerebellar tract and of cerebellar cortical Purkinje cells in initial sinusoidal cycles. The intensity and phase advance of the discharge in dorsal spinocerebellar tract neurons is altered little, but these features are usually increased in Purkinje cells during initial stretches compared to continuous cycling. In terms of overall motor control, these findings are compatible with concepts of movement control, modulated by the cerebellum, in which the discharge of antagonist motor neurons is regulated in concert with that of agonist muscles upon initiation and termination of movement.     

The effects of starvation and 20-hydroxy-ecdysone on feeding and pupariation of early 3rd-instar Calliphora vicina larvae

At a larval body weight of 15–16 mg, young 3rd-instar Calliphora vicina larvae pass from the obligatory feeding phase, when starvation results in cessation of growth or development, to the facultative feeding phase when starvation leads to subsequent pupariation. During normal larval development the transition from obligatory to facultative feeding phase appears to be independent of 20-hydroxyecdysone but absence (or insufficiency) of the hormone is the likely reason why starved larvae in the obligatory feeding phase fail to pupariate. 20-Hydroxy-ecdysone induces termination of feeding by normal, facultative phase larvae and is the probable signal for larval commitment to pupariation.     

Three phase model of the processive motor protein kinesin

Kinesin is a stepping motor that successively produces forward and backward 8-nm steps along microtubules. Under physiological conditions, the steps powering kinesin's motility are biased in one direction and drive various biological motile processes. So far, the physical mechanism underlying the unidirectional bias of the kinesin is not fully understood. Recently, Martin Bier have provided a stepper model [Martin Bier, 2003, Processive motor protein as an overdamped Brownian stepper, Phys. Rev. Lett. 91, 148104], in which the stepping cycle of kinesin includes two distinguished phases: (i) a power stroke phase and (ii) a ratcheted diffusion phase which is characterized as a "random diffusional search". At saturating ATP level, this model can fit the experimental results accurately. In this paper, we'll provide a modified Brownian stepper model, in which the dependence of ATP concentration is considered. In our model, the stepping cycle of kinesin is distinguished into three phases: an ATP-binding phase, a power stroke phase and a ratcheted diffusion phase. This modified model can reconstruct most of the experimental results accurately.     

Aberrant Fat Metabolism in Caenorhabditis elegans Mutants with Defects in the Defecation Motor Program

The molecular mechanisms by which dietary fatty acids are absorbed by the intestine, and the way in which the process is regulated are poorly understood. In a genetic screen for mutations affecting fat accumulation in the intestine of Caenorhabditis elegans, nematode worms, we have isolated mutations in the aex-5 gene, which encodes a Kex2/subtilisin-family, Ca²⁺-sensitive proprotein convertase known to be required for maturation of certain neuropeptides, and for a discrete step in an ultradian rhythmic phenomenon called the defecation motor program. We demonstrate that aex-5 mutants have markedly lower steady-state levels of fat in the intestine, and that this defect is associated with a significant reduction in the rate at which labeled fatty acid derivatives are taken up from the intestinal lumen. Other mutations affecting the defecation motor program also affect steady-state levels of triglycerides, suggesting that the program is required per se for the proper accumulation of neutral lipids. Our results suggest that an important function of the defecation motor program in C. elegans is to promote the uptake of an important class of dietary nutrients. They also imply that modulation of the program might be one way in which worms adjust nutrient uptake in response to altered metabolic status.     

Phase-sensitive fluorescence spectroscopy: A new method to resolve fluorescence lifetimes or emission spectra of components in a mixture of fluorophores

A novel phase fluorometric method is described which permits direct recording of individual emission spectra from a mixture of two flourescent compounds. Additionally, the lifetimes of each component may be determined by examination of the phase-sensitive fluorescence spectra. The method utilizes phase-sensitive detection of the sinusoidally modulated emission from a phase fluorometer. Resolution of the individual emission spectra in the mixture requires different fluorescence lifetimes for each components. Determination of the individual lifetime requires knowledge of the steady-state emission spectra of the components. Use of low-frequency (≈ 10 Hz) cross-correlated signals eliminates the need for high-frequency frequency (≈10⁶ Hz) phase-sensitive detection. A mixture of 2-p-toluidinyl-6-naphthalenesulfonic acid (TNS) and 6-propionyl-2-(dimethylamino)naphthalene (PRODAN) was used to demonstrate the possibility of phase resolution of fluorophore mixture and to confirm theoretical predictions. A mixture of dibenzo[a,h]anthracene and dibenzo[c,g]carbazole was used to demonstrate that phase resolution is possible for spectra which overlap strongly and which are highly structured. In addition, the possibility of using phase-sensitive emission spectra for the resolution of excited-state reactions was demonstrated with anthracene and its diethylaniline exciplex. From a sample whose steady-state emission displayed both components we directly recorded the emission spectrum of anthracene monomer and the exciplex. For all these samples the dependence of the individual intensities on the phase angle of the detector agreed precisely with that expected on the basis of the individual fluorescence lifetimes. The detector phase angles chosen for suppression of each component in the mixture also agreed with the measured lifetimes. Thus, phase-sensitive fluorescence spectra can reveal individual spectral distributions or lifetimes. This method will be useful in the analysis fluorescence emissions which frequently occur from proteins, membranes and other biological samples.     

Fundamental Constraints on the Abundances of Chemotaxis Proteins

Flagellated bacteria, such as Escherichia coli, perform directed motion in gradients of concentration of attractants and repellents in a process called chemotaxis. The E. coli chemotaxis signaling pathway is a model for signal transduction, but it has unique features. We demonstrate that the need for fast signaling necessitates high abundances of the proteins involved in this pathway. We show that further constraints on the abundances of chemotaxis proteins arise from the requirements of self-assembly both of flagellar motors and of chemoreceptor arrays. All these constraints are specific to chemotaxis, and published data confirm that chemotaxis proteins tend to be more highly expressed than their homologs in other pathways. Employing a chemotaxis pathway model, we show that the gain of the pathway at the level of the response regulator CheY increases with overall chemotaxis protein abundances. This may explain why, at least in one E. coli strain, the abundance of all chemotaxis proteins is higher in media with lower nutrient content. We also demonstrate that the E. coli chemotaxis pathway is particularly robust to abundance variations of the motor protein FliM.     

Short-term enhancement of motor neuron synaptic exocytosis during early aging extends lifespan in Caenorhabditis elegans

Age-related mobility decline is often associated with negative physical and psychological outcomes, such as frailty, in the elderly population. In C. elegans, during the early stage of the aging process, a progressive deficit of synaptic exocytosis in the motor neurons results in a functional decline at the neuromuscular junctions, which eventually leads to degeneration of both neurons and muscles. This age-dependent functional decline can be ameliorated by pharmacological interventions, such as arecoline, a muscarinic AChR agonist known to promote synaptic exocytosis at the neuromuscular junctions. In this study, we found that a short-term treatment of arecoline during the early stage of aging, when the NMJ functional decline begins, not only slows muscle tissue aging, but also extends lifespan in C. elegans. We have also demonstrated that arecoline acts on the GAR-2/PLCβ pathway in the motor neurons to increases longevity. Together, our findings suggest that synaptic transmission in aging motor neurons may serve as a potential target for pharmacological interventions to promote both health span and lifespan, when applied at the early stage aging.Impact statementThe functional decline of motor activity is a common feature in almost all aging animals that leads to frailty, loss of independence, injury, and even death in the elderly population. Thus, understanding the molecular mechanism that drives the initial stage of this functional decline and developing strategies to increase human healthspan and even lifespan by targeting this process would be of great interests to the field. In this study, we found that by precisely targeting the motor neurons to potentiate its synaptic releases either genetically or pharmacologically, we can not only delay the functional aging at NMJs but also slow the rate of aging at the organismal level. Most importantly, we have demonstrated that a critical window of time, that is the early stage of NMJs functional decline, is required for the beneficial effects. A short-term treatment within this time period is sufficient to extend the animals’ lifespan.     

Formation of germichrysone by tissue cultures of Cassia torosa: Induction of secondary metabolism in the lag phase

A callus culture of Cassia torosa which produced germichrysone, an octaketide hydroanthracene, in high yield was established on Murashige-Skoog's medium containing IAA (3 ppm) and benzyladenine (0.1 ppm). In six-week-old callus culture the main pigment was pinselin and the germichrysone content was markedly decreased. This may indicate that pinselin was formed from germichrysone. A shake culture was established in liquid Murashige-Skoogs medium containing IAA and benzyladenine. The course of germichrysone production in relation to growth was investigated with the shake culture and the production of germichrysone was found to possess two maxima. The first maximum was observed in lag phase and the second coincided with active growth. Incorporation experiments with [1-¹⁴C] acetate clearly demonstrated that the secondary metabolism was induced in the lag phase.     

Steady-state,Pre-steady-state,and Single-turnover Kinetic Measurement for DNA Glycosylase Activity

Human 8-oxoguanine DNA glycosylase (OGG1) excises the mutagenic oxidative DNA lesion 8-oxo-7,8-dihydroguanine (8-oxoG) from DNA. Kinetic characterization of OGG1 is undertaken to measure the rates of 8-oxoG excision and product release. When the OGG1 concentration is lower than substrate DNA, time courses of product formation are biphasic; a rapid exponential phase (i.e. burst) of product formation is followed by a linear steady-state phase. The initial burst of product formation corresponds to the concentration of enzyme properly engaged on the substrate, and the burst amplitude depends on the concentration of enzyme. The first-order rate constant of the burst corresponds to the intrinsic rate of 8-oxoG excision and the slower steady-state rate measures the rate of product release (product DNA dissociation rate constant, k_off). Here, we describe steady-state, pre-steady-state, and single-turnover approaches to isolate and measure specific steps during OGG1 catalytic cycling. A fluorescent labeled lesion-containing oligonucleotide and purified OGG1 are used to facilitate precise kinetic measurements. Since low enzyme concentrations are used to make steady-state measurements, manual mixing of reagents and quenching of the reaction can be performed to ascertain the steady-state rate (k_off). Additionally, extrapolation of the steady-state rate to a point on the ordinate at zero time indicates that a burst of product formation occurred during the first turnover (i.e. y-intercept is positive). The first-order rate constant of the exponential burst phase can be measured using a rapid mixing and quenching technique that examines the amount of product formed at short time intervals (<1 sec) before the steady-state phase and corresponds to the rate of 8-oxoG excision (i.e. chemistry). The chemical step can also be measured using a single-turnover approach where catalytic cycling is prevented by saturating substrate DNA with enzyme (E>S). These approaches can measure elementary rate constants that influence the efficiency of removal of a DNA lesion.     

Nonequilibrium Steady-State Differences in Partial Pressure of CO2 and in Concentration of Weak Acids and Bases between Blood and Tissue

Several investigators have demonstrated that under conditions where little or no gas exchange occurs across the alveolar capillary membrane the PCO₂ is higher in the alveolus than in the mixed venous blood, whereas there are no PO₂ differences. Gurtner et al. have explained the ΔPCO₂ by a model in which H⁺ dissociation of proteins due to an electrical field caused by a negatively charged capillary wall (Wien effect) sets up an intracapillary PCO₂ difference between wall and bulk phase which is maintained by blood flow. The model is not specific for CO₂ and predicts that weak acids should be concentrated in a manner similar to CO₂ whereas weak bases should be relatively excluded from the alveolar space. Measurements of the steady-state distribution of the uncharged forms of the weak acids 5,5-dimethyloxyazoladinedione (DMO) and barbital and of the weak base tris(hydroxymethyl)aminomethane (THAM) between mixed venous blood and a fluid-filled lobe of lung were made in living dogs. The results agree fairly well with the predicted values.     

Two-phase strategy of controlling motor coordination determined by task performance optimality

A quantitative model of optimal coordination between hand transport and grip aperture has been derived in our previous studies of reach-to-grasp movements without utilizing explicit knowledge of the optimality criterion or motor plant dynamics. The model’s utility for experimental data analysis has been demonstrated. Here we show how to generalize this model for a broad class of reaching-type, goal-directed movements. The model allows for measuring the variability of motor coordination and studying its dependence on movement phase. The experimentally found characteristics of that dependence imply that execution noise is low and does not affect motor coordination significantly. From those characteristics it is inferred that the cost of neural computations required for information acquisition and processing is included in the criterion of task performance optimality as a function of precision demand for state estimation and decision making. The precision demand is an additional optimized control variable that regulates the amount of neurocomputational resources activated dynamically. It is shown that an optimal control strategy in this case comprises two different phases. During the initial phase, the cost of neural computations is significantly reduced at the expense of reducing the demand for their precision, which results in speed-accuracy tradeoff violation and significant inter-trial variability of motor coordination. During the final phase, neural computations and thus motor coordination are considerably more precise to reduce the cost of errors in making a contact with the target object. The generality of the optimal coordination model and the two-phase control strategy is illustrated on several diverse examples.