An expanded view of the protein folding landscape of PDZ domains

Most protein domains fold in an apparently co-operative and two-state manner with only the native and denatured states significantly populated at any experimental condition. However, the protein folding energy landscape is often rugged and different transition states may be rate limiting for the folding reaction under different conditions, as seen for the PDZ protein domain family. We have here analyzed the folding kinetics of two PDZ domains and found that a previously undetected third transition state is rate limiting under conditions that stabilize the native state relative to the denatured state. In light of these results, we have re-analyzed previous folding data on PDZ domains and present a unified folding mechanism with three distinct transition states separated by two high-energy intermediates. Our data show that sequence composition tunes the relative stabilities of folding transition states within the PDZ family, while the overall mechanism is determined by topology. This model captures the kinetic folding mechanism of all PDZ domains studied to date.     

Folding and Misfolding of a PDZ Tandem Repeat

Although the vast majority of the human proteome is represented by multi-domain proteins, the study of multi-domain folding and misfolding is a relatively poorly explored field. The protein Whirlin is a multi-domain scaffolding protein expressed in the inner ear. It is characterized by the presence of tandem repeats of PDZ domains. The first two PDZ domains of Whirlin (PDZ1 and PDZ2 – namely P1P2) are structurally close and separated by a disordered short linker. We recently described the folding mechanism of the P1P2 tandem. The difference in thermodynamic stability of the two domains allowed us to selectively unfold one or both PDZ domains and to pinpoint the accumulation of a misfolded intermediate, which we demonstrated to retain physiological binding activity. In this work, we provide an extensive characterization of the folding and unfolding of P1P2. Based on the observed data, we describe an integrated kinetic analysis that satisfactorily fits the experiments and provides a valuable model to interpret multi-domain folding. The experimental and analytical approaches described in this study may be of general interest for the interpretation of complex multi-domain protein folding kinetics.     

Morphogenesis of a protein: folding pathways and the energy landscape

Current knowledge on the reaction whereby a protein acquires its native three-dimensional structure was obtained by and large through characterization of the folding mechanism of simple systems. Given the multiplicity of amino acid sequences and unique folds, it is not so easy, however, to draw general rules by comparing folding pathways of different proteins. In fact, quantitative comparison may be jeopardized not only because of the vast repertoire of sequences but also in view of a multiplicity of structures of the native and denatured states. We have tackled the problem of the relationships between the sequence information and the folding pathway of a protein, using a combination of kinetics, protein engineering and computational methods, applied to relatively simple systems. Our strategy has been to investigate the folding mechanism determinants using two complementary approaches, i.e. (i) the study of members of the same family characterized by a common fold, but substantial differences in amino acid sequence, or (ii) heteromorphic pairs characterized by largely identical sequences but with different folds. We discuss some recent data on protein-folding mechanisms by presenting experiments on different members of the PDZ domain family and their circularly permuted variants. Characterization of the energetics and structures of intermediates and TSs (transition states), obtained by Φ-value analysis and restrained MD (molecular dynamics) simulations, provides a glimpse of the malleability of the dynamic states and of the role of the topology of the native states and of the denatured states in dictating folding and misfolding pathways.     

The Folding Mechanism of BBL: Plasticity of Transition-State Structure Observed within an Ultrafast Folding Protein Family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Φ-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation-condensation through structurally very similar, polarized transition states. Here, we present a Φ-value analysis of the family member BBL and found that it also folds by a nucleation-condensation mechanism. The mean Φ values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Φ values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.     

The Role of Backbone Hydrogen Bonds in the Transition State for Protein Folding of a PDZ Domain

Backbone hydrogen bonds are important for the structure and stability of proteins. However, since conventional site-directed mutagenesis cannot be applied to perturb the backbone, the contribution of these hydrogen bonds in protein folding and stability has been assessed only for a very limited set of small proteins. We have here investigated effects of five amide-to-ester mutations in the backbone of a PDZ domain, a 90-residue globular protein domain, to probe the influence of hydrogen bonds in a β-sheet for folding and stability. The amide-to-ester mutation removes NH-mediated hydrogen bonds and destabilizes hydrogen bonds formed by the carbonyl oxygen. The overall stability of the PDZ domain generally decreased for all amide-to-ester mutants due to an increase in the unfolding rate constant. For this particular region of the PDZ domain, it is therefore clear that native hydrogen bonds are formed after crossing of the rate-limiting barrier for folding. Moreover, three of the five amide-to-ester mutants displayed an increase in the folding rate constant suggesting that the hydrogen bonds are involved in non-native interactions in the transition state for folding.     

