Cooperative folding in a multi-domain protein

Most protein domains are found in multi-domain proteins, yet most studies of protein folding have concentrated on small, single-domain proteins or on isolated domains from larger proteins. Spectrin domains are small (106 amino acid residues), independently folding domains consisting of three long alpha-helices. They are found in multi-domain proteins with a number of spectrin domains in tandem array. Structural studies have shown that in these arrays the last helix of one domain forms a continuous helix with the first helix of the following domain. It has been demonstrated that a number of spectrin domains are stabilised by their neighbours. Here we investigate the molecular basis for cooperativity between adjacent spectrin domains 16 and 17 from chicken brain alpha-spectrin (R16 and R17). We show that whereas the proteins unfold as a single cooperative unit at 25 degrees C, cooperativity is lost at higher temperatures and in the presence of stabilising salts. Mutations in the linker region also cause the cooperativity to be lost. However, the cooperativity does not rely on specific interactions in the linker region alone. Most mutations in the R17 domain cause a decrease in cooperativity, whereas proteins with mutations in the R16 domain still fold cooperatively. We propose a mechanism for this behaviour.     

Understanding the role of domain–domain linkers in the spatial orientation of domains in multi-domain proteins

Inter-domain linkers (IDLs)’ bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain–domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.     

Domain interdependence in the biosynthetic assembly of CFTR

The dimerization of their two nucleotide binding domains (NBDs) in a so-called "nucleotide-sandwich" is the hallmark of ATP cassette binding (ABC) proteins and the basis of their catalytic activities. The major disease-causing mutation in the cystic fibrosis transmembrane conductance regulator (CFTR or ABCC7), deletion of Phe508 in NBD1, does not grossly alter the structure of that domain but prevents conformational maturation of the whole CFTR protein, possibly by disrupting the native interaction between NBD1 and NBD2. However, the role of inter-domain interactions in CFTR folding has been brought into question by a recent report that all CFTR domains fold independently. Here we show that in addition to domain folding, correct inter-domain assembly is essential to form a stable unit that satisfies endoplasmic reticulum (ER) quality control. N-terminal domains depend on their more C-terminal neighbors, most essentially the second membrane-spanning domain (MSD2) but significantly, not NBD2. Wild-type C-terminal truncation constructs, completely devoid of NBD2 are transported out of the ER and to the cell surface where they form characteristic CFTR chloride channels with low open probability. The DeltaNBD2 wild-type protein matures and has similar stability as its full-length counterpart. Therefore, the catalytically crucial inter-NBD associations are not required to satisfy ER quality control mechanisms. The DeltaF508 mutation arrests the maturation of DeltaNBD2 just as it does full-length CFTR, indicating that DeltaF508 perturbs other portions of the molecule in addition to NBD2. We find that the mutation prevents formation of a compact MSD1, reflected in its susceptibility to protease digestion. This perturbation of MSD1 may in turn prevent its normal integration with MSD2. The dispensability of NBD2 in the folding of more N-terminal domains stands in contrast to the known hypersensitivity to proteolysis of NBD2 in the DeltaF508 protein.     

Comparing domain interactions within antibody Fabs with kappa and lambda light chains

IgG antibodies are multi-domain proteins with complex inter-domain interactions. Human IgG heavy chains (HCs) associate with light chains (LCs) of the κ or λ isotype to form mature antibodies capable of binding antigen. The HC/LC interaction involves 4 domains: VH and CH1 from the HC and VL and CL from the LC. Human Fabs with κ LCs have been well characterized for their unfolding behaviors and demonstrate a significant level of cooperativity and stabilization when all 4 domains are intact. Very little is known regarding the thermodynamic properties of human Fabs with λ LCs. Here, we dissect the domain contributions to Fab stability for both κ and λ LC-containing Fabs. We find the cooperativity of unfolding between the constant domains, CH1/Cλ, and variable domains, VH/Vλ, within λ LC-containing Fabs is significantly weaker than that of κ LC-containing Fabs. The data suggests there may not be an evolutionary necessity for strong variable/constant domain cooperativity within λ LC-containing Fabs. After investigating the biophysical properties of Fabs with mismatched variable and constant domain subunits (e.g., VH/Vκ paired with CH1/Cλ or T cell receptor Cα/Cβ), the major role of the constant domains for both κ- and λ-containing Fabs may be to reduce the hydrophobic exposure at the VH/VL interface. Even though Fabs with these non-native pairings were thermodynamically less stable, they secreted well from mammalian cells as well behaved monodisperse proteins, which was in contrast to what was observed with the VH/Vκ and VH/Vλ scFvs that secreted as a mixture of monomer and aggregates.     

