The fitness contribution of an allele at one genetic site may depend on the states of other sites, a phenomenon known as epistasis. Epistasis can profoundly influence the process of evolution in populations under selection, and  shape the course of protein evolution across divergent species. Whereas epistasis among adaptive substitutions has been the subject of extensive study, relatively little is known about epistasis under purifying selection. Here we use mechanistic models of thermodynamic stability in a ligand-binding protein to explore computationally the structure of epistatic interactions among substitutions that fix in protein sequences under purifying selection. We find that the selection coefficients of mutations that are nearly neutral when they fix are highly conditional on the presence of preceding mutations. In addition, substitutions which are initially neutral become increasingly entrenched over time due to antagonistic epistasis with subsequent substitutions. Our evolutionary model includes insertions and deletions, as well as point mutations, which allows us to quantify epistasis between these classes of mutations, and also to study the evolution of protein length. We find that protein length remains largely constant over time, because indels are more deleterious than point mutations. Our results imply that, even under purifying selection, protein sequence evolution is highly contingent on history and  it cannot be predicted by the phenotypic effects of mutations introduced into the wildtype sequencealone .
The fitness contribution of an allele at one genetic site may depend on alleles at other sites, a phenomenon known as epistasis. Epistasis can profoundly influence the process of evolution in populations under selection, and can shape the course of protein evolution across divergent species. Whereas epistasis between adaptive substitutions has been the subject of extensive study, relatively little is known about epistasis under purifying selection. Here we use mechanistic models of thermodynamic stability in a ligand-binding protein to explore the structure of epistatic interactions between substitutions that fix in protein sequences under purifying selection. We find that the selection coefficients of mutations that are nearly-neutral when they fix are highly contingent on the presence of preceding mutations. Conversely, mutations that are nearly-neutral when they fix are subsequently entrenched due to epistasis with later substitutions. Our evolutionary model includes insertions and deletions, as well as point mutations, and so it allows us to quantify epistasis within each of these classes of mutations, and also to study the evolution of protein length. We find that protein length remains largely constant over time, because indels are more deleterious than point mutations. Our results imply that, even under purifying selection, protein sequence evolution is highly contingent on history and so it cannot be predicted by the phenotypic effects of mutations assayed in the wild-type sequence .
R
the states of
alleles at
A
can
R
among
between
D
computationally
R
among
between
R
nearly neutral
nearly-neutral
R
conditional
contingent
R
In addition, substitutions which are initially neutral become increasingly entrenched over time due to antagonistic epistasis with subsequent
Conversely, mutations that are nearly-neutral when they fix are subsequently entrenched due to epistasis with later
R
which
and so it
R
between
within each of
A
so
R
introduced into the wildtype sequencealone
assayed in the wild-type sequence