Supramodular nature of GRIP1 revealed by the structure of its PDZ12 tandem in complex with the carboxyl tail of Fras1

The scaffold protein GRIP1 (glutamate receptor interacting protein 1) binds to and regulates both the trafficking and membrane organization of a large number of transmembrane proteins. Mutation of GRIP1 in mice displays essentially the same phenotype of the mutations of Fras1 or Frem2, which are the animal models of the human genetic disorder Fraser syndrome. However, the molecular basis governing the interaction between GRIP1 and Fras1/Frem2 is unknown. Here, we show that interaction between Fras1 and GRIP1 requires the first two PDZ domains (PDZ1 and PDZ2) to be connected in tandem, as the folding of PDZ1 strictly depends on the covalent attachment of PDZ2. The crystal structure of GRIP1 PDZ12 in complex with the Fras1 C-terminal peptide reveals that the PDZ12 tandem forms a supramodule in which only the peptide-binding groove of PDZ1 is bound with the Fras1 peptide. The GRIP1 PDZ12/Fras1 peptide complex not only provides a mechanistic explanation of the link between GRIP1 and the Fraser syndrome but may also serve as a foundation for searching for potential mutations in GRIP1 that could lead to the Fraser syndrome.     

Different Members of a Simple Three-Helix Bundle Protein Family Have Very Different Folding Rate Constants and Fold by Different Mechanisms

The 15th, 16th, and 17th repeats of chicken brain α-spectrin (R15, R16, and R17, respectively) are very similar in terms of structure and stability. However, R15 folds and unfolds 3 orders of magnitude faster than R16 and R17. This is unexpected. The rate-limiting transition state for R15 folding is investigated using protein engineering methods (Φ-value analysis) and compared with previously completed analyses of R16 and R17. Characterisation of many mutants suggests that all three proteins have similar complexity in the folding landscape. The early rate-limiting transition states of the three domains are similar in terms of overall structure, but there are significant differences in the patterns of Φ-values. R15 apparently folds via a nucleation-condensation mechanism, which involves concomitant folding and packing of the A- and C-helices, establishing the correct topology. R16 and R17 fold via a more framework-like mechanism, which may impede the search to find the correct packing of the helices, providing a possible explanation for the fast folding of R15.     

The binding of the PDZ tandem of syntenin to target proteins

PDZ domains are among the most abundant protein modules in the known genomes. Their main function is to provide scaffolds for membrane-associated protein complexes by binding to the cytosolic, C-terminal fragments of receptors, channels, and other integral membrane proteins. Here, using both heteronuclear NMR and single crystal X-ray diffraction, we show how peptides with different sequences, including those corresponding to the C-termini of syndecan, neurexin, and ephrin B, can simultaneously bind to both PDZ domains of the scaffolding protein syntenin. The PDZ2 domain binds these peptides in the canonical fashion, and an induced fit mechanism allows for the accommodation of a range of side chains in the P(0) and P(-)(2) positions. However, binding to the PDZ1 domain requires that the target peptide assume a noncanonical conformation. These data help explain how syntenin, and perhaps other PDZ-containing proteins, may preferentially bind to dimeric and clustered targets, and provide a mechanistic explanation for the previously reported cooperative ligand binding by syntenin's two PDZ domains.     

The Folding of a Family of Three-Helix Bundle Proteins: Spectrin R15 Has a Robust Folding Nucleus,Unlike Its Homologous Neighbours

