Denaturation of human Cu/Zn superoxide dismutase by guanidine hydrochloride: a dynamic fluorescence study

Insight Into Factors Influencing the Aggregation Process in Wild-Type and P66R Mutant SOD1: Computational and Spectroscopic Approaches

Disturbances in metal ion homeostasis associated with amyotrophic lateral sclerosis (ALS) have been described for several years, but the exact mechanism of involvement is not well understood. To elucidate the role of metalation in superoxide dismutase (SOD1) misfolding and aggregation, we comprehensively characterized the structural features (apo/holo forms) of WT-SOD1 and P66R mutant in loop IV. Using computational and experimental methodologies, we assessed the physicochemical properties of these variants and their correlation with protein aggregation at the molecular level. Modifications in apo-SOD1 compared to holo-SOD1 were more pronounced in flexibility, stability, hydrophobicity, and intramolecular interactions, as indicated by molecular dynamics simulations. The enzymatic activities of holo/apo-WT SOD1 were 1.30 and 1.88-fold of the holo/apo P66R mutant, respectively. Under amyloid-inducing conditions, decreased ANS fluorescence intensity in the apo-form relative to the holo-form suggested pre-fibrillar species and amyloid aggregate growth due to occluded hydrophobic pockets. FTIR spectroscopy revealed that apo-WT-SOD1 and apo-P66R exhibited a mixture of parallel and intermolecular β-sheet structures, indicative of aggregation propensity. Aggregate species were identified using TEM, Congo red staining, and ThT/ANS fluorescence spectroscopy. Thermodynamic analyses with GdnHCl demonstrated that metal deficit, mutation, and intramolecular disulfide bond reduction are essential for initiating SOD1 misfolding and aggregation. These disruptions destabilize the dimer-monomer equilibrium, promoting dimer dissociation into monomers and decreasing the thermodynamic stability of SOD1 variants, thus facilitating amyloid/amorphous aggregate formation. Our findings offer novel insights into protein aggregation mechanisms in disease pathology and highlight potential therapeutic strategies against toxic protein aggregation, including SOD1.

Mutation/metal deficiency in the "electrostatic loop" enhanced aggregation process in apo/holo SOD1 variants: implications for ALS diseases

Despite the many mechanisms it has created to prevent unfolding and aggregation of proteins, many diseases are caused by abnormal folding of proteins, which are called misfolding diseases. During this process, proteins undergo structural changes and become stable, insoluble beta-sheet aggregates called amyloid fibrils. Mutations/disruptions in metal ion homeostasis in the ALS-associated metalloenzyme superoxide dismutase (SOD1) reduce conformational stability, consistent with the protein aggregation hypothesis for neurodegenerative diseases. However, the exact mechanism of involvement is not well understood. Hence, to understand the role of mutation/ metal deficiency in SOD1 misfolding and aggregation, we investigated the effects of apo/holo SOD1 variants on structural properties using biophysical/experimental techniques. The MD results support the idea that the mutation/metal deficiency can lead to a change in conformation. The increased content of β-sheet structures in apo/holo SOD1 variants can be attributed to the aggregation tendency, which was confirmed by FTIR spectroscopy and dictionary of secondary structure in proteins (DSSP) results. Thermodynamic studies of GdnHCl showed that metal deficiency/mutation/intramolecular S-S reduction together are required to initiate misfolding/aggregation of SOD1. The results showed that apo/holo SOD1 variants under destabilizing conditions induced amyloid aggregates at physiological pH, which were detected by ThT/ANS fluorescence, as well as further confirmation of amyloid/amorphous species by TEM. This study confirms that mutations in the electrostatic loop of SOD1 lead to structural abnormalities, including changes in hydrophobicity, reduced disulfide bonds, and an increased propensity for protein denaturation. This process facilitates the formation of amyloid/amorphous aggregates ALS-associated.

Structural motifs in the early metallation steps of Zn(II) and Cd(II) binding to apo-metallothionein 1a

Many proteins require a metal cofactor to function and these metals are often involved in the protein folding process. The protein metallothionein (MT) has a dynamic structure capable of binding to a variety of metals with different stoichiometries. The most well-understood structure is the seven-metal, two domain structure formed upon metallation using Zn(II) or Cd(II). However, the partially metallated states and the pathways to form these clusters are less well-understood, although it is known that the pathways are pH dependent. Using stopped flow methods, it is shown that the metallation rates of the less cooperative Zn(II) binding pathway is much more impacted by low pH conditions that that of the more cooperative Cd(II) binding pathway. Electrospray ionization mass spectrometry (ESI-MS) methods reveal specific mixtures of bridging and terminally bound MxSy structures form in the first few metallation steps. Using a combination of methods, the data show that the result of unfolding this intrinsically disordered apo-MT structure using guanidinium chloride is that the formation of preliminary bridging structures that form in the first few metallation steps is impeded. The data show that more terminally bound structures form. Our conclusion is that the compact conformation of the native apo-MT at physiological pH allows for rapid formation of complex metal-thiolate structures with high affinity that provides protection from oxidation, a function that is suppressed upon unfolding. Overall, these results highlight both the importance of the apo-MT structure in the metallation pathway, but also the differences in Zn(II) and Cd(II) binding under different conditions.

CuATSM improves motor function and extends survival but is not tolerated at a high dose in SOD1G93A mice with a C57BL/6 background

