S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Testing demographic models of effective population size.

Due to practical difficulties in obtaining direct genetic estimates of effective sizes, conservation biologists have to rely on so-called 'demographic models' which combine life-history and mating-system parameters with F-statistics in order to produce indirect estimates of effective sizes. However, for the same practical reasons that prevent direct genetic estimates, the accuracy of demographic models is difficult to evaluate. Here we use individual-based, genetically explicit computer simulations in order to investigate the accuracy of two such demographic models aimed at investigating the hierarchical structure of populations. We show that, by and large, these models provide good estimates under a wide range of mating systems and dispersal patterns. However, one of the models should be avoided whenever the focal species' breeding system approaches monogamy with no sex bias in dispersal or when a substructure within social groups is suspected because effective sizes may then be strongly overestimated. The timing during the life cycle at which F-statistics are evaluated is also of crucial importance and attention should be paid to it when designing field sampling since different demographic models assume different timings. Our study shows that individual-based, genetically explicit models provide a promising way of evaluating the accuracy of demographic models of effective size and delineate their field of applicability.

Mixed mating in androdioecious Mercurialis annua inferred using progeny arrays and diploid-acting microsatellite loci in a hexaploid background.

Background and Aims The frequency at which males can be maintained with hermaphrodites in androdioecious populations is predicted to depend on the selfing rate, because self-fertilization by hermaphrodites reduces prospective siring opportunities for males. In particular, high selfing rates by hermaphrodites are expected to exclude males from a population. Here, the first estimates are provided of the mating system from two wild hexaploid populations of the androdioecious European wind-pollinated plant M. annua with contrasting male frequencies.Methods Four diploid microsatellite loci were used to genotype 19-20 progeny arrays from two populations of M. annua, one with males and one without. Mating-system parameters were estimated using the program MLTR.Key Results Both populations had similar, intermediate outcrossing rates (t(m) = 0.64 and 0.52 for the population with and without males, respectively). The population without males showed a lower level of correlated paternity and biparental inbreeding and higher allelic richness and gene diversity than the population with males.Conclusions The results demonstrate the utility of new diploid microsatellite loci for mating system analysis in a hexaploid plant. It would appear that androdioecious M. annua has a mixed-mating system in the wild, an uncommon finding for wind-pollinated species. This study sets a foundation for future research to assess the relative importance of the sexual system, plant-density variation and stochastic processes for the regulation of male frequencies in M. annua over space and time.

Dustiness and deagglomeration testing: interlaboratory comparison of systems for nanoparticle powders

Different types of aerosolization and deagglomeration testing systems exist for studying the properties of nanomaterial powders and their aerosols. However, results are dependent on the specific methods used. In order to have well-characterized aerosols, we require a better understanding of how system parameters and testing conditions influence the properties of the aerosols generated. In the present study, four experimental setups delivering different aerosolization energies were used to test the resultant aerosols of two distinct nanomaterials (hydrophobic and hydrophilic TiO2). The reproducibility of results within each system was good. However, the number concentrations and size distributions of the aerosols created varied across the four systems; for number concentrations, e.g., from 10(3) to 10(6) #/cm(3). Moreover, distinct differences were also observed between the two materials with different surface coatings. The article discusses how system characteristics and other pertinent conditions modify the test results. We propose using air velocity as a suitable proxy for estimating energy input levels in aerosolization systems. The information derived from this work will be especially useful for establishing standard operating procedures for testing nanopowders, as well as for estimating their release rates under different energy input conditions, which is relevant for occupational exposure.

How useful is satellite positioning system (GPS) to track gait parameters? A review.