Three homologous spectrin domains have remarkably different folding characteristics. We have previously shown that the slow-folding R16 and R17 spectrin domains can be altered to resemble the fast folding R15, in terms of speed of folding (and unfolding), landscape roughness and folding mechanism, simply by substituting five residues in the core. Here we show that, by contrast, R15 cannot be engineered to resemble R16 and R17. It is possible to engineer a slow-folding version of R15, but our analysis shows that this protein neither has a rougher energy landscape nor does change its folding mechanism. Quite remarkably, R15 appears to be a rare example of a protein with a folding nucleus that does not change in position or in size when its folding nucleus is disrupted. Thus, while two members of this protein family are remarkably plastic, the third has apparently a restricted folding landscape.     

Kinetic partitioning mechanism governs the folding of the third FnIII domain of tenascin-C: evidence at the single-molecule level

Statistical mechanics and molecular dynamics simulations proposed that the folding of proteins can follow multiple parallel pathways on a rugged energy landscape from unfolded state en route to their folded native states. Kinetic partitioning mechanism is one of the possible mechanisms underlying such complex folding dynamics. Here, we use single-molecule atomic force microscopy technique to directly probe the multiplicity of the folding pathways of the third fibronectin type III domain from the extracellular matrix protein tenascin-C (TNfn3). By stretching individual (TNfn3)₈ molecules, we forced TNfn3 domains to undergo mechanical unfolding and refolding cycles, allowing us to directly observe the folding pathways of TNfn3. We found that, after being mechanically unraveled and then relaxed to zero force, TNfn3 follows multiple parallel pathways to fold into their native states. The majority of TNfn3 fold into the native state in a simple two-state fashion, while a small percentage of TNfn3 were found to be trapped into kinetically stable folding intermediate states with well-defined three-dimensional structures. Furthermore, the folding of TNfn3 was also influenced by its neighboring TNfn3 domains. Complex misfolded states of TNfn3 were observed, possibly due to the formation of domain-swapped dimeric structures. Our studies revealed the ruggedness of the folding energy landscape of TNfn3 and provided direct experimental evidence that the folding dynamics of TNfn3 are governed by the kinetic partitioning mechanism. Our results demonstrated the unique capability of single-molecule AFM to probe the folding dynamics of proteins at the single-molecule level.     

Folding of tandem-linked domains

One of the factors, which influences protein folding in vivo, is a linkage of protein domains into multidomain tandems. However, relatively little is known about the impact of domain connectivity on protein folding mechanisms. In this article, we use coarse grained models of proteins to explore folding of tandem-linked domains (TLD). We found TLD folding to follow two scenarios. In the first, the tandem connectivity produces relatively minor impact on folding and the mechanisms of folding of tandem-linked and single domains remain similar. The second scenario involves qualitative changes in folding mechanism because of tandem linkage. As a result, protein domains, which fold via two-state mechanism as single isolated domains, may form new stable intermediates when inserted into tandems. The new intermediates are created by topological constraints imposed by the linkers between domains. In both cases tandem linkage slows down folding. We propose that the impact of tandem connectivity can be minimized, if the terminal secondary structure elements (SSEs) are flexible. In particular, two factors appear to facilitate TLD folding: (1) the interactions between terminal SSE are poorly ordered in the folding transition state, whereas nonterminal SSE are better structured, (2) the interactions between terminal SSE are weak in the native state. We apply these findings to wild-type proteins by examining experimental phi-value data and by performing all-atom molecular dynamics simulations. We show that immunoglobulin-like domains appear to utilize the factors, which minimize the impact of tandem connectivity on their folding. Several single domain proteins, which are likely to misfold in tandems, are also identified.     

Interdomain chaperoning between PSD-95, Dlg, and Zo-1 (PDZ) domains of glutamate receptor-interacting proteins.

The multiple PSD-95, Dlg, and Zo-1 (PDZ) domain protein, glutamate receptor-interacting protein (GRIP), is involved in the clustering and trafficking of the alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionate receptor by directly binding to the cytoplasmic tail of the receptor's GluR2 subunit. Both the forth and fifth PDZ domains (PDZ4 and PDZ5) of GRIP are required for effective binding to the receptor. Using NMR and circular dichroism spectroscopic techniques, we show that PDZ5 is completely unstructured in solution. Freshly prepared PDZ4 is largely folded, but the domain can spontaneously unfold. Neither PDZ4 nor PDZ5 binds to GluR2 in solution. Unexpectedly, when PDZ4 and PDZ5 are covalently connected (i.e. PDZ45), both PDZ domains become well folded and stable in solution. The covalent linkage of the two PDZ domains is essential for proper folding of the tandem PDZ domains and its effective binding to GluR2. The interdomain chaperoning effect observed in the PDZ domains of GRIP represents a previously uncharacterized function of PDZ domains.     