A conserved interdomain interaction is a determinant of folding cooperativity in the GST fold

The canonical glutathione transferase (GST) fold found in many monomeric and dimeric proteins consists of two domains that differ in structure and conformational dynamics. However, no evidence exists that the two domains unfold/fold independently at equilibrium, indicating the significance of interdomain interactions in governing cooperativity between domains. Bioinformatics analyses indicate the interdomain interface of the GST fold is large, predominantly hydrophobic with a high packing density explaining cooperative interdomain behavior. Structural alignments reveal a topologically conserved lock-and-key interaction across the domain interface in which a bulky hydrophobic residue ("key") protrudes from the surface of the N-domain and inserts into a pocket ("lock") in the C-domain. To better understand the molecular basis for the contribution of interdomain interactions toward cooperativity within the GST fold in the absence of any influence from quaternary interactions, studies were done with two monomeric GST proteins: Escherichia coli Grx2 (EcGrx2) and human CLIC1 (hCLIC1). Replacing the methionine "key" residue with alanine is structurally nondisruptive, whereas it significantly diminishes the folding cooperativity of both proteins. The loss in cooperativity between domains in the mutants is reflected by a change in the equilibrium folding mechanism from a wild-type two-state process to a three-state process, populating a stable folding intermediate.     

Domain rearrangements in protein evolution

Most eukaryotic proteins are multi-domain proteins that are created from fusions of genes, deletions and internal repetitions. An investigation of such evolutionary events requires a method to find the domain architecture from which each protein originates. Therefore, we defined a novel measure, domain distance, which is calculated as the number of domains that differ between two domain architectures. Using this measure the evolutionary events that distinguish a protein from its closest ancestor have been studied and it was found that indels are more common than internal repetition and that the exchange of a domain is rare. Indels and repetitions are common at both the N and C-terminals while they are rare between domains. The evolution of the majority of multi-domain proteins can be explained by the stepwise insertions of single domains, with the exception of repeats that sometimes are duplicated several domains in tandem. We show that domain distances agree with sequence similarity and semantic similarity based on gene ontology annotations. In addition, we demonstrate the use of the domain distance measure to build evolutionary trees. Finally, the evolution of multi-domain proteins is exemplified by a closer study of the evolution of two protein families, non-receptor tyrosine kinases and RhoGEFs.     

An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding–unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.     

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two‐state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non‐cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier‐less “downhill” folding, as well as for continuous “uphill” unfolding transitions, indicate that gradual non‐cooperative processes may be ubiquitous features on the free energy landscape of protein folding.     

Domain motions in GroEL upon binding of an oligopeptide

GroEL assists protein folding by preventing the interaction of partially folded molecules with other non-native proteins. It binds them, sequesters them, and then releases them so that they can fold in an ATP-driven cycle. Previous studies have also shown that protein substrates, GroES, and oligopeptides bind to partially overlapped sites on the apical domain surfaces of GroEL. In this study, we have determined the crystal structure at 3.0A resolution of a symmetric (GroEL-peptide)(14) complex. The binding of each of these small 12 amino acid residue peptides to GroEL involves interactions between three adjacent apical domains of GroEL. Each peptide interacts primarily with a single GroEL subunit. Residues R231 and R268 from adjacent subunits isolate each substrate-binding pocket, and prevent bound substrates from sliding into adjacent binding pockets. As a consequence of peptide binding, domains rotate and inter-domain interactions are greatly enhanced. The direction of rotation of the apical domain of each GroEL subunit is opposite to that of its intermediate domain. Viewed from outside, the apical domains rotate clockwise within one GroEL ring, while the ATP-induced apical domain rotation is counter-clockwise.     