Over the last century, numerous techniques have been developed to analyze the movement of humans while walking and running. The combined use of kinematics and kinetics methods, mainly based on high speed video analysis and forceplate, have permitted a comprehensive description of locomotion process in terms of energetics and biomechanics. While the different phases of a single gait cycle are well understood, there is an increasing interest to know how the neuro-motor system controls gait form stride to stride. Indeed, it was observed that neurodegenerative diseases and aging could impact gait stability and gait parameters steadiness. From both clinical and fundamental research perspectives, there is therefore a need to develop techniques to accurately track gait parameters stride-by-stride over a long period with minimal constraints to patients. In this context, high accuracy satellite positioning can provide an alternative tool to monitor outdoor walking. Indeed, the high-end GPS receivers provide centimeter accuracy positioning with 5-20 Hz sampling rate: this allows the stride-by-stride assessment of a number of basic gait parameters--such as walking speed, step length and step frequency--that can be tracked over several thousand consecutive strides in free-living conditions. Furthermore, long-range correlations and fractal-like pattern was observed in those time series. As compared to other classical methods, GPS seems a promising technology in the field of gait variability analysis. However, relative high complexity and expensiveness--combined with a usability which requires further improvement--remain obstacles to the full development of the GPS technology in human applications.

Establishing an association between a peri-operative perfusion score system (PerfSCORE) and post-operative patient morbidity/mortality during CPB cardiac surgery.

BACKGROUND: To date, there is no quality assurance program that correlates patient outcome to perfusion service provided during cardiopulmonary bypass (CPB). A score was devised, incorporating objective parameters that would reflect the likelihood to influence patient outcome. The purpose was to create a new method for evaluating the quality of care the perfusionist provides during CPB procedures and to deduce whether it predicts patient morbidity and mortality. METHODS: We analysed 295 consecutive elective patients. We chose 10 parameters: fluid balance, blood transfused, Hct, ACT, PaO2, PaCO2, pH, BE, potassium and CPB time. Distribution analysis was performed using the Shapiro-Wilcoxon test. This made up the PerfSCORE and we tried to find a correlation to mortality rate, patient stay in the ICU and length of mechanical ventilation. Univariate analysis (UA) using linear regression was established for each parameter. Statistical significance was established when p < 0.05. Multivariate analysis (MA) was performed with the same parameters. RESULTS: The mean age was 63.8 +/- 12.6 years with 70% males. There were 180 CABG, 88 valves, and 27 combined CABG/valve procedures. The PerfSCORE of 6.6 +/- 2.4 (0-20), mortality of 2.7% (8/295), CPB time 100 +/- 41 min (19-313), ICU stay 52 +/- 62 hrs (7-564) and mechanical ventilation of 10.5 +/- 14.8 hrs (0-564) was calculated. CPB time, fluid balance, PaO2, PerfSCORE and blood transfused were significantly correlated to mortality (UA, p < 0.05). Also, CPB time, blood transfused and PaO2 were parameters predicting mortality (MA, p < 0.01). Only pH was significantly correlated for predicting ICU stay (UA). Ultrafiltration (UF) and CPB time were significantly correlated (UA, p < 0.01) while UF (p < 0.05) was the only parameter predicting mechanical ventilation duration (MA). CONCLUSIONS: CPB time, blood transfused and PaO2 are independent risk factors of mortality. Fluid balance, blood transfusion, PaO2, PerfSCORE and CPB time are independent parameters for predicting morbidity. PerfSCORE is a quality of perfusion measure that objectively quantifies perfusion performance.

Oral cancer treatments and adherence: medication event monitoring system assessment for capecitabine

Background: Oncological treatments are traditionally administered via intravenous injection by qualified personnel. Oral formulas which are developing rapidly are preferred by patients and facilitate administration however they may increase non-adherence. In this study 4 common oral chemotherapeutics are given to 50 patients, who are still in the process of inclusion, divided into 4 groups. The aim is to evaluate adherence and offer these patients interdisciplinary support with the joint help of doctors and pharmacists. We present here the results for capecitabine. Materials and Methods: The final goal is to evaluate adhesion in 50 patients split into 4 groups according to oral treatments (letrozole/exemestane, imatinib/sunitinib, capecitabine and temozolomide) using persistence and quality of execution as parameters. These parameters are evaluated using a medication event monitoring system (MEMS®) in addition to routine oncological visits and semi-structured interviews. Patients were monitored for the entire duration of treatment up to a maximum of 1 year. Patient satisfaction was assessed at the end of the monitoring period using a standardized questionary. Results: Capecitabine group included 2 women and 8 men with a median age of 55 years (range: 36−77 years) monitored for an average duration of 100 days (range: 5-210 days). Persistence was 98% and quality of execution 95%. 5 patients underwent cyclic treatment (2 out of 3 weeks) and 5 patients continuous treatment. Toxicities higher than grade 1 were grade 2−3 hand-foot syndrome in 1 patient and grade 3 acute coronary syndrome in 1 patient both without impact on adherence. Patients were satisfied with the interviews undergone during the study (57% useful, 28% very useful, 15% useless) and successfully integrated the MEMS® in their daily lives (57% very easily, 43% easily) according to the results obtained by questionary at the end of the monitoring period. Conclusion: Persistence and quality of execution observed in our Capecitabine group of patients were excellent and better than expected compared to previously published studies. The interdisciplinary approach allowed us to better identify and help patients with toxicities to maintain adherence. Overall patients were satisfied with the global interdisciplinary follow-up. With longer follow up better evaluation of our method and its impact will be possible. Interpretation of the results of patients in the other groups of this ongoing trial will provide us information for a more detailed analysis.