The folding pathway of T4 lysozyme: the high-resolution structure and folding of a hidden intermediate

Folding intermediates have been detected and characterized for many proteins. However, their structures at atomic resolution have only been determined for two small single domain proteins: Rd-apocytochrome b(562) and engrailed homeo domain. T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. An early native-state hydrogen exchange experiment identified an intermediate with the C-terminal domain folded and the N-terminal domain unfolded. We have used a native-state hydrogen exchange-directed protein engineering approach to populate this intermediate and demonstrated that it is on the folding pathway and exists after the rate-limiting step. Here, we determined its high-resolution structure and the backbone dynamics by multi-dimensional NMR methods. We also characterized the folding behavior of the intermediate using stopped-flow fluorescence, protein engineering, and native-state hydrogen exchange. Unlike the folding intermediates of the two single-domain proteins, which have many non-native side-chain interactions, the structure of the hidden folding intermediate of T4 lysozyme is largely native-like. It folds like many small single domain proteins. These results have implications for understanding the folding mechanism and evolution of multi-domain proteins.     

Structure-Based Derivation of Protein Folding Intermediates and Energies from Optical Tweezers

Optical tweezers (OTs) measure the force-dependent time-resolved extension of a single macromolecule tethered between two trapped beads. From this measurement, it is possible to determine the folding intermediates, energies, and kinetics of the macromolecule. Previous data analysis generally has used the extension as a reaction coordinate to characterize the observed folding transitions. Despite its convenience, the extension poorly describes folding in the absence of force. Here, we chose the contour length of the unfolded polypeptide as a reaction coordinate and modeled the extensions of protein structures along their predicted folding pathways based on high-resolution structures of the proteins in their native states. We included the extension in our model to calculate the total extensions, energies, and transition rates of the proteins as a function of force. We fit these calculations to the corresponding experimental measurements and obtained the best-fit conformations and energies of proteins in different folding states. We applied our method to analyze single-molecule trajectories of two representative protein complexes responsible for membrane fusion, the HIV-1 glycoprotein 41 and the synaptic SNARE proteins, which involved transitions between two and five states, respectively. Nonlinear fitting of the model to the experimental data revealed the structures of folding intermediates and transition states and their associated energies. Our results demonstrate that the contour length is a useful reaction coordinate to characterize protein folding and that intrinsic extensions of protein structures should be taken into account to properly derive the conformations and energies of protein folding intermediates from single-molecule manipulation experiments.     

Modular organization of the PDZ domains in the human discs-large protein suggests a mechanism for coupling PDZ domain-binding proteins to ATP and the membrane cytoskeleton

