Dynamics of the head-neck complex in response to the trunk horizontal vibration: modeling and identification

Although many studies exist concerning the influence of seat vibration on the head in the seated human body, the dynamic response of the head-neck complex (HNC) to the trunk vibration has not been well investigated. Little quantitative knowledge exists about viscoelastic parameters of the neck. In this study, the dynamics of the HNC is identified when it is exposed to the trunk horizontal (fore-and-aft) vibration. The frequency response functions between the HNC angular velocity and the trunk horizontal acceleration, corresponding to four volunteers, are obtained in the frequency range of 0.5 Hz to 10 Hz. A fourth-order mathematical model, derived by considering a double-inverted-pendulum model for the HNC, is designed to simulate the dynamic response of the HNC to the trunk horizontal vibration. The frequency domain identification method is used to determine the coefficients of the mathematical model of the HNC. Good agreement has been obtained between experimental and simulation results. This indicates that the system, similar to the designed fourth-order model, has mainly two resonance frequencies. The viscoelastic parameters of the neck, including the spring and damping coefficients, are then obtained by use of the optimization method.     

Identification of the head-neck complex in response to trunk horizontal vibration

A method is proposed for identifying the head-neck complex (HNC) in the seated human body when it is exposed to the trunk horizontal (fore-and-aft) vibration. It is assumed that the HNC only has the anteroposterior (flexion/extension) motion in the sagittal plane. An electrohydraulic vibrator is used as a source of vibration. To generate the trunk horizontal vibration, the trunk of the seated subject is fixed to the seatback. The subjects are exposed to the random vibration at a magnitude of 1.60 ms^-2 rms (root-mean-square) for 50 s. The coherence and frequency response function are then obtained in the frequency range 0.5–3 Hz. The results show that the HNC behavior is quasilinear with a resonance frequency between 1 and 1.4 Hz. Accordingly, a two-dimensional single-inverted pendulum is considered as a model for the HNC. The frequency domain identification method is then used to estimate the unknown parameters, including the HNC viscoelastic and inertia parameters. The model is examined in a time domain using the random vibration. Good agreement is obtained between experimental and simulation results, indicating the reliability of the proposed method.     

Dynamics of the linuron hydrolase libA gene pool size in response to linuron application and environmental perturbations in agricultural soil and on-farm biopurification systems

libA, a gene encoding a novel type of linuron hydrolase, was recently identified in the linuron-mineralizing Variovorax sp. strain SRS16. In order to assess the contribution of libA to linuron degradation in environmental settings, libA abundance was monitored in response to the application of linuron and to environmental perturbations in agricultural soil microcosms and microcosms simulating the matrix of on-farm biopurification systems. libA numbers were measured by real-time PCR and linked to reported data of Variovorax community composition and linuron mineralization capacity. In the soil microcosms and one biopurification system setup, libA numbers responded to the application of linuron and environmental changes in congruency with the modulation of linuron mineralization capacity and the occurrence of a particular Variovorax phylotype (phylotype A). However, in another biopurification system setup, no such correlations were found. Our data suggest that in the simulated environmental settings, the occurrence of libA can be linked to the linuron mineralization capacity and that libA is primarily hosted by Variovorax phylotype A strains. However, the results also suggest that, apart from libA, other, as-yet-unknown isofunctional genes play an important role in linuron mineralization in the environment.     

Local and global cycles in an acarine predator-prey system: a frequency domain analysis

Techniques of time series analysis are applied to the dynamics of the phytophagous miteTetranychus urticae Koch and its predatorTyphlodromus occidentalis Nesbitt in six experimental mini-orchards, sampled weekly during two years (van de Klashorst et al., 1992). Autocovariance and crosscovariance functions characterize local and global behaviour in the time domain. Spectral density function and cross amplitude spectrum provide important information of the system's behaviour in the frequency domain. Because of a strikingly different behaviour, the total record was cut into two distinct periods. During the first period, all local systems oscillate with a frequency of four cycles per year and all in phase, resulting in a strongly periodical global behaviour at a high mean level. Since the mean local power spectra coincide with the global power spectrum over a wide range of frequencies, it is concluded that the total system of six orchards is homogeneous. During the second period, all local power spectra are mutually different, more or less smoothed with no apparent peak pointing to a periodical component. The resulting global power spectrum is almost flat at a low mean level. The system is heterogeneous with a quasi-stable global behaviour.From the cross-amplitude spectrum it became clear that the dynamics of the two species within the orchards remained strongly coupled over a wide range of frequencies.Since the experimental circumstances had not been intentionally changed, the origin of the drastic change in behaviour could not be identified. If high rates of dispersal have been the synchronizer during the first period, it is not clear why asynchrony suddenly occured in the second period.A comparison of the overall record with the results of Nachman's (1987) stochastic simulation model suggests that the change could possibly be of stochastic origin.     

