A mathematical model of integrin-mediated haptotactic cell migration

Haptotactic cell migration, a directed response to gradients of cell—extracellular matrix adhesion, is an important process in a number of biological phenomena such as wound healing and tumour cell invasion. Previously, mathematical models of haptotaxis have been developed on the premise that cells migrate in response to gradients in the density of the extracellular matrix. In this paper, we develop a novel mathematical model of haptotaxis which includes the adhesion receptors known as integrins and a description of their functional activation, local recruitment and protrusion as part of lamellipodia. Through the inclusion of integrins, the modelled cell matter is able to respond to a true gradient of cell–matrix adhesion, represented by functionally active integrins. We also show that previous matrix-mediated models are in fact a subset of the novel integrin-mediated models, characterised by specific choices of diffusion and haptotaxis coefficients in their model equations. Numerical solutions suggest the existence of travelling waves of cell migration that are confirmed via a phase plane analysis of a simplified model.

New approximation for free surface flow of groundwater: capillarity correction

An existing capillarity correction for free surface groundwater flow as modelled by the Boussinesq equation is re-investigated. Existing solutions, based on the shallow flow expansion, have considered only the zeroth-order approximation. Here, a second-order capillarity correction to tide-induced watertable fluctuations in a coastal aquifer adjacent to a sloping beach is derived. A new definition of the capillarity correction is proposed for small capillary fringes, and a simplified solution is derived. Comparisons of the two models show that the simplified model can be used in most cases. The significant effects of higher-order capillarity corrections on tidal fluctuations in a sloping beach are also demonstrated. (c) 2004 Elsevier Ltd. All rights reserved.

Novel strategies in feedforward adaptation to a position-dependent perturbation

To investigate the control mechanisms used in adapting to position-dependent forces, subjects performed 150 horizontal reaching movements over 25 cm in the presence of a position-dependent parabolic force field (PF). The PF acted only over the first 10 cm of the movement. On every fifth trial, a virtual mechanical guide (double wall) constrained subjects to move along a straight-line path between the start and target positions. Its purpose was to register lateral force to track formation of an internal model of the force field, and to look for evidence of possible alternative adaptive strategies. The force field produced a force to the right, which initially caused subjects to deviate in that direction. They reacted by producing deviations to the left, into the force field, as early as the second trial. Further adaptation resulted in rapid exponential reduction of kinematic error in the latter portion of the movement, where the greatest perturbation to the handpath was initially observed, whereas there was little modification of the handpath in the region where the PF was active. Significant force directed to counteract the PF was measured on the first guided trial, and was modified during the first half of the learning set. The total force impulse in the region of the PF increased throughout the learning trials, but it always remained less than that produced by the PF. The force profile did not resemble a mirror image of the PF in that it tended to be more trapezoidal than parabolic in shape. As in previous studies of force-field adaptation, we found that changes in muscle activation involved a general increase in the activity of all muscles, which increased arm stiffness, and selectively-greater increases in the activation of muscles which counteracted the PF. With training, activation was exponentially reduced, albeit more slowly than kinematic error. Progressive changes in kinematics and EMG occurred predominantly in the region of the workspace beyond the force field. We suggest that constraints on muscle mechanics limit the ability of the central nervous system to employ an inverse dynamics model to nullify impulse-like forces by generating mirror-image forces. Consequently, subjects adopted a strategy of slightly overcompensating for the first half of the force field, then allowing the force field to push them in the opposite direction. Muscle activity patterns in the region beyond the boundary of the force field were subsequently adjusted because of the relatively-slow response of the second-order mechanics of muscle impedance to the force impulse.

Effects of structural variation in xyloglucan polymers on interactions with bacterial cellulose

