Protein Stability and Folding Kinetics in the Nucleus and Endoplasmic Reticulum of Eucaryotic Cells

We measure the stability and folding relaxation rate of phosphoglycerate kinase (PGK) Förster resonance energy transfer (FRET) constructs localized in the nucleus or in the endoplasmic reticulum (ER) of eukaryotic cells. PGK has a more compact native state in the cellular compartments than in aqueous solution. Its native FRET signature is similar to that previously observed in a carbohydrate-crowding matrix, consistent with crowding being responsible for the compact native state of PGK in the cell. PGK folds through multiple states in vitro, but its folding kinetics is more two-state-like in the ER, so the folding mechanism can be modified by intracellular compartments. The nucleus increases PGK stability and folding rate over the cytoplasm and ER, even though the density of crowders in the nucleus is no greater than in the ER or cytoplasm. Nuclear folding kinetics (and to a lesser extent, thermodynamics) vary less from cell to cell than in the cytoplasm or ER, indicating a more homogeneous crowding and chemical environment in the nucleus.     

The Extracellular Protein VlsE Is Destabilized Inside Cells

We use U2OS cells as in vivo “test tubes” to study how the same cytoplasmic environment has opposite effects on the stability of two different proteins. Protein folding stability and kinetics were compared by fast relaxation imaging, which combines a temperature jump with fluorescence microscopy of FRET (Förster resonance energy transfer)-labeled proteins. While the stability of the cytoplasmic enzyme PGK (phosphoglycerate kinase) increases in cells, the stability of the cell surface antigen VlsE, which presumably did not evolve for stability inside cells, decreases. VlsE folding also slows down more than PGK folding in cells, relative to their respective aqueous buffer kinetics. Our FRET measurements provide evidence that VlsE is more compact inside cells than in aqueous buffer. Two kinetically distinct protein populations exist inside cells, making a connection with previous in vitro crowding studies. In addition, we confirm previous studies showing that VlsE is stabilized by 150 mg/mL of the carbohydrate crowder Ficoll, even though it is destabilized in the cytoplasm relative to aqueous buffer. We propose two mechanisms for the observed destabilization of VlsE in U2OS cells: long-range interactions competing with crowding or shape-dependent crowding favoring more compact states inside the cell over the elongated aqueous buffer native state.     

Quantification of Excluded Volume Effects on the Folding Landscape of Pseudomonas aeruginosa Apoazurin In Vitro

Proteins fold and function inside cells that are crowded with macromolecules. Here, we address the role of the resulting excluded volume effects by in vitro spectroscopic studies of Pseudomonas aeruginosa apoazurin stability (thermal and chemical perturbations) and folding kinetics (chemical perturbation) as a function of increasing levels of crowding agents dextran (sizes 20, 40, and 70 kDa) and Ficoll 70. We find that excluded volume theory derived by Minton quantitatively captures the experimental effects when crowding agents are modeled as arrays of rods. This finding demonstrates that synthetic crowding agents are useful for studies of excluded volume effects. Moreover, thermal and chemical perturbations result in free energy effects by the presence of crowding agents that are identical, which shows that the unfolded state is energetically the same regardless of method of unfolding. This also underscores the two-state approximation for apoazurin’s unfolding reaction and suggests that thermal and chemical unfolding experiments can be used in an interchangeable way. Finally, we observe increased folding speed and invariant unfolding speed for apoazurin in the presence of macromolecular crowding agents, a result that points to unfolded-state perturbations. Although the absolute magnitude of excluded volume effects on apoazurin is only on the order of 1–3 kJ/mol, differences of this scale may be biologically significant.     

Plasmin Triggers a Switch-Like Decrease in Thrombospondin-Dependent Activation of TGF-β1

