Cooperative Folding in a Multi-domain Protein

Thermodynamic Coupling in a Consensus-Designed Spectrin Repeat Protein

Cooperativity is a central feature to protein folding and is important for protein design. Repeat proteins are good systems for quantifying the thermodynamic basis of cooperativity. Analysis of repeat proteins composed of identical consensus repeats show that repeats strongly drive the folding of their neighbors through extensive tertiary contacts. Here, we use the consensus approach to quantify the cooperativity of folding of spectrin repeat arrays. Spectrin repeats are unique among tandem repeat proteins in that they share an elongated α-helix that spans neighboring repeats. We generate a consensus spectrin repeat sequence and show that this sequence is structured by CD and NMR spectroscopy, and is considerably more stable than extant spectrin repeats. By generating pairs of consensus spectrin repeats, we find tandem repeats to be further stabilized, demonstrating cooperative stabilization by neighboring repeats. Using an Ising model to analyze single- and tandem spectrin repeat unfolding, we find that the consensus stability increase results from intrinsic but not interfacial stabilization. By introducing mutations and insertions at the boundary between consensus repeats, we find that cooperativity is driven primarily by helical propagation; to a lesser extent, helix propagation also stabilizes partly folded states where one of two repeats is unfolded.

Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication

The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.

Inter-Domain Repulsion of Dumbbell-Shaped Calmodulin during Electrospray Ionization Revealed by Molecular Dynamics Simulations

The mechanisms whereby protein ions are released from nanodroplets at the liquid-gas interface have continued to be controversial since electrospray ionization (ESI) mass spectrometry was widely applied in biomolecular structure analysis in solution. Several viable pathways have been proposed and verified for single-domain proteins. However, the ESI mechanism of multi-domain proteins with more complicated and flexible structures remains unclear. Herein, dumbbell-shaped calmodulin was chosen as a multi-domain protein model to perform molecular dynamics simulations to investigate the structural evolution during the ESI process. For [Ca₄CAM], the protein followed the classical charge residue model. As the inter-domain electrostatic repulsion increased, the droplet was found to split into two sub-droplets, while stronger-repulsive apo-calmodulin unfolded during the early evaporation stage. We designated this novel ESI mechanism as the domain repulsion model, which provides new mechanistic insights into further exploration of proteins containing more domains. Our results suggest that greater attention should be paid to the effect of domain-domain interactions on structure retention during liquid-gas interface transfer when mass spectrometry is used as the developing technique in gas phase structural biology.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

First Principles Calculation of Protein-Protein Dimer Affinities of ALS-Associated SOD1 Mutants

Cu,Zn superoxide dismutase (SOD1) is a 32 kDa homodimer that converts toxic oxygen radicals in neurons to less harmful species. The dimerization of SOD1 is essential to the stability of the protein. Monomerization increases the likelihood of SOD1 misfolding into conformations associated with aggregation, cellular toxicity, and neuronal death in familial amyotrophic lateral sclerosis (fALS). The ubiquity of disease-associated mutations throughout the primary sequence of SOD1 suggests an important role of physicochemical processes, including monomerization of SOD1, in the pathology of the disease. Herein, we use a first-principles statistical mechanics method to systematically calculate the free energy of dimer binding for SOD1 using molecular dynamics, which involves sequentially computing conformational, orientational, and separation distance contributions to the binding free energy. We consider the effects of two ALS-associated mutations in SOD1 protein on dimer stability, A4V and D101N, as well as the role of metal binding and disulfide bond formation. We find that the penalty for dimer formation arising from the conformational entropy of disordered loops in SOD1 is significantly larger than that for other protein-protein interactions previously considered. In the case of the disulfide-reduced protein, this leads to a bound complex whose formation is energetically disfavored. Somewhat surprisingly, the loop free energy penalty upon dimerization is still significant for the holoprotein, despite the increased structural order induced by the bound metal cations. This resulted in a surprisingly modest increase in dimer binding free energy of only about 1.5 kcal/mol upon metalation of the protein, suggesting that the most significant stabilizing effects of metalation are on folding stability rather than dimer binding stability. The mutant A4V has an unstable dimer due to weakened monomer-monomer interactions, which are manifested in the calculation by a separation free energy surface with a lower barrier. The mutant D101N has a stable dimer partially due to an unusually rigid β-barrel in the free monomer. D101N also exhibits anticooperativity in loop folding upon dimerization. These computational calculations are, to our knowledge, the most quantitatively accurate calculations of dimer binding stability in SOD1 to date.

Slowest-first protein translation scheme: Structural asymmetry and co-translational folding

