Driven fragmentation of granular gases

The dynamics of homogeneously heated granular gases which fragment due to particle collisions is analyzed. We introduce a kinetic model which accounts for correlations induced at the grain collisions and analyze both the kinetics and relevant distribution functions these systems develop. The work combines analytical and numerical studies based on direct simulation Monte Carlo calculations. A broad family of fragmentation probabilities is considered, and its implications for the system kinetics are discussed. We show that generically these driven materials evolve asymptotically into a dynamical scaling regime. If the fragmentation probability tends to a constant, the grain number diverges at a finite time, leading to a shattering singularity. If the fragmentation probability vanishes, then the number of grains grows monotonously as a power law. We consider different homogeneous thermostats and show that the kinetics of these systems depends weakly on both the grain inelasticity and driving. We observe that fragmentation plays a relevant role in the shape of the velocity distribution of the particles. When the fragmentation is driven by local stochastic events, the longvelocity tail is essentially exponential independently of the heating frequency and the breaking rule. However, for a Lowe-Andersen thermostat, numerical evidence strongly supports the conjecture that the scaled velocity distribution follows a generalized exponential behavior f (c)~exp (−cⁿ), with n ≈1.2, regarding less the fragmentation mechanisms

Generalized food-deceptive orchid species flower earlier and occur at lower altitudes than rewarding ones

Aims Food-deceptive pollination, in which plants do not offer any food reward to their pollinators, is common within the Orchidaceae. As food-deceptive orchids are poorer competitors for pollinator visitation than rewarding orchids, their occurrence in a given habitat may be more constrained than that of rewarding orchids. In particular, the success of deceptive orchids strongly relies on several biotic factors such as interactions with co-flowering rewarding species and pollinators, which may vary with altitude and over time. Our study compares generalized food-deceptive (i.e. excluding sexually deceptive) and rewarding orchids to test whether (i) deceptive orchids flower earlier compared to their rewarding counterparts and whether (ii) the relative occurrence of deceptive orchids decreases with increasing altitude. Methods To compare the flowering phenology of rewarding and deceptive orchids, we analysed data compiled from the literature at the species level over the occidental Palaearctic area. Since flowering phenology can be constrained by the latitudinal distribution of the species and by their phylogenetic relationships, we accounted for these factors in our analysis. To compare the altitudinal distribution of rewarding and deceptive orchids, we used field observations made over the entire Swiss territory and over two Swiss mountain ranges. Important Findings We found that deceptive orchid species start flowering earlier than rewarding orchids do, which is in accordance with the hypotheses of exploitation of naive pollinators and/or avoidance of competition with rewarding co-occurring species. Also, the relative frequency of deceptive orchids decreases with altitude, suggesting that deception may be less profitable at high compared to low altitude.

IL-6 activated integrated BATF/IRF4 functions in lymphocytes are T-bet-independent and reversed by subcutaneous immunotherapy.

IL-6 plays a central role in supporting pathological TH2 and TH17 cell development and inhibiting the protective T regulatory cells in allergic asthma. TH17 cells have been demonstrated to regulate allergic asthma in general and T-bet-deficiency-induced asthma in particular. Here we found an inverse correlation between T-bet and Il-6 mRNA expression in asthmatic children. Moreover, experimental subcutaneous immunotherapy (SIT) in T-bet((-/-)) mice inhibited IL-6, IL-21R and lung TH17 cells in a setting of asthma. Finally, local delivery of an anti-IL-6R antibody in T-bet((-/-)) mice resulted in the resolution of this allergic trait. Noteworthy, BATF, crucial for the immunoglobulin-class-switch and TH2,TH17 development, was found down-regulated in the lungs of T-bet((-/-)) mice after SIT and after treatment with anti-IL-6R antibody, indicating a critical role of IL-6 in controlling BATF/IRF4 integrated functions in TH2, TH17 cells and B cells also in a T-bet independent fashion in allergic asthma.

Mitotic functions for SNAP45, a subunit of the small nuclear RNA-activating protein complex SNAPc.

