Kinetic Control in Protein Folding for Light Chain Amyloidosis and the Differential Effects of Somatic Mutations

Light chain amyloidosis is a devastating disease where immunoglobulin light chains form amyloid fibrils, resulting in organ dysfunction and death. Previous studies have shown a direct correlation between the protein thermodynamic stability and the propensity for amyloid formation for some proteins involved in light chain amyloidosis. Here we investigate the effect of somatic mutations on protein stability and in vitro fibril formation of single and double restorative mutants of the protein AL-103 compared to the wild-type germline control protein. A scan rate dependence and hysteresis in the thermal unfolding and refolding was observed for all proteins. This indicates that the unfolding/refolding reaction is kinetically determined with different kinetic constants for unfolding and refolding even though the process remains experimentally reversible. Our structural analysis of AL-103 and AL-103 delP95aIns suggests a kinetic coupling of the unfolding/refolding process with cis–trans prolyl isomerization. Our data reveal that the deletion of proline 95a (AL-103 delP95aIns), which removes the trans–cis di-proline motif present in the patient protein AL-103, results in a dramatic increment in the thermodynamic stability and a significant delay in fibril formation kinetics with respect to AL-103. Fibril formation is pH dependent; all proteins form fibrils at pH 2; reactions become slower and more stochastic as the pH increases up to pH 7. Based on these results, we propose that, in addition to thermodynamic stability, kinetic stability (possibly influenced by the presence of cis proline 95a) plays a major role in the AL-103 amyloid fibril formation process.     

Fibrillogenesis of apomyoglobin facilitated by aggregation sequence of yeast Sup35 in various regions

To examine the effect of aggregation sequence QGGYQQQYNP from yeast Sup35 on fibril formation of sperm whale apomyoglobin (apoMb), we constructed several mutants via substitution. Urea-induced unfolding of apoMb confirms that the substitution of the aggregation sequence does not significantly affect the stability of the mutants compared to wild type (WT) at pH 4.2. Under this condition, however, despite the difference in rate most apoMb mutants form fibrils more readily than WT with distinct morphology. These results suggest that the aggregation sequence facilitates fibril assembly of apoMb at acidic pH in vitro and this facilitation depends on the regions replaced.     

Hydrophobic Core Mutations Associated with Cataract Development in Mice Destabilize Human ��D-Crystallin

The human eye lens is composed of fiber cells packed with crystallins up to 450 mg/ml. Human γD-crystallin (HγD-Crys) is a monomeric, two-domain protein of the lens central nucleus. Both domains of this long lived protein have double Greek key β-sheet folds with well packed hydrophobic cores. Three mutations resulting in amino acid substitutions in the γ-crystallin buried cores (two in the N-terminal domain (N-td) and one in the C-terminal domain (C-td)) cause early onset cataract in mice, presumably an aggregated state of the mutant crystallins. It has not been possible to identify the aggregating precursor within lens tissues. To compare in vivo cataract-forming phenotypes with in vitro unfolding and aggregation of γ-crystallins, mouse mutant substitutions were introduced into HγD-Crys. The mutant proteins L5S, V75D, and I90F were expressed and purified from Escherichia coli. WT HγD-Crys unfolds in vitro through a three-state pathway, exhibiting an intermediate with the N-td unfolded and the C-td native-like. L5S and V75D in the N-td also displayed three-state unfolding transitions, with the first transition, unfolding of the N-td, shifted to significantly lower denaturant concentrations. I90F destabilized the C-td, shifting the overall unfolding transition to lower denaturant concentrations. During thermal denaturation, the mutant proteins exhibited lowered thermal stability compared with WT. Kinetic unfolding experiments showed that the N-tds of L5S and V75D unfolded faster than WT. I90F was globally destabilized and unfolded more rapidly. These results support models of cataract formation in which generation of partially unfolded species are precursors to the aggregated cataractous states responsible for light scattering.     

Experimental comparison of energy landscape features of ubiquitin family proteins

