Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

Effect of floor type on the performance, physiological and behavioural responses of finishing beef steers

peer-reviewedBackground:The study objective was to investigate the effect of bare concrete slats (Control), two types of mats [(Easyfix mats (mat 1) and Irish Custom Extruder mats (mat 2)] fitted on top of concrete slats, and wood-chip to simulate deep bedding (wood-chip placed on top of a plastic membrane overlying the concrete slats) on performance, physiological and behavioral responses of finishing beef steers. One-hundred and forty-four finishing steers (503 kg; standard deviation 51.8 kg) were randomly assigned according to their breed (124 Continental cross and 20 Holstein–Friesian) and body weight to one of four treatments for 148 days. All steers were subjected to the same weighing, blood sampling (jugular venipuncture), dirt and hoof scoring pre study (day 0) and on days 23, 45, 65, 86, 107, 128 and 148 of the study. Cameras were fitted over each pen for 72 h recording over five periods and subsequent 10 min sampling scans were analysed. Results: Live weight gain and carcass characteristics were similar among treatments. The number of lesions on the hooves of the animals was greater (P < 0.05) on mats 1 and 2 and wood-chip treatments compared with the animals on the slats. Dirt scores were similar for the mat and slat treatments while the wood-chip treatment had greater dirt scores. Animals housed on either slats or wood-chip had similar lying times. The percent of animals lying was greater for animals housed on mat 1 and mat 2 compared with those housed on concrete slats and wood chips. Physiological variables showed no significant difference among treatments. Conclusions: In this exploratory study, the performance or welfare of steers was not adversely affected by slats, differing mat types or wood-chip as underfoot material

A randomized phase II trial of mitoxantrone, estramustine and vinorelbine or bcl-2 modulation with 13-cis retinoic acid, interferon and paclitaxel in patients with metastatic castrate-resistant prostate cancer: ECOG 3899

<p>Abstract</p> <p>Background</p> <p>To test the hypothesis that modulation of Bcl-2 with 13-cis retinoic acid (CRA)/interferon-alpha2b (IFN) with paclitaxel (TAX), or mitoxantrone, estramustine and vinorelbine (MEV) will have clinical activity in men with metastatic castrate-resistant prostate cancer (CRPC).</p> <p>Methods</p> <p>70 patients were treated with either MEV (Arm A) in a 3-week cycle or CRA/IFN/TAX with an 8-week cycle (Arm B). Patients were assessed for response, toxicity, quality of life (QOL), and the effect of treatment on Bcl-2 levels in peripheral blood mononuclear cells (PBMC).</p> <p>Results</p> <p>The PSA response rates were 50% and 23%, measurable disease response rates (CR+PR) 14% and 15%, and median overall survival 19.4 months and 13.9 months on Arm A and Arm B respectively. Transient grade 4 neutropenia occurred in 18 and 2 patients, and grade 3 to 4 thrombosis in 7 patients and 1 patient in Arm A and Arm B respectively. Patients on Arm B reported a clinically significant decline in QOL between baseline and week 9/10 (.71 s.d.), and a significantly lower level of QOL than Arm A (p = 0.01). As hypothesized, Bcl-2 levels decreased with CRA/IFN therapy only in Arm B (p = 0.03).</p> <p>Conclusions</p> <p>Treatment with MEV was well tolerated and demonstrated clinical activity in patients with CRPC. Given the adverse effect of CRA/IFN/TAX on QOL, the study of other novel agents that target Bcl-2 family proteins is warranted. The feasibility of measuring Bcl-2 protein in a cooperative group setting is hypothesis generating and supports further study as a marker for Bcl-2 targeted therapy.</p> <p>Trial Registration</p> <p><b>Clinical Trials Registration number</b>: CDR0000067865</p

Modeling the differentiation of A- and C-type baroreceptor firing patterns

The baroreceptor neurons serve as the primary transducers of blood pressure for the autonomic nervous system and are thus critical in enabling the body to respond effectively to changes in blood pressure. These neurons can be separated into two types (A and C) based on the myelination of their axons and their distinct firing patterns elicited in response to specific pressure stimuli. This study has developed a comprehensive model of the afferent baroreceptor discharge built on physiological knowledge of arterial wall mechanics, firing rate responses to controlled pressure stimuli, and ion channel dynamics within the baroreceptor neurons. With this model, we were able to predict firing rates observed in previously published experiments in both A- and C-type neurons. These results were obtained by adjusting model parameters determining the maximal ion-channel conductances. The observed variation in the model parameters are hypothesized to correspond to physiological differences between A- and C-type neurons. In agreement with published experimental observations, our simulations suggest that a twofold lower potassium conductance in C-type neurons is responsible for the observed sustained basal firing, whereas a tenfold higher mechanosensitive conductance is responsible for the greater firing rate observed in A-type neurons. A better understanding of the difference between the two neuron types can potentially be used to gain more insight into the underlying pathophysiology facilitating development of targeted interventions improving baroreflex function in diseased individuals, e.g. in patients with autonomic failure, a syndrome that is difficult to diagnose in terms of its pathophysiology.Comment: Keywords: Baroreflex model, mechanosensitivity, A- and C-type afferent baroreceptors, biophysical model, computational mode

Interaction among apoptosis-associated sequence variants and joint effects on aggressive prostate cancer

<p>Abstract</p> <p>Background</p> <p>Molecular and epidemiological evidence demonstrate that altered gene expression and single nucleotide polymorphisms in the apoptotic pathway are linked to many cancers. Yet, few studies emphasize the interaction of variant apoptotic genes and their joint modifying effects on prostate cancer (PCA) outcomes. An exhaustive assessment of all the possible two-, three- and four-way gene-gene interactions is computationally burdensome. This statistical conundrum stems from the prohibitive amount of data needed to account for multiple hypothesis testing.</p> <p>Methods</p> <p>To address this issue, we systematically prioritized and evaluated individual effects and complex interactions among 172 apoptotic SNPs in relation to PCA risk and aggressive disease (i.e., Gleason score ≥ 7 and tumor stages III/IV). Single and joint modifying effects on PCA outcomes among European-American men were analyzed using statistical epistasis networks coupled with multi-factor dimensionality reduction (SEN-guided MDR). The case-control study design included 1,175 incident PCA cases and 1,111 controls from the prostate, lung, colo-rectal, and ovarian (PLCO) cancer screening trial. Moreover, a subset analysis of PCA cases consisted of 688 aggressive and 488 non-aggressive PCA cases. SNP profiles were obtained using the NCI Cancer Genetic Markers of Susceptibility (CGEMS) data portal. Main effects were assessed using logistic regression (LR) models. Prior to modeling interactions, SEN was used to pre-process our genetic data. SEN used network science to reduce our analysis from > 36 million to < 13,000 SNP interactions. Interactions were visualized, evaluated, and validated using entropy-based MDR. All parametric and non-parametric models were adjusted for age, family history of PCA, and multiple hypothesis testing.</p> <p>Results</p> <p>Following LR modeling, eleven and thirteen sequence variants were associated with PCA risk and aggressive disease, respectively. However, none of these markers remained significant after we adjusted for multiple comparisons. Nevertheless, we detected a modest synergistic interaction between <it>AKT3 rs2125230-PRKCQ rs571715 </it>and disease aggressiveness using SEN-guided MDR (p = 0.011).</p> <p>Conclusions</p> <p>In summary, entropy-based SEN-guided MDR facilitated the logical prioritization and evaluation of apoptotic SNPs in relation to aggressive PCA. The suggestive interaction between <it>AKT3-PRKCQ </it>and aggressive PCA requires further validation using independent observational studies.</p

Habitat Composition and Connectivity Predicts Bat Presence and Activity at Foraging Sites in a Large UK Conurbation

Background: Urbanization is characterized by high levels of sealed land-cover, and small, geometrically complex, fragmented land-use patches. The extent and density of urbanized land-use is increasing, with implications for habitat quality, connectivity and city ecology. Little is known about densification thresholds for urban ecosystem function, and the response of mammals, nocturnal and cryptic taxa are poorly studied in this respect. Bats (Chiroptera) are sensitive to changing urban form at a species, guild and community level, so are ideal model organisms for analyses of this nature. Methodology/Principal Findings: We surveyed bats around urban ponds in the West Midlands conurbation, United Kingdom (UK). Sites were stratified between five urban land classes, representing a gradient of built land-cover at the 1 km 2 scale. Models for bat presence and activity were developed using land-cover and land-use data from multiple radii around each pond. Structural connectivity of tree networks was used as an indicator of the functional connectivity between habitats. All species were sensitive to measures of urban density. Some were also sensitive to landscape composition and structural connectivity at different spatial scales. These results represent new findings for an urban area. The activity of Pipistrellus pipistrellus (Schreber 1774) exhibited a non-linear relationship with the area of built land-cover, being much reduced beyond the threshold of,60 % built surface. The presence of tree networks appears to mitigate the negative effects of urbanization for this species

Genomewide association study for onset age in Parkinson disease

<p>Abstract</p> <p>Background</p> <p>Age at onset in Parkinson disease (PD) is a highly heritable quantitative trait for which a significant genetic influence is supported by multiple segregation analyses. Because genes associated with onset age may represent invaluable therapeutic targets to delay the disease, we sought to identify such genetic modifiers using a genomewide association study in familial PD. There have been previous genomewide association studies (GWAS) to identify genes influencing PD susceptibility, but this is the first to identify genes contributing to the variation in onset age.</p> <p>Methods</p> <p>Initial analyses were performed using genotypes generated with the Illumina HumanCNV370Duo array in a sample of 857 unrelated, familial PD cases. Subsequently, a meta-analysis of imputed SNPs was performed combining the familial PD data with that from a previous GWAS of 440 idiopathic PD cases. The SNPs from the meta-analysis with the lowest p-values and consistency in the direction of effect for onset age were then genotyped in a replication sample of 747 idiopathic PD cases from the Parkinson Institute Biobank of Milan, Italy.</p> <p>Results</p> <p>Meta-analysis across the three studies detected consistent association (p < 1 × 10^-5) with five SNPs, none of which reached genomewide significance. On chromosome 11, the SNP with the lowest p-value (rs10767971; p = 5.4 × 10^-7) lies between the genes <it>QSER1 </it>and <it>PRRG4</it>. Near the PARK3 linkage region on chromosome 2p13, association was observed with a SNP (rs7577851; p = 8.7 × 10^-6) which lies in an intron of the <it>AAK1 </it>gene. This gene is closely related to <it>GAK</it>, identified as a possible PD susceptibility gene in the GWAS of the familial PD cases.</p> <p>Conclusion</p> <p>Taken together, these results suggest an influence of genes involved in endocytosis and lysosomal sorting in PD pathogenesis.</p

Alanine Zipper-Like Coiled-Coil Domains Are Necessary for Homotypic Dimerization of Plant GAGA-Factors in the Nucleus and Nucleolus

GAGA-motif binding proteins control transcriptional activation or repression of homeotic genes. Interestingly, there are no sequence similarities between animal and plant proteins. Plant BBR/BPC-proteins can be classified into two distinct groups: Previous studies have elaborated on group I members only and so little is known about group II proteins. Here, we focused on the initial characterization of AtBPC6, a group II protein from Arabidopsis thaliana. Comparison of orthologous BBR/BPC sequences disclosed two conserved signatures besides the DNA binding domain. A first peptide signature is essential and sufficient to target AtBPC6-GFP to the nucleus and nucleolus. A second domain is predicted to form a zipper-like coiled-coil structure. This novel type of domain is similar to Leucine zippers, but contains invariant alanine residues with a heptad spacing of 7 amino acids. By yeast-2-hybrid and BiFC-assays we could show that this Alanine zipper domain is essential for homotypic dimerization of group II proteins in vivo. Interhelical salt bridges and charge-stabilized hydrogen bonds between acidic and basic residues of the two monomers are predicted to form an interaction domain, which does not follow the classical knobs-into-holes zipper model. FRET-FLIM analysis of GFP/RFP-hybrid fusion proteins validates the formation of parallel dimers in planta. Sequence comparison uncovered that this type of domain is not restricted to BBR/BPC proteins, but is found in all kingdoms

The first Dutch SDHB founder deletion in paraganglioma – pheochromocytoma patients

Contains fulltext : 81280.pdf (publisher's version ) (Open Access)BACKGROUND: Germline mutations of the tumor suppressor genes SDHB, SDHC and SDHD play a major role in hereditary paraganglioma and pheochromocytoma. These three genes encode subunits of succinate dehydrogenase (SDH), the mitochondrial tricarboxylic acid cycle enzyme and complex II component of the electron transport chain. The majority of variants of the SDH genes are missense and nonsense mutations. To date few large deletions of the SDH genes have been described. METHODS: We carried out gene deletion scanning using MLPA in 126 patients negative for point mutations in the SDH genes. We then proceeded to the molecular characterization of deletions, mapping breakpoints in each patient and used haplotype analysis to determine whether the deletions are due to a mutation hotspot or if a common haplotype indicated a single founder mutation. RESULTS: A novel deletion of exon 3 of the SDHB gene was identified in nine apparently unrelated Dutch patients. An identical 7905 bp deletion, c.201-4429_287-933del, was found in all patients, resulting in a frameshift and a predicted truncated protein, p.Cys68HisfsX21. Haplotype analysis demonstrated a common haplotype at the SDHB locus. Index patients presented with pheochromocytoma, extra-adrenal PGL and HN-PGL. A lack of family history was seen in seven of the nine cases. CONCLUSION: The identical exon 3 deletions and common haplotype in nine patients indicates that this mutation is the first Dutch SDHB founder mutation. The predominantly non-familial presentation of these patients strongly suggests reduced penetrance. In this small series HN-PGL occurs as frequently as pheochromocytoma and extra-adrenal PGL