Effective size of two feral domestic cat populations (Felis catus L): effect of the mating system.

A variety of behavioural traits have substantial effects on the gene dynamics and genetic structure of local populations. The mating system is a plastic trait that varies with environmental conditions in the domestic cat (Felis catus) allowing an intraspecific comparison of the impact of this feature on genetic characteristics of the population. To assess the potential effect of the heterogenity of males' contribution to the next generation on variance effective size, we applied the ecological approach of Nunney & Elam (1994) based upon a demographic and behavioural study, and the genetic 'temporal methods' of Waples (1989) and Berthier et al. (2002) using microsatellite markers. The two cat populations studied were nearly closed, similar in size and survival parameters, but differed in their mating system. Immigration appeared extremely restricted in both cases due to environmental and social constraints. As expected, the ratio of effective size to census number (Ne/N) was higher in the promiscuous cat population (harmonic mean = 42%) than in the polygynous one (33%), when Ne was calculated from the ecological method. Only the genetic results based on Waples' estimator were consistent with the ecological results, but failed to evidence an effect of the mating system. Results based on the estimation of Berthier et al. (2002) were extremely variable, with Ne sometimes exceeding census size. Such low reliability in the genetic results should retain attention for conservation purposes.

The appropriate use of neurostimulation of the spinal cord and peripheral nervous system for the treatment of chronic pain and ischemic diseases: the neuromodulation appropriateness consensus committee.

INTRODUCTION: The Neuromodulation Appropriateness Consensus Committee (NACC) of the International Neuromodulation Society (INS) evaluated evidence regarding the safety and efficacy of neurostimulation to treat chronic pain, chronic critical limb ischemia, and refractory angina and recommended appropriate clinical applications. METHODS: The NACC used literature reviews, expert opinion, clinical experience, and individual research. Authors consulted the Practice Parameters for the Use of Spinal Cord Stimulation in the Treatment of Neuropathic Pain (2006), systematic reviews (1984 to 2013), and prospective and randomized controlled trials (2005 to 2013) identified through PubMed, EMBASE, and Google Scholar. RESULTS: Neurostimulation is relatively safe because of its minimally invasive and reversible characteristics. Comparison with medical management is difficult, as patients considered for neurostimulation have failed conservative management. Unlike alternative therapies, neurostimulation is not associated with medication-related side effects and has enduring effect. Device-related complications are not uncommon; however, the incidence is becoming less frequent as technology progresses and surgical skills improve. Randomized controlled studies support the efficacy of spinal cord stimulation in treating failed back surgery syndrome and complex regional pain syndrome. Similar studies of neurostimulation for peripheral neuropathic pain, postamputation pain, postherpetic neuralgia, and other causes of nerve injury are needed. International guidelines recommend spinal cord stimulation to treat refractory angina; other indications, such as congestive heart failure, are being investigated. CONCLUSIONS: Appropriate neurostimulation is safe and effective in some chronic pain conditions. Technological refinements and clinical evidence will continue to expand its use. The NACC seeks to facilitate the efficacy and safety of neurostimulation.

