GSHMC: an efficient method for molecular simulation

The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.

A Gaussian-mixture ensemble transform filter

We generalize the popular ensemble Kalman filter to an ensemble transform filter, in which the prior distribution can take the form of a Gaussian mixture or a Gaussian kernel density estimator. The design of the filter is based on a continuous formulation of the Bayesian filter analysis step. We call the new filter algorithm the ensemble Gaussian-mixture filter (EGMF). The EGMF is implemented for three simple test problems (Brownian dynamics in one dimension, Langevin dynamics in two dimensions and the three-dimensional Lorenz-63 model). It is demonstrated that the EGMF is capable of tracking systems with non-Gaussian uni- and multimodal ensemble distributions. Copyright © 2011 Royal Meteorological Society

A dynamical systems framework for intermittent data assimilation

We consider the problem of discrete time filtering (intermittent data assimilation) for differential equation models and discuss methods for its numerical approximation. The focus is on methods based on ensemble/particle techniques and on the ensemble Kalman filter technique in particular. We summarize as well as extend recent work on continuous ensemble Kalman filter formulations, which provide a concise dynamical systems formulation of the combined dynamics-assimilation problem. Possible extensions to fully nonlinear ensemble/particle based filters are also outlined using the framework of optimal transportation theory.

A nonparametric ensemble transform method for Bayesian inference

Many applications, such as intermittent data assimilation, lead to a recursive application of Bayesian inference within a Monte Carlo context. Popular data assimilation algorithms include sequential Monte Carlo methods and ensemble Kalman filters (EnKFs). These methods differ in the way Bayesian inference is implemented. Sequential Monte Carlo methods rely on importance sampling combined with a resampling step, while EnKFs utilize a linear transformation of Monte Carlo samples based on the classic Kalman filter. While EnKFs have proven to be quite robust even for small ensemble sizes, they are not consistent since their derivation relies on a linear regression ansatz. In this paper, we propose another transform method, which does not rely on any a priori assumptions on the underlying prior and posterior distributions. The new method is based on solving an optimal transportation problem for discrete random variables. © 2013, Society for Industrial and Applied Mathematics

Ensemble transform Kalman-Bucy filters

Two recent works have adapted the Kalman–Bucy filter into an ensemble setting. In the first formulation, the ensemble of perturbations is updated by the solution of an ordinary differential equation (ODE) in pseudo-time, while the mean is updated as in the standard Kalman filter. In the second formulation, the full ensemble is updated in the analysis step as the solution of single set of ODEs in pseudo-time. Neither requires matrix inversions except for the frequently diagonal observation error covariance. We analyse the behaviour of the ODEs involved in these formulations. We demonstrate that they stiffen for large magnitudes of the ratio of background error to observational error variance, and that using the integration scheme proposed in both formulations can lead to failure. A numerical integration scheme that is both stable and is not computationally expensive is proposed. We develop transform-based alternatives for these Bucy-type approaches so that the integrations are computed in ensemble space where the variables are weights (of dimension equal to the ensemble size) rather than model variables. Finally, the performance of our ensemble transform Kalman–Bucy implementations is evaluated using three models: the 3-variable Lorenz 1963 model, the 40-variable Lorenz 1996 model, and a medium complexity atmospheric general circulation model known as SPEEDY. The results from all three models are encouraging and warrant further exploration of these assimilation techniques.

The time-course of morphological processing in a second language

We report findings from psycholinguistic experiments investigating the detailed timing of processing morphologically complex words by proficient adult second (L2) language learners of English in comparison to adult native (L1) speakers of English. The first study employed the masked priming technique to investigate -ed forms with a group of advanced Arabic-speaking learners of English. The results replicate previously found L1/L2 differences in morphological priming, even though in the present experiment an extra temporal delay was offered after the presentation of the prime words. The second study examined the timing of constraints against inflected forms inside derived words in English using the eye-movement monitoring technique and an additional acceptability judgment task with highly advanced Dutch L2 learners of English in comparison to adult L1 English controls. Whilst offline the L2 learners performed native-like, the eye-movement data showed that their online processing was not affected by the morphological constraint against regular plurals inside derived words in the same way as in native speakers. Taken together, these findings indicate that L2 learners are not just slower than native speakers in processing morphologically complex words, but that the L2 comprehension system employs real-time grammatical analysis (in this case, morphological information) less than the L1 system.

The short-chain fatty acid acetate reduces appetite via a central homeostatic mechanism