The Modular Organization of Domain Structures: Insights into Protein–Protein Binding

Domains are the building blocks of proteins and play a crucial role in protein–protein interactions. Here, we propose a new approach for the analysis and prediction of domain–domain interfaces. Our method, which relies on the representation of domains as residue-interacting networks, finds an optimal decomposition of domain structures into modules. The resulting modules comprise highly cooperative residues, which exhibit few connections with other modules. We found that non-overlapping binding sites in a domain, involved in different domain–domain interactions, are generally contained in different modules. This observation indicates that our modular decomposition is able to separate protein domains into regions with specialized functions. Our results show that modules with high modularity values identify binding site regions, demonstrating the predictive character of modularity. Furthermore, the combination of modularity with other characteristics, such as sequence conservation or surface patches, was found to improve our predictions. In an attempt to give a physical interpretation to the modular architecture of domains, we analyzed in detail six examples of protein domains with available experimental binding data. The modular configuration of the TEM1-β-lactamase binding site illustrates the energetic independence of hotspots located in different modules and the cooperativity of those sited within the same modules. The energetic and structural cooperativity between intramodular residues is also clearly shown in the example of the chymotrypsin inhibitor, where non–binding site residues have a synergistic effect on binding. Interestingly, the binding site of the T cell receptor β chain variable domain 2.1 is contained in one module, which includes structurally distant hot regions displaying positive cooperativity. These findings support the idea that modules possess certain functional and energetic independence. A modular organization of binding sites confers robustness and flexibility to the performance of the functional activity, and facilitates the evolution of protein interactions.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Conformational transitions in adenylate kinase. Allosteric communication reduces misligation

Large conformational changes in the LID and NMP domains of adenylate kinase (AKE) are known to be key to ligand binding and catalysis, yet the order of binding events and domain motion is not well understood. Combining the multiple available structures for AKE with the energy landscape theory for protein folding, a theoretical model was developed for allostery, order of binding events, and efficient catalysis. Coarse-grained models and nonlinear normal mode analysis were used to infer that intrinsic structural fluctuations dominate LID motion, whereas ligand-protein interactions and cracking (local unfolding) are more important during NMP motion. In addition, LID-NMP domain interactions are indispensable for efficient catalysis. LID domain motion precedes NMP domain motion, during both opening and closing. These findings provide a mechanistic explanation for the observed 1:1:1 correspondence between LID domain closure, NMP domain closure, and substrate turnover. This catalytic cycle has likely evolved to reduce misligation, and thus inhibition, of AKE. The separation of allosteric motion into intrinsic structural fluctuations and ligand-induced contributions can be generalized to further our understanding of allosteric transitions in other proteins.     

Complex folding kinetics of a multidomain protein

Spectrin domains are three-helix bundles, commonly found in large tandem arrays. Equilibrium studies have shown that spectrin domains are significantly stabilized by their neighbors. In this work we show that domain:domain interactions can also have profound effects on their kinetic behavior. We have studied the folding of a tandem pair of spectrin domains (R1617) using a combination of single- and double-jump stopped flow experiments (monitoring folding by both circular dichroism and fluorescence). Mutant proteins were also used to investigate the complex folding kinetics. We find that, although the domains fold and unfold individually, there is a single rate-determining step for both folding and unfolding of the protein. This is consistent with the equilibrium observation of cooperative folding of the entire two-domain protein. The results may have important biological implications. Not only will the protein fold more efficiently during cotranslational folding, but the ability of the multidomain protein to withstand thermal unfolding in the cell will be dramatically increased. This study suggests that caution has to be exercised when extrapolating from single domains to larger proteins with a number of independently folding modules arranged in tandem. The multidomain protein spectrin is certainly more than "the sum of its parts".     

Two immunoglobulin tandem proteins with a linking β-strand reveal unexpected differences in cooperativity and folding pathways