The synthetic copper-containing compound, CuATSM, has emerged as one of the most promising drug candidates developed for the treatment of amyotrophic lateral sclerosis (ALS). Multiple studies have reported CuATSM treatment provides therapeutic efficacy in various mouse models of ALS without any observable adverse effects. Moreover, recent results from an open label clinical study suggested that daily oral dosing with CuATSM slows disease progression in patients with both sporadic and familial ALS, providing encouraging support for CuATSM in the treatment of ALS. Here, we assessed CuATSM in high copy SOD1G93A mice on the congenic C57BL/6 background, treating at 100 mg/kg/day by gavage, starting at 70 days of age. This dose in this specific model has not been assessed previously. Unexpectedly, we report a subset of mice initially administered CuATSM exhibited signs of clinical toxicity, that necessitated euthanasia in extremis after 3-51 days of treatment. Following a 1-week washout period, the remaining mice resumed treatment at the reduced dose of 60 mg/kg/day. At this revised dose, treatment with CuATSM slowed disease progression and increased survival relative to vehicle-treated littermates. This work provides the first evidence that CuATSM produces positive disease-modifying outcomes in high copy SOD1G93A mice on a congenic C57BL/6 background. Furthermore, results from the 100 mg/kg/day phase of the study support dose escalation determination of tolerability as a prudent step when assessing treatments in previously unassessed models or genetic backgrounds.

Utilizing Ion Mobility-Mass Spectrometry to Investigate the Unfolding Pathway of Cu/Zn Superoxide Dismutase

Native mass spectrometry has emerged as a powerful tool for structural biology as it enables the evaluation of molecules as they occur in their physiological conditions. Ion mobility spectrometry-mass spectrometry (IMS-MS) has shown essential in these analyses as it allows the measurement of the shape of a molecule, denoted as its collision cross section (CCS), and mass. The structural information garnered from native IMS-MS provides insight into the tertiary and quaternary structure of proteins and can be used to validate NMR or crystallographic X-ray structures. Additionally, due to the rapid nature (millisecond measurements) and ability of IMS-MS to analyze heterogeneous solutions, it can be used to address structural questions not possible with traditional structural approaches. Herein, we applied multiple solution conditions to systematically denature bovine Cu/Zn-superoxide dismutase (SOD1) and assess its unfolding pathway from the holo-dimer to the holo-monomer, single-metal monomer, and apo-monomer. Additionally, we compared and noted 1–2% agreement between CCS values from both drift tube IMS and trapped IMS for the SOD1 holo-monomer and holo-dimer. The observed CCS values were in excellent agreement with computational CCS values predicted from the homo-dimer crystal structure, showcasing the ability to use both IMS-MS platforms to provide valuable structural information for molecular modeling of protein interactions and structural assessments.

Chaperonin-Based Biolayer Interferometry To Assess the Kinetic Stability of Metastable, Aggregation-Prone Proteins

Stabilizing the folded state of metastable and/or aggregation-prone proteins through exogenous ligand binding is an appealing strategy for decreasing disease pathologies caused by protein folding defects or deleterious kinetic transitions. Current methods of examining binding of a ligand to these marginally stable native states are limited because protein aggregation typically interferes with analysis. Here, we describe a rapid method for assessing the kinetic stability of folded proteins and monitoring the effects of ligand stabilization for both intrinsically stable proteins (monomers, oligomers, and multidomain proteins) and metastable proteins (e.g., low Tm) that uses a new GroEL chaperonin-based biolayer interferometry (BLI) denaturant pulse platform. A kinetically controlled denaturation isotherm is generated by exposing a target protein, immobilized on a BLI biosensor, to increasing denaturant concentrations (urea or GuHCl) in a pulsatile manner to induce partial or complete unfolding of the attached protein population. Following the rapid removal of the denaturant, the extent of hydrophobic unfolded/partially folded species that remains is detected by an increased level of GroEL binding. Because this kinetic denaturant pulse is brief, the amplitude of binding of GroEL to the immobilized protein depends on the duration of the exposure to the denaturant, the concentration of the denaturant, wash times, and the underlying protein unfolding-refolding kinetics; fixing all other parameters and plotting the GroEL binding amplitude versus denaturant pulse concentration result in a kinetically controlled denaturation isotherm. When folding osmolytes or stabilizing ligands are added to the immobilized target proteins before and during the denaturant pulse, the diminished population of unfolded/partially folded protein manifests as a decreased level of GroEL binding and/or a marked shift in these kinetically controlled denaturation profiles to higher denaturant concentrations. This particular platform approach can be used to identify small molecules and/or solution conditions that can stabilize or destabilize thermally stable proteins, multidomain proteins, oligomeric proteins, and, most importantly, aggregation-prone metastable proteins.

Electrostatic effect of the ribosomal surface on nascent polypeptide dynamics

The crucial molecular events accompanying protein folding in the cell are still largely unexplored. As nascent polypeptides emerge from the ribosomal exit tunnel, they come in close proximity with the highly negatively charged ribosomal surface. How is the nascent polypeptide influenced by the ribosomal surface? We address this question via the intrinsically disordered protein PIR and a number of its variably charged mutants. Two different populations are identified: one is highly spatially biased, and the other is highly dynamic. The more negatively charged nascent polypeptides emerging from the ribosome are richer in the extremely dynamic population. Hence, nascent proteins with a net negative charge are less likely to interact with the ribosome. Surprisingly, the amplitude of the local motions of the highly dynamic population is much wider than that of disordered polypeptides under physiological conditions, implying that proximity to the ribosomal surface enhances the molecular flexibility of a subpopulation of the nascent protein, much like a denaturing agent would. This effect could be important for a proper structural channeling of the nascent protein and the prevention of cotranslational kinetic trapping. Interestingly, a significant population of the highly spatially biased nascent chain, probably interacting extensively with the ribosome, is present even for very negatively charged nascent proteins. This "sticking" effect likely serves to protect nascent proteins (e.g., from cotranslational aggregation). In all, our results highlight the influence of the ribosome in nascent protein dynamics and show that the ribosome's function in protein biogenesis extends well beyond catalysis of peptide bond formation.

Slow histidine H/D exchange protocol for thermodynamic analysis of protein folding and stability using mass spectrometry