Increased intake of dietary carbohydrate that is fermented in the colon by the microbiota has been reported to decrease body weight, although the mechanism remains unclear. Here we use in vivo11C-acetate and PET-CT scanning to show that colonic acetate crosses the blood–brain barrier and is taken up by the brain. Intraperitoneal acetate results in appetite suppression and hypothalamic neuronal activation patterning. We also show that acetate administration is associated with activation of acetyl-CoA carboxylase and changes in the expression profiles of regulatory neuropeptides that favour appetite suppression. Furthermore, we demonstrate through 13C high-resolution magic-angle-spinning that 13C acetate from fermentation of 13C-labelled carbohydrate in the colon increases hypothalamic 13C acetate above baseline levels. Hypothalamic 13C acetate regionally increases the 13C labelling of the glutamate–glutamine and GABA neuroglial cycles, with hypothalamic 13C lactate reaching higher levels than the ‘remaining brain’. These observations suggest that acetate has a direct role in central appetite regulation.

Density of insect-pollinated grassland plants decreases with increasing surrounding land-use intensity

Pollinator declines have raised concerns about the persistence of plant species that depend on insect pollination, in particular by bees, for their reproduction. The impact of pollinator declines remains unknown for species-rich plant communities found in temperate seminatural grasslands. We investigated effects of land-use intensity in the surrounding landscape on the distribution of plant traits related to insect pollination in 239 European seminatural grasslands. Increasing arable land use in the surrounding landscape consistently reduced the density of plants depending on bee and insect pollination. Similarly, the relative abundance of bee-pollination-dependent plants increased with higher proportions of non-arable agricultural land (e.g. permanent grassland). This was paralleled by an overall increase in bee abundance and diversity. By isolating the impact of the surrounding landscape from effects of local habitat quality, we show for the first time that grassland plants dependent on insect pollination are particularly susceptible to increasing land-use intensity in the landscape.

Rare variants in LRRK1 and Parkinson's disease

Approximately 20 % of individuals with Parkinson's disease (PD) report a positive family history. Yet, a large portion of causal and disease-modifying variants is still unknown. We used exome sequencing in two affected individuals from a family with late-onset PD to identify 15 potentially causal variants. Segregation analysis and frequency assessment in 862 PD cases and 1,014 ethnically matched controls highlighted variants in EEF1D and LRRK1 as the best candidates. Mutation screening of the coding regions of these genes in 862 cases and 1,014 controls revealed several novel non-synonymous variants in both genes in cases and controls. An in silico multi-model bioinformatics analysis was used to prioritize identified variants in LRRK1 for functional follow- up. However, protein expression, subcellular localization, and cell viability were not affected by the identified variants. Although it has yet to be proven conclusively that variants in LRRK1 are indeed causative of PD, our data strengthen a possible role for LRRK1 in addition to LRRK2 in the genetic underpinnings of PD but, at the same time, highlight the difficulties encountered in the study of rare variants identified by next-generation sequencing in diseases with autosomal dominant or complex patterns of inheritance.

Negative online consumer reviews: effects of different responses

Purpose When consumers buy online, they are often confronted with consumer reviews. A negative consumer review on an online shopping website may keep consumers from buying the product. Therefore, negative online consumer reviews are a serious problem for brands. This paper aims to investigate the effects of different response options to a negative consumer review. Design/methodology/approach In an online experiment of 446 participants different response options towards a negative consumer review on an online shopping website are examined. The experimental data is analysed with simple linear regression models using product purchase intentions as the outcome variable. Findings The results indicate that a positive customer review counteracts a negative consumer review more effectively than a positive brand response, whereas brand strength moderates this relationship. Including a reference to an independent, trusted source in a brand or a customer response is only a limited strategy for increasing the effectiveness of a response. Research limitations/implications Additional research in other product categories and with other subjects than students is suggested to validate the findings. In future research, multiple degrees of the phrasing’s strength of the reference could be used. Practical implications Assuming high quality products, brands should encourage their customers to write reviews. Strong brands can also reassure consumers by responding whereas weak brands cannot. Originality/value This research contributes to the online consumer reviews literature with new insights about the role of brand strength and referencing to an independent, trusted source.

Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5