Small ubiquitin-related modifiers (SUMO1 and SUMO2) are ubiquitin family proteins, structurally similar to ubiquitin, differing in terms of their amino acid sequence and functions. Therefore, they provide a great platform for investigating sequence-structure-stability-function relationship. Here, we used chemical denaturation in comparing the folding-unfolding pathways of the SUMO proteins with their structural homologue ubiquitin (UF45W-pseudo wild-type [WT] tryptophan variant) with structurally analogous tryptophan mutations (SUMO1 [S1F66W], SUMO2 [S2F62W]). Equilibrium denaturation studies report that ubiquitin is the most stable protein among the three. The observed denaturant-dependent folding rates of SUMOs are much lower than ubiquitin and primarily exhibit a two-state folding pathway unlike ubiquitin, which has a kinetic folding intermediate. We hypothesize that, as SUMO proteins start off as slow folders, they avoid stabilizing their folding intermediates and the presence of which might further slow-down their folding rates. The denaturant-dependent unfolding of ubiquitin is the fastest, followed by SUMO2, and slowest for SUMO1. However, the spontaneous unfolding rate constant is the lowest for ubiquitin (~40 times), and similar for SUMOs. This correlation between thermodynamic stability and kinetic stability is achieved by having different unfolding transition state positions with respect to the solvent-accessible surface area, as quantified by the Tanford β _u values: ubiquitin (0.42) > SUMO2 (0.20) > SUMO1 (0.16). The results presented here highlight the unique energy landscape features which help in optimizing the folding-unfolding rates within a structurally homologous protein family.     

Trans-dimerization of JAM-A regulates Rap2 and is mediated by a domain that is distinct from the cis-dimerization interface

Junctional adhesion molecule-A (JAM-A) is a tight junction–associated signaling protein that regulates epithelial cell proliferation, migration, and barrier function. JAM-A dimerization on a common cell surface (in cis) has been shown to regulate cell migration, and evidence suggests that JAM-A may form homodimers between cells (in trans). Indeed, transfection experiments revealed accumulation of JAM-A at sites between transfected cells, which was lost in cells expressing cis- or predicted trans-dimerization null mutants. Of importance, microspheres coated with JAM-A containing alanine substitutions to residues 43NNP45 (NNP-JAM-A) within the predicted trans-dimerization site did not aggregate. In contrast, beads coated with cis-null JAM-A demonstrated enhanced clustering similar to that observed with wild-type (WT) JAM-A. In addition, atomic force microscopy revealed decreased association forces in NNP-JAM-A compared with WT and cis-null JAM-A. Assessment of effects of JAM-A dimerization on cell signaling revealed that expression of trans- but not cis-null JAM-A mutants decreased Rap2 activity. Furthermore, confluent cells, which enable trans-dimerization, had enhanced Rap2 activity. Taken together, these results suggest that trans-dimerization of JAM-A occurs at a unique site and with different affinity compared with dimerization in cis. Trans-dimerization of JAM-A may thus act as a barrier-inducing molecular switch that is activated when cells become confluent.     

Contributions of hydrophobic domain interface interactions to the folding and stability of human gammaD-crystallin

Human gammaD-crystallin (HgammaD-Crys) is a monomeric eye lens protein composed of two highly homologous beta-sheet domains. The domains interact through interdomain side chain contacts forming two structurally distinct regions, a central hydrophobic cluster and peripheral residues. The hydrophobic cluster contains Met43, Phe56, and Ile81 from the N-terminal domain (N-td) and Val132, Leu145, and Val170 from the C-terminal domain (C-td). Equilibrium unfolding/refolding of wild-type HgammaD-Crys in guanidine hydrochloride (GuHCl) was best fit to a three-state model with transition midpoints of 2.2 and 2.8 M GuHCl. The two transitions likely corresponded to sequential unfolding/refolding of the N-td and the C-td. Previous kinetic experiments revealed that the C-td refolds more rapidly than the N-td. We constructed alanine substitutions of the hydrophobic interface residues to analyze their roles in folding and stability. After purification from E. coli, all mutant proteins adopted a native-like structure similar to wild type. The mutants F56A, I81A, V132A, and L145A had a destabilized N-td, causing greater population of the single folded domain intermediate. Compared to wild type, these mutants also had reduced rates for productive refolding of the N-td but not the C-td. These data suggest a refolding pathway where the domain interface residues of the refolded C-td act as a nucleating center for refolding of the N-td. Specificity of domain interface interactions is likely important for preventing incorrect associations in the high protein concentrations of the lens nucleus.     

Molecular characterization of the human COQ5 C-methyltransferase in coenzyme Q10 biosynthesis

