Epistasis shapes the fitness landscape of an allosteric specificity switch

Epistasis is a major determinant in the emergence of novel protein function. In allosteric proteins, direct interactions between inducer-binding mutations propagate through the allosteric network, manifesting as epistasis at the level of biological function. Elucidating this relationship between local interactions and their global effects is essential to understanding evolution of allosteric proteins. We integrate computational design, structural and biophysical analysis to characterize the emergence of novel inducer specificity in an allosteric transcription factor. Adaptive landscapes of different inducers of the designed mutant show that a few strong epistatic interactions constrain the number of viable sequence pathways, revealing ridges in the fitness landscape leading to new specificity. The structure of the designed mutant shows that a striking change in inducer orientation still retains allosteric function. Comparing biophysical and functional properties suggests a nonlinear relationship between inducer binding affinity and allostery. Our results highlight the functional and evolutionary complexity of allosteric proteins.

Prediction of allosteric sites and signalling: insights from benchmarking datasets

Allostery is a pervasive mechanism which regulates the activity of proteins in living systems through binding of a molecule at a distant site from the orthosteric site of the protein. The universality of allosteric regulation complemented by the benefits of highly specific, potentially non-toxic and protein activity modulating allosteric drugs makes uncovering allosteric sites on proteins invaluable for drug discovery. However, there are few computational methods to effectively predict them. Bond-to-bond propensity analysis, a recently developed method, has successfully predicted allosteric sites for a diverse group of proteins with only the knowledge of the orthosteric sites and the corresponding ligands in 19 of 20 cases. The method is based on an energy-weighted atomistic protein graph and allows for computationally highly efficient analysis in atomistic detail. We here extended the analysis onto 432 structures of 146 proteins from two existing benchmarking datasets for allosteric proteins: ASBench and CASBench. We further refined the metrics to account for the cumulative effect of residues with high propensities and the crucial residues in a given site with two additional measures. The allosteric site is recovered for 95/113 proteins (99/118 structures) from ASBench and 32/33 proteins (304/314 structures) from CASBench, with the only a priori knowledge being the orthosteric site residues. Knowing the orthosteric ligands of the protein, the allosteric site is identified for 32/33 proteins (308/314 structures) from CASBench.

On the Emergence of Orientational Order in Folded Proteins with Implications for Allostery

The beautiful structures of single- and multi-domain proteins are clearly ordered in some fashion but cannot be readily classified using group theory methods that are successfully used to describe periodic crystals. For this reason, protein structures are considered to be aperiodic, and may have evolved this way for functional purposes, especially in instances that require a combination of softness and rigidity within the same molecule. By analyzing the solved protein structures, we show that orientational symmetry is broken in the aperiodic arrangement of the secondary structure elements (SSEs), which we deduce by calculating the nematic order parameter, P2. We find that the folded structures are nematic droplets with a broad distribution of P2. We argue that a non-zero value of P2, leads to an arrangement of the SSEs that can resist external forces, which is a requirement for allosteric proteins. Such proteins, which resist mechanical forces in some regions while being flexible in others, transmit signals from one region of the protein to another (action at a distance) in response to binding of ligands (oxygen, ATP, or other small molecules).

On the Emergence of Orientational Order in Folded Proteins with Implications for Allostery

The beautiful structures of single and multi-domain proteins are clearly ordered in some fashion but cannot be readily classified using group theory methods that are successfully used to describe periodic crystals. For this reason, protein structures are considered to be aperiodic, and may have evolved this way for functional purposes, especially in instances that require a combination of softness and rigidity within the same molecule. By analyzing the solved protein structures, we show that orientational symmetry is broken in the aperiodic arrangement of the secondary structural elements (SSEs), which we deduce by calculating the nematic order parameter, $P_2$. We find that the folded structures are nematic droplets with a broad distribution of $P_2$. We argue that non-zero values of $P_2$, leads to an arrangement of the SSEs that can resist mechanical forces, which is a requirement for allosteric proteins. Such proteins, which resist mechanical forces in some regions while being flexible in others, transmit signals from one region of the protein to another (action at a distance) in response to binding of ligands (oxygen, ATP or other small molecules).

On the emergence of orientational order in folded proteins with implications for allostery

The beautiful structures of single and multi-domain proteins are clearly ordered in some fashion but cannot be readily classified using group theory methods that are successfully used to describe periodic crystals. For this reason, protein structures are considered to be aperiodic, and may have evolved this way for functional purposes, especially in instances that require a combination of softness and rigidity within the same molecule. By analyzing the solved protein structures, we show that orientational symmetry is broken in the aperiodic arrangement of the secondary structural elements (SSEs), which we deduce by calculating the nematic order parameter, P2. We find that the folded structures are nematic droplets with a broad distribution of P2. We argue that non-zero values of P2, leads to an arrangement of the SSEs that can resist external stresses forces, which is a requirement for allosteric proteins. Such proteins, which resist mechanical forces in some regions while being flexible in others, transmit signals from one region of the protein to another (action at a distance) in response to binding of ligands (oxygen, ATP or other small molecules).

Epistasis shapes the fitness landscape of an allosteric specificity switch

Epistasis is a major determinant in the emergence of novel protein function. In allosteric proteins, direct interactions between inducer-binding mutations propagate through the allosteric network, manifesting as epistasis at the level of biological function. Elucidating this relationship between local interactions and their global effects is essential to understanding evolution of allosteric proteins. We integrate computational design, structural and biophysical analysis to characterize the emergence of novel inducer specificity in an allosteric transcription factor. Adaptive landscapes of different inducers of the engineered mutant show that a few strong epistatic interactions constrain the number of viable sequence pathways, revealing ridges in the fitness landscape leading to new specificity. Crystallographic evidence shows a single mutation drives specificity by reshaping the binding pocket. Comparison of biophysical and functional landscapes emphasizes the nonlinear relationship between local inducer affinity and global function (allostery). Our results highlight the functional and evolutionary complexity of allosteric proteins.

Mechanics of allostery: contrasting the induced fit and population shift scenarios

In allosteric proteins, binding a ligand can affect function at a distant location, for example by changing the binding affinity of a substrate at the active site. The induced fit and population shift models, which differ by the assumed number of stable configurations, explain such cooperative binding from a thermodynamic viewpoint. Yet, understanding what mechanical principles constrain these models remains a challenge. Here we provide an empirical study on 34 proteins supporting the idea that allosteric conformational change generally occurs along a soft elastic mode presenting extended regions of high shear. We argue, based on a detailed analysis of how the energy profile along such a mode depends on binding, that in the induced fit scenario there is an optimal stiffness for cooperative binding, where N is the number of residues involved in the allosteric response. We find that the population shift scenario is more robust to mutation affecting stiffness, as binding becomes more and more cooperative with stiffness up to the same characteristic value , beyond which cooperativity saturates instead of decaying. We confirm numerically these findings in a non-linear mechanical model. Dynamical considerations suggest that a stiffness of order is favorable in that scenario as well, supporting that for proper function proteins must evolve a functional elastic mode that is softer as their size increases. In consistency with this view, we find a significant anticorrelation between the stiffness of the allosteric response and protein size in our data set.