Described here is a mass spectrometry-based protocol to study the thermodynamic stability of proteins and protein-ligand complexes using the chemical denaturant dependence of the slow H/D exchange reaction of the imidazole C(2) proton in histidine side chains. The protocol is developed using several model protein systems including: ribonuclease (Rnase) A, myoglobin, bovine carbonic anhydrase (BCA) II, hemoglobin (Hb), and the hemoglobin-haptoglobin (Hb-Hp) protein complex. Folding free energies consistent with those previously determined by other more conventional techniques were obtained for the two-state folding proteins, Rnase A and myoglobin. The protocol successfully detected a previously observed partially unfolded intermediate stabilized in the BCA II folding/unfolding reaction, and it could be used to generate a K(d) value of 0.24 nM for the Hb-Hp complex. The compatibility of the protocol with conventional mass spectrometry-based proteomic sample preparation and analysis methods was also demonstrated in an experiment in which the protocol was used to detect the binding of zinc to superoxide dismutase in the yeast cell lysate sample. The yeast cell sample analyses also helped define the scope of the technique, which requires the presence of globally protected histidine residues in a protein's three-dimensional structure for successful application.

FALS mutations in Cu, Zn superoxide dismutase destabilize the dimer and increase dimer dissociation propensity: a large-scale thermodynamic analysis

Mutations in the dimeric enzyme Cu, Zn superoxide dismutase (SOD1) leading to its aggregation are implicated in the toxicity in familial amyotrophic lateral sclerosis (FALS). We and others have previously shown that aggregation occurs by a pathway involving dimer dissociation, metal-loss from monomers and multimeric assembly of apo-SOD1 monomers. We postulate that FALS mutations cause enhanced aggregation by affecting one or more steps in the pathway, and computationally test this postulate for 75 known mis-sense FALS mutants of SOD1. Based on an extensive thermodynamic analysis of the stability of apo-dimer and apo-monomer forms of these mutants, we classify the mutations into the following groups: 70 out of 75 mutations in SOD1 lead to (i) decreased dimer stability, and/or (ii) increased dimer dissociation, compared to wild type, and four mutations lead to (iii) decreased monomer stability compared to wild type. Our results suggest that enhanced aggregation of SOD1 in FALS occurs due to an increased population of mutant SOD1 apo-monomers compared to wild type. The dissociation of multimeric proteins induced by diverse mutations may be a common theme in several human diseases.

Structure, folding, and misfolding of Cu,Zn superoxide dismutase in amyotrophic lateral sclerosis

Fourteen years after the discovery that mutations in Cu, Zn superoxide dismutase (SOD1) cause a subset of familial amyotrophic lateral sclerosis (fALS), the mechanism by which mutant SOD1 exerts toxicity remains unknown. The two principle hypotheses are (a) oxidative damage stemming from aberrant SOD1 redox chemistry, and (b) misfolding of the mutant protein. Here we review the structure and function of wild-type SOD1, as well as the changes to the structure and function in mutant SOD1. The relative merits of the two hypotheses are compared and a common unifying principle is outlined. Lastly, the potential for therapies targeting SOD1 misfolding is discussed.

Mapping the folding free energy surface for metal-free human Cu,Zn superoxide dismutase

Mutations at many different sites in the gene encoding human Cu,Zn superoxide dismutase (SOD) are known to be causative agents in amyotrophic lateral sclerosis (ALS). One explanation for the molecular basis of this pathology is the aggregation of marginally soluble, partially structured states whose populations are enhanced in the protein variants. As a benchmark for testing this hypothesis, the equilibrium and kinetic properties of the reversible folding reaction of a metal-free variant of SOD were investigated. Reversibility was achieved by replacing the two non-essential cysteine residues with non-oxidizable analogs, C6A/C111S, to produce apo-AS-SOD. The metal-free pseudo-wild-type protein is folded and dimeric in the absence of chemical denaturants, and its equilibrium folding behavior is well described by an apparent two-state mechanism involving the unfolded monomer and the native dimer. The apparent free energy of folding in the absence of denaturant and at standard state is -20.37(+/- 1.04) kcal (mol dimer)(-1). A global analysis of circular dichroism kinetic traces for both unfolding and refolding reactions, combined with results from small angle X-ray scattering and time-resolved fluorescence anisotropy measurements, supports a sequential mechanism involving the unfolded monomer, a folded monomeric intermediate, and the native dimer. The rate-limiting monomer folding reaction is followed by a near diffusion-limited self-association reaction to form the native dimer. The relative population of the folded monomeric intermediate is predicted not to exceed 0.5% at micromolar concentrations of protein under equilibrium and both strongly unfolding and refolding conditions for metal-free pseudo-wild-type SOD.

Folding of Cu/Zn superoxide dismutase suggests structural hotspots for gain of neurotoxic function in ALS: parallels to precursors in amyloid disease

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease linked to misfolding of the ubiquitous enzyme Cu/Zn superoxide dismutase (SOD). In contrast to other protein-misfolding disorders with similar neuropathogenesis, ALS is not always associated with the in vivo deposition of protein aggregates. Thus, under the assumption that all protein-misfolding disorders share at primary level a similar disease mechanism, ALS constitutes an interesting disease model for identifying the yet-mysterious precursor states from which the cytotoxic pathway emerges. In this study, we have mapped out the conformational repertoire of the apoSOD monomer through analysis of its folding behavior. The results allow us to target the regions of the SOD structure that are most susceptible to unfolding locally under physiological conditions, leading to the exposure of structurally promiscuous interfaces that are normally hidden in the protein's interior. The structure of this putative ALS precursor is strikingly similar to those implicated in amyloid disease.

Mechanism and Thermodynamics of Guanidinium Chloride-induced Denaturation of ALS-associated Mutant Cu,Zn Superoxide Dismutases

