Abstract:BackgroundProtein sites evolve at different rates due to functional and biophysical constraints. It is usually considered that the main structural determinant of a site’s rate of evolution is its Relative Solvent Accessibility (RSA). However, a recent comparative study has shown that the main structural determinant is the site’s Local Packing Density (LPD). LPD is related with dynamical flexibility, which has also been shown to correlate with sequence variability. Our purpose is to investigate the mechanism that connects a site’s LPD with its rate of evolution.ResultsWe consider two models: an empirical Flexibility Model and a mechanistic Stress Model. The Flexibility Model postulates a linear increase of site-specific rate of evolution with dynamical flexibility. The Stress Model, introduced here, models mutations as random perturbations of the protein’s potential energy landscape, for which we use simple Elastic Network Models (ENMs). To account for natural selection we assume a single active conformation and use basic statistical physics to derive a linear relationship between site-specific evolutionary rates and the local stress of the mutant’s active conformation.We compare both models on a large and diverse dataset of enzymes. In a protein-by-protein study we found that the Stress Model outperforms the Flexibility Model for most proteins. Pooling all proteins together we show that the Stress Model is strongly supported by the total weight of evidence. Moreover, it accounts for the observed nonlinear dependence of sequence variability on flexibility. Finally, when mutational stress is controlled for, there is very little remaining correlation between sequence variability and dynamical flexibility.ConclusionsWe developed a mechanistic Stress Model of evolution according to which the rate of evolution of a site is predicted to depend linearly on the local mutational stress of the active conformation. Such local stress is proportional to LPD, so that this model explains the relationship between LPD and evolutionary rate. Moreover, the model also accounts for the nonlinear dependence between evolutionary rate and dynamical flexibility.

On the variation of structural divergence among residues in enzyme evolution

Functional Sites Induce Long-Range Evolutionary Constraints in Enzymes

Conformational variation in enzyme catalysis: A structural study on catalytic residues

Evolution is not uniform along coding sequences

The divergence time of protein structures modelled by Markov matrices and its relation to the divergence of sequences

Conformational dynamics and enzyme evolution

Beyond Stability Constraints: A Biophysical Model of Enzyme Evolution with Selection on Stability and Activity

Protein structure prediction from sequence variation

Dissecting the roles of local packing density and longer-range effects in protein sequence evolution

Evolution of sparsity and modularity in a model of protein allostery

Beyond Thermodynamic Constraints: Evolutionary Sampling Generates Realistic Protein Sequence Variation

Evolutionary Conservation of Mechanical Strain Distributions in Functional Transitions of Protein Structures

A new formulation of protein evolutionary models that account for structural constraints

Substitution models of protein evolution with selection on enzymatic activity

The role of evolutionary selection in the dynamics of protein structure evolution

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

The Evolution of Protein Structures and Structural Ensembles Under Functional Constraint

Paradigms of convergent evolution in enzymes

A mechanistic stress model of protein evolution accounts for site-specific evolutionary rates and their relationship with packing density and flexibility

Inter-residue, inter-protein and inter-family coevolution: bridging the scales