The human homologue (hDIg) of the Drosophila discs-large tumor suppressor (DIg) is a multidomain protein consisting of a carboxyl- terminal guanylate kinase-like domain, an SH3 domain, and three slightly divergent copies of the PDZ (DHR/GLGF) domain. Here have examined the structural organization of the three PDZ domains of hDIg using a combination of protease digestion and in vitro binding measurements. Our results show that the PDZ domains are organized into two conformationally stable modules one (PDZ, consisting of PDZ domains 1 and 2, and the other (PDZ) corresponding to the third PDZ domain. Using amino acid sequencing and mass spectrometry, we determined the boundaries of the PDZ domains after digestion with endoproteinase Asp- N, trypsin, and alpha-chymotrypsin. The purified PDZ1+2, but not the PDZ3 domain, contains a high affinity binding site for the cytoplasmic domain of Shaker-type K+ channels. Similarly, we demonstrate that the PDZ1+2 domain can also specifically bind to ATP. Furthermore, we provide evidence for an in vivo interaction between hDIg and protein 4.1 and show that the hDIg protein contains a single high affinity protein 4.1-binding site that is not located within the PDZ domains. The results suggest a mechanism by which PDZ domain-binding proteins may be coupled to ATP and the membrane cytoskeleton via hDlg.     

The kinetics of PDZ domain-ligand interactions and implications for the binding mechanism

PDZ domains are protein adapter modules present in a few hundred human proteins. They play important roles in scaffolding and signal transduction. PDZ domains usually bind to the C termini of their target proteins. To assess the binding mechanism of this interaction we have performed the first in-solution kinetic study for PDZ domains and peptides corresponding to target ligands. Both PDZ3 from postsynaptic density protein 95 and PDZ2 from protein tyrosine phosphatase L1 bind their respective target peptides through an apparent A + B --> A.B mechanism without rate-limiting conformational changes. But a mutant with a fluorescent probe (Trp) outside of the binding pocket suggests that slight changes in the structure take place upon binding in protein tyrosine phosphatase-L1 PDZ2. For PDZ3 from postsynaptic density protein 95 the pH dependence of the binding reaction is consistent with a one-step mechanism with one titratable group. The salt dependence of the interaction shows that the formation of electrostatic interactions is rate-limiting for the association reaction but not for dissociation of the complex.     

Evidence for ezrin-radixin-moesin-binding phosphoprotein 50 (EBP50) self-association through PDZ-PDZ interactions

Ezrin-radixin-moesin (ERM)-binding phosphoprotein 50 (EBP50) is a versatile membrane-cytoskeleton linking protein that binds to the COOH-tail of specific integral membrane proteins through its two PDZ domains. These EBP50 binding interactions have been implicated in sequestering interactive sets of proteins into common microdomains, regulating the activity of interacting proteins, and modulating membrane protein trafficking. With only two PDZ domains, it is unclear how EBP50 forms multiprotein complexes. Other PDZ proteins increase their breadth and diversity of protein interactions through oligomerization. Hypothesizing that EBP50 self-associates to amplify its functional capacity, far-Western blotting of cholangiocyte epithelial cell proteins with EBP50 fusion protein revealed that EBP50 binds to a 50-kDa protein. Far-Western blotting of EBP50 isolated by two-dimensional gel electrophoresis or immunoprecipitation demonstrates that the 50-kDa binding partner is itself EBP50. Further, co-transfection/co-precipitation studies show the self-association can occur in an intracellular environment. In vitro analysis of the EBP50-EBP50 binding interaction indicates it is both saturable and of relatively high affinity. Analysis of truncated EBP50 proteins indicates EBP50 self-association is mediated through its PDZ domains. The ability to self-associate provides a mechanism for EBP50 to expand its capacity to form multiprotein complexes and regulate membrane transport events.     

Structural and functional analysis of the PDZ domains of human HtrA1 and HtrA3

High-temperature requirement A (HtrA) and its homologs contain a serine protease domain followed by one or two PDZ domains. Bacterial HtrA proteins and the mitochondrial protein HtrA2/Omi maintain cell function by acting as both molecular chaperones and proteases to manage misfolded proteins. The biological roles of the mammalian family members HtrA1 and HtrA3 are less clear. We report a detailed structural and functional analysis of the PDZ domains of human HtrA1 and HtrA3 using peptide libraries and affinity assays to define specificity, structural studies to view the molecular details of ligand recognition, and alanine scanning mutagenesis to investigate the energetic contributions of individual residues to ligand binding. In common with HtrA2/Omi, we show that the PDZ domains of HtrA1 and HtrA3 recognize hydrophobic polypeptides, and while C-terminal sequences are preferred, internal sequences are also recognized. However, the details of the interactions differ, as different domains rely on interactions with different residues within the ligand to achieve high affinity binding. The results suggest that mammalian HtrA PDZ domains interact with a broad range of hydrophobic binding partners. This promiscuous specificity resembles that of bacterial HtrA family members and suggests a similar function for recognizing misfolded polypeptides with exposed hydrophobic sequences. Our results support a common activation mechanism for the HtrA family, whereby hydrophobic peptides bind to the PDZ domain and induce conformational changes that activate the protease. Such a mechanism is well suited to proteases evolved for the recognition and degradation of misfolded proteins.