Mutations in human copper zinc superoxide dismutase (hSOD) that are associated with amyotrophic lateral sclerosis (ALS) have been proposed to destabilize the protein and thereby enhance toxic protein aggregation. In previous studies, denaturation of metallated (holo) hSODs was found to be irreversible, and complicated by the formation of intermolecular disulfide bonds. Here, ALS-associated mutations (E100G, G93A, G85R and A4V) are introduced into a pseudo wild-type background containing no free cysteine residues. The guanidinium chloride-induced denaturation of the holo proteins is generally found to be highly reversible (except for A4V, which tended to aggregate), enabling quantitative analysis of the effects of the mutations on protein stability. Denaturation and renaturation curves were monitored by tryptophan fluorescence, circular dichroism, enzyme activity, chemical cross-linking and analytical sedimentation, as a function of equilibration time and protein concentration. There is strong kinetic hysteresis, with curves requiring exceptionally long times (many days for pseudo wild-type) to reach equilibrium, and evidence for the formation of kinetic and equilibrium intermediate(s), which are more highly populated at lower protein concentrations. The effects of metal dissociation were included in the data fitting. The full protein concentration dependence is best described using a three-state model involving metallated native dimer, metallated monomeric intermediate and unfolded monomers with no bound metals; however, at high protein concentrations the unfolding approaches a two-state transition with metal binding to both the native dimers and unfolded monomers. We show that the E100G, G93A and G85R mutations decrease overall protein stability, largely by decreasing monomer stability with little effect on dimer dissociation. Comparison of the chemical denaturation data with ALS disease characteristics suggests that aggregation of some mutant hSOD may occur through increased population of partially folded states that are less stable than the monomeric intermediate and accessed from the destabilized holo protein.

Sequence and structural determinants of Cu, Zn superoxide dismutase aggregation

Diverse point mutations in the enzyme Cu, Zn superoxide dismutase (SOD1) are linked to its aggregation in the familial form of the disease amyotrophic lateral sclerosis. The disease-associated mutations are known to destabilize the protein, but the structural basis of the aggregation of the destabilized protein and the structure of aggregates are not well understood. Here, we investigate in silico the sequence and structural determinants of SOD1 aggregation: (1) We identify sequence fragments in SOD1 that have a high aggregation propensity, using only the sequence of SOD1, and (2) we perform molecular dynamics simulations of the SOD1 dimer folding and misfolding. In both cases, we identify identical regions of the protein as having high propensity to form intermolecular interactions. These regions correspond to the N- and C-termini, and two crossover loops and two beta-strands in the Greek-key native fold of SOD1. Our results suggest that the high aggregation propensity of mutant SOD1 may result from a synergy of two factors: the presence of highly amyloidogenic sequence fragments ("hot spots"), and the presence of these fragments in regions of the protein that are structurally most likely to form intermolecular contacts under destabilizing conditions. Therefore, we postulate that the balance between the self-association of aggregation-prone sequences and the specific structural context of these sequences in the native state determines the aggregation propensity of proteins.

Dissociation of Human Copper-Zinc Superoxide Dismutase Dimers Using Chaotrope and Reductant

The dissociation of apo- and metal-bound human copper-zinc superoxide dismutase (SOD1) dimers induced by the chaotrope guanidine hydrochloride (GdnHCl) or the reductant Tris(2-carboxyethyl)phosphine (TCEP) has been analyzed using analytical ultracentrifugation. Global fitting of sedimentation equilibrium data under native solution conditions (without GdnHCl or TCEP) demonstrate that both the apo- and metal-bound forms of SOD1 are stable dimers. Sedimentation velocity experiments show that apo-SOD1 dimers dissociate cooperatively over the range 0.5-1.0 M GdnHCl. In contrast, metal-bound SOD1 dimers possess a more compact shape and dissociate at significantly higher GdnHCl concentrations (2.0-3.0 M). Reduction of the intrasubunit disulfide bond within each SOD1 subunit by 5-10 mM TCEP promotes dissociation of apo-SOD1 dimers, whereas the metal-bound enzyme remains a stable dimer under these conditions. The Cys-57 --> Ser mutant of SOD1, a protein incapable of forming the intrasubunit disulfide bond, sediments as a monomer in the absence of metal ions and as a dimer when metals are bound. Taken together, these data indicate that the stability imparted to the human SOD1 dimer by metal binding and the formation of the intrasubunit disulfide bond are mediated by independent molecular mechanisms. By combining the sedimentation data with previous crystallographic results, a molecular explanation is provided for the existence of different SOD1 macromolecular shapes and multiple SOD1 dimeric species with different stabilities.

Protein topology determines binding mechanism

Protein recognition and binding, which result in either transient or long-lived complexes, play a fundamental role in many biological functions, but sometimes also result in pathologic aggregates. We use a simplified simulation model to survey a range of systems where two highly flexible protein chains form a homodimer. In all cases, this model, which corresponds to a perfectly funneled energy landscape for folding and binding, reproduces the macroscopic experimental observations on whether folding and binding are coupled in one step or whether intermediates occur. Owing to the minimal frustration principle, we find that, as in the case of protein folding, the native topology is the major factor that governs the choice of binding mechanism. Even when the monomer is stable on its own, binding sometimes occurs fastest through unfolded intermediates, thus showing the speedup envisioned in the fly-casting scenario for molecular recognition.

Folding of Cu, Zn Superoxide Dismutase and Familial Amyotrophic Lateral Sclerosis

Cu, Zn superoxide dismutase (SOD1) has been implicated in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis (ALS). It has been suggested that mutant mediated SOD1 misfolding/aggregation is an integral part of the pathology of ALS. We study the folding thermodynamics and kinetics of SOD1 using a hybrid molecular dynamics approach. We reproduce the experimentally observed SOD1 folding thermodynamics and find that the residues which contribute the most to SOD1 thermal stability are also crucial for apparent two-state folding kinetics. Surprisingly, we find that these residues are located on the surface of the protein and not in the hydrophobic core. Mutations in some of the identified residues are found in patients with the disease. We argue that the identified residues may play an important role in aggregation. To further characterize the folding of SOD1, we study the role of cysteine residues in folding and find that non-native disulfide bond formation may significantly alter SOD1 folding dynamics and aggregation propensity.

The recovery of dipolar relaxation times from fluorescence decays as a tool to probe local dynamics in single tryptophan proteins