Proteins are translated from the N to the C terminus, raising the basic question of how this innate directionality affects their evolution. To explore this question, we analyze 16,200 structures from the Protein Data Bank (PDB). We find remarkable enrichment of α helices at the C terminus and β strands at the N terminus. Furthermore, this α-β asymmetry correlates with sequence length and contact order, both determinants of folding rate, hinting at possible links to co-translational folding (CTF). Hence, we propose the "slowest-first" scheme, whereby protein sequences evolved structural asymmetry to accelerate CTF: the slowest of the cooperatively folding segments are positioned near the N terminus so they have more time to fold during translation. A phenomenological model predicts that CTF can be accelerated by asymmetry in folding rate, up to double the rate, when folding time is commensurate with translation time; analysis of the PDB predicts that structural asymmetry is indeed maximal in this regime. This correspondence is greater in prokaryotes, which generally require faster protein production. Altogether, this indicates that accelerating CTF is a substantial evolutionary force whose interplay with stability and functionality is encoded in secondary structure asymmetry.

Analysis of Tandem Repeat Protein Folding Using Nearest-Neighbor Models

Cooperativity is a hallmark of protein folding, but the thermodynamic origins of cooperativity are difficult to quantify. Tandem repeat proteins provide a unique experimental system to quantify cooperativity due to their internal symmetry and their tolerance of deletion, extension, and in some cases fragmentation into single repeats. Analysis of repeat proteins of different lengths with nearest-neighbor Ising models provides values for repeat folding ([Formula: see text]) and inter-repeat coupling (ΔG_i-1,i). In this article, we review the architecture of repeat proteins and classify them in terms of ΔG_i and ΔG_i-1,i; this classification scheme groups repeat proteins according to their degree of cooperativity. We then present various statistical thermodynamic models, based on the 1D-Ising model, for analysis of different classes of repeat proteins. We use these models to analyze data for highly and moderately cooperative and noncooperative repeat proteins and relate their fitted parameters to overall structural features.

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.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

A Change in Domain Cooperativity Drives the Function of Calnuc

With the increasing incidence of neurodegenerative disorders, there is an urgent need to understand the protein folding process. Examining the folding process of multidomain proteins remains a prime challenge, as their complex conformational dynamics make them highly susceptible to misfolding and/or aggregation. The presence of multiple domains in a protein can lead to interaction between the partially folded domains, thereby driving misfolding and/or aggregation. Calnuc is one such multidomain protein for which Ca²⁺ binding plays a pivotal role in governing its structural dynamics and stability and, presumably, in directing its interactions with other proteins. We demonstrate differential structural dynamics between the Ca²⁺-free and Ca²⁺-bound forms of calnuc. In the absence of Ca²⁺, full-length calnuc displays equilibrium structural transitions with four intermediate states, reporting a sum of the behavioral properties of its individual domains. Fragment-based studies illustrate the sequential events of structure adoption proceeding in the following order: EF domain followed by the NT and LZ domains in the apo state. On the other hand, Ca²⁺ binding increases domain cooperativity and enables the protein to fold as a single unit. Single-tryptophan mutant proteins, designed in a domain-dependent manner, confirm an increase in the number of interdomain interactions in the Ca²⁺-bound form as compared to the Ca²⁺-free state of the protein, thereby providing insight into its folding process. The attenuated domain crosstalk in apo-calnuc is likely to influence and regulate its physiologically important intermolecular interactions.

Cotranslational folding cooperativity of contiguous domains of α-spectrin

Proteins synthesized in the cell can begin to fold during translation before the entire polypeptide has been produced, which may be particularly relevant to the folding of multidomain proteins. Here, we study the cotranslational folding of adjacent domains from the cytoskeletal protein α-spectrin using force profile analysis (FPA). Specifically, we investigate how the cotranslational folding behavior of the R15 and R16 domains are affected by their neighboring R14 and R16, and R15 and R17 domains, respectively. Our results show that the domains impact each other's folding in distinct ways that may be important for the efficient assembly of α-spectrin, and may reduce its dependence on chaperones. Furthermore, we directly relate the experimentally observed yield of full-length protein in the FPA assay to the force exerted by the folding protein in piconewtons. By combining pulse-chase experiments to measure the rate at which the arrested protein is converted into full-length protein with a Bell model of force-induced rupture, we estimate that the R16 domain exerts a maximal force on the nascent chain of ∼15 pN during cotranslational folding.

Unraveling the Mechanical Unfolding Pathways of a Multidomain Protein: Phosphoglycerate Kinase

