Generalized Rosenfeld scalings for tracer diffusivities in not-so-simple fluids: Mixtures and soft particles

Rosenfeld [Phys. Rev. A 15, 2545 (1977)] noticed that casting transport coefficients of simple monatomic, equilibrium fluids in specific dimensionless forms makes them approximately single-valued functions of excess entropy. This has predictive value because, while the transport coefficients of dense fluids are difficult to estimate from first principles, excess entropy can often be accurately predicted from liquid-state theory. Here, we use molecular simulations to investigate whether Rosenfeld's observation is a special case of a more general scaling law relating mobility of particles in mixtures to excess entropy. Specifically, we study tracer diffusivities, static structure, and thermodynamic properties of a variety of one- and two-component model fluid systems with either additive or non-additive interactions of the hard-sphere or Gaussian-core form. The results of the simulations demonstrate that the effects of mixture concentration and composition, particle-size asymmetry and additivity, and strength of the interparticle interactions in these fluids are consistent with an empirical scaling law relating the excess entropy to a new dimensionless (generalized Rosenfeld) form of tracer diffusivity, which we introduce here. The dimensionless form of the tracer diffusivity follows from knowledge of the intermolecular potential and the transport / thermodynamic behavior of fluids in the dilute limit. The generalized Rosenfeld scaling requires less information, and provides more accurate predictions, than either Enskog theory or scalings based on the pair-correlation contribution to the excess entropy. As we show, however, it also suffers from some limitations, especially for systems that exhibit significant decoupling of individual component tracer diffusivities.Comment: 15 pages, 10 figure

CodonTest: Modeling Amino Acid Substitution Preferences in Coding Sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes

Modeling HIV-1 Drug Resistance as Episodic Directional Selection

The evolution of substitutions conferring drug resistance to HIV-1 is both episodic, occurring when patients are on antiretroviral therapy, and strongly directional, with site-specific resistant residues increasing in frequency over time. While methods exist to detect episodic diversifying selection and continuous directional selection, no evolutionary model combining these two properties has been proposed. We present two models of episodic directional selection (MEDS and EDEPS) which allow the a priori specification of lineages expected to have undergone directional selection. The models infer the sites and target residues that were likely subject to directional selection, using either codon or protein sequences. Compared to its null model of episodic diversifying selection, MEDS provides a superior fit to most sites known to be involved in drug resistance, and neither one test for episodic diversifying selection nor another for constant directional selection are able to detect as many true positives as MEDS and EDEPS while maintaining acceptable levels of false positives. This suggests that episodic directional selection is a better description of the process driving the evolution of drug resistance

HIV-Specific Probabilistic Models of Protein Evolution

Comparative sequence analyses, including such fundamental bioinformatics techniques as similarity searching, sequence alignment and phylogenetic inference, have become a mainstay for researchers studying type 1 Human Immunodeficiency Virus (HIV-1) genome structure and evolution. Implicit in comparative analyses is an underlying model of evolution, and the chosen model can significantly affect the results. In general, evolutionary models describe the probabilities of replacing one amino acid character with another over a period of time. Most widely used evolutionary models for protein sequences have been derived from curated alignments of hundreds of proteins, usually based on mammalian genomes. It is unclear to what extent these empirical models are generalizable to a very different organism, such as HIV-1–the most extensively sequenced organism in existence. We developed a maximum likelihood model fitting procedure to a collection of HIV-1 alignments sampled from different viral genes, and inferred two empirical substitution models, suitable for describing between-and within-host evolution. Our procedure pools the information from multiple sequence alignments, and provided software implementation can be run efficiently in parallel on a computer cluster. We describe how the inferred substitution models can be used to generate scoring matrices suitable for alignment and similarity searches. Our models had a consistently superior fit relative to the best existing models and to parameter-rich data-driven models when benchmarked on independent HIV-1 alignments, demonstrating evolutionary biases in amino-acid substitution that are unique to HIV, and that are not captured by the existing models. The scoring matrices derived from the models showed a marked difference from common amino-acid scoring matrices. The use of an appropriate evolutionary model recovered a known viral transmission history, whereas a poorly chosen model introduced phylogenetic error. We argue that our model derivation procedure is immediately applicable to other organisms with extensive sequence data available, such as Hepatitis C and Influenza A viruses

Investigating the running abilities of Tyrannosaurus rex using stress-constrained multibody dynamic analysis

The running ability of Tyrannosaurus rex has been intensively studied due to its relevance to interpretations of feeding behaviour and the biomechanics of scaling in giant predatory dinosaurs. Different studies using differing methodologies have produced a very wide range of top speed estimates and there is therefore a need to develop techniques that can improve these predictions. Here we present a new approach that combines two separate biomechanical techniques (multibody dynamic analysis and skeletal stress analysis) to demonstrate that true running gaits would probably lead to unacceptably high skeletal loads in T. rex. Combining these two approaches reduces the high-level of uncertainty in previous predictions associated with unknown soft tissue parameters in dinosaurs, and demonstrates that the relatively long limb segments of T. rex—long argued to indicate competent running ability—would actually have mechanically limited this species to walking gaits. Being limited to walking speeds contradicts arguments of high-speed pursuit predation for the largest bipedal dinosaurs like T. rex, and demonstrates the power of multiphysics approaches for locomotor reconstructions of extinct animals