Transforming growth factor-β1 (TGF-β1) is a potent regulator of extracellular matrix production, wound healing, differentiation, and immune response, and is implicated in the progression of fibrotic diseases and cancer. Extracellular activation of TGF-β1 from its latent form provides spatiotemporal control over TGF-β1 signaling, but the current understanding of TGF-β1 activation does not emphasize cross talk between activators. Plasmin (PLS) and thrombospondin-1 (TSP1) have been studied individually as activators of TGF-β1, and in this work we used a systems-level approach with mathematical modeling and in vitro experiments to study the interplay between PLS and TSP1 in TGF-β1 activation. Simulations and steady-state analysis predicted a switch-like bistable transition between two levels of active TGF-β1, with an inverse correlation between PLS and TSP1. In particular, the model predicted that increasing PLS breaks a TSP1-TGF-β1 positive feedback loop and causes an unexpected net decrease in TGF-β1 activation. To test these predictions in vitro, we treated rat hepatocytes and hepatic stellate cells with PLS, which caused proteolytic cleavage of TSP1 and decreased activation of TGF-β1. The TGF-β1 activation levels showed a cooperative dose response, and a test of hysteresis in the cocultured cells validated that TGF-β1 activation is bistable. We conclude that switch-like behavior arises from natural competition between two distinct modes of TGF-β1 activation: a TSP1-mediated mode of high activation and a PLS-mediated mode of low activation. This switch suggests an explanation for the unexpected effects of the plasminogen activation system on TGF-β1 in fibrotic diseases in vivo, as well as novel prognostic and therapeutic approaches for diseases with TGF-β dysregulation.     

Low temperature study of myoglobin-ligand rebinding kinetics with Mössbauer spectroscopy

We have studied the recombination kinetics of carboxymyoglobin (after photodissociation of the CO ligand) by Mössbauer spectroscopy for temperatures in the range 4.2 – 60 K. The observed kinetics display non-exponential behaviour which was monitored over periods of a few days. It is shown that the time dependence of the kinetics can be reduced to a single universal function of the temperature-dependent variable (t/τ _1/2(T)) ^β(T) . The half-decay time τ _1/2(T) and the scaling parameter β(T) are analysed for the presence of tunneling effects. The non-Arrhenius temperature dependence of the half-decay time below 60 K is interpreted as activated tunneling in models with an Eckart barrier or a fluctuating barrier.     

Flow-Induced β-Hairpin Folding of the Glycoprotein Ibα β-Switch

Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces β-hairpin folding of the glycoprotein Ibα β-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-μs molecular dynamics simulations, we investigate the kinetics of flow-induced β-hairpin folding. Simulations sampling different flow velocities reveal that under flow, β-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded β-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the β-switch has two stable conformations imprinted on it, namely, loop and hairpin—with flow inducing a transition between the two.     

Structural Rearrangements Linked to Global Folding Pathways of the Azoarcus Group I Ribozyme

Stable RNAs must fold into specific three-dimensional structures to be biologically active, yet many RNAs form metastable structures that compete with the native state. Our previous time-resolved footprinting experiments showed that Azoarcus group I ribozyme forms its tertiary structure rapidly (τ < 30 ms) without becoming significantly trapped in kinetic intermediates. Here, we use stopped-flow fluorescence spectroscopy to probe the global folding kinetics of a ribozyme containing 2-aminopurine in the loop of P9. The modified ribozyme was catalytically active and exhibited two equilibrium folding transitions centered at 0.3 and 1.6 mM Mg²⁺, consistent with previous results. Stopped-flow fluorescence revealed four kinetic folding transitions with observed rate constants of 100, 34, 1, and 0.1 s^− 1 at 37 °C. From comparison with time-resolved Fe(II)-ethylenediaminetetraacetic acid footprinting of the modified ribozyme under the same conditions, these folding transitions were assigned to formation of the I_C intermediate, tertiary folding and docking of the nicked P9 tetraloop, reorganization of the P3 pseudoknot, and refolding of nonnative conformers, respectively. The footprinting results show that 50-60% of the modified ribozyme folds in less than 30 ms, while the rest of the RNA population undergoes slow structural rearrangements that control the global folding rate. The results show how small perturbations to the structure of the RNA, such as a nick in P9, populate kinetic folding intermediates that are not observed in the natural ribozyme.     

