Joint PDF modelling of turbulent flow and dispersion in an urban street canyon

The joint probability density function (PDF) of turbulent velocity and concentration of a passive scalar in an urban street canyon is computed using a newly developed particle-in-cell Monte Carlo method. Compared to moment closures, the PDF methodology provides the full one-point one-time PDF of the underlying fields containing all higher moments and correlations. The small-scale mixing of the scalar released from a concentrated source at the street level is modelled by the interaction by exchange with the conditional mean (IECM) model, with a micro-mixing time scale designed for geometrically complex settings. The boundary layer along no-slip walls (building sides and tops) is fully resolved using an elliptic relaxation technique, which captures the high anisotropy and inhomogeneity of the Reynolds stress tensor in these regions. A less computationally intensive technique based on wall functions to represent boundary layers and its effect on the solution are also explored. The calculated statistics are compared to experimental data and large-eddy simulation. The present work can be considered as the first example of computation of the full joint PDF of velocity and a transported passive scalar in an urban setting. The methodology proves successful in providing high level statistical information on the turbulence and pollutant concentration fields in complex urban scenarios.Comment: Accepted in Boundary-Layer Meteorology, Feb. 19, 200

FITBAR: a web tool for the robust prediction of prokaryotic regulons

<p>Abstract</p> <p>Background</p> <p>The binding of regulatory proteins to their specific DNA targets determines the accurate expression of the neighboring genes. The <it>in silico </it>prediction of new binding sites in completely sequenced genomes is a key aspect in the deeper understanding of gene regulatory networks. Several algorithms have been described to discriminate against false-positives in the prediction of new binding targets; however none of them has been implemented so far to assist the detection of binding sites at the genomic scale.</p> <p>Results</p> <p>FITBAR (Fast Investigation Tool for Bacterial and Archaeal Regulons) is a web service designed to identify new protein binding sites on fully sequenced prokaryotic genomes. This tool consists in a workbench where the significance of the predictions can be compared using different statistical methods, a feature not found in existing resources. The Local Markov Model and the Compound Importance Sampling algorithms have been implemented to compute the P-value of newly discovered binding sites. In addition, FITBAR provides two optimized genomic scanning algorithms using either log-odds or entropy-weighted position-specific scoring matrices. Other significant features include the production of a detailed genomic context map for each detected binding site and the export of the search results in spreadsheet and portable document formats. FITBAR discovery of a high affinity <it>Escherichia coli </it>NagC binding site was validated experimentally <it>in vitro </it>as well as <it>in vivo </it>and published.</p> <p>Conclusions</p> <p>FITBAR was developed in order to allow fast, accurate and statistically robust predictions of prokaryotic regulons. This feature constitutes the main advantage of this web tool over other matrix search programs and does not impair its performance. The web service is available at <url>http://archaea.u-psud.fr/fitbar</url>.</p

The p53HMM algorithm: using profile hidden markov models to detect p53-responsive genes

<p>Abstract</p> <p>Background</p> <p>A computational method (called p53HMM) is presented that utilizes Profile Hidden Markov Models (PHMMs) to estimate the relative binding affinities of putative p53 response elements (REs), both p53 single-sites and cluster-sites. These models incorporate a novel "Corresponded Baum-Welch" training algorithm that provides increased predictive power by exploiting the redundancy of information found in the repeated, palindromic p53-binding motif. The predictive accuracy of these new models are compared against other predictive models, including position specific score matrices (PSSMs, or weight matrices). We also present a new dynamic acceptance threshold, dependent upon a putative binding site's distance from the Transcription Start Site (TSS) and its estimated binding affinity. This new criteria for classifying putative p53-binding sites increases predictive accuracy by reducing the false positive rate.</p> <p>Results</p> <p>Training a Profile Hidden Markov Model with corresponding positions matching a combined-palindromic p53-binding motif creates the best p53-RE predictive model. The p53HMM algorithm is available on-line: <url>http://tools.csb.ias.edu</url></p> <p>Conclusion</p> <p>Using Profile Hidden Markov Models with training methods that exploit the redundant information of the homotetramer p53 binding site provides better predictive models than weight matrices (PSSMs). These methods may also boost performance when applied to other transcription factor binding sites.</p

Hand use predicts the structure of representations in sensorimotor cortex.

Fine finger movements are controlled by the population activity of neurons in the hand area of primary motor cortex. Experiments using microstimulation and single-neuron electrophysiology suggest that this area represents coordinated multi-joint, rather than single-finger movements. However, the principle by which these representations are organized remains unclear. We analyzed activity patterns during individuated finger movements using functional magnetic resonance imaging (fMRI). Although the spatial layout of finger-specific activity patterns was variable across participants, the relative similarity between any pair of activity patterns was well preserved. This invariant organization was better explained by the correlation structure of everyday hand movements than by correlated muscle activity. This also generalized to an experiment using complex multi-finger movements. Finally, the organizational structure correlated with patterns of involuntary co-contracted finger movements for high-force presses. Together, our results suggest that hand use shapes the relative arrangement of finger-specific activity patterns in sensory-motor cortex

Discriminative motif discovery in DNA and protein sequences using the DEME algorithm