A cellulose/xyloglucan framework is considered to form the basis for the mechanical properties of primary plant cell walls and hence to have a major influence on the biomechanical properties of growing, fleshy plant tissues. In this study, structural variants of xyloglucan have been investigated as components of composites with bacterial cellulose as a simplified model for the cellulose/xyloglucan framework of primary plant cell walls. Evidence for molecular binding to cellulose with perturbation of cellulose crystallinity was found for all xyloglucan types. High molecular mass samples gave homogeneous centimeter-scale composites with extensive cross-linking of cellulose with xyloglucan. Lower molecular mass xyloglucans gave heterogeneous composites having a range of microscopic structures with little, if any, cross-linking. Xyloglucans with reduced levels of galactose substitution had evidence of self-association, competitive with cellulose binding. At comparable molecular mass, fucose substitution resulted in a modest promotion of microscopic features characteristic of primary cell walls. Taken together, the data are evidence that galactose substitution of the xyloglucan core structure is a major determinant of cellulose composite formation and properties, with additional fucose substitution acting as a secondary modulator. These conclusions are consistent with reported structural and mechanical properties of Arabidopsis mutants lacking specific facose and/or galactose residues.

Measuring and predicting the dynamics of linear monodisperse entangled polymers in rapid flow through an abrupt contraction. A small angle neutron scattering study

Small-angle neutron scattering measurements on a series of monodisperse linear entangled polystyrene melts in nonlinear flow through an abrupt 4:1 contraction have been made. Clear signatures of melt deformation and subsequent relaxation can be observed in the scattering patterns, which were taken along the centerline. These data are compared with the predictions of a recently derived molecular theory. Two levels of molecular theory are used: a detailed equation describing the evolution of molecular structure over all length scales relevant to the scattering data and a simplified version of the model, which is suitable for finite element computations. The velocity field for the complex melt flow is computed using the simplified model and scattering predictions are made by feeding these flow histories into the detailed model. The modeling quantitatively captures the full scattering intensity patterns over a broad range of data with independent variation of position within the contraction geometry, bulk flow rate and melt molecular weight. The study provides a strong, quantitative validation of current theoretical ideas concerning the microscopic dynamics of entangled polymers which builds upon existing comparisons with nonlinear mechanical stress data. Furthermore, we are able to confirm the appreciable length scale dependence of relaxation in polymer melts and highlight some wider implications of this phenomenon.

Investigation of factors related to employment outcome following traumatic brain injury: a critical review and conceptual model

Purpose : The purpose of this article is to critically review the literature to examine factors that are most consistently related to employment outcome following traumatic brain injury (TBI), with a particular focus on metacognitive skills. It also aims to develop a conceptual model of factors related to employment outcome. Method : The first stage of the review considered 85 studies published between 1980 and December 2003 which investigated factors associated with employment outcome following TBI. English-language studies were identified through searches of Medline and PsycINFO, as well as manual searches of journals and reference lists. The studies were evaluated and rated by two independent raters (Kappa = 0.835) according to the quality of their methodology based upon nine criteria. Fifty studies met the criteria for inclusion in the second stage of the review, which examined the relationship between a broad range of variables and employment outcome. Results : The factors most consistently associated with employment outcome included pre-injury occupational status, functional status at discharge, global cognitive functioning, perceptual ability, executive functioning, involvement in vocational rehabilitation services and emotional status. Conclusions : A conceptual model is presented which emphasises the importance of metacognitive, emotional and social environment factors for improving employment outcome.

Earthquake statistics in a Block Slider Model and a fully dynamic Fault Model

We examine the event statistics obtained from two differing simplified models for earthquake faults. The first model is a reproduction of the Block-Slider model of Carlson et al. (1991), a model often employed in seismicity studies. The second model is an elastodynamic fault model based upon the Lattice Solid Model (LSM) of Mora and Place (1994). We performed simulations in which the fault length was varied in each model and generated synthetic catalogs of event sizes and times. From these catalogs, we constructed interval event size distributions and inter-event time distributions. The larger, localised events in the Block-Slider model displayed the same scaling behaviour as events in the LSM however the distribution of inter-event times was markedly different. The analysis of both event size and inter-event time statistics is an effective method for comparative studies of differing simplified models for earthquake faults.

Links between tree species, symbiotic fungal diversity and ecosystem functioning in simplified tropical ecosystems