Synchrotron x-ray microdiffraction reveals intrinsic structural features of amyloid deposits in situ

Amyloidoses are increasingly recognized as a major public health concern in Western countries. All amyloidoses share common morphological, structural, and tinctorial properties. These consist of staining by specific dyes, a fibrillar aspect in electron microscopy and a typical cross-β folding in x-ray diffraction patterns. Most studies that aim at deciphering the amyloid structure rely on fibers generated in vitro or extracted from tissues using protocols that may modify their intrinsic structure. Therefore, the fine details of the in situ architecture of the deposits remain unknown. Here, we present to our knowledge the first data obtained on ex vivo human renal tissue sections using x-ray microdiffraction. The typical cross-β features from fixed paraffin-embedded samples are similar to those formed in vitro or extracted from tissues. Moreover, the fiber orientation maps obtained across glomerular sections reveal an intrinsic texture that is correlated with the glomerulus morphology. These results are of the highest importance to understanding the formation of amyloid deposits and are thus expected to trigger new incentives for tissue investigation. Moreover, the access to intrinsic structural parameters such as fiber size and orientation using synchrotron x-ray microdiffraction, could provide valuable information concerning in situ mechanisms and deposit formation with potential benefits for diagnostic and therapeutic purposes.     

Hierarchic finite level energy landscape model: to describe the refolding kinetics of phosphoglycerate kinase

One of the most intriguing predictions of energy landscape models is the existence of non-exponential protein folding kinetics caused by hierarchical structures in the landscapes. Here we provide the strongest evidence so far of such hierarchy and determine the time constants and weights of the kinetic components of the suggested hierarchic energy landscape. To our knowledge, the idea of hierarchical folding energy barriers has never been tested over such a broad timescale. Refolding of yeast phosphoglycerate kinase was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 ms to 15 min. The strategy to build a model that describes folding of yeast phosphoglycerate kinase was to start from the simplest paradigm and modify it stepwise to the necessary minimal extent after repeated comparisons with the experiments. We made no a priori assumptions about the folding landscape. The result was a hierarchic finite level landscape model that quantitatively describes the refolding of yeast phosphoglycerate kinase from 1 ms to 15 min. The early steps of the folding process happen in the upper region of the landscape, where the surface has a hierarchic structure. This leads to stretched kinetics in the early phase of the folding. The lower region of the energy landscape is dominated by a trap that reflects the accumulation of molten globule intermediate state. From this intermediate, the protein can reach the global energy minimum corresponding to the native state through a cross-barrier folding step.     

Karyopherin-Centric Control of Nuclear Pores Based on Molecular Occupancy and Kinetic Analysis of Multivalent Binding with FG Nucleoporins

Intrinsically disordered Phe-Gly nucleoporins (FG Nups) within nuclear pore complexes exert multivalent interactions with transport receptors (Karyopherins (Kaps)) that orchestrate nucleocytoplasmic transport. Current FG-centric views reason that selective Kap translocation is promoted by alterations in the barrier-like FG Nup conformations. However, the strong binding of Kaps with the FG Nups due to avidity contradicts rapid Kap translocation in vivo. Here, using surface plasmon resonance, we innovate a means to correlate in situ mechanistic (molecular occupancy and conformational changes) with equilibrium (binding affinity) and kinetic (multivalent binding kinetics) aspects of Karyopherinβ1 (Kapβ1) binding to four different FG Nups. A general feature of the FxFG domains of Nup214, Nup62, and Nup153 is their capacity to extend and accommodate large numbers of Kapβ1 molecules at physiological Kapβ1 concentrations. A notable exception is the GLFG domain of Nup98, which forms a partially penetrable cohesive layer. Interestingly, we find that a slowly exchanging Kapβ1 phase forms an integral constituent within the FG Nups that coexists with a fast phase, which dominates transport kinetics due to limited binding with the pre-occupied FG Nups at physiological Kapβ1 concentrations. Altogether, our data reveal an emergent Kap-centric barrier mechanism that may underlie mechanistic and kinetic control in the nuclear pore complex.     

