The Theory of Single Domain Allostery
There are times when you need to flog your own work. I recently published a paper  that provides a concrete model of single-domain allostery. It also provides a clean computational model of interior sidechain-sidechain interactions. I think it's my best work yet, and I will now flog the shit out of it for you.
1: Bosco K. Ho and David A. Agard. 2010. Conserved Tertiary Couplings Stabilize Elements in the PDZ fold, leading to Characteristic Patterns of Domain Conformational Flexibility. Protein Science 19:398-411.
What is allostery and why should you care about it?
From the greek word for other (allo), allostery refers to how events occurring at one site on a protein induces changes at another site, far, far away. Allostery was originally conceived to explain cooperative binding in oligomeric proteins made up of identical subunits. Binding sites are typically found at domain interfaces, where the binding of ligand at one site forces changes in the domain interface that lead to quaternary structure rearrangements. The rearrangements propagate to the empty binding sites on the other domain interfaces, resulting in better (cooperative) binding and hence, allostery. The concept of allostery now means any long-distance effect in a protein system.
As such, allostery is the foundation stone of information processing in the cell. Let's say we have two pathways, pathway A and pathway B, which are at first separate. We then select protein A from pathway A, and protein B from pathway B. Then we design an allosteric adaptor protein C where binding protein A to protein C will induce better binding for protein B at another site on protein C. The allostery in the adaptor protein C now intermingles pathway A to pathway B.
Allostery is thus the key to identifying pathway interactions in protein structures and finding a theory of allostery in protein systems is one of the major challenges in computational structural biology. The problem has split in all sorts of sub-problems such as oligomeric proteins and nucleotide-protein systems. But for me, the most challenging sub-problem is the allostery of single domains.
Singe domain allostery: allostery through dynamics
The surprising thing is that even a single domain can undergo allostery. What makes this difficult to understand? Well, if binding a ligand transmit changes to another site, then, in the absence of quaternary structure, these changes must tunnel through the body of the protein. Here is the rub: how do you model dynamics through the body of a protein?
One hypothesis for how single domain allostery might work is the Cooper-Dryden model , which argues that allostery occurs through flexibility in a protein. The idea is that large pieces of a single domain, including the binding site, are intrinsically flexible. Subsequent binding of a ligand will rigidify, not just the binding site, but other parts of the protein. These newly rigidified sites can serve as allosteric binding site.
2: Cooper A, Dryden DTF (1984) Allostery without conformational change—a plausible model. Eur Biophys J Biophys Lett 11:103–109.
However, the Cooper-Dryden model is only a generic thermodynamic argument. What is lacking is a way of calculating the intrinsic flexibility of a protein from the crystal structure, and showing that this explains allostery in a real protein system.
The same PDZ fold gives rise to different dynamics
The poster child of single-domain allostery is the PDZ domain. A shockingly small domain of ~60 amino acids, PDZ domains are ubiquitous scaffolding proteins. The PDZ domains have a well-defined binding groove that mainly binds C-terminii of other proteins or short peptides. However, there are at least 2 or three other binding sites in the PDZ fold. A quick survey of these interactions show that indeed, these little buggers act as allosteric adaptors.
Nevertheless, no one had actually catalogued the differences in allosteric changes across the PDZ domains. I found 5 PDZ domains with crystal structures that show distinctly different allosteric response. The differences can be seen in the mobility of the α-helices in the different PDZ domains upon ligand-binding, apo structures in magenta, ligand in blue, and ligand-bound structure in light blue:
As per the Cooper-Dryden model, we can attribute allostery in the PDZ domain to the intrinsic dynamics of the α-helices. Upon binding, the α-helices rigidifies creating new rigid surfaces, which are available for binding to other ligands, leading to allostery.
Here's the challenge: we have five PDZ domains (all with the same fold) with different dynamics in the α-helices. Is there a way to calculate the conformational flexibility from the apo crystal structures of the PDZ domains?
A computational model of sidechain-sidechain interactions
Well, of course there is a solution to the computational problem of single domain flexibility, otherwise I would be leading you down the garden path. Clearly, all 5 PDZ domains have the same fold, so major differences must lie in the sidechain-sidechain interactions. Now one might argue that long-time molecular-dynamic (MD) simulations of the 5 structures will tell you about the flexibility, but there are problems with this. First, it's expensive to run these simulations, but second, even if you could identify the flexible regions, there is currently no trajectory analysis method that allows you to clearly deduce sidechain-sidechain interactions.