Three-dimensional response of a lumped parameter head-neck model due to impact and impulsive loading

A numerical procedure developed previously for predicting sagittal-plane motion of the human head-neck system due to impact and impulsive loading has been extended to three dimensions. In both situations, a lumped parameter approach is employed, but the current model lumps the mechanical response of each intervertebral joint into a single force-deformation relation evaluated from mechanical properties assembled by various investigators. Computations were performed to obtain the response of the model to a two-dimensional case of flexion whiplash, to one three-dimensional case of side impact to the skull and to another involving base acceleration normal to the sagittal plane. Agreement of the kinematic variables with the results of both the previous two-dimensional analysis and experimental data from a volunteer run is satisfactory, but somewhat poorer correspondence was found for the three-dimensional predictions upon comparison with data obtained from a physical model and from a volunteer when subjected to the prescribed loading. The differences in response are attributed to higher stiffness of facet separation of the model relative both to the structure and the volunteer, to insufficient damping, as well as to substantial differences in the mechanical deformation characteristic of the components of the prototypes and the numerical model.     

Reflex delays of the trunk muscles in response to postural perturbations: A reliability study

It has been reported that altered neuromuscular control of the trunk is associated with lower back pain. In this context reflex delays of the trunk muscles have often been assessed but the reliability of the tests has not been well established. The aim of this study was to test the reliability of measuring reflex delays of the trunk muscles after two types of postural perturbations. 24 Healthy subjects participated in the intra-session study and 13 of them repeated the test protocol within 1–3 weeks, to determine inter-session reliability. Postural reflex delays to unexpected loading and unloading of the arms were assessed in a standing unrestrained position. Each subject performed 40 trials of each test in order to evaluate muscle responses of 5 trunk muscles using surface electromyography. Overall reliability increased with higher number of the averaged trials. Good intra-session (ICC_3,1>0.75) and moderate (ICC_3,1>0.60) inter-session reliability were reached in most of the monitored trunk muscles. Within the performed number of trials we did not observe any significant systematic intra- or inter-session bias effect. Averaging a higher number of consecutive trials would be recommended in future research and clinical practice.     

Identification of small scale biochemical networks based on general type system perturbations

New technologies enable acquisition of large data-sets containing genomic, proteomic and metabolic information that describe the state of a cell. These data-sets call for systematic methods enabling relevant information about the inner workings of the cell to be extracted. One important issue at hand is the understanding of the functional interactions between genes, proteins and metabolites. We here present a method for identifying the dynamic interactions between biochemical components within the cell, in the vicinity of a steady-state. Key features of the proposed method are that it can deal with data obtained under perturbations of any system parameter, not only concentrations of specific components, and that the direct effect of the perturbations does not need to be known. This is important as concentration perturbations are often difficult to perform in biochemical systems and the specific effects of general type perturbations are usually highly uncertain, or unknown. The basis of the method is a linear least-squares estimation, using time-series measurements of concentrations and expression profiles, in which system states and parameter perturbations are estimated simultaneously. An important side-effect of also employing estimation of the parameter perturbations is that knowledge of the system's steady-state concentrations, or activities, is not required and that deviations from steady-state prior to the perturbation can be dealt with. Time derivatives are computed using a zero-order hold discretization, shown to yield significant improvements over the widely used Euler approximation. We also show how network interactions with dynamics that are too fast to be captured within the available sampling time can be determined and excluded from the network identification. Known and unknown moiety conservation relationships can be processed in the same manner. The method requires that the number of samples equals at least the number of network components and, hence, is at present restricted to relatively small-scale networks. We demonstrate herein the performance of the method on two small-scale in silico genetic networks.     

The response of oscillating glycolysis to perturbations in the NADH/NAD system: a comparison between experiments and a computer model.

The glycolytic pathway is described by a set of coupled non linear differential equations of first order with respect to time. The individual terms of these equations consist of enzyme velocities assuming a steady state hypothesis for the enzymatic forms. These are specified and the system is solved numerically. Oscillations are explained by interaction of PFK with the adenylate system. The conditions for the occurrence of oscillations are tested in a series of computer runs. The phase relations between intermediates of the model agree with those found in yeast cells. As an application of the model the disturbation of oscillations by the addition of acetaldehyde is simulated. The predictions of the model agree with experimental results.     

The mechanical behavior of the human forearm in response to transient perturbations

