Thermodynamics of Neutral Protein Evolution

Naturally evolving proteins gradually accumulate mutations while continuing to fold to thermodynamically stable native structures. This process of neutral protein evolution is an important mode of genetic change, and forms the basis for the molecular clock. Here we present a mathematical theory that predicts the number of accumulated mutations, the index of dispersion, and the distribution of stabilities in an evolving protein population from knowledge of the stability effects (ddG values) for single mutations. Our theory quantitatively describes how neutral evolution leads to marginally stable proteins, and provides formulae for calculating how fluctuations in stability cause an overdispersion of the molecular clock. It also shows that the structural influences on the rate of sequence evolution that have been observed in earlier simulations can be calculated using only the single-mutation ddG values. We consider both the case when the product of the population size and mutation rate is small and the case when this product is large, and show that in the latter case proteins evolve excess mutational robustness that is manifested by extra stability and increases the rate of sequence evolution. Our basic method is to treat protein evolution as a Markov process constrained by a minimal requirement for stable folding, enabling an evolutionary description of the proteins solely in terms of the experimentally measureable ddG values. All of our theoretical predictions are confirmed by simulations with model lattice proteins. Our work provides a mathematical foundation for understanding how protein biophysics helps shape the process of evolution

Strong magnetoresistance induced by long-range disorder

We calculate the semiclassical magnetoresistivity $\rho_{xx}(B)$ of non-interacting fermions in two dimensions moving in a weak and smoothly varying random potential or random magnetic field. We demonstrate that in a broad range of magnetic fields the non-Markovian character of the transport leads to a strong positive magnetoresistance. The effect is especially pronounced in the case of a random magnetic field where $\rho_{xx}(B)$ becomes parametrically much larger than its B=0 value.Comment: REVTEX, 4 pages, 2 eps figure

Population genetics of translational robustness

Recent work has shown that expression level is the main predictor of a gene’s evolutionary rate, and that more highly expressed genes evolve slower. A possible explanation for this observation is selection for proteins which fold properly despite mistranslation, in short selection for translational robustness. Translational robustness leads to the somewhat paradoxical prediction that highly expressed genes are extremely tolerant to missense substitutions but nevertheless evolve very slowly. Here, we study a simple theoretical model of translational robustness that allows us to gain analytic insight into how this paradoxical behavior arises.Comment: 32 pages, 4 figures, Genetics in pres

Visual adaptation to convexity in macaque area V4

Aftereffects are perceptual illusions caused by visual adaptation to one or more stimulus attribute, such as orientation, motion, or shape. Neurophysiological studies seeking to understand the basis of visual adaptation have observed firing rate reduction and changes in tuning of stimulus-selective neurons following periods of prolonged visual stimulation. In the domain of shape, recent psychophysical work has shown that adaptation to a convex pattern induces a subsequently seen rectangle to appear slightly concave. In the present study, we investigate the possible contribution of V4 neurons of rhesus monkeys, which are thought to be involved in the coding of convexity, to shape-specific adaptation. Visually responsive neurons were monitored during the brief presentation of simple shapes varying in their convexity level. Each test presentation was preceded by either a blank period or several seconds of adaptation to a convex or concave stimulus, presented in two different sizes. Adaptation consistently shifted the tuning of neurons away from the convex or concave adapter, including shifting response to the neutral rectangle in the direction of the opposite convexity. This repulsive shift resembled the known perceptual distortion associated with adaptation to such stimuli. In addition, adaptation caused a nonspecific response decrease, as well as a specific decrease for repeated stimuli. The latter effects were observed whether or not the adapting and test stimuli matched closely in their size. Taken together, these results provide evidence for shape-specific adaptation of neurons in area V4, which may contribute to the perception of the convexity aftereffect

On quasilinear parabolic evolution equations in weighted Lp-spaces II

Our study of abstract quasi-linear parabolic problems in time-weighted L_p-spaces, begun in [17], is extended in this paper to include singular lower order terms, while keeping low initial regularity. The results are applied to reaction-diffusion problems, including Maxwell-Stefan diffusion, and to geometric evolution equations like the surface-diffusion flow or the Willmore flow. The method presented here will be applicable to other parabolic systems, including free boundary problems.Comment: 21 page

Zero-frequency anomaly in quasiclassical ac transport: Memory effects in a two-dimensional metal with a long-range random potential or random magnetic field

We study the low-frequency behavior of the {\it ac} conductivity $\sigma(\omega)$ of a two-dimensional fermion gas subject to a smooth random potential (RP) or random magnetic field (RMF). We find a non-analytic $\propto|\omega|$ correction to ${\rm Re} \sigma$, which corresponds to a $1/t^2$ long-time tail in the velocity correlation function. This contribution is induced by return processes neglected in Boltzmann transport theory. The prefactor of this $|\omega|$-term is positive and proportional to $(d/l)^2$ for RP, while it is of opposite sign and proportional to $d/l$ in the weak RMF case, where $l$ is the mean free path and $d$ the disorder correlation length. This non-analytic correction also exists in the strong RMF regime, when the transport is of a percolating nature. The analytical results are supported and complemented by numerical simulations.Comment: 12 pages, RevTeX, 7 figure