Simulations of cross-amyloid aggregation of amyloid-β and islet amyloid polypeptide fragments

Amyloid-β (Aβ) and islet amyloid polypeptide (IAPP) are small peptides, classified as amyloids, that have the potential to self-assemble and form cytotoxic species, such as small soluble oligomers and large insoluble fibrils. The formation of Aβ aggregates facilitates the progression of Alzheimer’s disease (AD), while IAPP aggregates induce pancreatic β-cell apoptosis, leading to exacerbation of type 2 diabetes (T2D). Cross-amyloid interactions between Aβ and IAPP have been described both in vivo and in vitro, implying the role of Aβ or IAPP as modulators of cytotoxic self-aggregation of each species, and suggesting that Aβ-IAPP interactions are a potential molecular link between AD and T2D. Using molecular dynamics (MD) simulations, “hotspot” regions of the two peptides were studied to understand the formation of hexamers in a heterogeneous and homogeneous peptide-containing environment. Systems of only Aβ_(16–22) peptides formed antiparallel, β-barrel-like structures, while systems of only IAPP_(20–29) peptides formed stacked, parallel β-sheets and had relatively unstable aggregation structures after 2 μs of simulation time. Systems containing both Aβ and IAPP (1:1 ratio) hexamers showed antiparallel, β-barrel-like structures, with an interdigitated arrangement of Aβ_(16–22) and IAPP_(20–29). These β-barrel structures have features of cytotoxic amyloid species identified in previous literature. Ultimately, this work seeks to provide atomistic insight into both the mechanism behind cross-amyloid interactions and structural morphologies of these toxic amyloid species.     

Macromolecular crowding modulates α-synuclein amyloid fiber growth

The crowdedness of living cells, hundreds of milligrams per milliliter of macromolecules, may affect protein folding, function, and misfolding. Still, such processes are most often studied in dilute solutions in vitro. To assess consequences of the in vivo milieu, we here investigated the effects of macromolecular crowding on the amyloid fiber formation reaction of α-synuclein, the amyloidogenic protein in Parkinson’s disease. For this, we performed spectroscopic experiments probing individual steps of the reaction as a function of the macromolecular crowding agent Ficoll70, which is an inert sucrose-based polymer that provides excluded-volume effects. The experiments were performed at neutral pH at quiescent conditions to avoid artifacts due to shaking and glass beads (typical conditions for α-synuclein), using amyloid fiber seeds to initiate reactions. We find that both primary nucleation and fiber elongation steps during α-synuclein amyloid formation are accelerated by the presence of 140 and 280 mg/mL Ficoll70. Moreover, in the presence of Ficoll70 at neutral pH, secondary nucleation appears favored, resulting in faster overall α-synuclein amyloid formation. In contrast, sucrose, a small-molecule osmolyte and building block of Ficoll70, slowed down α-synuclein amyloid formation. The ability of cell environments to modulate reaction kinetics to a large extent, such as severalfold faster individual steps in α-synuclein amyloid formation, is an important consideration for biochemical reactions in living systems.     

The effect of Zn2+ on Euphausia superba arginine kinase: Unfolding and aggregation studies

Arginine kinase plays an important role in the cellular energy metabolism of invertebrates. We investigated the effects of Zn²⁺ on the enzymatic activity and unfolding and aggregation of Euphausia superba arginine kinase (ESAK). Zn²⁺ inhibited the activity of ESAK (IC₅₀ = 0.027 ± 0.002 mM) following first-order kinetics consistent with the transition from a mono-phasic to a bi-phasic reaction. Double-reciprocal Lineweaver–Burk plots indicated that Zn²⁺ induced non-competitive inhibition of arginine and ATP. Circular dichroism spectra and spectrofluorometry results showed that Zn²⁺ induced secondary and tertiary structural changes in ESAK with exposure of hydrophobic surfaces and directly induced ESAK aggregation. The addition of osmolytes such as glycine and proline successfully blocked ESAK aggregation, recovering the conformation and activity of ESAK. Our study demonstrates the effect of Zn²⁺ on ESAK enzymatic function and folding and unfolding mechanisms, and might provide important insights into other metabolic enzymes of invertebrates in extreme climatic marine environments.     