The small nuclear RNA-activating protein complex SNAP(c) is required for transcription of small nuclear RNA genes and binds to a proximal sequence element in their promoters. SNAP(c) contains five types of subunits stably associated with each other. Here we show that one of these polypeptides, SNAP45, also known as PTF delta, localizes to centrosomes during parts of mitosis, as well as to the spindle midzone during anaphase and the mid-body during telophase. Consistent with localization to these mitotic structures, both down- and up-regulation of SNAP45 lead to a G(2)/M arrest with cells displaying abnormal mitotic structures. In contrast, down-regulation of SNAP190, another SNAP(c) subunit, leads to an accumulation of cells with a G(0)/G(1) DNA content. These results are consistent with the proposal that SNAP45 plays two roles in the cell, one as a subunit of the transcription factor SNAP(c) and another as a factor required for proper mitotic progression.

New gene functions in megakaryopoiesis and platelet formation.

Platelets are the second most abundant cell type in blood and are essential for maintaining haemostasis. Their count and volume are tightly controlled within narrow physiological ranges, but there is only limited understanding of the molecular processes controlling both traits. Here we carried out a high-powered meta-analysis of genome-wide association studies (GWAS) in up to 66,867 individuals of European ancestry, followed by extensive biological and functional assessment. We identified 68 genomic loci reliably associated with platelet count and volume mapping to established and putative novel regulators of megakaryopoiesis and platelet formation. These genes show megakaryocyte-specific gene expression patterns and extensive network connectivity. Using gene silencing in Danio rerio and Drosophila melanogaster, we identified 11 of the genes as novel regulators of blood cell formation. Taken together, our findings advance understanding of novel gene functions controlling fate-determining events during megakaryopoiesis and platelet formation, providing a new example of successful translation of GWAS to function.

Negative Fukui functions: new insights based on electronegativity equalization

The occurrence of negative values for Fukui functions was studied through the electronegativity equalization method. Using algebraic relations between Fukui functions and different other conceptual DFT quantities on the one hand and the hardness matrix on the other hand, expressions were obtained for Fukui functions for several archetypical small molecules. Based on EEM calculations for large molecular sets, no negative Fukui functions were found

Critical thoughts on computing atom condensed Fukui functions

Different procedures to obtain atom condensed Fukui functions are described. It is shown how the resulting values may differ depending on the exact approach to atom condensed Fukui functions. The condensed Fukui function can be computed using either the fragment of molecular response approach or the response of molecular fragment approach. The two approaches are nonequivalent; only the latter approach corresponds in general with a population difference expression. The Mulliken approach does not depend on the approach taken but has some computational drawbacks. The different resulting expressions are tested for a wide set of molecules. In practice one must make seemingly arbitrary choices about how to compute condensed Fukui functions, which suggests questioning the role of these indicators in conceptual density-functional theory

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

Review article: Intestinal epithelia and barrier functions.

The mucosal epithelia of the digestive tract acts as a selective barrier, permeable to ions, small molecules and macromolecules. These epithelial cells aid the digestion of food and absorption of nutrients. They contribute to the protection against pathogens and undergo continuous cell renewal which facilitates the elimination of damaged cells. Both innate and adaptive defence mechanisms protect the gastrointestinal-mucosal surfaces against pathogens. Interaction of microorganisms with epithelial cells triggers a host response by activating specific transcription factors which control the expression of chemokines and cytokines. This host response is characterized by the recruitment of macrophages and neutrophils at the site of infection. Disruption of epithelial signalling pathways that recruit migratory immune cells results in a chronic inflammatory response. The adaptive defence mechanism relies on the collaboration of epithelial cells (resident sampling system) with antigen-presenting and lymphoid cells (migratory sampling system); in order to obtain samples of foreign antigen, these samples must be transported across the barriers without affecting the integrity of the barrier. These sampling systems are regulated by both environmental and host factors. Fates of the antigen may differ depending on the way in which they cross the epithelial barrier, i.e. via interaction with motile dendritic cells or epithelial M cells in the follicle-associated epithelium.