Coq5 catalyzes the only C-methylation involved in the biosynthesis of coenzyme Q (Q or ubiquinone) in humans and yeast Saccharomyces cerevisiae. As one of eleven polypeptides required for Q production in yeast, Coq5 has also been shown to assemble with the multi-subunit complex termed the CoQ-synthome. In humans, mutations in several COQ genes cause primary Q deficiency, and a decrease in Q biosynthesis is associated with mitochondrial, cardiovascular, kidney and neurodegenerative diseases. In this study, we characterize the human COQ5 polypeptide and examine its complementation of yeast coq5 point and null mutants. We show that human COQ5 RNA is expressed in all tissues and that the COQ5 polypeptide is associated with the mitochondrial inner membrane on the matrix side. Previous work in yeast has shown that point mutations within or adjacent to conserved COQ5 methyltransferase motifs result in a loss of Coq5 function but not Coq5 steady state levels. Here, we show that stabilization of the CoQ-synthome within coq5 point mutants or by over-expression of COQ8 in coq5 null mutants permits the human COQ5 homolog to partially restore coq5 mutant growth on respiratory media and Q₆ content. Immunoblotting against the human COQ5 polypeptide in isolated yeast mitochondria shows that the human Coq5 polypeptide migrates in two-dimensional blue-native/SDS-PAGE at the same high molecular mass as other yeast Coq proteins. The results presented suggest that human and Escherichia coli Coq5 homologs expressed in yeast retain C-methyltransferase activity but are capable of rescuing the coq5 yeast mutants only when the CoQ-synthome is assembled.     

Ubiquitin-metallothionein fusion protein expression in yeast. A genetic approach for analysis of ubiquitin functions

We have established a Saccharomyces cerevisiae genetic system that expresses the fusion protein ubiquitin-metallothionein. We have evaluated the effects of amino-terminal ubiquitination of metallothionein on the stability and function of metallothionein. The fusion protein of wild type ubiquitin and metallothionein was rapidly processed in vivo to release free ubiquitin and metallothionein. Site-directed mutants of ubiquitin-metallothionein expressed in yeast were used to study the specificity of the (alpha-NH2-ubiquitin) protein endopeptidases. The data suggest that amino-terminal ubiquitination is not a signal for the proteolysis of yeast metallothionein in yeast. We have also discovered that expression of selected ubiquitin mutants blocked the growth of yeast. The data suggest that in addition to its function as a proteolytic signal, ubiquitination of proteins plays multiple roles in the cell.     

Structural stability of GlcV, the nucleotide binding domain of the glucose ABC transporter of Sulfolobus solfataricus

GlcV is the nucleotide binding domain of the ABC-type glucose transporter of the hyperthermoacidophile Sulfolobus solfataricus. GlcV consists of two domains, an N-terminal domain containing the typical nucleotide binding-fold and a C-terminal β-barrel domain with unknown function. The unfolding and structural stability of the wild-type (wt) protein and three mutants that are blocked at different steps in the ATP hydrolytic cycle were studied. The G144A mutant is unable to dimerize, while the E166A and E166Q mutants are defective in ATP hydrolysis and dimer dissociation. Unfolding of the wt GlcV and G144A GlcV occurred with a single transition, whereas the E166A and E166Q mutants showed a second transition at a higher melting temperature indicating an increased stability of the ABCα/β subdomain. The structural stability of GlcV was increased in the presence of nucleotides suggesting that the transition corresponds to the unfolding of the NBD domain. Unfolding of the C-terminal domain appears to occur at temperatures above the unfolding of the NBD which coincides with the aggregation of the protein. Analysis of the domain organization of GlcV by trypsin digestion demonstrates cleavage of the NBD domain into three fragments, while nucleotides protect against proteolysis. The cleaved GlcV protein retained the ability to bind nucleotides and to dimerize. These data indicate that the wt GlcV NBD domain unfolds as a single domain protein, and that its stability is modified by mutations in the glutamate after the Walker B motif and by nucleotide binding.     

Structural Insights into the Folding Defects of Oncogenic pVHL Lead to Correction of Its Function In Vitro

Loss of function mutations in the von Hippel-Lindau (pVHL) tumor suppressor protein are tumorigenic. In silico analysis of the structure and folding of WT pVHL identified in its core an aromatic tetrahedron, essential for stabilizing the protein. The mutations disrupt the aromatic tetrahedron, leading to misfolding of pVHL. Using biophysical methods we confirmed the in silico predictions, demonstrating that mutant pVHL proteins have lower stability than the WT, distort the core domain and as a result reduce the ability of the protein to bind its target HIF-1α. Using bacterial pVHL-EGFP based assay we screened for osmolytes capable of restoring folding of mutant pVHL. Among them, Arginine was the most effective and was verified by in vitro assays as a potent re-folder of pVHL. This resulted in functional restoration of the mutant proteins to the level of the WT.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

Kinetic stability of Cu/Zn superoxide dismutase is dependent on its metal ligands: implications for ALS