Induced Secondary Structure and Polymorphism in an Intrinsically Disordered Structural Linker of the CNS: Solid-State NMR and FTIR Spectroscopy of Myelin Basic Protein Bound to Actin

The 18.5 kDa isoform of myelin basic protein (MBP) is a peripheral membrane protein that maintains the structural integrity of the myelin sheath of the central nervous system by conjoining the cytoplasmic leaflets of oligodendrocytes and by linking the myelin membrane to the underlying cytoskeleton whose assembly it strongly promotes. It is a multifunctional, intrinsically disordered protein that behaves primarily as a structural stabilizer, but with elements of a transient or induced secondary structure that represent binding sites for calmodulin or SH3-domain-containing proteins, inter alia. In this study we used solid-state NMR (SSNMR) and Fourier transform infrared (FTIR) spectroscopy to study the conformation of 18.5 kDa MBP in association with actin microfilaments and bundles. FTIR spectroscopy of fully ¹³C,¹⁵N-labeled MBP complexed with unlabeled F-actin showed induced folding of both protein partners, viz., some increase in β-sheet content in actin, and increases in both α-helix and β-sheet content in MBP, albeit with considerable extended structure remaining. Solid-state NMR spectroscopy revealed that MBP in MBP-actin assemblies is structurally heterogeneous but gains ordered secondary structure elements (both α-helical and β-sheet), particularly in the terminal fragments and in a central immunodominant epitope. The overall conformational polymorphism of MBP is consistent with its in vivo roles as both a linker (membranes and cytoskeleton) and a putative signaling hub.     

Beta-Amyloid Oligomers Activate Apoptotic BAK Pore for Cytochrome c Release

In Alzheimer’s disease, cytochrome c-dependent apoptosis is a crucial pathway in neuronal cell death. Although beta-amyloid (Aβ) oligomers are known to be the neurotoxins responsible for neuronal cell death, the underlying mechanisms remain largely elusive. Here, we report that the oligomeric form of synthetic Aβ of 42 amino acids elicits death of HT-22 cells. But, when expression of a bcl-2 family protein BAK is suppressed by siRNA, Aβ oligomer-induced cell death was reduced. Furthermore, significant reduction of cytochrome c release was observed with mitochondria isolated from BAK siRNA-treated HT-22 cells. Our in vitro experiments demonstrate that Aβ oligomers bind to BAK on the membrane and induce apoptotic BAK pores and cytochrome c release. Thus, the results suggest that Aβ oligomers function as apoptotic ligands and hijack the intrinsic apoptotic pathway to cause unintended neuronal cell death.     

Nonnative protein polymers: structure, morphology, and relation to nucleation and growth

Thermally induced aggregates of α-chymotrypsinogen A and bovine granulocyte-colony stimulating factor in acidic solutions were characterized by a combination of static and dynamic light scattering, spectroscopy, transmission electron microscopy, and monomer loss kinetics. The resulting soluble, high-molecular weight aggregates (∼10³-10⁵ kDa) are linear, semiflexible polymer chains that do not appreciably associate with one another under the conditions at which they were formed, with classic power-law scaling of the radius of gyration and hydrodynamic radius with weight-average molecular weight (M_w). Aggregates in both systems are composed of nonnative monomers with elevated levels of β-sheet secondary structure, and bind thioflavine T. In general, the aggregate size distributions showed low polydispersity by light scattering. Together with the inverse scaling of M_w with protein concentration, the results clearly indicate that aggregation proceeds via nucleated (chain) polymerization. For α-chymotrypsinogen A, the scaling behavior is combined with the kinetics of aggregation to deduce separate values for the characteristic timescales for nucleation (τ_n) and growth (τ_g), as well as the stoichiometry of the nucleus (x). The analysis illustrates a general procedure to noninvasively and quantitatively determine τ_n, τ_g, and x for soluble (chain polymer) aggregates, as well as the relationship between τ_n/τ_g and aggregate M_w.