Static and dynamic components of mechanical impedance of human forearm were evaluated by applying two kinds of perturbations: 1) large viscoelastic loads, and 2) small pseudo-random perturbations. When the task involved the active resistance of the perturbations, both stiffness and viscosity increased relatively to their values in the passive task, the increment in stiffness being larger than that in viscosity. The time course of such changes was investigated during the transition between the two operating points defined by the instructions do not resist and resist the applied perturbations. The changes in stiffness and viscosity were relatively slow, those in the latter lagging behind those in the former.     

Uncovering interactions in the frequency domain

Oscillatory activity plays a critical role in regulating biological processes at levels ranging from subcellular, cellular, and network to the whole organism, and often involves a large number of interacting elements. We shed light on this issue by introducing a novel approach called partial Granger causality to reliably reveal interaction patterns in multivariate data with exogenous inputs and latent variables in the frequency domain. The method is extensively tested with toy models, and successfully applied to experimental datasets, including (1) gene microarray data of HeLa cell cycle; (2) in vivo multi-electrode array (MEA) local field potentials (LFPs) recorded from the inferotemporal cortex of a sheep; and (3) in vivo LFPs recorded from distributed sites in the right hemisphere of a macaque monkey.     

Canonical sensitivities: a useful tool to deal with large perturbations in metabolic network modeling

The dynamic range of metabolic models can be extended to deal with large perturbations by introducing the related concepts of "generalized" kinetic order and "canonical" sensitivities. Generalized kinetic orders are built as a well-defined non linear combination of the canonical sensitivities coefficients, which in turn are obtained by a least-squares regression on central composite factorial design data. In a such way, the whole domain of the operating variables is mapped without need to determine locally neither the first nor the second order model derivatives. The method was validated through numerical simulations, its predictions being compared with those coming from a Michaelis-Menten formalism taken as reference. In parallel, two variants of the Power-law formalism (S-system, least-squares GMA) also were tested. The canonical sensitivities method produced the widest range to predict metabolite concentrations and metabolic fluxes at the steady states. In addition, the variation pattern for the logarithmic gains and for the characteristic eigenvalues have been accurately determined from a unique overall model, being both required to make realistic analysis in metabolic engineering. The achieved information also can be expressed in terms of those typical coefficients derived from the Metabolic Control Analysis (MCA). Even if current first order Power-law or MCA formalisms were used, the canonical sensitivities approach provides a significant advantage, since complete sets of homologous, accurate, locally valid metabolic coefficients can be simultaneously recovered from the array proposed, being representative of the whole range of the operating variables instead of a unique nominal condition as is usual.     

Lipid phase perturbations and the unfolded protein response

Recent studies of the consequences of ganglioside accumulation in lysosomal storage disease and free cholesterol accumulation in cell membranes in atherosclerosis suggest an unexpected link between perturbation of the endoplasmic reticulum membrane's lipid phase, induction of the unfolded protein response, and cell death.     

Genome-wide gene expression profiles in response to plastid division perturbations

Plastids are vital organelles involved in important metabolic functions that directly affect plant growth and development. Plastids divide by binary fission involving the coordination of numerous protein components. A tight control of the plastid division process ensures that: there is a full plastid complement during and after cell division, specialized cell types have optimal plastid numbers; the division rate is modulated in response to stress, metabolic fluxes and developmental status. However, how this control is exerted by the host nucleus is unclear. Here, we report a genome-wide microarray analysis of three accumulation and replication of chloroplasts (arc) mutants that show a spectrum of altered plastid division characteristics. To ensure a comprehensive data set, we selected arc3, arc5 and arc11 because they harbour mutations in protein components of both the stromal and cytosolic division machinery, are of different evolutionary origin and display different phenotypic severities in terms of chloroplast number, size and volume. We show that a surprisingly low number of genes are affected by altered plastid division status, but that the affected genes encode proteins important for a variety of fundamental plant processes.     

Noninvasive methodology for analyzing the response of the system to electromagnetic fields of industrial frequency

Characteristics of the main biochemical components in surface structures and whole cells of unicellular green algae D. tertiolecta and T. viridis before and after the influence of an electromagnetic field (EMF) of industrial frequency were obtained by the method of noninvasive analysis in nonpolarized light. We showed fundamental differences in response to the influence of EMF among cells of algae cultures, which manifest themselves in changes in the main biochemical components and their concentration gradients. These changes led to the transformation of metabolism.     

Evaluation of viscoelastic property inhomogeneities in soft tissues exposed to low frequency perturbations

This work presents the methods of recognition of inhomogeneities of tissue shear viscoelastic properties using partial data on the internal displacements in an object exposed to low-frequency perturbation. An approach to detect tissue inhomogeneities using the single displacement component is presented.     

