Allostery describes altered protein function at one site due to a perturbation at another site. One mechanism of allostery involves correlated motions, which can occur even in the absence of substantial conformational change. We present a novel method, "MutInf", to identify statistically significant correlated motions from equilibrium molecular dynamics simulations. Our approach analyzes both backbone and sidechain motions using internal coordinates to account for the gearlike twists that can take place even in the absence of the large conformational changes typical of traditional allosteric proteins. We quantify correlated motions using a mutual information metric, which we extend to incorporate data from multiple short simulations and to filter out correlations that are not statistically significant. Applying our approach to uncover mechanisms of cooperative small molecule binding in human interleukin-2, we identify clusters of correlated residues from 50 ns of molecular dynamics simulations. Interestingly, two of the clusters with the strongest correlations highlight known cooperative small-molecule binding sites and show substantial correlations between these sites. These cooperative binding sites on interleukin-2 are correlated not only through the hydrophobic core of the protein but also through a dynamic polar network of hydrogen bonding and electrostatic interactions. Since this approach identifies correlated conformations in an unbiased, statistically robust manner, it should be a useful tool for finding novel or "orphan" allosteric sites in proteins of biological and therapeutic importance.

The Ras-specific guanine nucleotide-exchange factors Son of sevenless (Sos) and Ras guanine nucleotide-releasing factor 1 (RasGRF1) transduce extracellular stimuli into Ras activation by catalyzing the exchange of Ras-bound GDP for GTP. A truncated form of RasGRF1 containing only the core catalytic Cdc25 domain is sufficient for stimulating Ras nucleotide exchange, whereas the isolated Cdc25 domain of Sos is inactive. At a site distal to the catalytic site, nucleotide-bound Ras binds to Sos, making contacts with the Cdc25 domain and with a Ras exchanger motif (Rem) domain. This allosteric Ras binding stimulates nucleotide exchange by Sos, but the mechanism by which this stimulation occurs has not been defined. We present a crystal structure of the Rem and Cdc25 domains of Sos determined at 2.0-Å resolution in the absence of Ras. Differences between this structure and that of Sos bound to two Ras molecules show that allosteric activation of Sos by Ras occurs through a rotation of the Rem domain that is coupled to a rotation of a helical hairpin at the Sos catalytic site. This motion relieves steric occlusion of the catalytic site, allowing substrate Ras binding and nucleotide exchange. A structure of the isolated RasGRF1 Cdc25 domain determined at 2.2-Å resolution, combined with computational analyses, suggests that the Cdc25 domain of RasGRF1 is able to maintain an active conformation in isolation because the helical hairpin has strengthened interactions with the Cdc25 domain core. These results indicate that RasGRF1 lacks the allosteric activation switch that is crucial for Sos activity.Cdc25 ͉ nucleotide-exchange factor ͉ crystal structure ͉ Ras exchanger motif (Rem) domain ͉ Ras guanine nucleotide-releasing factor 1

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.

No abstract