Phosphoglycerate kinase (PGK) is a highly conserved enzyme that is crucial for glycolysis. PGK is a monomeric protein composed of two similar domains and has been the focus of many studies for investigating interdomain interactions within the native state and during folding. Previous studies used traditional biophysical methods (such as circular dichroism, tryptophan fluorescence, and NMR) to measure signals over a large ensemble of molecules, which made it difficult to observe transient changes in stability or structure during unfolding and refolding of single molecules. Here, we unfold single molecules of PGK using atomic force spectroscopy and steered molecular dynamic computer simulations to examine the conformational dynamics of PGK during its unfolding process. Our results show that after the initial forced separation of its domains, yeast PGK (yPGK) does not follow a single mechanical unfolding pathway; instead, it stochastically follows two distinct pathways: unfolding from the N-terminal domain or unfolding from the C-terminal domain. The truncated yPGK N-terminal domain unfolds via a transient intermediate, whereas the structurally similar isolated C-terminal domain has no detectable intermediates throughout its mechanical unfolding process. The N-terminal domain in the full-length yPGK displays a strong unfolding intermediate 13% of the time, whereas the truncated domain (yPGKNT) transitions through the intermediate 81% of the time. This effect indicates that the mechanical properties of yPGK cannot be simply deduced from the mechanical properties of its constituents. We also find that Escherichia coli PGK is significantly less mechanically stable as compared to yPGK, contrary to bulk unfolding measurements. Our results support the growing body of observations that the folding behavior of multidomain proteins is difficult to predict based solely on the studies of isolated domains.

Modulating long-range energetics via helix stabilization: A case study using T4 lysozyme

Cooperative protein folding requires distant regions of a protein to interact and provide mutual stabilization. The mechanism of this long-distance coupling remains poorly understood. Here, we use T4 lysozyme (T4L*) as a model to investigate long-range communications across two subdomains of a globular protein. T4L* is composed of two structurally distinct subdomains, although it behaves in a two-state manner at equilibrium. The subdomains of T4L* are connected via two topological connections: the N-terminal helix that is structurally part of the C-terminal subdomain (the A-helix) and a long helix that spans both subdomains (the C-helix). To understand the role that the C-helix plays in cooperative folding, we analyzed a circularly permuted version of T4L* (CP13*), whose subdomains are connected only by the C-helix. We demonstrate that when isolated as individual fragments, both subdomains of CP13* can fold autonomously into marginally stable conformations. The energetics of the N-terminal subdomain depend on the formation of a salt bridge known to be important for stability in the full-length protein. We show that the energetic contribution of the salt bridge to the stability of the N-terminal fragment increases when the C-helix is stabilized, such as occurs upon folding of the C-terminal subdomain. These results suggest a model where long-range energetic coupling is mediated by helix stabilization and not specific tertiary interactions.

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1 Å·ps^-1) and constant forces (700 and 900 pN) using molecular dynamics (MD) simulations at T = 300 K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.

Protein Assembly and Building Blocks: Beyond the Limits of the LEGO Brick Metaphor

Proteins, like other biomolecules, have a modular and hierarchical structure. Various building blocks are used to construct proteins of high structural complexity and diverse functionality. In multidomain proteins, for example, domains are fused to each other in different combinations to achieve different functions. Although the LEGO brick metaphor is justified as a means of simplifying the complexity of three-dimensional protein structures, several fundamental properties (such as allostery or the induced-fit mechanism) make deviation from it necessary to respect the plasticity, softness, and cross-talk that are essential to protein function. In this work, we illustrate recently reported protein behavior in multidomain proteins that deviates from the LEGO brick analogy. While earlier studies showed that a protein domain is often unaffected by being fused to another domain or becomes more stable following the formation of a new interface between the tethered domains, destabilization due to tethering has been reported for several systems. We illustrate that tethering may sometimes result in a multidomain protein behaving as "less than the sum of its parts". We survey these cases for which structure additivity does not guarantee thermodynamic additivity. Protein destabilization due to fusion to other domains may be linked in some cases to biological function and should be taken into account when designing large assemblies.

Nonspecific yet decisive: Ubiquitination can affect the native-state dynamics of the modified protein

Ubiquitination is one of the most common post-translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site-specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48-linked tetra-ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half-life through proteasomal degradation.

Nucleophosmin contains amyloidogenic regions that are able to form toxic aggregates under physiological conditions