<p>Abstract</p> <p>Background</p> <p>Motif discovery aims to detect short, highly conserved patterns in a collection of unaligned DNA or protein sequences. Discriminative motif finding algorithms aim to increase the sensitivity and selectivity of motif discovery by utilizing a second set of sequences, and searching only for patterns that can differentiate the two sets of sequences. Potential applications of discriminative motif discovery include discovering transcription factor binding site motifs in ChIP-chip data and finding protein motifs involved in thermal stability using sets of orthologous proteins from thermophilic and mesophilic organisms.</p> <p>Results</p> <p>We describe DEME, a discriminative motif discovery algorithm for use with protein and DNA sequences. Input to DEME is two sets of sequences; a "positive" set and a "negative" set. DEME represents motifs using a probabilistic model, and uses a novel combination of global and local search to find the motif that optimally discriminates between the two sets of sequences. DEME is unique among discriminative motif finders in that it uses an informative Bayesian prior on protein motif columns, allowing it to incorporate prior knowledge of residue characteristics. We also introduce four, synthetic, discriminative motif discovery problems that are designed for evaluating discriminative motif finders in various biologically motivated contexts. We test DEME using these synthetic problems and on two biological problems: finding yeast transcription factor binding motifs in ChIP-chip data, and finding motifs that discriminate between groups of thermophilic and mesophilic orthologous proteins.</p> <p>Conclusion</p> <p>Using artificial data, we show that DEME is more effective than a non-discriminative approach when there are "decoy" motifs or when a variant of the motif is present in the "negative" sequences. With real data, we show that DEME is as good, but not better than non-discriminative algorithms at discovering yeast transcription factor binding motifs. We also show that DEME can find highly informative thermal-stability protein motifs. Binaries for the stand-alone program DEME is free for academic use and is available at <url>http://bioinformatics.org.au/deme/</url></p

Analysis of Variance in Neuroreceptor Ligand Imaging Studies

Radioligand positron emission tomography (PET) with dual scan paradigms can provide valuable insight into changes in synaptic neurotransmitter concentration due to experimental manipulation. The residual t-test has been utilized to improve the sensitivity of the t-test in PET studies. However, no further development of statistical tests using residuals has been proposed so far to be applied in cases when there are more than two conditions. Here, we propose the residual f-test, a one-way analysis of variance (ANOVA), and examine its feasibility using simulated [11C]raclopride PET data. We also re-visit data from our previously published [11C]raclopride PET study, in which 10 individuals underwent three PET scans under different conditions. We found that the residual f-test is superior in terms of sensitivity than the conventional f-test while still controlling for type 1 error. The test will therefore allow us to reliably test hypotheses in the smaller sample sizes often used in explorative PET studies

Ranking insertion, deletion and nonsense mutations based on their effect on genetic information

<p>Abstract</p> <p>Background</p> <p>Genetic variations contribute to normal phenotypic differences as well as diseases, and new sequencing technologies are greatly increasing the capacity to identify these variations. Given the large number of variations now being discovered, computational methods to prioritize the functional importance of genetic variations are of growing interest. Thus far, the focus of computational tools has been mainly on the prediction of the effects of amino acid changing single nucleotide polymorphisms (SNPs) and little attention has been paid to indels or nonsense SNPs that result in premature stop codons.</p> <p>Results</p> <p>We propose computational methods to rank insertion-deletion mutations in the coding as well as non-coding regions and nonsense mutations. We rank these variations by measuring the extent of their effect on biological function, based on the assumption that evolutionary conservation reflects function. Using sequence data from budding yeast and human, we show that variations which that we predict to have larger effects segregate at significantly lower allele frequencies, and occur less frequently than expected by chance, indicating stronger purifying selection. Furthermore, we find that insertions, deletions and premature stop codons associated with disease in the human have significantly larger predicted effects than those not associated with disease. Interestingly, the large-effect mutations associated with disease show a similar distribution of predicted effects to that expected for completely random mutations.</p> <p>Conclusions</p> <p>This demonstrates that the evolutionary conservation context of the sequences that harbour insertions, deletions and nonsense mutations can be used to predict and rank the effects of the mutations.</p

Gene Expression Patterns of Dengue Virus-Infected Children from Nicaragua Reveal a Distinct Signature of Increased Metabolism

Dengue is a widespread viral disease for which over 3 billion people are at risk. There are no drug treatments or vaccines available for this disease. It is also difficult for physicians to predict which patients are at highest risk for the severe manifestations known as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS). We used genome-wide transcriptional profiling analysis to study peripheral blood responses to dengue among patients from Nicaragua. We found that patients with severe manifestations involving shock had very different transcriptional profiles from dengue patients with mild and moderate illness. We then compared our results with other microarray experiments on dengue patients available from public databases and confirmed that dengue is often associated with large changes to the metabolic processes within cells. This approach could identify prognostic markers for severe dengue as well as provide a better understanding of the pathophysiology associated with different grades of disease severity

Steroid receptor coactivator-1 modulates the function of Pomc neurons and energy homeostasis

Hypothalamic neurons expressing the anorectic peptide Pro-opiomelanocortin (Pomc) regulate food intake and body weight. Here, we show that Steroid Receptor Coactivator-1 (SRC-1) interacts with a target of leptin receptor activation, phosphorylated STAT3, to potentiate Pomc transcription. Deletion of SRC-1 in Pomc neurons in mice attenuates their depolarization by leptin, decreases Pomc expression and increases food intake leading to high-fat diet-induced obesity. In humans, fifteen rare heterozygous variants in SRC-1 found in severely obese individuals impair leptin-mediated Pomc reporter activity in cells, whilst four variants found in non-obese controls do not. In a knock-in mouse model of a loss of function human variant (SRC-1L1376P), leptin-induced depolarization of Pomc neurons and Pomc expression are significantly reduced, and food intake and body weight are increased. In summary, we demonstrate that SRC-1 modulates the function of hypothalamic Pomc neurons, and suggest that targeting SRC-1 may represent a useful therapeutic strategy for weight loss.Peer reviewe