Over 100 mutants of the enzyme Cu/Zn superoxide dismutase (SOD) have been implicated in the neurodegenerative disease familial amyotrophic lateral sclerosis (FALS). Growing evidence suggests that the aggregation of SOD mutants may play a causative role in FALS and that aberrant copper chemistry, decreased thermodynamic stability, and decreased affinity for metals may contribute independently or synergistically to this process. Since the loss of the copper and zinc ions significantly decreases the thermodynamic stability of SOD, it is expected that this would also decrease its kinetic stability, thereby facilitating partial or global unfolding transitions that may lead to misfolding and aggregation. Here we used wild-type (WT) SOD and five FALS-related mutants (G37R, H46R, G85R, D90A, and L144F) to show that the metals contribute significantly to the kinetic stability of the protein, with demetalated (apo) SOD showing acid-induced unfolding rates about 60-fold greater than the metalated (holo) protein. However, the unfolding rates of SOD WT and mutants were similar to each other in both the holo and apo states, indicating that regardless of the effect of mutation on thermodynamic stability, the kinetic barrier toward SOD unfolding is dependent on the presence of metals. Thus, these results suggest that pathogenic SOD mutations that do not significantly alter the stability of the protein may still lead to SOD aggregation by compromising its ability to bind or retain its metals and thereby decrease its kinetic stability. Furthermore, the mutant-like decrease in the kinetic stability of apo WT SOD raises the possibility that the loss of metals in WT SOD may be involved in nonfamilial forms of ALS.     

Improving the thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 using a computationally aided method

Good protein thermostability is very important for the protein application. In this report, we propose a strategy which contained a prediction method to select residues related to protein thermal stability, but not related to protein function, and an experiment method to screen the mutants with enhanced thermostability. The prediction strategy was based on the calculated site evolutionary entropy and unfolding free energy difference between the mutant and wild-type (WT) methyl parathion hydrolase enzyme from Ochrobactrum sp. M231 [Ochr-methyl parathion hydrolase (MPH)]. As a result, seven amino acid sites within Ochr-MPH were selected and used to construct seven saturation mutagenesis libraries. The results of screening these libraries indicated that six sites could result in mutated enzymes exhibiting better thermal stability than the WT enzyme. A stepwise evolutionary approach was designed to combine these selected mutants and a mutant with four point mutations (S274Q/T183E/K197L/S192M) was selected. The T _m and T ₅₀ of the mutant enzyme were 11.7 and 10.2 °C higher, respectively, than that of the WT enzyme. The success of this design methodology for Ochr-MPH suggests that it was an efficient strategy for enhancing protein thermostability and suitable for protein engineering.     

Exploring sequence/folding space: folding studies on multiple hydrophobic core mutants of ubiquitin

The stability, dynamic, and structural properties of ubiquitin and two multiple hydrophobic core mutants were studied. One of the mutants (U4) has seven substitutions in the hydrophobic core (M1L, I3L, V5I, I13F, L15V, V17M, and V26L). On average, its side chains are larger than the wild-type, and it can thus be thought of as having an overpacked core. The other mutant (U7) has two substitutions (I3V and I13V). On average, it has smaller side chains than the wild-type, and it can therefore be considered to be underpacked. The three proteins are well-folded and show similar backbone dynamics (T(1), T(2), and HNOE values), indicating that the regular secondary structure extends over the same residue ranges. The crystallographic structure of U4 was determined. The final R(factor) and R(free) are 0.198 and 0.248, respectively, at 2.18 A resolution. The structure of U4 is very similar to wild-type ubiquitin. Remarkably, there are almost no changes in the positions of the C(alpha) atoms along the entire backbone, and the hydrogen-bonding network is maintained. The mutations of the hydrophobic core are accommodated by small movements of side chains in the core of mutated and nonmutated residues. Unfolding and refolding kinetic studies revealed that U4 unfolds with the highest rates; however, its refolding rate constants are very similar to those of the wild-type protein. Conversely, U7 seems to be the most destabilized protein; its refolding rate constant is smaller than the other two proteins. This was confirmed by stopped-flow techniques and by H/D exchange methodologies. This work illustrates the possibility of repacking the hydrophobic core of small proteins and has important implications in the de novo design of stable proteins.     

Role of the Deubiquitylating Enzyme DmUsp5 in Coupling Ubiquitin Equilibrium to Development and Apoptosis in Drosophila melanogaster