Nucleophosmin (NPM)-1 is a multifunctional protein involved in a variety of biologic processes and has been implicated in the pathogenesis of several human malignancies. To gain insight into the role of isolated fragments in NPM1 activities, we dissected the C-terminal domain (CTD) into its helical fragments. In this study, we observed the unexpected structural behavior of the peptide fragment corresponding to helix (H)2 (residues 264-277). This peptide has a strong tendency to form amyloidlike assemblies endowed with fibrillar morphology and β-sheet structure, under physiologic conditions, as shown by circular dichroism, thioflavin T, and Congo red binding assays; dynamic light scattering; and atomic force microscopy. The aggregates are also toxic to neuroblastoma cells, as determined using 3-(4;5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide reduction and Ca(2+) influx assays. We also found that the extension of the H2 sequence beyond its N terminus, comprising the connecting loop with H1, delayed aggregation and its associated cytotoxicity, suggesting that contiguous regions of H2 have a protective role in preventing aggregation. Our findings and those in the literature suggest that the helical structures present in the CTD are important in preventing harmful aggregation. These findings could elucidate the pathogenesis of acute myeloid leukemia (AML) caused by NPM1 mutants. Because the CTD is not properly folded in these mutants, we hypothesize that the aggregation propensity of this NPM1 region is involved in the pathogenesis of AML. Preliminary assays on NPM1-Cter-MutA, the most frequent AML-CTD mutation, revealed its significant propensity for aggregation. Thus, the aggregation phenomena should be seriously considered in studies aimed at unveiling the molecular mechanisms of this pathology.

Probing conformational stability and dynamics of erythroid and nonerythroid spectrin: effects of urea and guanidine hydrochloride

We have studied the conformational stability of the two homologous membrane skeletal proteins, the erythroid and non-erythroid spectrins, in their dimeric and tetrameric forms respectively during unfolding in the presence of urea and guanidine hydrochloride (GuHCl). Fluorescence and circular dichroism (CD) spectroscopy have been used to study the changes of intrinsic tryptophan fluorescence, anisotropy, far UV-CD and extrinsic fluorescence of bound 1-anilinonapthalene-8-sulfonic acid (ANS). Chemical unfolding of both proteins were reversible and could be described as a two state transition. The folded erythroid spectrin and non-erythroid spectrin were directly converted to unfolded monomer without formation of any intermediate. Fluorescence quenching, anisotropy, ANS binding and dynamic light scattering data suggest that in presence of low concentrations of the denaturants (up-to 1M) hydrogen bonding network and van der Waals interaction play a role inducing changes in quaternary as well as tertiary structures without complete dissociation of the subunits. This is the first report of two large worm like, multi-domain proteins obeying twofold rule which is commonly found in small globular proteins. The free energy of stabilization (ΔGuH20) for the dimeric spectrin has been 20 kcal/mol lesser than the tetrameric from.

Significantly enhanced heme retention ability of myoglobin engineered to mimic the third covalent linkage by nonaxial histidine to heme (vinyl) in synechocystis hemoglobin

Heme proteins, which reversibly bind oxygen and display a particular fold originally identified in myoglobin (Mb), characterize the "hemoglobin (Hb) superfamily." The long known and widely investigated Hb superfamily, however, has been enriched by the discovery and investigation of new classes and members. Truncated Hbs typify such novel classes and exhibit a distinct two-on-two α-helical fold. The truncated Hb from the freshwater cyanobacterium Synechocystis exhibits hexacoordinate heme chemistry and bears an unusual covalent bond between the nonaxial His(117) and a heme porphyrin 2-vinyl atom, which remains tightly associated with the globin unlike any other. It seems to be the most stable Hb known to date, and His(117) is the dominant force holding the heme. Mutations of amino acid residues in the vicinity did not influence this covalent linkage. Introduction of a nonaxial His into sperm whale Mb at the topologically equivalent position and in close proximity to vinyl group significantly increased the heme stability of this prototype globin. Reversed phase chromatography, electrospray ionization-MS, and MALDI-TOF analyses confirmed the presence of covalent linkage in Mb I107H. The Mb mutant with the engineered covalent linkage was stable to denaturants and exhibited ligand binding and auto-oxidation rates similar to the wild type protein. This indeed is a novel finding and provides a new perspective to the evolution of Hbs. The successful attempt at engineering heme stability holds promise for the production of stable Hb-based blood substitute.

Mechanism of Assembly of the Non-Covalent Spectrin Tetramerization Domain from Intrinsically Disordered Partners

Interdomain interactions of spectrin are critical for maintenance of the erythrocyte cytoskeleton. In particular, “head-to-head” dimerization occurs when the intrinsically disordered C-terminal tail of β-spectrin binds the N-terminal tail of α-spectrin, folding to form the “spectrin tetramer domain”. This non-covalent three-helix bundle domain is homologous in structure and sequence to previously studied spectrin domains. We find that this tetramer domain is surprisingly kinetically stable. Using a protein engineering Φ-value analysis to probe the mechanism of formation of this tetramer domain, we infer that the domain folds by the docking of the intrinsically disordered β-spectrin tail onto the more structured α-spectrin tail.

Understanding pathogenic single-nucleotide polymorphisms in multidomain proteins--studies of isolated domains are not enough

Studying the effects of pathogenic mutations is more complex in multidomain proteins when compared with single domains: mutations occurring at domain boundaries may have a large effect on a neighbouring domain that will not be detected in a single-domain system. To demonstrate this, we present a study that utilizes well-characterized model protein domains from human spectrin to investigate the effect of disease-and non-disease-causing single point mutations occurring at the boundaries of human spectrin repeats. Our results show that mutations in the single domains have no clear correlation with stability and disease; however, when studied in a tandem model system, the disease-causing mutations are shown to disrupt stabilizing interactions that exist between domains. This results in a much larger decrease in stability than would otherwise have been predicted, and demonstrates the importance of studying such mutations in the correct protein context.

Dynamic control of electron transfers in diflavin reductases

Diflavin reductases are essential proteins capable of splitting the two-electron flux from reduced pyridine nucleotides to a variety of one electron acceptors. The primary sequence of diflavin reductases shows a conserved domain organization harboring two catalytic domains bound to the FAD and FMN flavins sandwiched by one or several non-catalytic domains. The catalytic domains are analogous to existing globular proteins: the FMN domain is analogous to flavodoxins while the FAD domain resembles ferredoxin reductases. The first structural determination of one member of the diflavin reductases family raised some questions about the architecture of the enzyme during catalysis: both FMN and FAD were in perfect position for interflavin transfers but the steric hindrance of the FAD domain rapidly prompted more complex hypotheses on the possible mechanisms for the electron transfer from FMN to external acceptors. Hypotheses of domain reorganization during catalysis in the context of the different members of this family were given by many groups during the past twenty years. This review will address the recent advances in various structural approaches that have highlighted specific dynamic features of diflavin reductases.

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.

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.

Stabilization of the spectrin-like domains of nesprin-1α by the evolutionarily conserved "adaptive" domain

Nesprins are located at the outer and inner membranes of the nuclear envelope and help link the cytoskeleton to the nucleoskeleton. Nesprin-1α, located at the inner nuclear membrane, binds to A-type lamins and emerin and has homology to spectrin-repeat proteins. However, the mechanical and thermodynamic properties of the spectrin-like repeats (SLRs) of nesprin-1α and the potential structural contributions of the unique central domain were untested. In other spectrin superfamily proteins, tandem spectrin-repeat domains undergo cooperatively coupled folding and unfolding. We hypothesized that the large central domain, which interrupts SLRs and is conserved in other nesprin isoforms, might confer unique structural properties. To test this model we measured the thermal unfolding of nesprin-1α fragments using circular dichroism and dynamic light scattering. The SLRs in nesprin-1α were found to have structural and thermodynamic properties typical of spectrins. The central domain had relatively little secondary structure as an isolated fragment, but significantly stabilized larger SLR-containing molecules by increasing their overall helicity, thermal stability and cooperativity of folding. We suggest this domain, now termed the 'adaptive' domain (AD), also strengthens dimerization and inhibits unfolding. Further engineering of the isolated AD, and AD-containing nesprin molecules, may yield new information about the higher-order association of cooperative protein motifs.

The HD-exchange motions of ribosomal protein S6 are insensitive to reversal of the protein-folding pathway

An increasing number of protein structures are found to encompass multiple folding nuclei, allowing their structures to be formed by several competing pathways. A typical example is the ribosomal protein S6, which comprises two folding nuclei (sigma1 and sigma2) defining two competing pathways in the folding energy landscape: sigma1 --> sigma2 and sigma2 --> sigma1. The balance between the two pathways, and thus the order of folding events, is easily controlled by circular permutation. In this study, we make use of this ability to manipulate the folding pathway to demonstrate that the dynamic motions of the S6 structure are independent of how the protein folds. The HD-exchange protection factors remain the same upon complete reversal of the folding order. The phenomenon arises because the HD-exchange motions and the high-energy excitations controlling the folding pathway occur at separated free-energy levels: the Boltzmann distribution of unproductive unfolding attempts samples all unfolding channels in parallel, even those that end up in excessively high barriers. Accordingly, the findings provide a simple rationale for how to interpret native-state dynamics without the need to invoke fluctuations off the normal unfolding reaction coordinate.

Reversible unfolding of dimeric phosphofructokinase-2 from Escherichia coli reveals a dominant role of inter-subunit contacts for stability

Escherichia coli phosphofructokinase-2 (Pfk-2) is a homodimer whose subunits consist of a large domain and an additional beta-sheet that provides the interfacial contacts between the subunits, creating a beta-barrel flattened-like structure with the adjacent subunit's beta-sheet. To determine how the structural organization of Pfk-2 determines its stability, the reversible unfolding of the enzyme was characterized under equilibrium conditions by enzymatic activity, circular dichroism, fluorescence and hydrodynamic measurements. Pfk-2 undergoes a cooperative unfolding/dissociation process with the accumulation of an expanded and unstructured monomeric intermediate with a marginal stability and a large solvent accessibility with respect to the native dimer.

A Two-amino Acid Mutation Encountered in Duchenne Muscular Dystrophy Decreases Stability of the Rod Domain 23 (R23) Spectrin-like Repeat of Dystrophin

Lack of functional dystrophin causes severe Duchenne muscular dystrophy. The subsarcolemmal location of dystrophin, as well as its association with both cytoskeleton and membrane, suggests a role in the mechanical regulation of muscular membrane stress. In particular, phenotype rescue in a Duchenne muscular dystrophy mice model has shown that some parts of the central rod domain of dystrophin, constituted by 24 spectrin-like repeats, are essential. In this study, we made use of rare missense pathogenic mutations in the dystrophin gene and analyzed the biochemical properties of the isolated repeat 23 bearing single or double mutations E2910V and N2912D found in muscle dystrophy with severity grading. No dramatic effect on secondary and tertiary structure of the repeat was found in mutants compared with wild type as revealed by circular dichroism and NMR. Thermal and chemical unfolding data from circular dichroism and tryptophan fluorescence show significant decrease of stability for the mutants, and stopped-flow spectroscopy shows decreased refolding rates. The most deleterious single mutation is the N2912D replacement, although we observe additive effects of the two mutations on repeat stability. Based on three-dimensional structures built by homology molecular modeling, we discuss the modifications of the mutation-induced repeat stability. We conclude that the main forces involved in repeat stability are electrostatic inter-helix interactions that are disrupted following mutations. This study represents the first analysis at the protein level of the consequences of missense mutations in the human dystrophin rod domain. Our results suggest that it may participate in mechanical weakening of dystrophin-deficient muscle.

Studying the folding of multidomain proteins

There have been relatively few detailed comprehensive studies of the folding of protein domains (or modules) in the context of their natural covalently linked neighbors. This is despite the fact that a significant proportion of the proteome consists of multidomain proteins. In this review we highlight some key experimental investigations of the folding of multidomain proteins to draw attention to the difficulties that can arise in analyzing such systems. The evidence suggests that interdomain interactions can significantly affect stability, folding, and unfolding rates. However, preliminary studies suggest that folding pathways are unaffected-to this extent domains can be truly considered to be independent folding units. Nonetheless, it is clear that interactions between domains cannot be ignored, in particular when considering the effects of mutations.

Rational design of a conformation-switchable Ca2+- and Tb3+-binding protein without the use of multiple coupled metal-binding sites

Ca2+, as a messenger of signal transduction, regulates numerous target molecules via Ca2+-induced conformational changes. Investigation into the determinants for Ca2+-induced conformational change is often impeded by cooperativity between multiple metal-binding sites or protein oligomerization in naturally occurring proteins. To dissect the relative contributions of key determinants for Ca2+-dependent conformational changes, we report the design of a single-site Ca2+-binding protein (CD2.trigger) created by altering charged residues at an electrostatically sensitive location on the surface of the host protein rat Cluster of Differentiation 2 (CD2).CD2.trigger binds to Tb3+ and Ca2+ with dissociation constants of 0.3 +/- 0.1 and 90 +/- 25 microM, respectively. This protein is largely unfolded in the absence of metal ions at physiological pH, but Tb3+ or Ca2+ binding results in folding of the native-like conformation. Neutralization of the charged coordination residues, either by mutation or protonation, similarly induces folding of the protein. The control of a major conformational change by a single Ca2+ ion, achieved on a protein designed without reliance on sequence similarity to known Ca2+-dependent proteins and coupled metal-binding sites, represents an important step in the design of trigger proteins.

Cooperativity, connectivity, and folding pathways of multidomain proteins

Multidomain proteins are ubiquitous in both prokaryotic and eukaryotic proteomes. Study on protein folding, however, has concentrated more on the isolated single domains of proteins, and there have been relatively few systematic studies on the effects of domain-domain interactions on folding. We here discuss this issue by examining human gammaD-crystallin, spore coat protein S, and a tandem array of the R16 and R17 domains of spectrin as example proteins by using a structure-based model of folding. The calculated results consistently explain the experimental data on folding pathways and effects of mutational perturbations, supporting the view that the connectivity of two domains and the distribution of domain-domain interactions in the native conformation are factors to determine kinetic and equilibrium properties of cooperative folding.

Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.

Navigating the RNA folding landscape

Analysis of individual RNA folding reactions reveals that, as in proteins, cooperative interactions selectively drive RNA toward its biologically active, native conformation. This new work establishes a platform for future investigations of the physical principles underlying the assembly of large RNA enzymes.

Direct measurement of tertiary contact cooperativity in RNA folding

All structured biological macromolecules must overcome the thermodynamic folding problem to populate a unique functional state among a vast ensemble of unfolded and alternate conformations. The exploration of cooperativity in protein folding has helped reveal and distinguish the underlying mechanistic solutions to this folding problem. Analogous dissections of RNA tertiary stability remain elusive, however, despite the central biological importance of folded RNA molecules and the potential to reveal fundamental properties of structured macromolecules via comparisons of protein and RNA folding. We report a direct quantitative measure of tertiary contact cooperativity in a folded RNA. We precisely measured the stability of an independently folding P4-P6 domain from the Tetrahymena thermophila group I intron by single molecule fluorescence resonance energy transfer (smFRET). Using wild-type and mutant RNAs, we found that cooperativity between the two tertiary contacts enhances P4-P6 stability by 3.2 +/- 0.2 kcal/mol.

The folding pathway of a single domain in a multidomain protein is not affected by its neighbouring domain

Domains are the structural, functional, and evolutionary components of proteins. Most folding studies to date have concentrated on the folding of single domains, but more than 70% of human proteins contain more than one domain, and interdomain interactions can affect both the stability and the folding kinetics. Whether the folding pathway is altered by interdomain interactions is not yet known. Here we investigated the effect of a folded neighbouring domain on the folding pathway of spectrin R16 (the 16th α-helical repeat from chicken brain α-spectrin) by using the two-domain construct R1516. The R16 folds faster and unfolds more slowly in the presence of its folded neighbour R15 (the 15th α-helical repeat from chicken brain α-spectrin). An extensive Φ-value analysis of the R16 domain in R1516 was completed to compare the transition state of the R16 domain alone with that of the R16 domain in a multidomain construct. The results indicate that the folding pathways are the same. This result validates the current approach of breaking up larger proteins into domains for the study of protein folding pathways.

Mechanical biochemistry of proteins one molecule at a time

The activity of proteins and their complexes often involves the conversion of chemical energy (stored or supplied) into mechanical work through conformational changes. Mechanical forces are also crucial for the regulation of the structure and function of cells and tissues. Thus, the shape of eukaryotic cells (and by extension, that of the multicellular organisms they form) is the result of cycles of mechanosensing, mechanotransduction, and mechanoresponse. Recently developed single-molecule atomic force microscopy techniques can be used to manipulate single molecules, both in real time and under physiological conditions, and are ideally suited to directly quantify the forces involved in both intra- and intermolecular protein interactions. In combination with molecular biology and computer simulations, these techniques have been applied to characterize the unfolding and refolding reactions in a variety of proteins. Single-molecule mechanical techniques are providing fundamental information on the structure and function of proteins and are becoming an indispensable tool to understand how these molecules fold and work.

Nonassembled human chorionic gonadotropin subunits and alphaalpha-homodimers use fast-track processing in the secretory pathway in contrast to alphabeta-heterodimers

In multimeric glycoproteins, like glycoprotein hormones, mutual subunit interactions are required for correct folding, assembly, and transport in the secretory pathway. However, character and time course of these interactions need further elucidation. The influence of the glycoprotein hormone alpha-subunit (GPHalpha) on the folding of the human chorionic gonadotropin (hCG) beta-subunit (hCGbeta) in hCG alphabeta-heterodimers was investigated in [(35)S]Met/Cys-labeled JEG-3 cells. Completeness of disulfide bridge formation during the time course of folding was estimated by labeling with [(3)H]N-ethylmaleinimide of free thiol groups not yet consumed. Subunit association took place between immature hCGbeta (high (3)H/(35)S ratio) and almost completely folded GPHalpha. Analysis revealed a highly dynamic maturation process comprising of at least eight main hCGbeta folding intermediates (molecular masses from 107 to 28 kDa) that could be micro-preparatively isolated and characterized. These hCGbeta variants developed while being associated with GPHalpha. The 107-kDa variant was identified as a complex with calnexin. In contrast to hCG alphabeta-heterodimers, free nonassociated hCGbeta, free large GPHalpha, and GPHalphaalpha homodimers showed a fast-track-like processing in the secretory pathway. At 10 min before hCG secretion, sialylation of these variants had already been completed in the late Golgi, whereas hCG alphabeta-heterodimers had still not arrived medial Golgi. This shows that the GPHalpha in the hCG alphabeta-heterodimers decelerates the maturation of the hCGbeta portion in the heterodimer complex. This results in a postponed approval of hCG alphabeta-heterodimers by the endoplasmic reticulum quality control unlike GPHalphaalpha homodimers, free hCGbeta, and GPHalpha subunits.

Distinguishing specific and nonspecific interdomain interactions in multidomain proteins

Multidomain proteins account for over two-thirds of the eukaryotic genome. Although there have been extensive studies into the biophysical properties of isolated domains, few have investigated how the domains interact. Spectrin is a well-characterized multidomain protein with domains linked in tandem array by contiguous helices. Several of these domains have been shown to be stabilized by their neighbors. Until now, this stabilization has been attributed to specific interactions between the natural neighbors, however we have recently observed that nonnatural neighboring domains can also induce a significant amount of stabilization. Here we investigate this nonnative stabilizing effect. We created spectrin-titin domain pairs of both spectrin R16 and R17 with a single titin I27 domain at either the N- or the C-terminus and found that spectrin domains are significantly stabilized, through slowed unfolding, by nonnative interactions at the C-terminus only. Of particular importance, we show that specific interactions between natural folded neighbors at either terminus confer even greater stability by additionally increasing the folding rate constants. We demonstrate that it is possible to distinguish between natural stabilizing interactions and nonspecific stabilizing effects through examination of the kinetics of well chosen mutant proteins. This work adds to the complexity of studying multidomain proteins.

The folding and evolution of multidomain proteins

Analyses of genomes show that more than 70% of eukaryotic proteins are composed of multiple domains. However, most studies of protein folding focus on individual domains and do not consider how interactions between domains might affect folding. Here, we address this by analysing the three-dimensional structures of multidomain proteins that have been characterized experimentally and observe that where the interface is small and loosely packed, or unstructured, the folding of the domains is independent. Furthermore, recent studies indicate that multidomain proteins have evolved mechanisms to minimize the problems of interdomain misfolding.

Mechanical unfolding of proteins: insights into biology, structure and folding

Since user-friendly atomic force microscopes came onto the market a few years ago, scientists have explored the response of many proteins to applied force. This field has now matured beyond the phenomenological with exciting recent developments, particularly with regards to research into biological questions. For example, detailed mechanistic studies have suggested how mechanically active proteins perform their functions. Also, in vitro forced unfolding has been compared with in vivo protein import and degradation. Additionally, investigations have been carried out that probe the relationship between protein structure and response to applied force, an area that has benefited significantly from synergy between experiments and simulations. Finally, recent technological developments offer exciting new avenues for experimental studies.

Malleability of protein folding pathways: a simple reason for complex behaviour

Although the structures of native proteins are generally unique, the pathways by which they form are often free to vary. Some proteins fold by a multitude of different pathways, whereas others seem restricted to only one choice. An explanation for this variation in folding behaviour has recently emerged from studies of transition state changes: the number of accessible pathways is linked to the number of nucleation motifs contained within the native topology. We refer to these nucleation motifs as 'foldons', as they approach the size of an independent cooperative unit. Thus, with respect to pathway malleability and the composition of the folding funnel, proteins can be seen as modular assemblies of competing foldons. For the split beta-alpha-beta fold, these foldons are two-strand-helix motifs coupled by spatial overlap.

Spectrin organization and dynamics: new insights

Spectrin is the major constituent protein of the erythrocyte cytoskeleton which forms a filamentous network on the cytoplasmic face of the membrane by providing a scaffold for a variety of proteins. In this review, several aspects of spectrin organization are highlighted, particularly with respect to its ability to bind hydrophobic ligands and its interaction with membrane surfaces. The characteristic binding of the fluorescent hydrophobic probes Prodan and pyrene to spectrin, which allows an estimation of the polarity of the hydrophobic probe binding site, is illustrated. In addition, the contribution of uniquely localized and conserved tryptophan residues in the 'spectrin repeats' in these processes is discussed. A functional implication of the presence of hydrophobic binding sites in spectrin is its recently discovered chaperone-like activity. Interestingly, spectrin exhibits residual structural integrity even after denaturation which could be considered as a hallmark of cytoskeletal proteins. Future research could provide useful information about the possible role played by spectrin in cellular physiology in healthy and diseased states.

Thermal stabilities of brain spectrin and the constituent repeats of subunits

The different genes that encode mammalian spectrins give rise to proteins differing in their apparent stiffness. To explore this, we have compared the thermal stabilities of the structural repeats of brain spectrin subunits (alphaII and betaII) with those of erythrocyte spectrin (alphaI and betaI). The unfolding transition midpoints (T(m)) of the 36 alphaII- and betaII-spectrin repeats extend between 24 and 82 degrees C, with an average higher by some 10 degrees C than that of the alphaI- and betaI-spectrin repeats. This difference is reflected in the T(m) values of the intact brain and erythrocyte spectrins. Two of three tandem-repeat constructs from brain spectrin exhibited strong cooperative coupling, with elevation of the T(m) of the less stable partner corresponding to coupling free energies of approximately -4.4 and -3.5 kcal/mol. The third tandem-repeat construct, by contrast, exhibited negligible cooperativity. Tandem-repeat mutants, in which a part of the "linker" helix that connects the two domains was replaced with a corresponding helical segment from erythroid spectrin, showed only minor perturbation of the thermal melting profiles, without breakdown of cooperativity. Thus, the linker regions, which tolerate few point mutations without loss of cooperative function, have evidently evolved to permit conformational coupling in specified regions. The greater structural stability of the repeats in alphaII- and betaII-spectrin may account, at least in part, for the higher rigidity of brain compared to erythrocyte spectrin.

Spectrin domains lose cooperativity in forced unfolding

Spectrin is a multidomain cytoskeletal protein, the component three-helix bundle domains are expected to experience mechanical force in vivo. In thermodynamic and kinetic studies, neighboring domains of chicken brain alpha-spectrin R16 and R17 have been shown to behave cooperatively. Is this cooperativity maintained under force? The effect of force on these spectrin domains was investigated using atomic force microscopy. The response of the individual domains to force was compared to that of the tandem repeat R1617. Importantly, nonhelical linkers (all-beta immunoglobulin domains) were used to avoid formation of nonnative helical linkers. We show that, in contrast to previous studies on spectrin repeats, only 3% of R1617 unfolding events gave an increase in contour length consistent with cooperative two-domain unfolding events. Furthermore, the unfolding forces for R1617 were the same as those for the unfolding of R16 or R17 alone. This is a strong indication that the cooperative unfolding behavior observed in the stopped-flow studies is absent between these spectrin domains when force is acting as a denaturant. Our evidence suggests that the rare double unfolding events result from misfolding between adjacent repeats. We suggest that this switch from cooperative to independent behavior allows multidomain proteins to maintain integrity under applied force.