The dipolar relaxation process induced by the excitation of the single tryptophan residue of four proteins (staphylococcal nuclease, ribonuclease-T1, phosphofructokinase, and superoxide dismutase) has been studied by dynamic fluorescence measurements. A new algorithm taking into account the relaxation effect has been applied to the fluorescence decay function obtained by phase-shift and demodulation data. This approach only requires that fluorescence be collected through the whole emission spectrum, avoiding the time-consuming determination of the data at different emission wavelengths, as usual with time-resolved emission spectroscopy. The results nicely match those reported in the literature for staphylococcal nuclease and ribonuclease-T1, demonstrating the validity of the model. Furthermore, this new methodology provides an alternative explanation for the complex decay of phosphofructokinase and human superoxide dismutase suggesting the presence of a relaxation process even in proteins that lack a lifetime-dependent spectral shift. These findings may have important implications on the analysis of small-scale protein dynamics, since dielectric relaxation directly probes a local structural change around the excited state of tryptophan.

Superoxide dismutase folding/unfolding pathway: role of the metal ions in modulating structural and dynamical features

The unfolding and refolding of a monomeric mutant of copper-zinc superoxide dismutase was investigated by NMR spectroscopy in the copper-reduced form and by using guanidinium chloride as denaturing agent. It is found that the protein gives rise to a series of intermediates at low guanidinium concentration and to a globular unfolded state at a guanidinium concentration higher than 3.5 M, which closely resembles a random coil structure, but with a high degree of compactness. At 1.9 M guanidinium, the intermediate and unfolded forms are present in equimolar ratio. The intermediate states show changes in the 1H and 15N chemical shifts with respect to the native protein. The perturbations on the signals occur at different GdmCl concentrations for different regions of the protein. The residues affected first are located in the loops and in the beta3 strand, followed by changes in the sheet formed by beta4, beta5, beta7, beta8 strands. The transition into the unfolded structure implies the detachment of the metal ions from the native coordination sites, even if non-specific interactions with the metal ions remain. R(1), R(2), [1H]-15N NOE, and CLEANEX experiments provide information on the mobility at the various stages showing how protein rigidity is lost during unfolding. The whole process is reversible. The oxidized species behaves in a similar way. The apo protein shows formation of 50% of the unfolded species at a guanidinium concentration of 0.4 M, thus demonstrating the importance of metal ions with respect to the unfolding process and protein structure stability. Hints to understand the whole folding process are obtained and discussed.

Common denominator of Cu/Zn superoxide dismutase mutants associated with amyotrophic lateral sclerosis: decreased stability of the apo state

More than 100 point mutations of the superoxide scavenger Cu/Zn superoxide dismutase (SOD; EC ) have been associated with the neurodegenerative disease amyotrophic lateral sclerosis (ALS). However, these mutations are scattered throughout the protein and provide no clear functional or structural clues to the underlying disease mechanism. Therefore, we undertook to look for folding-related defects by comparing the unfolding behavior of five ALS-associated mutants with distinct structural characteristics: A4V at the interface between the N and C termini, C6F in the hydrophobic core, D90A at the protein surface, and G93A and G93C, which decrease backbone flexibility. With the exception of the disruptive replacements A4V and C6F, the mutations only marginally affect the stability of the native protein, yet all mutants share a pronounced destabilization of the metal-free apo state: the higher the stability loss, the lower the mean survival time for ALS patients carrying the mutation. Thus organism-level pathology may be directly related to the properties of the immature state of a protein rather than to those of the native species.

Familial amyotrophic lateral sclerosis-associated mutations decrease the thermal stability of distinctly metallated species of human copper/zinc superoxide dismutase

We report the thermal stability of wild type (WT) and 14 different variants of human copper/zinc superoxide dismutase (SOD1) associated with familial amyotrophic lateral sclerosis (FALS). Multiple endothermic unfolding transitions were observed by differential scanning calorimetry for partially metallated SOD1 enzymes isolated from a baculovirus system. We correlated the metal ion contents of SOD1 variants with the occurrence of distinct melting transitions. Altered thermal stability upon reduction of copper with dithionite identified transitions resulting from the unfolding of copper-containing SOD1 species. We demonstrated that copper or zinc binding to a subset of "WT-like" FALS mutants (A4V, L38V, G41S, G72S, D76Y, D90A, G93A, and E133Delta) conferred a similar degree of incremental stabilization as did metal ion binding to WT SOD1. However, these mutants were all destabilized by approximately 1-6 degrees C compared with the corresponding WT SOD1 species. Most of the "metal binding region" FALS mutants (H46R, G85R, D124V, D125H, and S134N) exhibited transitions that probably resulted from unfolding of metal-free species at approximately 4-12 degrees C below the observed melting of the least stable WT species. We conclude that decreased conformational stability shared by all of these mutant SOD1s may contribute to SOD1 toxicity in FALS.

Conformational stability of dimeric and monomeric forms of the C-terminal domain of human immunodeficiency virus-1 capsid protein

The unfolding equilibrium of the C-terminal domain of human immunodeficiency virus-1 (HIV-1) capsid protein has been analyzed by circular dichroism and fluorescence spectroscopy. The results for the dimeric, natural domain are consistent with a three-state model (N(2)<-->2I<-->2U). The dimer (N(2)) dissociates and partially unfolds in a coupled cooperative process, into a monomeric intermediate (I) of very low conformational stability. This intermediate, which is the only significantly populated form at low (1 microM) protein concentrations, fully preserves the secondary structure but has lost part of the tertiary (intramonomer) interactions found in the dimer. In a second transition, the intermediate cooperatively unfolds into denatured monomer (U). The latter process is the equivalent of a two-state unfolding transition observed for a monomeric domain in which Trp184 at the dimer interface had been truncated to Ala. A highly conserved, disulfide-bonded cysteine, but not the disulfide bond itself, and three conserved residues within the major homology region of the retroviral capsid are important for the conformational stability of the monomer. All these residues are involved also in the association process, despite being located far away from the dimerization interface. It is proposed that dimerization of the C-terminal domain of the HIV-1 capsid protein involves induced-fit recognition, and the conformational reorganization also improves substantially the low intrinsic stability of each monomeric half.

PABMB Lecture. Protein dynamics, folding and misfolding: from basic physical chemistry to human conformational diseases