We studied the relationships among plant and arbuscular mycorrhizal (AM) fungal diversity, and their effects on ecosystem function, in a series of replicate tropical forestry plots in the La Selva Biological Station, Costa Rica. Forestry plots were 12 yr old and were either monocultures of three tree species, or polycultures of the tree species with two additional understory species. Relationships among the AM fungal spore community, host species, plant community diversity and ecosystem phosphorus-use efficiency (PUE) and net primary productivity (NPP) were assessed. Analysis of the relative abundance of AM fungal spores found that host tree species had a significant effect on the AM fungal community, as did host plant community diversity (monocultures vs polycultures). The Shannon diversity index of the AM fungal spore community differed significantly among the three host tree species, but was not significantly different between monoculture and polyculture plots. Over all the plots, significant positive relationships were found between AM fungal diversity and ecosystem NPP, and between AM fungal community evenness and PUE. Relative abundance of two of the dominant AM fungal species also showed significant correlations with NPP and PUE. We conclude that the AM fungal community composition in tropical forests is sensitive to host species, and provide evidence supporting the hypothesis that the diversity of AM fungi in tropical forests and ecosystem NPP covaries.

Degrees of freedom and motor planning in purposive movement

This paper presents empirical evidence suggesting that healthy humans can perform a two degree of freedom visuo-motor pursuit tracking task with the same response time delay as a one degree of freedom task. In contrast, the time delay of the response is influenced markedly by the nature of the motor synergy required to produce it. We suggest a conceptual account of this evidence based on adaptive model theory, which combines theories of intermittency from psychology and adaptive optimal control from engineering. The intermittent response planning stage has a fixed period. It possesses multiple optimal trajectory generators such that multiple degrees of freedom can be planned concurrently, without requiring an increase in the planning period. In tasks which require unfamiliar motor synergies, or are deemed to be incompatible, internal adaptive models representing movement dynamics are inaccurate. This means that the actual response which is produced will deviate from the one which is planned. For a given target-response discrepancy, corrective response trajectories of longer duration are planned, consistent with the principle of speed-accuracy trade-off. Compared to familiar or compatible tasks, this results in a longer response time delay and reduced accuracy. From the standpoint of the intermittency approach, the findings of this study help make possible a more integral and predictive account of purposive action. (c) 2005 Elsevier B.V. All rights reserved.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

A multiresolution terrain model for efficient visualization query processing

Multiresolution Triangular Mesh (MTM) models are widely used to improve the performance of large terrain visualization by replacing the original model with a simplified one. MTM models, which consist of both original and simplified data, are commonly stored in spatial database systems due to their size. The relatively slow access speed of disks makes data retrieval the bottleneck of such terrain visualization systems. Existing spatial access methods proposed to address this problem rely on main-memory MTM models, which leads to significant overhead during query processing. In this paper, we approach the problem from a new perspective and propose a novel MTM called direct mesh that is designed specifically for secondary storage. It supports available indexing methods natively and requires no modification to MTM structure. Experiment results, which are based on two real-world data sets, show an average performance improvement of 5-10 times over the existing methods.

A neuro-mechanical model for interpersonal coordination

The present study investigates the coordination between two people oscillating handheld pendulums, with a special emphasis on the influence of the mechanical properties of the effector systems involved. The first part of the study is an experiment in which eight pairs of participants are asked to coordinate the oscillation of their pendulum with the other participant's in an in-phase or antiphase fashion. Two types of pendulums, A and B, having different resonance frequencies (Freq A=0.98 Hz and Freq B=0.64 Hz), were used in different experimental combinations. Results confirm that the preferred frequencies produced by participants while manipulating each pendulum individually were close to the resonance frequencies of the pendulums. In their attempt to synchronize with one another, participants met at common frequencies that were influenced by the mechanical properties of the two pendulums involved. In agreement with previous studies, both the variability of the behavior and the shift in the intended relative phase were found to depend on the task-effector asymmetry, i.e., the difference between the mechanical properties of the effector systems involved. In the second part of the study, we propose a model to account for these results. The model consists of two cross-coupled neuro-mechanical units, each composed of a neural oscillator driving a wrist-pendulum system. Taken individually, each unit reproduced the natural tendency of the participants to freely oscillate a pendulum close to its resonance frequency. When cross-coupled through the vision of the pendulum of the other unit, the two units entrain each other and meet at a common frequency influenced by the mechanical properties of the two pendulums involved. The ability of the proposed model to address the other effects observed as a function of the different conditions of the pendulum and intended mode of coordination is discussed.

