Rejoinder to 'Performance of different synchronisation measures in real data: A case study on electroencephalograhic signals

We agree with Duckrow and Albano [Phys. Rev. E 67, 063901 (2003)] and Quian Quiroga et al. [Phys. Rev. E 67, 063902 (2003)] that mutual information (MI) is a useful measure of dependence for electroencephalogram (EEG) data, but we show that the improvement seen in the performance of MI on extracting dependence trends from EEG is more dependent on the type of MI estimator rather than any embedding technique used. In an independent study we conducted in search for an optimal MI estimator, and in particular for EEG applications, we examined the performance of a number of MI estimators on the data set used by Quian Quiroga et al. in their original study, where the performance of different dependence measures on real data was investigated [Phys. Rev. E 65, 041903 (2002)]. We show that for EEG applications the best performance among the investigated estimators is achieved by k-nearest neighbors, which supports the conjecture by Quian Quiroga et al. in Phys. Rev. E 67, 063902 (2003) that the nearest neighbor estimator is the most precise method for estimating MI.

Comment on "Performance of different synchronization measures in real data: A case study on electroencephalographic signals"

We agree with Duckrow and Albano [Phys. Rev. E 67, 063901 (2003)] and Quian Quiroga [Phys. Rev. E 67, 063902 (2003)] that mutual information (MI) is a useful measure of dependence for electroencephalogram (EEG) data, but we show that the improvement seen in the performance of MI on extracting dependence trends from EEG is more dependent on the type of MI estimator rather than any embedding technique used. In an independent study we conducted in search for an optimal MI estimator, and in particular for EEG applications, we examined the performance of a number of MI estimators on the data set used by Quian Quiroga in their original study, where the performance of different dependence measures on real data was investigated [Phys. Rev. E 65, 041903 (2002)]. We show that for EEG applications the best performance among the investigated estimators is achieved by k-nearest neighbors, which supports the conjecture by Quian Quiroga in Phys. Rev. E 67, 063902 (2003) that the nearest neighbor estimator is the most precise method for estimating MI.

Effect of changing the amount and type of fat and carbohydrate on insulin sensitivity and cardiovascular risk: the RISCK (Reading, Imperial, Surrey, Cambridge, and Kings) trial

Background: Insulin sensitivity (Si) is improved by weight loss and exercise, but the effects of the replacement of saturated fatty acids (SFAs) with monounsaturated fatty acids (MUFAs) or carbohydrates of high glycemic index (HGI) or low glycemic index (LGI) are uncertain. Objective: We conducted a dietary intervention trial to study these effects in participants at risk of developing metabolic syndrome. Design: We conducted a 5-center, parallel design, randomized controlled trial [RISCK (Reading, Imperial, Surrey, Cambridge, and Kings)]. The primary and secondary outcomes were changes in Si (measured by using an intravenous glucose tolerance test) and cardiovascular risk factors. Measurements were made after 4 wk of a high-SFA and HGI (HS/HGI) diet and after a 24-wk intervention with HS/HGI (reference), high-MUFA and HGI (HM/HGI), HM and LGI (HM/LGI), low-fat and HGI (LF/HGI), and LF and LGI (LF/LGI) diets. Results: We analyzed data for 548 of 720 participants who were randomly assigned to treatment. The median Si was 2.7 × 10−4 mL · μU−1 · min−1 (interquartile range: 2.0, 4.2 × 10−4 mL · μU−1 · min−1), and unadjusted mean percentage changes (95% CIs) after 24 wk treatment (P = 0.13) were as follows: for the HS/HGI group, −4% (−12.7%, 5.3%); for the HM/HGI group, 2.1% (−5.8%, 10.7%); for the HM/LGI group, −3.5% (−10.6%, 4.3%); for the LF/HGI group, −8.6% (−15.4%, −1.1%); and for the LF/LGI group, 9.9% (2.4%, 18.0%). Total cholesterol (TC), LDL cholesterol, and apolipoprotein B concentrations decreased with SFA reduction. Decreases in TC and LDL-cholesterol concentrations were greater with LGI. Fat reduction lowered HDL cholesterol and apolipoprotein A1 and B concentrations. Conclusions: This study did not support the hypothesis that isoenergetic replacement of SFAs with MUFAs or carbohydrates has a favorable effect on Si. Lowering GI enhanced reductions in TC and LDL-cholesterol concentrations in subjects, with tentative evidence of improvements in Si in the LF-treatment group. This trial was registered at clinicaltrials.gov as ISRCTN29111298.

Gene expression profiling of human hibernating myocardium: Increased expression of B-type natriuretic peptide and proenkephalin in hypocontractile vs normally-contracting regions of the heart

