Partially unfolded forms and non-two-state folding of a beta-sandwich: FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

The case for defined protein folding pathways

We consider the differences between the many-pathway protein folding model derived from theoretical energy landscape considerations and the defined-pathway model derived from experiment. A basic tenet of the energy landscape model is that proteins fold through many heterogeneous pathways by way of amino acid-level dynamics biased toward selecting native-like interactions. The many pathways imagined in the model are not observed in the structure-formation stage of folding by experiments that would have found them, but they have now been detected and characterized for one protein in the initial prenucleation stage. Analysis presented here shows that these many microscopic trajectories are not distinct in any functionally significant way, and they have neither the structural information nor the biased energetics needed to select native vs. nonnative interactions during folding. The opposed defined-pathway model stems from experimental results that show that proteins are assemblies of small cooperative units called foldons and that a number of proteins fold in a reproducible pathway one foldon unit at a time. Thus, the same foldon interactions that encode the native structure of any given protein also naturally encode its particular foldon-based folding pathway, and they collectively sum to produce the energy bias toward native interactions that is necessary for efficient folding. Available information suggests that quantized native structure and stepwise folding coevolved in ancient repeat proteins and were retained as a functional pair due to their utility for solving the difficult protein folding problem.

Mechanism of Protein Denaturation: Partial Unfolding of the P22 Coat Protein I-Domain by Urea Binding

The I-domain is an insertion domain of the bacteriophage P22 coat protein that drives rapid folding and accounts for over half of the stability of the full-length protein. We sought to determine the role of hydrogen bonds (H-bonds) in the unfolding of the I-domain by examining (3)JNC' couplings transmitted through H-bonds, the temperature and urea-concentration dependence of (1)HN and (15)N chemical shifts, and native-state hydrogen exchange at urea concentrations where the domain is predominantly folded. The native-state hydrogen-exchange data suggest that the six-stranded β-barrel core of the I-domain is more stable against unfolding than a smaller subdomain comprised of a short α-helix and three-stranded β-sheet. H-bonds, separately determined from solvent protection and (3)JNC' H-bond couplings, are identified with an accuracy of 90% by (1)HN temperature coefficients. The accuracy is improved to 95% when (15)N temperature coefficients are also included. In contrast, the urea dependence of (1)HN and (15)N chemical shifts is unrelated to H-bonding. The protein segments with the largest chemical-shift changes in the presence of urea show curved or sigmoidal titration curves suggestive of direct urea binding. Nuclear Overhauser effects to urea for these segments are also consistent with specific urea-binding sites in the I-domain. Taken together, the results support a mechanism of urea unfolding in which denaturant binds to distinct sites in the I-domain. Disordered segments bind urea more readily than regions in stable secondary structure. The locations of the putative urea-binding sites correlate with the lower stability of the structure against solvent exchange, suggesting that partial unfolding of the structure is related to urea accessibility.

Phosphorylation in the catalytic cleft stabilizes and attracts domains of a phosphohexomutase

Phosphorylation can modulate the activities of enzymes. The phosphoryl donor in the catalytic cleft of α-D-phosphohexomutases is transiently dephosphorylated while the reaction intermediate completes a 180° reorientation within the cleft. The phosphorylated form of 52 kDa bacterial phosphomannomutase/phosphoglucomutase is less accessible to dye or protease, more stable to chemical denaturation, and widely stabilized against NMR-detected hydrogen exchange across the core of domain 3 to juxtaposed domain 4 (each by ≥ 1.3 kcal/mol) and parts of domains 1 and 2. However, phosphorylation accelerates hydrogen exchange in specific regions of domains 1 and 2, including a metal-binding residue in the active site. Electrostatic field lines reveal attraction across the catalytic cleft between phosphorylated Ser-108 and domain 4, but repulsion when Ser-108 is dephosphorylated. Molecular dynamics (MD) simulated the dephosphorylated form to be expanded due to enhanced rotational freedom of domain 4. The contacts and fluctuations of the MD trajectories enabled correct simulation of more than 80% of sites that undergo either protection or deprotection from hydrogen exchange due to phosphorylation. Electrostatic attraction in the phosphorylated enzyme accounts for 1) domain 4 drawing closer to domains 1 and 3; 2) decreased accessibility; and 3) increased stability within these domains. The electrostriction due to phosphorylation may help capture substrate, whereas the opening of the cleft upon transient dephosphorylation allows rotation of the intermediate. The long-range effects of phosphorylation on hydrogen exchange parallel reports on protein kinases, suggesting a conceptual link among these multidomain, phosphoryl transfer enzymes.

Mechanistic insights into phosphoprotein-binding FHA domains

[Structure: see text]. FHA domains are protein modules that switch signals in diverse biological pathways by monitoring the phosphorylation of threonine residues of target proteins. As part of the effort to gain insight into cellular avoidance of cancer, FHA domains involved in the cellular response to DNA damage have been especially well-characterized. The complete protein where the FHA domain resides and the interaction partners determine the nature of the signaling. Thus, a key biochemical question is how do FHA domains pick out their partners from among thousands of alternatives in the cell? This Account discusses the structure, affinity, and specificity of FHA domains and the formation of their functional structure. Although FHA domains share sequence identity at only five loop residues, they all fold into a beta-sandwich of two beta-sheets. The conserved arginine and serine of the recognition loops recognize the phosphorylation of the threonine targeted. Side chains emanating from loops that join beta-strand 4 with 5, 6 with 7, or 10 with 11 make specific contacts with amino acids of the ligand that tailor sequence preferences. Many FHA domains choose a partner in extended conformation, somewhat according to the residue three after the phosphothreonine in sequence (pT + 3 position). One group of FHA domains chooses a short carboxylate-containing side chain at pT + 3. Another group chooses a long, branched aliphatic side chain. A third group prefers other hydrophobic or uncharged polar side chains at pT + 3. However, another FHA domain instead chooses on the basis of pT - 2, pT - 3, and pT + 1 positions. An FHA domain from a marker of human cancer instead chooses a much longer protein fragment that adds a beta-strand to its beta-sheet and that presents hydrophobic residues from a novel helix to the usual recognition surface. This novel recognition site and more remote sites for the binding of other types of protein partners were predicted for the entire family of FHA domains by a bioinformatics approach. The phosphopeptide-dependent dynamics of an FHA domain, SH2 domain, and PTB domain suggest a common theme: rigid, preformed binding surfaces support van der Waals contacts that provide favorable binding enthalpy. Despite the lack of pronounced conformational changes in FHA domains linked to binding events, more subtle adjustments may be possible. In the one FHA domain tested, phosphothreonine peptide binding is accompanied by increased flexibility just outside the binding site and increased rigidity across the beta-sandwich. The folding of the same FHA domain progresses through near-native intermediates that stabilize the recognition loops in the center of the phosphoprotein-binding surface; this may promote rigidity in the interface and affinity for targets phosphorylated on threonine.

Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.