The study of the folding of single domains, in the context of their multidomain environment, is important because more than 70% of eukaryotic proteins are composed of multiple domains. The structures of the tandem immunoglobulin (Ig) domain pairs A164-A165 and A168-A169, from the A-band of the giant muscle protein titin, reveal that they form tightly associated domain arrangements, connected by a continuous β-strand. We investigate the thermodynamic and kinetic properties of these tandem domain pairs. While A164-A165 apparently behaves as a single cooperative unit at equilibrium, unfolding without the accumulation of a large population of intermediates, domains in A168-A169 behave independently. Although A169 appears to be stabilized in the tandem protein, we show that this is due to nonspecific stabilization by extension. We elucidate the folding and unfolding pathways of both tandem pairs and show that cooperativity in A164-A165 is a manifestation of the relative refolding and unfolding rate constants of each individual domain. We infer that the differences between the two tandem pairs result from a different pattern of interactions at the domain/domain interface.     

De novo folding of GFP fusion proteins: high efficiency in eukaryotes but not in bacteria

Eukaryotic genomes encode a considerably higher fraction of multi-domain proteins than their prokaryotic counterparts. It has been postulated that efficient co-translational and sequential domain folding has facilitated the explosive evolution of multi-domain proteins in eukaryotes by the recombination of pre-existent domains. Here, we tested whether eukaryotes and bacteria differ generally in the folding efficiency of multi-domain proteins generated by domain recombination. To this end, we compared the folding behavior of a series of recombinant proteins comprised of green fluorescent protein (GFP) fused to four different robustly folding proteins through six different linkers upon expression in Escherichia coli and the yeast Saccharomyces cerevisiae. We found that, unlike yeast, bacteria are remarkably inefficient at folding these fusion proteins, even at comparable levels of expression. In vitro and in vivo folding experiments demonstrate that the GFP domain imposes significant constraints on de novo folding of its fusion partners in bacteria, consistent with a largely post-translational folding mechanism. This behavior may result from an interference of GFP with adjacent domains during folding due to the particular topology of the beta-barrel GFP structure. By following the accumulation of enzymatic activity, we found that the rate of appearance of correctly folded fusion protein per ribosome is indeed considerably higher in yeast than in bacteria.     

Trigger factor lacking the PPIase domain can enhance the folding of eukaryotic multi-domain proteins in Escherichia coli

Recombinant expression of eukaryotic proteins in bacteria often results in misfolding and aggregation. The ribosome-binding Trigger factor (TF) is the first molecular chaperone that interacts with nascent polypeptide chains in bacteria. Here we show that mutant TF lacking the PPIase domain (TFNC) is more efficient than wild-type TF in enhancing the folding yield of multi-domain proteins such as firefly luciferase. We find that TFNC has a shorter residence time on nascent chains, thus facilitating co-translational folding. By delaying folding relative to translation, the PPIase domain may increase the propensity of misfolding for certain eukaryotic proteins that rely on a mechanism of co-translational, domain-wise folding.     

UNC-45B Chaperone: The Role of its Domains in the Interaction with the Myosin Motor Domain

The proper folding of many proteins can only be achieved by interaction with molecular chaperones. The molecular chaperone UNC-45B is required for the folding of striated muscle myosin II. However, the precise mechanism by which it contributes to proper folding of the myosin head remains unclear. UNC-45B contains three domains: an N-terminal TPR domain known to bind Hsp90, a Central domain of unknown function, and a C-terminal UCS domain known to interact with the myosin head. Here we used fluorescence titrations methods, dynamic light scattering, and single-molecule atomic force microscopy (AFM) unfolding/refolding techniques to study the interactions of the UCS and Central domains with the myosin motor domain. We found that both the UCS and the Central domains bind to the myosin motor domain. Our data show that the domains bind to distinct subsites on the myosin head, suggesting distinct roles in forming the myosin−UNC-45B complex. To determine the chaperone activity of the UCS and Central domains, we used two different methods: 1), prevention of misfolding using single-molecule AFM, and 2), prevention of aggregation using dynamic light scattering. Using the first method, we found that the UCS domain is sufficient to prevent misfolding of a titin mechanical reporter. Application of the second method showed that the UCS domain but not the Central domain prevents the thermal aggregation of the myosin motor domain. We conclude that while both the UCS and the Central domains bind the myosin head with high affinity, only the UCS domain displays chaperone activity.