A greater understanding of the molecular basis of hibernating myocardium may assist in identifying those patients who would most benefit from revascularization. Paired heart biopsies were taken from hypocontractile and normally-contracting myocardium (identified by cardiovascular magnetic resonance) from 6 patients with chronic stable angina scheduled for bypass grafting. Gene expression profiles of hypocontractile and normally-contracting samples were compared using Affymetrix microarrays. The data for patients with confirmed hibernating myocardium were analysed separately and a different, though overlapping, set (up to 380) of genes was identified which may constitute a molecular fingerprint for hibernating myocardium. The expression of B-type natriuretic peptide (BNP) was increased in hypocontractile relative to normally-contracting myocardium. The expression of BNP correlated most closely with the expression of proenkephalin and follistatin 3, which may constitute additional heart failure markers. Our data illustrate differential gene expression in hypocontractile and/hibernating myocardium relative to normally-contracting myocardium within individual human hearts. Changes in expression of these genes, including increased relative expression of natriuretic and other factors, may constitute a molecular signature for hypocontractile and/or hibernating myocardium.

Assessing bee species richness in two Mediterranean communities: importance of habitat type and sampling techniques

The decline of bees has raised concerns regarding their conservation and the maintenance of ecosystem services they provide to bee-pollinated wild flowers and crops. Although the Mediterranean region is a hotspot for bee species richness, their status remains poorly studied. There is an urgent need for cost-effective, reliable, and unbiased sampling methods that give good bee species richness estimates. This study aims: (a) to assess bee species richness in two common Mediterranean habitat types: semi-natural scrub (phrygana) and managed olive groves; (b) to compare species richness in those systems to that of other biogeographic regions, and (c) to assess whether six different sampling methods (pan traps, variable and standardized transect walks, observation plots and trap nests), previously tested in other European biogeographic regions, are suitable in Mediterranean communities. Eight study sites, four per habitat type, were selected on the island of Lesvos, Greece. The species richness observed was high compared to other habitat types worldwide for which comparable data exist. Pan traps collected the highest proportion of the total bee species richness across all methods at the scale of a study site. Variable and standardized transect walks detected the highest total richness over all eight study sites. Trap nests and observation plots detected only a limited fraction of the bee species richness. To assess the total bee species richness in bee diversity hotspots, such as the studied habitats, we suggest a combination of transect walks conducted by trained bee collectors and pan trap sampling

Cost733cat: a database of weather and circulation type classifications

A new database of weather and circulation type catalogs is presented comprising 17 automated classification methods and five subjective classifications. It was compiled within COST Action 733 "Harmonisation and Applications of Weather Type Classifications for European regions" in order to evaluate different methods for weather and circulation type classification. This paper gives a technical description of the included methods using a new conceptual categorization for classification methods reflecting the strategy for the definition of types. Methods using predefined types include manual and threshold based classifications while methods producing types derived from the input data include those based on eigenvector techniques, leader algorithms and optimization algorithms. In order to allow direct comparisons between the methods, the circulation input data and the methods' configuration were harmonized for producing a subset of standard catalogs of the automated methods. The harmonization includes the data source, the climatic parameters used, the classification period as well as the spatial domain and the number of types. Frequency based characteristics of the resulting catalogs are presented, including variation of class sizes, persistence, seasonal and inter-annual variability as well as trends of the annual frequency time series. The methodological concept of the classifications is partly reflected by these properties of the resulting catalogs. It is shown that the types of subjective classifications compared to automated methods show higher persistence, inter-annual variation and long-term trends. Among the automated classifications optimization methods show a tendency for longer persistence and higher seasonal variation. However, it is also concluded that the distance metric used and the data preprocessing play at least an equally important role for the properties of the resulting classification compared to the algorithm used for type definition and assignment.

Evidence both L-type and non-L-type voltage-dependent calcium channels contribute to cerebral artery vasospasm following loss of NO in the rat