Semi-stochastic cell-level computational modeling of the immune system response to bacterial infections and the effects of antibiotics

A mathematical model for the immune system response to bacterial infections is proposed. The formalism is based on modeling the chemokine-determined transmigration of leukocytes from a venule through the venule walls and the subsequent in-tissue migration and engulfment of the pathogens that are responsible for the infection. The model is based on basic principles, such as Poiseuille blood flow through the venule, fundamental solutions of the diffusion–reaction equation for the concentration field of pathogen-released chemokines, linear chemotaxis of the leukocytes, random walk of pathogens, and stochastic processes for the death and division of pathogens. Thereby, a computationally tractable and, as far as we know, original framework has been obtained, which is used to incorporate the interaction of a substantial number of leukocytes and thereby to unravel the significance of biological processes and parameters regarding the immune system response. The developed model provides a neat way for visualization of the biophysical mechanism of the immune system response. The simulations indicate a weak correlation between the immune system response in terms of bacterial clearing time and the leukocyte stiffness, and a significant decrease in the clearing time with increasing in-blood leukocyte density, decreasing pathogen motility, and increasing venule wall transmissivity. Finally, the increase in the pathogen death rate and decrease in pathogen motility induce a decrease in the clearing time of the infection. The adjustment of the latter two quantities mimic the administration of antibiotics.     

Dependence of peripheral tremor on mechanical perturbations: a modeling study

The present study scrutinizes the popular view that tremors of central origin but not those of peripheral origin are largely resistant to mechanical perturbations. We explore the effects of perturbations in a well-established model of peripheral tremor and document that (a) tremor frequency can remain unchanged when spring or weight loads are added, (b) entrainment by external drives can be limited to drives of similar frequency, and (c) resetting of tremor phase by torque pulses can remain fractional. This resistance to mechanical perturbations arises in the model because peripheral neuromuscular dynamics act as a limit-cycle oscillator which, by its very nature, will absorb moderate changes to signals and parameters. We conclude from our study that resistance to mechanical perturbations is not an exclusive property of central tremors, but rather may also be found in peripheral tremors. Other criteria are therefore needed to distinguish between different origins of tremor. Received: 8 April 1998 / Accepted in revised form: 6 October 1998     

Sympathetic nervous system representation in time and frequency domain indices of heart rate variability

This study evaluated the contributions of sympathetic and parasympathetic modulation to heart rate variability during situations in which vagal and sympathetic tone predominated. In a placebo-controlled, randomized, double blind blockade study, six young healthy male individuals received propranolol (0.2 mg x kg(-1)), atropine (0.04 mg x kg(-1)), propranolol plus atropine, or placebo infusions over 4 days. Time-domain indices were calculated during 40 min of rest and 20 min of exercise at 70% of maximal exercise intensity. Spectrum analysis, using fast Fourier transformation, was also performed at rest and during the exercise. The time-domain indices standard deviation of R-R intervals, mean of the standard deviations of all R-R intervals for all 5-min segments, percentage of number of pairs of adjacent R-R intervals differing by more than 50 ms, and square root of the mean of the sum of squares of differences between adjacent R-R intervals were reduced after atropine and propranolol plus atropine. Propranolol alone caused no appreciable change in any of the time-domain indices. At rest, all spectrum components were similar after placebo and propranolol infusions, but following parasympathetic and double autonomic blockade there was a reduction in all components of the spectrum analysis, except for the low:high ratio. During exercise, partial and double blockade did not change significantly any of the spectrum components. Thus, time and frequency-domain indices of heart rate variability were able to detect vagal activity, but could not detect sympathetic activity. During exercise, spectrum analysis is not capable of evaluating autonomic modulation of heart rate.     

Color vision in the comb frequency domain

In 1982, Horace Barlow considered the question of human trichromacy in the context of information theory: according to the Sampling Theorem, three types of receptors covering the visible spectrum (400-700 nm) might be sufficient to reconstruct the color signal. Although Barlow was led to reject the direct application of the Sampling Theorem to explain color dimensionality, the theoretical framework offers a fresh point of view for analyzing the color system in conjunction with the physical characteristics of natural color signals. This review aims to illustrate that if the strict mathematical reconstruction (as implied by the Sampling Theorem) is replaced by a pragmatic approximation of color signals, then trichromacy, with its subsequent opponent-color process, could be regarded as an optimization of color constancy abilities in the spectral environment of primates. Higher dimension systems (tetrachromacy) found in other species can also serve the purpose of color constancy optimization in environments where color signals exhibit a finer spectral structure.