A new method for quantifying residue conservation and its applications to the protein folding nucleus

Kinetic versus thermodynamic control of mutational effects on protein homeostasis: A perspective from computational modeling and experiment

The effect of mutations in individual proteins on protein homeostasis, or "proteostasis," can in principle depend on the mutations' effects on the thermodynamics or kinetics of folding, or both. Here, we explore this issue using a computational model of in vivo protein folding that we call FoldEcoSlim. Our model predicts that kinetic versus thermodynamic control of mutational effects on proteostasis hinges on the relationship between how fast a protein's folding reaction reaches equilibrium and a critical time scale that characterizes the lifetime of a protein in its environment: for rapidly dividing bacteria, this time scale is that of cell division; for proteins that are produced in heterologous expression systems, this time scale is the amount of time before the protein is harvested; for proteins that are synthesized in and then exported from the eukaryotic endoplasmic reticulum, this time scale is that of protein secretion, and so forth. This prediction was validated experimentally by examining the expression yields of the wild type and several destabilized mutants of a model protein, the mouse ortholog of cellular retinoic acid-binding protein 1.

Codon substitution models based on residue similarity and their applications

Codon models are now widely used to draw evolutionary inferences from alignments of homologous sequence data. Incorporating physicochemical properties of amino acids into codon models, two novel codon substitution models describing the evolution of protein-coding DNA sequences are presented based on the similarity scores of amino acids. To describe substitutions between codons a continue-time Markov process is used. Transition/transversion rate bias and nonsynonymous codon usage bias are allowed in the models. In our implementation, the parameters are estimated by maximum-likelihood (ML) method as in previous studies. Furthermore, instantaneous mutations involving more than one nucleotide position of a codon are considered in the second model. Then the two suggested models are applied to five real data sets. The analytic results indicate that the new codon models considering physicochemical properties of amino acids can provide a better fit to the data comparing with existing codon models, and then produce more reliable estimates of certain biologically important measures than existing methods.

Robustness of the residue conservation score reflecting both frequencies and physicochemistries

Measuring residue conservation at aligned positions has many applications in biology. Recently, a new conservation score has been defined. Unlike the previous methods, the new approach considers both residue frequencies and physicochemistries. Specifically, it measures physicochemistries based on BLOSUM matrices disregarding the meaning of the entries in such matrices, which may involve the problem of log-log probability. In this paper we present a conservation measure that also reflects both frequencies and physicochemistries while considering the fact that the entries of BLOSUM matrices are already interpreted as log probability. When the supposed score is applied to 14 protein examples, the results show that these two conservation scores are equivalent aside from the different score ranges. The method is also used to score the functional sites of three protein families. Compared with the widely used entropy-based methods, the resulting scores are more robust and consistent in the sense that the functional sites are much more conserved because of functional constraints.