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.

A Study on Structural Flexibility in Escherichia Coli Genome

Dynamic Flexibility in the Escherichia Coli Genome

Sequence-dependent flexibility in promoter sequences.

Flexibility of nucleic acids: From DNA to RNA

The Relationship among Gene Expression, Folding Free Energy and Codon Usage Bias in Escherichia Coli

The Role of CopG Mediated DNA Bending on the Regulation of the Σ 54 -Dependent Promoters in E. Coli

Change in flexibility of DNA with binding ligands

Flexibility of $Α$-Helices: Results of a Statistical Analysis of Database Protein Structures

A Study on Conformations of Dna Chain under the Elastic Rod Model

Length scale dependent elasticity in DNA from coarse-grained and all-atom models

Searching for sequence features that control DNA flexibility

Combined Study of Gene Expression and Chromosome Three-Dimensional Structure in Escherichia coli During Growth Process

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

A comparative kinetic and thermodynamic perspective of the σ-competition model in Escherichia coli

Influence of disordered packing pattern on elastic modulus of single-stranded DNA film on substrate

Correlated local bending of DNA double helix and its effect on the cyclization of short DNA fragments

Theoretical Analysis of Twist/Bend Ratio and Mechanical Moduli of Bacterial Flagellar Hook and Filament

Structural attributes of nucleotide sequences in promoter regions of supercoiling-sensitive genes: How to relate microarray expression data with genomic sequences

Entropic elasticity of DNA with a permanent kink

Entropic elasticity and dynamics of the bacterial chromosome: A simulation study

Optimal self-assembly pathways towards colloidal lattices with tunable flexibility