I=3/2 KπK \pi Scattering in the Nonrelativisitic Quark Potential Model

We study $I=3/2$ elastic $K\pi$ scattering to Born order using nonrelativistic quark wavefunctions in a constituent-exchange model. This channel is ideal for the study of nonresonant meson-meson scattering amplitudes since s-channel resonances do not contribute significantly. Standard quark model parameters yield good agreement with the measured S- and P-wave phase shifts and with PCAC calculations of the scattering length. The P-wave phase shift is especially interesting because it is nonzero solely due to $SU(3)_f$ symmetry breaking effects, and is found to be in good agreement with experiment given conventional values for the strange and nonstrange constituent quark masses.Comment: 12 pages + 2 postscript figures, Revtex, MIT-CTP-210

Adaptive Evolution of the Myo6 Gene in Old World Fruit Bats (Family: Pteropodidae)

PMCID: PMC3631194This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Do contaminants originating from state-of-the-art treated wastewater impact the ecological quality of surface waters?

Since the 1980s, advances in wastewater treatment technology have led to considerably improved surface water quality in the urban areas of many high income countries. However, trace concentrations of organic wastewater-associated contaminants may still pose a key environmental hazard impairing the ecological quality of surface waters. To identify key impact factors, we analyzed the effects of a wide range of anthropogenic and environmental variables on the aquatic macroinvertebrate community. We assessed ecological water quality at 26 sampling sites in four urban German lowland river systems with a 0–100% load of state-of-the-art biological activated sludge treated wastewater. The chemical analysis suite comprised 12 organic contaminants (five phosphor organic flame retardants, two musk fragrances, bisphenol A, nonylphenol, octylphenol, diethyltoluamide, terbutryn), 16 polycyclic aromatic hydrocarbons, and 12 heavy metals. Non-metric multidimensional scaling identified organic contaminants that are mainly wastewater-associated (i.e., phosphor organic flame retardants, musk fragrances, and diethyltoluamide) as a major impact variable on macroinvertebrate species composition. The structural degradation of streams was also identified as a significant factor. Multiple linear regression models revealed a significant impact of organic contaminants on invertebrate populations, in particular on Ephemeroptera, Plecoptera, and Trichoptera species. Spearman rank correlation analyses confirmed wastewater-associated organic contaminants as the most significant variable negatively impacting the biodiversity of sensitive macroinvertebrate species. In addition to increased aquatic pollution with organic contaminants, a greater wastewater fraction was accompanied by a slight decrease in oxygen concentration and an increase in salinity. This study highlights the importance of reducing the wastewater-associated impact on surface waters. For aquatic ecosystems in urban areas this would lead to: (i) improvement of the ecological integrity, (ii) reduction of biodiversity loss, and (iii) faster achievement of objectives of legislative requirements, e.g., the European Water Framework Directive

Detection of low prevalence somatic mutations in solid tumors with ultra-deep targeted sequencing

Ultra-deep targeted sequencing (UDT-Seq) can identify subclonal somatic mutations in tumor samples. Early assays' limited breadth and depth restrict their clinical utility. Here, we target 71 kb of mutational hotspots in 42 cancer genes. We present novel methods enhancing both laboratory workflow and mutation detection. We evaluate UDT-Seq true sensitivity and specificity (> 94% and > 99%, respectively) for low prevalence mutations in a mixing experiment and demonstrate its utility using six tumor samples. With an improved performance when run on the Illumina Miseq, the UDT-Seq assay is well suited for clinical applications to guide therapy and study clonal selection in heterogeneous samples

Phylodynamic Reconstruction Reveals Norovirus GII.4 Epidemic Expansions and their Molecular Determinants

Noroviruses are the most common cause of viral gastroenteritis. An increase in the number of globally reported norovirus outbreaks was seen the past decade, especially for outbreaks caused by successive genogroup II genotype 4 (GII.4) variants. Whether this observed increase was due to an upswing in the number of infections, or to a surveillance artifact caused by heightened awareness and concomitant improved reporting, remained unclear. Therefore, we set out to study the population structure and changes thereof of GII.4 strains detected through systematic outbreak surveillance since the early 1990s. We collected 1383 partial polymerase and 194 full capsid GII.4 sequences. A Bayesian MCMC coalescent analysis revealed an increase in the number of GII.4 infections during the last decade. The GII.4 strains included in our analyses evolved at a rate of 4.3–9.0×10−3 mutations per site per year, and share a most recent common ancestor in the early 1980s. Determinants of adaptation in the capsid protein were studied using different maximum likelihood approaches to identify sites subject to diversifying or directional selection and sites that co-evolved. While a number of the computationally determined adaptively evolving sites were on the surface of the capsid and possible subject to immune selection, we also detected sites that were subject to constrained or compensatory evolution due to secondary RNA structures, relevant in virus-replication. We highlight codons that may prove useful in identifying emerging novel variants, and, using these, indicate that the novel 2008 variant is more likely to cause a future epidemic than the 2007 variant. While norovirus infections are generally mild and self-limiting, more severe outcomes of infection frequently occur in elderly and immunocompromized people, and no treatment is available. The observed pattern of continually emerging novel variants of GII.4, causing elevated numbers of infections, is therefore a cause for concern