Measurement of the dynamics in ski jumping using a wearable inertial sensor-based system.

Abstract Dynamics is a central aspect of ski jumping, particularly during take-off and stable flight. Currently, measurement systems able to measure ski jumping dynamics (e.g. 3D cameras, force plates) are complex and only available in few research centres worldwide. This study proposes a method to determine dynamics using a wearable inertial sensor-based system which can be used routinely on any ski jumping hill. The system automatically calculates characteristic dynamic parameters during take-off (position and velocity of the centre of mass perpendicular to the table, force acting on the centre of mass perpendicular to the table and somersault angular velocity) and stable flight (total aerodynamic force). Furthermore, the acceleration of the ski perpendicular to the table was quantified to characterise the skis lift at take-off. The system was tested with two groups of 11 athletes with different jump distances. The force acting on the centre of mass, acceleration of the ski perpendicular to the table, somersault angular velocity and total aerodynamic force were different between groups and correlated with the jump distances. Furthermore, all dynamic parameters were within the range of prior studies based on stationary measurement systems, except for the centre of mass mean force which was slightly lower.

A 2D/3D image analysis system to track fluorescently labeled structures in rod-shaped cells: application to measure spindle pole asymmetry during mitosis.

BACKGROUND: The yeast Schizosaccharomyces pombe is frequently used as a model for studying the cell cycle. The cells are rod-shaped and divide by medial fission. The process of cell division, or cytokinesis, is controlled by a network of signaling proteins called the Septation Initiation Network (SIN); SIN proteins associate with the SPBs during nuclear division (mitosis). Some SIN proteins associate with both SPBs early in mitosis, and then display strongly asymmetric signal intensity at the SPBs in late mitosis, just before cytokinesis. This asymmetry is thought to be important for correct regulation of SIN signaling, and coordination of cytokinesis and mitosis. In order to study the dynamics of organelles or large protein complexes such as the spindle pole body (SPB), which have been labeled with a fluorescent protein tag in living cells, a number of the image analysis problems must be solved; the cell outline must be detected automatically, and the position and signal intensity associated with the structures of interest within the cell must be determined. RESULTS: We present a new 2D and 3D image analysis system that permits versatile and robust analysis of motile, fluorescently labeled structures in rod-shaped cells. We have designed an image analysis system that we have implemented as a user-friendly software package allowing the fast and robust image-analysis of large numbers of rod-shaped cells. We have developed new robust algorithms, which we combined with existing methodologies to facilitate fast and accurate analysis. Our software permits the detection and segmentation of rod-shaped cells in either static or dynamic (i.e. time lapse) multi-channel images. It enables tracking of two structures (for example SPBs) in two different image channels. For 2D or 3D static images, the locations of the structures are identified, and then intensity values are extracted together with several quantitative parameters, such as length, width, cell orientation, background fluorescence and the distance between the structures of interest. Furthermore, two kinds of kymographs of the tracked structures can be established, one representing the migration with respect to their relative position, the other representing their individual trajectories inside the cell. This software package, called "RodCellJ", allowed us to analyze a large number of S. pombe cells to understand the rules that govern SIN protein asymmetry. CONCLUSIONS: "RodCell" is freely available to the community as a package of several ImageJ plugins to simultaneously analyze the behavior of a large number of rod-shaped cells in an extensive manner. The integration of different image-processing techniques in a single package, as well as the development of novel algorithms does not only allow to speed up the analysis with respect to the usage of existing tools, but also accounts for higher accuracy. Its utility was demonstrated on both 2D and 3D static and dynamic images to study the septation initiation network of the yeast Schizosaccharomyces pombe. More generally, it can be used in any kind of biological context where fluorescent-protein labeled structures need to be analyzed in rod-shaped cells. AVAILABILITY: RodCellJ is freely available under http://bigwww.epfl.ch/algorithms.html, (after acceptance of the publication).

Paralysis of ventilated newborn babies does not influence resistance of the total respiratory system.

