James O. Wrabl

The ensemble nature of allostery.

Allostery is the process by which biological macromolecules (mostly proteins) transmit the effect of binding at one site to another, often distal, functional site, allowing for regulation of activity. Recent experimental observations demonstrating that allostery can be facilitated by dynamic and intrinsically disordered proteins have resulted in a new paradigm for understanding allosteric mechanisms, which focuses on the conformational ensemble and the statistical nature of the interactions responsible for the transmission of information. Analysis of allosteric ensembles reveals a rich spectrum of regulatory strategies, as well as a framework to unify the description of allosteric mechanisms from different systems.

Structural and Energetic Basis of Allostery

Allostery is a biological phenomenon of fundamental importance in regulation and signaling, and efforts to understand this process have led to the development of numerous models. In spite of individual successes in understanding the structural determinants of allostery in well-documented systems, much less success has been achieved in identifying a set of quantitative and transferable ground rules that provide an understanding of how allostery works. Are there organizing principles that allow us to relate structurally different proteins, or are the determinants of allostery unique to each system? Using an ensemble-based model, we show that allosteric phenomena can be formulated in terms of conformational free energies of the cooperative elements in a protein and the coupling interactions between them. Interestingly, the resulting allosteric ground rules provide a framework to reconcile observations that challenge purely structural models of site-to-site coupling, including (a) allostery in the absence of pathways of structural distortions, (b) allostery in the absence of any structural change, and (c) the ability of allosteric ligands to act as agonists under some circumstances and antagonists under others. The ensemble view of allostery that emerges provides insights into the energetic prerequisites of site-to-site coupling and thus into how allostery works.

CASP5 Assessment of Fold Recognition Target Predictions

We present an overview of the fifth round of Critical Assessment of Protein Structure Prediction (CASP5) fold recognition category. Prediction models were evaluated by using six different structural measures and four different alignment measures, and these scores were compared to those assigned manually over a diverse subset of target domains. Scores were combined to compare overall performance of participating groups and to estimate rank significance. The methods used by a few groups outperformed all other methods in terms of the evaluated criteria and could be considered state‐of‐the‐art in structure prediction. We discuss a few examples of difficult fold recognition targets to highlight the progress of ab initio‐type methods on difficult structure analogs and the difficulties of predicting multidomain targets and selecting prediction models. We also compared the results of manual groups to those of automatic servers evaluated in parallel by CAFASP, showing that the top performing automated server structure predictions approached those of the best manual predictors. Proteins 2003;53:395–409.

The role of protein conformational fluctuations in allostery, function, and evolution

It is now well-known that proteins exist at equilibrium as ensembles of conformational states rather than as unique static structures. Here we review from an ensemble perspective important biological effects of such spontaneous fluctuations on protein allostery, function, and evolution. However, rather than present a thorough literature review on each subject, we focus instead on connecting these phenomena through the ensemble-based experimental, theoretical, and computational investigations from our laboratory over the past decade. Special emphasis is given to insights that run counter to some of the prevailing ideas that have emerged over the past 40 years of structural biology research. For instance, when proteins are viewed as conformational ensembles rather than as single structures, the commonly held notion of an allosteric pathway as an obligate series of individual structural distortions loses its meaning. Instead, allostery can result from energetic linkage between distal sites as one Boltzmann distribution of states transitions to another. Additionally, the emerging principles from this ensemble view of proteins have proven surprisingly useful in describing the role of intrinsic disorder in inter-domain communication, functional adaptation mediated by mutational control of fluctuations, and evolutionary conservation of the energetics of protein stability.

A model of the changes in denatured state structure underlying m value effects in staphylococcal nuclease

Hydrogen exchange kinetics were measured on the native states of wild type staphylococcal nuclease and four mutants with values of mGuHCl (defined as dΔG/d[guanidine hydrochloride]) ranging from 0.8 to 1.4 of the wild type value. Residues within the five-strand β-barrel of wild type and E75A and D77A, two mutants with reduced values of m GuHCl, were significantly more protected from exchange than expected on the basis of global stability as measured by fluorescence. In contrast, mutants V23A and M26G with elevated values of mGuHCl approach a flat profile of more or less constant protection independent of position in the structure. Differences in exchange protection between the C-terminus and the β-barrel region correlate with mGuHCl, suggesting that a residual barrel-like structure becomes more highly populated in the denatured states of m- mutants and less populated in m+ mutants. Variations in the population of such a molten globule-like structure would account for the large changes in solvent accessible surface area of the denatured state thought to underlie m value effects.

The Torsin-family AAA+ Protein OOC-5 Contains a Critical Disulfide Adjacent to Sensor-II That Couples Redox State to Nucleotide Binding

A subgroup of the AAA+ proteins that reside in the endoplasmic reticulum and the nuclear envelope including human torsinA, a protein mutated in hereditary dystonia, is called the torsin family of AAA+ proteins. A multiple-sequence alignment of this family with Hsp100 proteins of known structure reveals a conserved cysteine in the C-terminus of torsin proteins within the Sensor-II motif. A structural model predicts this cysteine to be a part of an intramolecular disulfide bond, suggesting that it may function as a redox sensor to regulate ATPase activity. In vitro experiments with OOC-5, a torsinA homolog from Caenorhabditis elegans, demonstrate that redox changes that reduce this disulfide bond affect the binding of ATP and ADP and cause an attendant local conformational change detected by limited proteolysis. Transgenic worms expressing an ooc-5 gene with cysteine-to-serine mutations that disrupt the disulfide bond have a very low embryo hatch rate compared with wild-type controls, indicating these two cysteines are essential for OOC-5 function. We propose that the Sensor-II in torsin family proteins is a redox-regulated sensor. This regulatory mechanism may be central to the function of OOC-5 and human torsinA.