Proteins exhibit a variety of motions ranging from amino acid side-chain rotations to the motions of large domains. Recognition of their conformational flexibility has led to the view that protein molecules undergo fast dynamic interconversion between different conformational substates. This proposal has received support from a wide variety of experimental techniques and from computer simulations of protein dynamics. More recently, studies of the subunit dissociation of oligomeric proteins induced by hydrostatic pressure have shown that the characteristic times for subunit exchange between oligomers and for interconversion between different conformations may be rather slow (hours or days). In such cases, proteins cannot be treated as an ensemble of rapidly interconverting conformational substates, but rather as a persistently heterogeneous population of different long-lived conformers. This is reminiscent of the deterministic behavior exhibited by macroscopic bodies, and may have important implications for our understanding of protein folding and biological functions. Here, we propose that the deterministic behavior of proteins may be closely related to the genesis of conformational diseases, a class of pathological conditions that includes transmissible spongiform encephalopathies, Alzheimer's disease and other amyloidosis.

Role of pyridoxal 5'-phosphate in the structural stabilization of O-acetylserine sulfhydrylase

Proteins belonging to the superfamily of pyridoxal 5'-phosphate-dependent enzymes are currently classified into three functional groups and five distinct structural fold types. The variation within this enzyme group creates an ideal system to investigate the relationships among amino acid sequences, folding pathways, and enzymatic functions. The number of known three-dimensional structures of pyridoxal 5'-phosphate-dependent enzymes is rapidly increasing, but only for relatively few have the folding mechanisms been characterized in detail. The dimeric O-acetylserine sulfhydrylase from Salmonella typhimurium belongs to the beta-family and fold type II group. Here we report the guanidine hydrochloride-induced unfolding of the apo- and holoprotein, investigated using a variety of spectroscopic techniques. Data from absorption, fluorescence, circular dichroism, (31)P nuclear magnetic resonance, time-resolved fluorescence anisotropy, and photon correlation spectroscopy indicate that the O-acetylserine sulfhydrylase undergoes extensive disruption of native secondary and tertiary structure before monomerization. Also, we have observed that the holo-O-acetylserine sulfhydrylase exhibits a greater conformational stability than the apoenzyme form. The data are discussed in light of the fact that the role of the coenzyme in structural stabilization varies among the pyridoxal 5'-phosphate-dependent enzymes and does not seem to be linked to the particular enzyme fold type.

Selective neoglycosylation increases the structural stability of vicilin, the 7S storage globulin from pea seeds

The effects of glycosylation on the stability and subunit interactions of vicilin, the major storage protein in pea seeds, were investigated. Glycosylated vicilin derivatives were prepared by alkylation of lysine epsilon-amino groups with various carbohydrates. Average modification levels of 13.4 +/- 3.0, 11.1 +/- 3.6, 7.5 +/- 4.2, and 4.7 +/- 0.3 moles of carbohydrate/mol of vicilin were obtained with glucose, galactose, galacturonic acid, and lactose, respectively. Nondenaturing polyacrylamide gel electrophoresis and size-exclusion chromatography indicated that the quaternary structure and hydrodynamic radius of vicilin were not affected by glycosylation at the levels used. We have previously shown that application of hydrostatic pressure causes dissociation of vicilin subunits [C. Pedrosa and S. T. Ferreira (1994) Biochemistry 33, 4046-4055]. Analysis of pressure dissociation data allowed determination of the Gibbs free energy change (deltaG(diss)) and molar volume change (deltaV(diss)) of dissociation of vicilin subunits. For unmodified vicilin, deltaG(diss) = 18.2 kcal/mol and deltaV(diss) = -102 ml/mol. Glycosylated vicilin derivatives were significantly stabilized against subunit dissociation, with deltaG(diss) of 19.4, 19.2, 20.6, and 22.1 kcal/mol for glucose, galactose, lactose, and galacturonic acid derivatives, respectively. No changes in deltaV(diss) were found for the glucose and galactose derivatives, whereas deltaV(diss) of -128 and -135 ml/mol, respectively, were found for the lactose and galacturonic acid derivatives. The glycosylated derivatives also appeared more resistant to unfolding by guanidine hydrochloride than unmodified vicilin. Intrinsic fluorescence lifetime measurements showed that glycosylation caused a significant increase in heterogeneity of the fluorescence decay, possibly reflecting increased conformational heterogeneity of glycosylated derivatives relative to unmodified vicilin. These results indicate that the stability and subunit interactions of vicilin may be modulated by mild, selective glycosylation at low modification levels, an effect that may be of interest in the study of other oligomeric proteins.

Role of the tertiary and quaternary structures in the stability of dimeric copper, zinc superoxide dismutases

The equilibrium unfolding process of human Cu,Zn superoxide dismutase has been quantitatively monitored through circular dichroism and fluorescence spectroscopy as a function of increasing guanidinium hydrochloride concentration. The process occurs through the formation of a monomeric intermediate species following a three-state transition equilibrium. Comparison with the stability of the prokaryotic Cu,Zn SOD from P. leiognathi shows that the eukaryotic enzyme is more stable than the prokaryotic enzyme by approximately 3 kcal/mol. This difference is due to the monomer-to-unfolded equilibrium, while the dimer-to-monomer equilibrium is comparable for the two enzymes despite their different intersubunit interactions. These results are confirmed by the unfolding of the copper-depleted derivatives. The Cu,Zn superoxide dismutase represents a good example of how evolution has found two independent quaternary assemblies maintaining the same dimer stability.

Evidence of stable monomeric species in the unfolding of Cu,Zn superoxide dismutase from Photobacterium leiognathi