Protein ubiquitylation is a dynamic process that affects the function and stability of proteins and controls essential cellular processes ranging from cell proliferation to cell death. This process is regulated through the balanced action of E3 ubiquitin ligases and deubiquitylating enzymes (DUB) which conjugate ubiquitins to, and remove them from target proteins, respectively. Our genetic analysis has revealed that the deubiquitylating enzyme DmUsp5 is required for maintenance of the ubiquitin equilibrium, cell survival and normal development in Drosophila. Loss of the DmUsp5 function leads to late larval lethality accompanied by the induction of apoptosis. Detailed analyses at a cellular level demonstrated that DmUsp5 mutants carry multiple abnormalities, including a drop in the free monoubiquitin level, the excessive accumulation of free polyubiquitins, polyubiquitylated proteins and subunits of the 26S proteasome. A shortage in free ubiquitins results in the induction of a ubiquitin stress response previously described only in the unicellular budding yeast. It is characterized by the induction of the proteasome-associated deubiquitylase DmUsp14 and sensitivity to cycloheximide. Removal of DmUsp5 also activates the pro-apoptotic machinery thereby resulting in widespread apoptosis, indicative of an anti-apoptotic role of DmUsp5. Collectively, the pleiotropic effects of a loss of DmUsp5 function can be explained in terms of the existence of a limited pool of free monoubiquitins which makes the ubiquitin-dependent processes mutually interdependent.     

Single-Molecule Studies on PolySUMO Proteins Reveal Their Mechanical Flexibility

Proteins with β-sandwich and β-grasp topologies are resistant to mechanical unfolding as shown by single-molecule force spectroscopy studies. Their high mechanical stability has generally been associated with the mechanical clamp geometry present at the termini. However, there is also evidence for the importance of interactions other than the mechanical clamp in providing mechanical stability, which needs to be tested thoroughly. Here, we report the mechanical unfolding properties of ubiquitin-like proteins (SUMO1 and SUMO2) and their comparison with those of ubiquitin. Although ubiquitin and SUMOs have similar size and structural topology, they differ in their sequences and structural contacts, making them ideal candidates to understand the variations in the mechanical stability of a given protein topology. We observe a two-state unfolding pathway for SUMO1 and SUMO2, similar to that of ubiquitin. Nevertheless, the unfolding forces of SUMO1 (∼130 pN) and SUMO2 (∼120 pN) are lower than that of ubiquitin (∼190 pN) at a pulling speed of 400 nm/s, indicating their lower mechanical stability. The mechanical stabilities of SUMO proteins and ubiquitin are well correlated with the number of interresidue contacts present in their structures. From pulling speed-dependent mechanical unfolding experiments and Monte Carlo simulations, we find that the unfolding potential widths of SUMO1 (∼0.51 nm) and SUMO2 (∼0.33 nm) are much larger than that of ubiquitin (∼0.19 nm), indicating that SUMO1 is six times and SUMO2 is three times mechanically more flexible than ubiquitin. These findings might also be important in understanding the functional differences between ubiquitin and SUMOs.     

The TREX1 C-terminal Region Controls Cellular Localization through Ubiquitination

TREX1 is an autonomous 3′-exonuclease that degrades DNA to prevent inappropriate immune activation. The TREX1 protein is composed of 314 amino acids; the N-terminal 242 amino acids contain the catalytic domain, and the C-terminal region (CTR) localizes TREX1 to the cytosolic compartment. In this study, we show that TREX1 modification by ubiquitination is controlled by a highly conserved sequence in the CTR to affect cellular localization. Transfection of TREX1 deletion constructs into human cells demonstrated that this sequence is required for ubiquitination at multiple lysine residues through a “non-canonical” ubiquitin linkage. A proteomic approach identified ubiquilin 1 as a TREX1 CTR-interacting protein, and this interaction was verified in vitro and in vivo. Cotransfection studies indicated that ubiquilin 1 localizes TREX1 to cytosolic punctate structures dependent upon the TREX1 CTR and lysines within the TREX1 catalytic core. Several TREX1 mutants linked to the autoimmune diseases Aicardi-Goutières syndrome and systemic lupus erythematosus that exhibit full catalytic function were tested for altered ubiquitin modification and cellular localization. Our data show that these catalytically competent disease-causing TREX1 mutants exhibit differential levels of ubiquitination relative to WT TREX1, suggesting a novel mechanism of dysfunction. Furthermore, these differentially ubiquitinated disease-causing mutants also exhibit altered ubiquilin 1 co-localization. Thus, TREX1 post-translational modification indicates an additional mechanism by which mutations disrupt TREX1 biology, leading to human autoimmune disease.