We recently found block of NO synthase in rat middle cerebral artery caused spasm, associated with depolarizing oscillations in membrane potential (Em) similar in form but faster in frequency (circa 1 Hz) to vasomotion. T-type voltage-gated Ca2+ channels contribute to cerebral myogenic tone and vasomotion, so we investigated the significance of T-type and other ion channels for membrane potential oscillations underlying arterial spasm. Smooth muscle cell membrane potential (Em) and tension were measured simultaneously in rat middle cerebral artery. NO synthase blockade caused temporally coupled depolarizing oscillations in cerebrovascular Em with associated vasoconstriction. Both events were accentuated by block of smooth muscle BKCa. Block of T-type channels or inhibition of Na+/K+-ATPase abolished the oscillations in Em and reduced vasoconstriction. Oscillations in Em were either attenuated or accentuated by reducing [Ca2+]o or block of KV, respectively. TRAM-34 attenuated oscillations in both Em and tone, apparently independent of effects against KCa3.1. Thus, rapid depolarizing oscillations in Em and tone observed after endothelial function has been disrupted reflect input from T-type calcium channels in addition to L-type channels, while other depolarizing currents appear to be unimportant. These data suggest that combined block of T and L-type channels may represent an effective approach to reverse cerebral vasospasm.

Data assimilation for state and parameter estimation: application to morphodynamic modelling

Data assimilation is predominantly used for state estimation; combining observational data with model predictions to produce an updated model state that most accurately approximates the true system state whilst keeping the model parameters fixed. This updated model state is then used to initiate the next model forecast. Even with perfect initial data, inaccurate model parameters will lead to the growth of prediction errors. To generate reliable forecasts we need good estimates of both the current system state and the model parameters. This paper presents research into data assimilation methods for morphodynamic model state and parameter estimation. First, we focus on state estimation and describe implementation of a three dimensional variational(3D-Var) data assimilation scheme in a simple 2D morphodynamic model of Morecambe Bay, UK. The assimilation of observations of bathymetry derived from SAR satellite imagery and a ship-borne survey is shown to significantly improve the predictive capability of the model over a 2 year run. Here, the model parameters are set by manual calibration; this is laborious and is found to produce different parameter values depending on the type and coverage of the validation dataset. The second part of this paper considers the problem of model parameter estimation in more detail. We explain how, by employing the technique of state augmentation, it is possible to use data assimilation to estimate uncertain model parameters concurrently with the model state. This approach removes inefficiencies associated with manual calibration and enables more effective use of observational data. We outline the development of a novel hybrid sequential 3D-Var data assimilation algorithm for joint state-parameter estimation and demonstrate its efficacy using an idealised 1D sediment transport model. The results of this study are extremely positive and suggest that there is great potential for the use of data assimilation-based state-parameter estimation in coastal morphodynamic modelling.

Western Chess: Zugzwang data

This data is derived from Eugene Nalimov's Depth-to-Mate Endgame Tables for Western Chess. While having the move is normally advantageous, there are positions where the side-to-move would have a better theoretical result if it were the other side to move. These are (Type A) 'zugzwang' positions where the 'obligation to act' is unwelcome. This data provides lists of all zugzwangs in sub-7-man chess, and summary data about those sets of zugzwangs including exemplar zugzwangs of maximum depth.

Effects of gelation temperature on the Mozzarella-type curd made from buffalo and cows’ milk. 1: Rheology & microstructure