Paralysis with pancuronium bromide is used in newborn infants to facilitate ventilatory support during respiratory failure. Changes in lung mechanics have been attributed to paralysis. The aim of this study was to examine whether or not paralysis per se has an influence on the passive respiratory mechanics, resistance (Rrs) and compliance (Crs) of the respiratory system in newborn infants. In 30 infants with acute respiratory failure, Rrs was measured during paralysis with pancuronium bromide and after stopping pancuronium bromide (group A). Rrs was also measured in an additional 10 ventilated infants in a reversed fashion (group B): Rrs was measured first in nonparalysed infants and then they were paralysed, mainly for diagnostic procedures, and the Rrs measurement repeated. As Rrs is highly dependent on lung volume, several parameters, that depend directly on lung volume were recorded: inspiratory oxygen fraction (FI,O2), arterial oxygen tension/alveolar oxygen tension (a/A) ratio and volume above functional residual capacity (FRC). In group A, the Rrs was not different during (0.236+/-0.09 cmH2O x s x mL(-1)) and after (0.237+/-0.07 cmH2O x s x mL(-1)) paralysis. Also, in group B, Rrs did not change (0.207+/-0.046 versus 0.221+/-0.046 cm x s x mL(-1) without versus with pancuronium bromide). FI,O2, a/A ratio and volume above FRC remained constant during paralysis. These data demonstrate that paralysis does not influence the resistance of the total respiratory system in ventilated term and preterm infants when measured at comparable lung volumes.

An Automatic Approach for Learning and Tuning Gaussian Interval Type-2 Fuzzy Membership Functions Applied to Lung CAD Classification System

The potential of type-2 fuzzy sets for managing high levels of uncertainty in the subjective knowledge of experts or of numerical information has focused on control and pattern classification systems in recent years. One of the main challenges in designing a type-2 fuzzy logic system is how to estimate the parameters of type-2 fuzzy membership function (T2MF) and the Footprint of Uncertainty (FOU) from imperfect and noisy datasets. This paper presents an automatic approach for learning and tuning Gaussian interval type-2 membership functions (IT2MFs) with application to multi-dimensional pattern classification problems. T2MFs and their FOUs are tuned according to the uncertainties in the training dataset by a combination of genetic algorithm (GA) and crossvalidation techniques. In our GA-based approach, the structure of the chromosome has fewer genes than other GA methods and chromosome initialization is more precise. The proposed approach addresses the application of the interval type-2 fuzzy logic system (IT2FLS) for the problem of nodule classification in a lung Computer Aided Detection (CAD) system. The designed IT2FLS is compared with its type-1 fuzzy logic system (T1FLS) counterpart. The results demonstrate that the IT2FLS outperforms the T1FLS by more than 30% in terms of classification accuracy.

GPS analysis of human locomotion: further evidence for long-range correlations in stride-to-stride fluctuations of gait parameters.

During free walking, gait is automatically adjusted to provide optimal mechanical output and minimal energy expenditure; gait parameters, such as cadence, fluctuate from one stride to the next around average values. It was described that this fluctuation exhibited long-range correlations and fractal-like patterns. In addition, it was suggested that these long-range correlations disappeared if the participant followed the beep of metronome to regulate his or her pace. Until now, these fractal fluctuations were only observed for stride interval, because no technique existed to adequately analyze an extended time of free walking. The aim of the present study was to measure walking speed (WS), step frequency (SF) and step length (SL) with high accuracy (<1 cm) satellite positioning method (global positioning system or GPS) in order to detect long-range correlations in the stride-to-stride fluctuations. Eight participants walked 30 min under free and constrained (metronome) conditions. Under free walking conditions, DFA (detrended fluctuation analysis) and surrogate data tests showed that the fluctuation of WS, SL and SF exhibited a fractal pattern (i.e., scaling exponent alpha: 0.5 < alpha < 1) in a large majority of participants (7/8). Under constrained conditions (metronome), SF fluctuations became significantly anti-correlated (alpha < 0.5) in all participants. However, the scaling exponent of SL and WS was not modified. We conclude that, when the walking pace is controlled by an auditory signal, the feedback loop between the planned movement (at supraspinal level) and the sensory inputs induces a continual shifting of SF around the mean (persistent anti-correlation), but with no effect on the fluctuation dynamics of the other parameters (SL, WS).