The equilibrium unfolding process of Photobacterium leiognathi Cu,Zn superoxide dismutase has been quantitatively monitored through circular dichroism (CD) and fluorescence spectroscopy, upon increasing the guanidinium hydrochloride concentration. The study has been undertaken for both the holo- and the copper-free derivative to work out the role of copper in protein stability. In both cases the unfolding was reversible. The denaturation curve derived from CD and fluorescence spectroscopy was not coincident, suggesting that the denaturation process occurs through a three-state model with formation of an intermediate monomeric species. The occurrence of an intermediate species has been unambiguously demonstrated following CD and steady-state fluorescence spectra of the enzyme at various concentrations in presence of a fixed amounts of guanidinium hydrochloride.

Apohorseradish peroxidase unfolding and refolding: intrinsic tryptophan fluorescence studies

The unfolding and refolding of apohorseradish peroxidase, as a function of guanidinium chloride concentration, were monitored by the intrinsic fluorescence intensity, polarization, and lifetime of the single tryptophan residue. The unfolding was reversible and characterized by at least three distinct stages-the intensity and lifetime data, for example, were both characterized by an initial increase followed by a decrease and then a plateau region. The lifetime data, in the absence and presence of guanidinium chloride, were heterogeneous and fit best to a model consisting of a major Gaussian distribution component and a minor, short discrete component. The observed increase in intensity in the initial stage of the unfolding process is attributed to the conversion of this short component into the longer, distributed component as the guanidinium chloride concentration increases. Our results clarify and amplify previous studies on the unfolding of apohorseradish peroxidase by guanidinium chloride.

Mechanism and evolution of protein dimerization

We have investigated the mechanism and the evolutionary pathway of protein dimerization through analysis of experimental structures of dimers. We propose that the evolution of dimers may have multiple pathways, including (1) formation of a functional dimer directly without going through an ancestor monomer, (2) formation of a stable monomer as an intermediate followed by mutations of its surface residues, and (3), a domain swapping mechanism, replacing one segment in a monomer by an equivalent segment from an identical chain in the dimer. Some of the dimers which are governed by a domain swapping mechanism may have evolved at an earlier stage of evolution via the second mechanism. Here, we follow the theory that the kinetic pathway reflects the evolutionary pathway. We analyze the structure-kinetics-evolution relationship for a collection of symmetric homodimers classified into three groups: (1) 14 dimers, which were referred to as domain swapping dimers in the literature; (2) nine 2-state dimers, which have no measurable intermediates in equilibrium denaturation; and (3), eight 3-state dimers, which have stable intermediates in equilibrium denaturation. The analysis consists of the following stages: (i) The dimer is divided into two structural units, which have twofold symmetry. Each unit contains a contiguous segment from one polypeptide chain of the dimer, and its complementary contiguous segment from the other chain. (ii) The division is repeated progressively, with different combinations of the two segments in each unit. (iii) The coefficient of compactness is calculated for the units in all divisions. The coefficients obtained for different cuttings of a dimer form a compactness profile. The profile probes the structural organization of the two chains in a dimer and the stability of the monomeric state. We describe the features of the compactness profiles in each of the three dimer groups. The profiles identify the swapping segments in domain swapping dimers, and can usually predict whether a dimer has domain swapping. The kinetics of dimerization indicates that some dimers which have been assigned in the literature as domain swapping cases, dimerize through the 2-state kinetics, rather than through swapping segments of performed monomers. The compactness profiles indicate a wide spectrum in the kinetics of dimerization: dimers having no intermediate stable monomers; dimers having an intermediate with a stable monomer structure; and dimers having an intermediate with a stable structure in part of the monomer. These correspond to the multiple evolutionary pathways for dimer formation. The evolutionary mechanisms proposed here for dimers are applicable to other oligomers as well.

Spectroscopic study of conformational changes in subdomain 1 of G-actin: influence of divalent cations

Temperature dependence of the fluorescence intensity and anisotropy decay of N-(iodoacetyl)-N'-(5-sulfo-1-naphthyl)ethylenediamine attached to Cys374 of actin monomer was investigated to characterize conformational differences between Ca- and Mg-G-actin. The fluorescence lifetime is longer in Mg-G-actin than that in Ca-G-actin in the temperature range of 5-34 degrees C. The width of the lifetime distribution is smaller by 30% in Mg-saturated actin monomer at 5 degrees C, and the difference becomes negligible above 30 degrees C. The semiangle of the cone within which the fluorophore can rotate is larger in Ca-G-actin at all temperatures. Electron paramagnetic resonance measurements on maleimide spin-labeled (on Cys374) monomer actin gave evidence that exchange of Ca2+ for Mg2+ induced a rapid decrease in the mobility of the label immediately after the addition of Mg2+. These results suggest that the C-terminal region of the monomer becomes more rigid as a result of the replacement of Ca2+ by Mg2+. The change can be related to the difference between the polymerization abilities of the two forms of G-actin.

Protein unfolding by peptidylarginine deiminase. Substrate specificity and structural relationships of the natural substrates trichohyalin and filaggrin

Peptidylarginine deiminases, which are commonly found in mammalian cells, catalyze the deimination of protein-bound arginine residues to citrullines. However, very little is known about their substrate requirements and the significance or consequences of this postsynthetic modification. We have explored this reaction in vitro with two known substrates filaggrin and trichohyalin. First, the degree and rate of modification of arginines to citrullines directly correlates with the structural order of the substrate. In filaggrin, which has little structural order, the reaction proceeded rapidly to >95% completion. However, in the highly alpha-helical protein trichohyalin, the reaction proceeded slowly to about 25% and could be forced to a maximum of about 65%. Second, the rate and degree of modification depends on the sequence location of the target arginines. Third, we show by gel electrophoresis, circular dichroism, and fluorescence spectroscopy that the reaction interferes with organized protein structure: the net formation of >/=10% citrulline results in protein denaturation. Cyanate modification of the lysines in model alpha-helix-rich proteins to homocitrullines also results in loss of organized structure. These data suggest that the ureido group on the citrulline formed by the peptidylarginine deiminase enzyme modification functions to unfold proteins due to decrease in net charge, loss of potential ionic bonds, and interference with H bonds.

Probing the structure and mobility of Pseudomonas aeruginosa azurin by circular dichroism and dynamic fluorescence anisotropy

