Ligand modulation of allosteric networks in an ancestral steroid receptor

Understanding the evolution of binding specificity, a heavily studied area of research, is key for determining how protein sequence changes alter function. Ligand-activation in the steroid receptor subfamily of transcription factors operates via a common allosteric mechanism which permits extant receptors to respond specifically to their cognate hormones. Here, we combine atomistic simulations with graph theory-based modeling of the inter-residue interactions within protein complexes to gain insight into how allostery drove selectivity in an ancestral receptor. An inactive ligand complex displays weakened allosteric communication, as quantified by suboptimal paths linking two functional surfaces. When function-switching mutations are incorporated, responses in allosteric networks are consistent with ligand activation status. Further analysis reveals residues that modulate features distinguishing active and inactive complexes, identifying a key, conserved residue that is crucial for activation in steroid receptors. We have identified a computational method using dynamic network analysis to probe the allosteric mechanisms driving the evolution of ligand specificity in hormone receptors, determining how residue substitutions altered allosteric networks to permit gain or loss of ligand response. These results may have general utility in elucidating how modern steroid receptors are activated by endogenous and xenobiotic molecules.