Revealing evolutionary constraints on proteins through sequence analysis

Statistical analysis of alignments of large numbers of protein sequences has revealed "sectors" of collectively coevolving amino acids in several protein families. Here, we show that selection acting on any functional property of a protein, represented by an additive trait, can give rise to such a sector. As an illustration of a selected trait, we consider the elastic energy of an important conformational change within an elastic network model, and we show that selection acting on this energy leads to correlations among residues. For this concrete example and more generally, we demonstrate that the main signature of functional sectors lies in the small-eigenvalue modes of the covariance matrix of the selected sequences. However, secondary signatures of these functional sectors also exist in the extensively-studied large-eigenvalue modes. Our simple, general model leads us to propose a principled method to identify functional sectors, along with the magnitudes of mutational effects, from sequence data. We further demonstrate the robustness of these functional sectors to various forms of selection, and the robustness of our approach to the identification of multiple selected traits.Comment: 37 pages, 28 figure

Exposure of the Hidden Anti-Ferromagnetism in Paramagnetic CdSe:Mn Nanocrystals

We present theoretical and experimental investigations of the magnetism of paramagnetic semiconductor CdSe:Mn nanocrystals and propose an efficient approach to the exposure and analysis of the underlying anti-ferromagnetic interactions between magnetic ions therein. A key advance made here is the build-up of an analysis method with the exploitation of group theory technique that allows us to distinguish the anti-ferromagnetic interactions between aggregative Mn2+ ions from the overall pronounced paramagnetism of magnetic ion doped semiconductor nanocrystals. By using the method, we clearly reveal and identify the signatures of anti-ferromagnetism from the measured temperature dependent magnetisms, and furthermore determine the average number of Mn2+ ions and the fraction of aggregative ones in the measured CdSe:Mn nanocrystals.Comment: 26 pages, 5 figure

A Bayesian measurement error model for two-channel cell-based RNAi data with replicates

RNA interference (RNAi) is an endogenous cellular process in which small double-stranded RNAs lead to the destruction of mRNAs with complementary nucleoside sequence. With the production of RNAi libraries, large-scale RNAi screening in human cells can be conducted to identify unknown genes involved in a biological pathway. One challenge researchers face is how to deal with the multiple testing issue and the related false positive rate (FDR) and false negative rate (FNR). This paper proposes a Bayesian hierarchical measurement error model for the analysis of data from a two-channel RNAi high-throughput experiment with replicates, in which both the activity of a particular biological pathway and cell viability are monitored and the goal is to identify short hair-pin RNAs (shRNAs) that affect the pathway activity without affecting cell activity. Simulation studies demonstrate the flexibility and robustness of the Bayesian method and the benefits of having replicates in the experiment. This method is illustrated through analyzing the data from a RNAi high-throughput screening that searches for cellular factors affecting HCV replication without affecting cell viability; comparisons of the results from this HCV study and some of those reported in the literature are included.Comment: Published in at http://dx.doi.org/10.1214/11-AOAS496 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org