The UV dynamic fluorescence and CD of several Pseudomonas aeruginosa azurins bearing single amino acid mutation have been studied. Two classes of mutants were examined. In the first class, two hydrophobic residues in the core of the protein, Ile 7 and Phe 110, nearest to the azurin single tryptophan Trp 48, were substituted by a serine (mutants 17S and F110S). In the second class, two residues in the outer sphere of the copper ligand field were changed, obtaining the following mutants: M44K, H35F, H35L, and H35Q. All these proteins showed two fluorescence lifetimes in the copper-containing form, but only one in the copper-free form. The lifetime of the latter derivatives was different from either those of the metal-bound samples, definitely ruling out the presence of apo-like species in the holo protein. Copper-free 17S and F110S showed a more complex fluorescence decay profile requiring a distribution of lifetimes rather than a single lifetime. Holo F110S was also better fitted, in the limit of confidence, with two distributions rather than a pair of lifetimes. Time-resolved anisotropy of these two mutants as well as of wild-type (wt) protein showed two components (rotational times for wt < or = 200 ps and 7 ns, respectively). These components were not affected significantly by copper removal in the case of wt protein. Instead, the short rotational component of the mutants dropped dramatically to values near zero, indicating a much greater mobility of the tryptophanyl residue in the mutant apo azurins. These data were supported by CD measurements showing a small effect of the copper presence in the region below 250 nm, i.e., in the secondary structure, but almost a collapse of the aromatic asymmetry at 270-295 nm related to a relaxation of the structural constraint around the tryptophan. Altogether these data show that copper does not play a structural role in wt azurin, whereas it is crucial in the stabilization of 17S and F110S mutants. Furthermore, although the metal site geometry is rigidly kept in wt apo-azurin, it regains the native form only in the presence of the metal in the "core" mutants. This finding is important for the theory of entatic states in metalloproteins (Williams RJP, 1995, Eur J Biochem 234:363-381).

A stochastic model for the sigmoidal behaviour of cooperative biological systems

A stochastic model for cooperative transitions in biological systems based on a Markov chain is proposed. This model requires only two parameters, the mean probability, p, and the coupling capacity, Deltap, which measure the probability of forming a new weak bond depending on the number of similar bonds already formed and it is also responsible for the transition. In this paper we show how the model works for a large number of identical molecules and how it can be useful for studying the noise around the centre of the transition where, increasing the degree of cooperativity, i.e. the number n in the well-known Hill equation, the width of the noise increases along with its fractal dimension. A simple relationship between the degree of cooperativity and the parameter Deltap is proposed, suggesting that the cooperativity of real biological transitions is related to the coupling capacity Deltap of the present model.

Conformational stability of dimeric proteins: quantitative studies by equilibrium denaturation

The conformational stability of dimeric globular proteins can be measured by equilibrium denaturation studies in solvents such as guanidine hydrochloride or urea. Many dimeric proteins denature with a 2-state equilibrium transition, whereas others have stable intermediates in the process. For those proteins showing a single transition of native dimer to denatured monomer, the conformational stabilities, delta Gu (H2O), range from 10 to 27 kcal/mol, which is significantly greater than the conformational stability found for monomeric proteins. The relative contribution of quaternary interactions to the overall stability of the dimer can be estimated by comparing delta Gu (H2O) from equilibrium denaturation studies to the free energy associated with simple dissociation in the absence of denaturant. In many cases the large stabilization energy of dimers is primarily due to the intersubunit interactions and thus gives a rationale for the formation of oligomers. The magnitude of the conformational stability is related to the size of the polypeptide in the subunit and depends upon the type of structure in the subunit interface. The practical use, interpretation, and utility of estimation of conformational stability of dimers by equilibrium denaturation methods are discussed.

Similarity of fluorescence lifetime distributions for single tryptophan proteins in the random coil state

The picosecond time-resolved fluorescence decay data of nine single-tryptophan (trp) proteins and two multi-trp proteins in their native and denatured states were analyzed by the maximum entropy method (MEM). In the denatured state (6 M guanidine hydrochloride) a majority of the single-trp proteins show bimodal (at 25 degrees C) and trimodal (at 85 degrees C) distributions with similar patterns and similar values for average lifetimes. In the native state of the proteins the lifetime distributions were bimodal or trimodal. These results (multimodal distributions) are contradictory to the unimodal Lorentzian distribution of lifetimes reported for some proteins in the native and denatured states. MEM analysis gives a unimodal distribution of lifetimes only when the signal-to-noise ratio is poor in the time-resolved fluorescence decay data. The unimodal distribution model is therefore not realistic for proteins in the native and denatured states. The fluorescence decay components of the bi- or trimodal distribution are associated with the rotamer structures of the indole moiety when the protein is in the random coil state.

Probing alpha-helical secondary structure at a specific site in model peptides via restriction of tryptophan side-chain rotamer conformation

The relationship between alpha-helical secondary structure and the fluorescence properties of an intrinsic tryptophan residue were investigated. A monomeric alpha-helix forming peptide and a dimeric coiled-coil forming peptide containing a central tryptophan residue were synthesized. The fluorescence parameters of the tryptophan residue were determined for these model systems at a range of fractional alpha-helical contents. The steady-state emission maximum was independent of the fractional alpha-helical content. A minimum of three exponential decay times was required to fully describe the time-resolved fluorescence data. Changes were observed in the decay times and more significantly, in their relative contributions that could be correlated with alpha-helix content. The results were also shown to be consistent with a model in which the decay times were independent of both alpha-helix content and emission wavelength. In this model the relative contributions of the decay time components were directly proportional to the alpha-helix content. Data were also analyzed according to a continuous distribution of exponential decay time model, employing global analysis techniques. The recovered distributions had "widths" that were both poorly defined and independent of peptide conformation. We propose that the three decay times are associated with the three ground-state chi 1 rotamers of the tryptophan residue and that the changes in the relative contributions of the decay times are the result of conformational constraints, imposed by the alpha-helical main-chain, on the chi 1 rotamer populations.