mu O-conotoxin MrVIB selectively blocks Na(v)1.8 sensory neuron specific sodium channels and chronic pain behavior without motor deficits

The tetroclotoxin-resistant voltage-gated sodium channel (VGSC) Na(v)1.8 is expressed predominantly by damage-sensing primary afferent nerves and is important for the development and maintenance of persistent pain states. Here we demonstrate that mu O-conotoxin MrVIB from Conus marmoreus displays substantial selectivity for Na(v)1.8 and inhibits pain behavior in models of persistent pain. In rat sensory neurons, submicromolar concentrations of MrVIB blocked tetroclotoxin-resistant current characteristic of Na(v)1.8 but not Na(v)1.9 or tetroclotoxin-sensitive VGSC currents. MrVIB blocked human Nav1.8 expressed in Xenopus oocytes with selectivity at least 10-fold greater than other VGSCs. In neuropathic and chronic inflammatory pain models, allodynia and hyperalgesia were both reduced by intrathecal infusion of MrVIB (0.03-3 nmol), whereas motor side effects occurred only at 30-fold higher doses. In contrast, the nonselective VGSC blocker lignocaine displayed no selectivity for allodynia and hyperalgesia versus motor side effects. The actions of MrVIB reveal that VGSC antagonists displaying selectivity toward Na(v)1.8 can alleviate chronic pain behavior with a greater therapeutic index than nonselective antagonists.

Motor unit number estimation - A Bayesian approach

All muscle contractions are dependent on the functioning of motor units. In diseases such as amyotrophic lateral sclerosis (ALS), progressive loss of motor units leads to gradual paralysis. A major difficulty in the search for a treatment for these diseases has been the lack of a reliable measure of disease progression. One possible measure would be an estimate of the number of surviving motor units. Despite over 30 years of motor unit number estimation (MUNE), all proposed methods have been met with practical and theoretical objections. Our aim is to develop a method of MUNE that overcomes these objections. We record the compound muscle action potential (CMAP) from a selected muscle in response to a graded electrical stimulation applied to the nerve. As the stimulus increases, the threshold of each motor unit is exceeded, and the size of the CMAP increases until a maximum response is obtained. However, the threshold potential required to excite an axon is not a precise value but fluctuates over a small range leading to probabilistic activation of motor units in response to a given stimulus. When the threshold ranges of motor units overlap, there may be alternation where the number of motor units that fire in response to the stimulus is variable. This means that increments in the value of the CMAP correspond to the firing of different combinations of motor units. At a fixed stimulus, variability in the CMAP, measured as variance, can be used to conduct MUNE using the "statistical" or the "Poisson" method. However, this method relies on the assumptions that the numbers of motor units that are firing probabilistically have the Poisson distribution and that all single motor unit action potentials (MUAP) have a fixed and identical size. These assumptions are not necessarily correct. We propose to develop a Bayesian statistical methodology to analyze electrophysiological data to provide an estimate of motor unit numbers. Our method of MUNE incorporates the variability of the threshold, the variability between and within single MUAPs, and baseline variability. Our model not only gives the most probable number of motor units but also provides information about both the population of units and individual units. We use Markov chain Monte Carlo to obtain information about the characteristics of individual motor units and about the population of motor units and the Bayesian information criterion for MUNE. We test our method of MUNE on three subjects. Our method provides a reproducible estimate for a patient with stable but severe ALS. In a serial study, we demonstrate a decline in the number of motor unit numbers with a patient with rapidly advancing disease. Finally, with our last patient, we show that our method has the capacity to estimate a larger number of motor units.