Gaps in structurally similar proteins: Towards improvement of multiple sequence alignment

An algorithm was developed to locally optimize gaps from the FSSP database. Over 2 million gaps were identified from all versus all FSSP structure comparisons, and datasets of non‐identical gaps and flanking regions comprising between 90,000 and 135,000 sequence fragments were extracted for statistical analysis. Relative to background frequencies, gaps were enriched in residue types with small side chains and high turn propensity (D, G, N, P, S), and were depleted in residue types with hydrophobic side chains (C, F, I, L, V, W, Y). In contrast, regions flanking a gap exhibited opposite trends in amino acid frequencies, i.e., enrichment in hydrophobic residues and a high degree of secondary structure. Log‐odds scores of residue type as a function of position in or around a gap were derived from the statistics. Three simple experiments demonstrated that these scores contained significant predictive information. First, regions where gaps were observed in single sequences taken from HOMSTRAD structure‐based multiple sequence alignments generally scored higher than regions where gaps were not observed. Second, given the correct pairwise‐aligned cores, the actual positions of gaps could be reproduced from sequence more accurately using the structurally‐derived statistics than by using random pairwise alignments. Finally, revision of the Clustal‐W residue‐specific gap opening parameters with this new information improved the agreement of Clustal‐W alignments with the structure‐based alignments. At least three applications for these results are envisioned: improvement of gap penalties in pairwise (or multiple) sequence alignment, prediction of regions of single sequences likely (or unlikely) to contain indels, and more accurate placement of gaps in automated pairwise structure alignment. Proteins 2003;53:000–000.

Thermodynamic environments in proteins: Fundamental determinants of fold specificity

To investigate the relationship between an amino acid sequence and its corresponding protein fold, a database of thermodynamic stability information was assembled as a function of residue type from 81 nonhomologous proteins. This information was obtained using the COREX algorithm, which computes an ensemble‐based description of the native state of proteins. Dissection of the COREX stability constant into its fundamental energetic components resulted in 12 thermodynamic environments describing the tertiary architecture of protein folds. Because of the observation that residue types partitioned unequally between these environments, it was hypothesized that thermodynamic environments contained energetic information that connected sequence to fold. To test the significance of this hypothesis, the thermodynamic stability information was incorporated into a three‐dimensional–to–one‐dimensional scoring matrix, and simple fold recognition experiments were performed in a manner such that information about the fold target was never included in the scoring. For 60 out of 81 fold targets, the correct sequence for the target scored in the top 5% of 3858 decoy sequences, with Z‐scores ranging from 1.76 to 12.23. Furthermore, a scoring matrix assembled from the residues of 40 nonhomologous all‐α proteins was used to thread sequences against 12 nonhomologous all‐β protein targets. In 10 of 12 cases, sequences known to adopt the native all‐β structure scored in the top 5% of 3858 decoy sequences, with Z‐scores ranging from 1.99 to 7.94. These results indicate that energetic information encoded by thermodynamic environments represents a fundamental property of proteins that underlies classifications based on secondary structure.

Grouping of amino acid types and extraction of amino acid properties from multiple sequence alignments using variance maximization

Understanding of amino acid type co‐occurrence in trusted multiple sequence alignments is a prerequisite for improved sequence alignment and remote homology detection algorithms. Two objective approaches were used to investigate co‐occurrence, both based on variance maximization of the weighted residue frequencies in columns taken from a large alignment database. The first approach discretely grouped amino acid types, and the second approach extracted orthogonal properties of amino acids using principal components analysis. The grouping results corresponded to amino acid physical properties such as side chain hydrophobicity, size, or backbone flexibility, and an optimal arrangement of approximately eight groups was observed. However, interpretation of the orthogonal properties was more complex. Although the principal components accounting for the largest variances exhibited modest correlations with hydrophobicity and conservation of glycine, in general principal components did not correspond to physical properties of amino acids. Although not intuitive, these amino acid mathematical properties were demonstrated to be robust and to improve local pairwise alignment accuracy, relative to 20 amino acid frequencies alone, for a simple test case. Proteins 2005.

Genetically tunable frustration controls allostery in an intrinsically disordered transcription factor

Intrinsically disordered proteins (IDPs) present a functional paradox because they lack stable tertiary structure, but nonetheless play a central role in signaling, utilizing a process known as allostery. Historically, allostery in structured proteins has been interpreted in terms of propagated structural changes that are induced by effector binding. Thus, it is not clear how IDPs, lacking such well-defined structures, can allosterically affect function. Here, we show a mechanism by which an IDP can allosterically control function by simultaneously tuning transcriptional activation and repression, using a novel strategy that relies on the principle of ‘energetic frustration’. We demonstrate that human glucocorticoid receptor tunes this signaling in vivo by producing translational isoforms differing only in the length of the disordered region, which modulates the degree of frustration. We expect this frustration-based model of allostery will prove to be generally important in explaining signaling in other IDPs.