Splitting statistical potentials into meaningful scoring functions: testing the prediction of near-native structures from decoy conformations

Background: Recent advances on high-throughput technologies have produced a vast amount of protein sequences, while the number of high-resolution structures has seen a limited increase. This has impelled the production of many strategies to built protein structures from its sequence, generating a considerable amount of alternative models. The selection of the closest model to the native conformation has thus become crucial for structure prediction. Several methods have been developed to score protein models by energies, knowledge-based potentials and combination of both.Results: Here, we present and demonstrate a theory to split the knowledge-based potentials in scoring terms biologically meaningful and to combine them in new scores to predict near-native structures. Our strategy allows circumventing the problem of defining the reference state. In this approach we give the proof for a simple and linear application that can be further improved by optimizing the combination of Zscores. Using the simplest composite score () we obtained predictions similar to state-of-the-art methods. Besides, our approach has the advantage of identifying the most relevant terms involved in the stability of the protein structure. Finally, we also use the composite Zscores to assess the conformation of models and to detect local errors.Conclusion: We have introduced a method to split knowledge-based potentials and to solve the problem of defining a reference state. The new scores have detected near-native structures as accurately as state-of-art methods and have been successful to identify wrongly modeled regions of many near-native conformations.

Evolutionary origin and functions of retrogene introns.

Retroposed genes (retrogenes) originate via the reverse transcription of mature messenger RNAs from parental source genes and are therefore usually devoid of introns. Here, we characterize a particular set of mammalian retrogenes that acquired introns upon their emergence and thus represent rare cases of intron gain in mammals. We find that although a few retrogenes evolved introns in their coding or 3' untranslated regions (untranslated region, UTR), most introns originated together with untranslated exons in the 5' flanking regions of the retrogene insertion site. They emerged either de novo or through fusions with 5' UTR exons of host genes into which the retrogenes inserted. Generally, retrogenes with introns display high transcription levels and show broader spatial expression patterns than other retrogenes. Our experimental expression analyses of individual intron-containing retrogenes show that 5' UTR introns may indeed promote higher expression levels, at least in part through encoded regulatory elements. By contrast, 3' UTR introns may lead to downregulation of expression levels via nonsense-mediated decay mechanisms. Notably, the majority of retrogenes with introns in their 5' flanks depend on distant, sometimes bidirectional CpG dinucleotide-enriched promoters for their expression that may be recruited from other genes in the genomic vicinity. We thus propose a scenario where the acquisition of new 5' exon-intron structures was directly linked to the recruitment of distant promoters by these retrogenes, a process potentially facilitated by the presence of proto-splice sites in the genomic vicinity of retrogene insertion sites. Thus, the primary role and selective benefit of new 5' introns (and UTR exons) was probably initially to span the often substantial distances to potent CpG promoters driving retrogene transcription. Later in evolution, these introns then obtained additional regulatory roles in fine tuning retrogene expression levels. Our study provides novel insights regarding mechanisms underlying the origin of new introns, the evolutionary relevance of intron gain, and the origin of new gene promoters.

Generalized low-density codes with BCH constituents for full-diversity near-outage performance

A new graph-based construction of generalized low density codes (GLD-Tanner) with binary BCH constituents is described. The proposed family of GLD codes is optimal on block erasure channels and quasi-optimal on block fading channels. Optimality is considered in the outage probability sense. Aclassical GLD code for ergodic channels (e.g., the AWGN channel,the i.i.d. Rayleigh fading channel, and the i.i.d. binary erasure channel) is built by connecting bitnodes and subcode nodes via a unique random edge permutation. In the proposed construction of full-diversity GLD codes (referred to as root GLD), bitnodes are divided into 4 classes, subcodes are divided into 2 classes, and finally both sides of the Tanner graph are linked via 4 random edge permutations. The study focuses on non-ergodic channels with two states and can be easily extended to channels with 3 states or more.