The rheology and microstructure of Mozzarella-type curds made from buffalo and cows’ milk were measured at gelation temperatures of 28, 34 and 39 °C after chymosin addition. The maximum curd strength (G′) was obtained at a gelation temperature of 34 °C in both types of bovine milk. The viscoelasticity (tan δ) of both curds was increased with increasing gelation temperature. The rennet coagulation time was reduced with increase of gelation temperature in both types of milk. Frequency sweep data (0.1–10Hz was recorded 90 min after chymosin addition, and both milk samples showed characteristics of weak viscoelastic gel systems. When both milk samples were subjected to shear stress to break the curd system at constant shear rate, 95 min after chymosin addition, the maximum yield stress was obtained at the gelation temperatures of 34 °C and 28 °C in buffalo and cows’ curd respectively. The cryo-SEM and CLSM techniques were used to observe the microstructure of Mozzarella-type curd. The porosity was measured using image J software. The cryo-SEM and CLSM micrographs showed that minimum porosity was observed at the gelation temperature of 34 °C in both types of milk. Buffalo curd showed minimum porosity at similar gelation temperature when compared to cows’ curd. This may be due to higher protein concentration in buffalo milk.

The role of type 1 and curli fimbriae of Shiga toxin-producing Escherichia coli in adherence to abiotic surfaces

Biofilm formation on abiotic surfaces may provide a source of microbial contamination and may also enhance microbial environmental survival. The role of fimbrial expression by Shiga toxin-producing Escherichia coli (STEC) in biofilm formation is poorly understood. This study aimed to investigate the role of STEC type 1 and curli fimbriae in adhesion to and biofilm formation on abiotic surfaces. None of 13 O157:H7 isolates expressed either fimbrial type whereas 11 of 13 and 5 of 13 non-O157 STEC elaborated type 1 fimbriae and curli fimbriae, respectively. Mutants made by allelic exchange of a diarrhoeal non-O157 STEC isolate, O128:H2 (E41509), unable to elaborate type 1 and curli fimbriae were made for adherence and biofilm assays. Elaboration of type 1 fimbriae was necessary for the adhesion to abiotic surfaces whereas curliation was associated with both adherence and subsequent biofilm formation. STEC O157:H7 adhered to thermanox and glass but poorly to polystyrene. Additionally, STEC O157:H7 failed to form biofilms. These data indicate that certain STEC isolates are able to form biofilms and that the elaboration of curli fimbriae may enhance biofilm formation leading to possible long-term survival and a potential source of human infection.

Towards data warehousing and mining of protein unfolding simulation data

OBJECTIVES: The prediction of protein structure and the precise understanding of protein folding and unfolding processes remains one of the greatest challenges in structural biology and bioinformatics. Computer simulations based on molecular dynamics (MD) are at the forefront of the effort to gain a deeper understanding of these complex processes. Currently, these MD simulations are usually on the order of tens of nanoseconds, generate a large amount of conformational data and are computationally expensive. More and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. METHODS: To adequately organize, manage, and analyze the data generated by unfolding simulation studies, we designed a data warehouse system that is embedded in a grid environment to facilitate the seamless sharing of available computer resources and thus enable many groups to share complex molecular dynamics simulations on a more regular basis. RESULTS: To gain insight into the conformational fluctuations and stability of the monomeric forms of the amyloidogenic protein transthyretin (TTR), molecular dynamics unfolding simulations of the monomer of human TTR have been conducted. Trajectory data and meta-data of the wild-type (WT) protein and the highly amyloidogenic variant L55P-TTR represent the test case for the data warehouse. CONCLUSIONS: Web and grid services, especially pre-defined data mining services that can run on or 'near' the data repository of the data warehouse, are likely to play a pivotal role in the analysis of molecular dynamics unfolding data.

Distributed classification for pocket data mining

Distributed and collaborative data stream mining in a mobile computing environment is referred to as Pocket Data Mining PDM. Large amounts of available data streams to which smart phones can subscribe to or sense, coupled with the increasing computational power of handheld devices motivates the development of PDM as a decision making system. This emerging area of study has shown to be feasible in an earlier study using technological enablers of mobile software agents and stream mining techniques [1]. A typical PDM process would start by having mobile agents roam the network to discover relevant data streams and resources. Then other (mobile) agents encapsulating stream mining techniques visit the relevant nodes in the network in order to build evolving data mining models. Finally, a third type of mobile agents roam the network consulting the mining agents for a final collaborative decision, when required by one or more users. In this paper, we propose the use of distributed Hoeffding trees and Naive Bayes classifers in the PDM framework over vertically partitioned data streams. Mobile policing, health monitoring and stock market analysis are among the possible applications of PDM. An extensive experimental study is reported showing the effectiveness of the collaborative data mining with the two classifers.

Carbon monoxide inhibits L-type Ca2+ channels via redox modulation of key cysteine residues by mitochondrial reactive oxygen species

Conditions of stress, such as myocardial infarction, stimulate up-regulation of heme oxygenase (HO-1) to provide cardioprotection. Here, we show that CO, a product of heme catabolism by HO-1, directly inhibits native rat cardiomyocyte L-type Ca2+ currents and the recombinant alpha1C subunit of the human cardiac L-type Ca2+ channel. CO (applied via a recognized CO donor molecule or as the dissolved gas) caused reversible, voltage-independent channel inhibition, which was dependent on the presence of a spliced insert in the cytoplasmic C-terminal region of the channel. Sequential molecular dissection and point mutagenesis identified three key cysteine residues within the proximal 31 amino acids of the splice insert required for CO sensitivity. CO-mediated inhibition was independent of nitric oxide and protein kinase G but was prevented by antioxidants and the reducing agent, dithiothreitol. Inhibition of NADPH oxidase and xanthine oxidase did not affect the inhibitory actions of CO. Instead, inhibitors of complex III (but not complex I) of the mitochondrial electron transport chain and a mitochondrially targeted antioxidant (Mito Q) fully prevented the effects of CO. Our data indicate that the cardioprotective effects of HO-1 activity may be attributable to an inhibitory action of CO on cardiac L-type Ca2+ channels. Inhibition arises from the ability of CO to promote generation of reactive oxygen species from complex III of mitochondria. This in turn leads to redox modulation of any or all of three critical cysteine residues in the channel's cytoplasmic C-terminal tail, resulting in channel inhibition.