The new Max-Planck-Institute Earth System Model (MPI-ESM) is used in the Coupled Model Intercomparison Project phase 5 (CMIP5) in a series of climate change experiments for either idealized CO2-only forcing or forcings based on observations and the Representative Concentration Pathway (RCP) scenarios. The paper gives an overview of the model configurations, experiments related forcings, and initialization procedures and presents results for the simulated changes in climate and carbon cycle. It is found that the climate feedback depends on the global warming and possibly the forcing history. The global warming from climatological 1850 conditions to 2080–2100 ranges from 1.5°C under the RCP2.6 scenario to 4.4°C under the RCP8.5 scenario. Over this range, the patterns of temperature and precipitation change are nearly independent of the global warming. The model shows a tendency to reduce the ocean heat uptake efficiency toward a warmer climate, and hence acceleration in warming in the later years. The precipitation sensitivity can be as high as 2.5% K−1 if the CO2 concentration is constant, or as small as 1.6% K−1, if the CO2 concentration is increasing. The oceanic uptake of anthropogenic carbon increases over time in all scenarios, being smallest in the experiment forced by RCP2.6 and largest in that for RCP8.5. The land also serves as a net carbon sink in all scenarios, predominantly in boreal regions. The strong tropical carbon sources found in the RCP2.6 and RCP8.5 experiments are almost absent in the RCP4.5 experiment, which can be explained by reforestation in the RCP4.5 scenario.

Morphological characterization of periodontium-derived human stem cells

The aim of this study has been to characterize adult human somatic periodontium-derived stem cells (PDSCS) isolated from human periodontium and to follow their differentiation after cell culture. PDSCS were isolated from human periodontal tissue and cultured as spheres in serum-free medium. After 10 days the primary spheres were dissociated and the secondary spheres sub-cultured for another 1-2 weeks. Cells from different time points were analyzed, and immunohistochemical and electron microscopic investigations carried out. Histological analysis showed differentiation of spheres deriving from the PDSCS with central production of extracellular matrix beginning 3 days after sub-culturing. Isolated PDSCS developed pseudopodia which contained actin. Tubulin was found in the central portion of the cells. Pseudopodia between different cells anastomosed, indicating intercellular transport. Immunostaining for osteopontin demonstrated a positive reaction in primary spheres and within extracellular matrix vesicles after sub-culturing. In cell culture under serum-free conditions human PDSCS form spheres which are capable of producing extracellular matrix. Further investigations have do be carried out to investigate the capability of these cells to differentiate into osteogenic progenitor cells.

Subdiffraction fluorescence imaging of biomolecular structure and distributions with quantum dots

We introduce semiconductor quantum dot-based fluorescence imaging with approximately 2-fold increased optical resolution in three dimensions as a method that allows both studying cellular structures and spatial organization of biomolecules in membranes and subcellular organelles. Target biomolecules are labelled with quantum dots via immunocytochemistry. The resolution enhancement is achieved by three-photon absorption of quantum dots and subsequent fluorescence emission from a higher-order excitonic state. Different from conventional multiphoton microscopy, this approach can be realized on any confocal microscope without the need for pulsed excitation light. We demonstrate quantum dot triexciton imaging (QDTI) of the microtubule network of U373 cells, 3D imaging of TNF receptor 2 on the plasma membrane of HeLa cells, and multicolor 3D imaging of mitochondrial cytochrome c oxidase and actin in COS-7 cells.

TNF-α influences the lateral dynamics of TNF receptor I in living cells

In mammalian cells, inflammation is mainly mediated by the binding of tumor necrosis factor alpha to tumor necrosis factor receptor 1. In this study, we investigated lateral dynamics of TNF-R1 before and after ligand binding using high-density single-particle tracking in combination with photoactivated localization microscopy. Our single-molecule data indicates the presence of tumor necrosis factor receptor 1 with different mobilities in the plasma membrane, suggesting different molecular organizations. Cholesterol depletion led to a decrease of slow receptor species and a strong increase in the average diffusion coefficient. Moreover, as a consequence of tumor necrosis factor-alpha treatment, the mean diffusion coefficient moderately increased while its distribution narrowed. Based on our observation, we propose a refined mechanism on the structural arrangement and activation of tumor necrosis factor receptor 1 in the plasma membrane.

Identification of novel cholesteatoma-related gene expression signatures using full-genome microarrays

This study demonstrates that the expression profile of cholesteatoma is similar to a metastatic tumour and chronically inflamed tissue. Based on the investigated profiles we present novel protein-protein interaction and signal transduction networks, which include cholesteatoma-regulated transcripts and may be of great value for drug targeting and therapy development.