Folding Nucleus Research Articles

Visualizing nucleation, condensation and propagation of β-tubulin folding in chaperonin TRiC.

The folding nucleus (FN) initiates protein folding and enables an efficient folding pathway. Here we directly visualize the tubulin FN consisting of a nonnative, partially assembled Rossmann fold, in the closed chamber of human chaperonin TRiC. Chaperonin TRiC interacts with non-natively folded secondary structural elements, stabilizing the nucleus for transition into its first native domain. Through progressive folding, the unfolded sequence goes through drastic spatial arrangement in the TRiC chamber to sample the conformational space, mediated by the highly dynamic CCT tails. The observed presence of individual nonnative secondary structures first in the nonnative FN and then around the incrementally folded native domains supports the hypothesis that tubulin folding in TRiC is a hierarchical process of nucleation, condensation and propagation in cooperation with TRiC subunits.

Proline Peptide Bond Isomerization in Ubiquitin Under Folding and Denaturing Conditions by Pressure-Jump NMR

Proline isomerization is widely recognized as a kinetic bottleneck in protein folding, amplified for proteins rich in Pro residues. We introduced repeated hydrostatic pressure jumps between native and pressure-denaturing conditions inside an NMR sample cell to study proline isomerization in the pressure-sensitized L50A ubiquitin mutant. Whereas in two unfolded heptapeptides, X-Pro peptide bonds isomerized ca 1.6-fold faster at 1 bar than at 2.5 kbar, for ubiquitin ca eight-fold faster isomerization was observed for Pro-38 and ca two-fold for Pro-19 and Pro-37 relative to rates measured in the pressure-denatured state. Activation energies for isomerization in pressure-denatured ubiquitin were close to literature values of 20 kcal/mole for denatured polypeptides but showed a substantial drop to 12.7 kcal/mole for Pro-38 at atmospheric pressure. For ubiquitin isomers with a cis E18-P19 peptide bond, the 1-bar NMR spectrum showed sharp resonances with near random coil chemical shifts for the C-terminal half of the protein, characteristic of an unfolded chain, while most of the N-terminal residues were invisible due to exchange broadening, pointing to a metastable partially folded state for this previously recognized ‘folding nucleus’. For cis-P37 isomers, a drop in pressure resulted in the rapid loss of nearly all unfolded-state NMR resonances, while the recovery of native state intensity revealed a slow component attributed to cis → trans isomerization of P37. This result implies that the NMR-invisible cis-P37 isomer adopts a molten globule state that encompasses the entire length of the ubiquitin chain, suggestive of a structure that mostly resembles the folded state.

Mechanistic Insights into the Folding Mechanism of Region V in Ice-Binding Protein Secreted by Marinomonas primoryensis Revealed by Single-Molecule Force Spectroscopy.

The Gram-negative bacteria Marinomonas primoryensis secrete an ice-binding protein (MpIBP), which is a vital bacterial adhesin facilitating the adaptation and survival of the bacteria in the harsh Antarctic environment. The C-terminal region of MpIBP, known as region V (RV), is the first domain to be exported into the Ca2+-rich extracellular environment and acts as a folding nucleus for the entire adhesin. However, the mechanisms underlying the secretion and folding of RV remain poorly understood. Here, we used optical tweezers (OT) to investigate the secretion and folding mechanisms of RV at the single-molecule level. In the absence of Ca2+, apo-RV remains unstructured, while Ca2+-bound RV folds into a mechanically stable structure. The folding of RV could occur via the formation of an intermediate state. Even though this folding intermediate is "hidden" during the folding process of wild type RV in vitro, it likely forms in vivo and plays an important role in facilitating protein secretion. Additionally, our results revealed that the N-terminal part of the RV can significantly stabilize its C-terminal structure. Our study paves the way for further investigations into the structure and functions of MpIBP that help bacteria survive in challenging environments.

Structural and functional analysis of actin point mutations leading to nemaline myopathy to elucidate their role in actin function.

In this work, we analyzed 78 mutations in the actin protein that cause the disease nemaline myopathy. We analyzed how these mutations are distributed in important regions of the actin molecule (folding nucleus, core of the filament, amyloidogenic regions, disordered regions, regions involved in interaction with other proteins). It was found that 54 mutations (43 residues) fall into the folding nucleus (Ф ≥ 0.5), 11 mutations (10 residues) into the filament core, 14 mutations into the amyloidogenic regions (11 residues), 14 mutations (9 residues) in the unstructured regions, and 24 mutations (22 residues) in regions involved in interaction with other proteins. It was also found that the occurrence of single mutations G44V, V45F, T68I, P72R, K338Iand S350L leads to the appearance of new amyloidogenic regions that are not present in native actin. The largest number of mutations (54 out of 78) occurs in the folding nucleus; these mutations are important for folding and therefore can affect the protein folding rate. We have shown that almost all of the considered mutations are associated with the structural characteristics of the actin molecule, and some of the residues we have considered have several important characteristics.

A short structural extension dictates the early stages of folding of a PDZ domain

PDZ domains are highly abundant protein-protein interaction modules in human. One of the most extensively characterized PDZ domain, the third PDZ domain from PSD-95 (PDZ3), contains an α-helical C-terminal extension that has a key role in the function of the domain. Here we compared the folding of PDZ3 with a truncated variant (PDZ3Δα3), lacking the additional helix, by means of the so-called Φ-value analysis, an experimental technique that allows inferring the structure of folding transition states. Experiments reveal subtle but detectable differences in the folding of PDZ3Δα3 versus PDZ3, as probed by structural characterization of the folding transition states. These differences appear more remarkable in the early stages of folding, with a detectable shift of the folding nucleus. The presented results allow demonstrating that the native state exerts a weak bias at the early stages of folding, which appear to be characterized by alternative pathways.

Characterization of early and late transition states of the folding pathway of a SH2 domain.

Albeit SH2 domains are abundant protein–protein interaction modules with fundamental roles in the regulation of several physiological and molecular pathways in the cell, the available information about the determinants of their thermodynamic stability and folding properties are still very limited. In this work, we provide a quantitative characterization of the folding pathway of the C‐terminal SH2 domain of SHP2, conducted through a combination of site‐directed mutagenesis and kinetic (un)folding experiments (Φ‐value analysis). The energetic profile of the folding reaction of the C‐SH2 domain is described by a three‐state mechanism characterized by the presence of two transition states and a high‐energy intermediate. The production of 29 site‐directed variants allowed us to calculate the degree of native‐like interactions occurring in the early and late events of the folding reaction. Data analysis highlights the presence of a hydrophobic folding nucleus surrounded by a lower degree of structure in the early events of folding, further consolidated as the reaction proceeds towards the native state. Interestingly, residues physically located in the functional region of the domain reported unusual Φ‐values, a hallmark of the presence of transient misfolding. We compared our results with previous ones obtained for the N‐terminal SH2 domain of SHP2. Notably, a conserved complex folding mechanism implying the presence of a folding intermediate arise from comparison, and the relative stability of such intermediate appears to be highly sequence dependent. Data are discussed under the light of previous works on SH2 domains.

Functionalization of a symmetric protein scaffold: Redundant folding nuclei and alternative oligomeric folding pathways.

Successful de novo protein design ideally targets specific folding kinetics, stability thermodynamics, and biochemical functionality, and the simultaneous achievement of all these criteria in a single step design is challenging. Protein design is potentially simplified by separating the problem into two steps: (a) an initial design of a protein "scaffold" having appropriate folding kinetics and stability thermodynamics, followed by (b) appropriate functional mutation-possibly involving insertion of a peptide functional "cassette." This stepwise approach can also separate the orthogonal effects of the "stability/function" and "foldability/function" tradeoffs commonly observed in protein design. If the scaffold is a protein architecture having an exact rotational symmetry, then there is the potential for redundant folding nuclei and multiple equivalent sites of functionalization; thereby enabling broader functional adaptation. We describe such a "scaffold" and functional "cassette" design strategy applied to a β-trefoil threefold symmetric architecture and a heparin ligand functionality. The results support the availability of redundant folding nuclei within this symmetric architecture, and also identify a minimal peptide cassette conferring heparin affinity. The results also identify an energy barrier of destabilization that switches the protein folding pathway from monomeric to trimeric, thereby identifying another potential advantage of symmetric protein architecture in de novo design.

Peptide Model of the Mutant Proinsulin Syndrome. I. Design and Clinical Correlation.

The mutant proinsulin syndrome is a monogenic cause of diabetes mellitus due to toxic misfolding of insulin’s biosynthetic precursor. Also designated mutant INS-gene induced diabetes of the young (MIDY), this syndrome defines molecular determinants of foldability in the endoplasmic reticulum (ER) of β-cells. Here, we describe a peptide model of a key proinsulin folding intermediate and variants containing representative clinical mutations; the latter perturb invariant core sites in native proinsulin (LeuB15→Pro, LeuA16→Pro, and PheB24→Ser). The studies exploited a 49-residue single-chain synthetic precursor (designated DesDi), previously shown to optimize in vitro efficiency of disulfide pairing. Parent and variant peptides contain a single disulfide bridge (cystine B19-A20) to provide a model of proinsulin’s first oxidative folding intermediate. The peptides were characterized by circular dichroism and redox stability in relation to effects of the mutations on (a) in vitro foldability of the corresponding insulin analogs and (b) ER stress induced in cell culture on expression of the corresponding variant proinsulins. Striking correlations were observed between peptide biophysical properties, degree of ER stress and age of diabetes onset (neonatal or adolescent). Our findings suggest that age of onset reflects the extent to which nascent structure is destabilized in proinsulin’s putative folding nucleus. We envisage that such peptide models will enable high-resolution structural studies of key folding determinants and in turn permit molecular dissection of phenotype-genotype relationships in this monogenic diabetes syndrome. Our companion study (next article in this issue) employs two-dimensional heteronuclear NMR spectroscopy to define site-specific perturbations in the variant peptides.

Structural characterization of the thermal unfolding pathway of human VEGFR1 D2 domain.

Folding stability is a crucial feature of protein evolution and is essential for protein functions. Thus, the comprehension of protein folding mechanisms represents an important complement to protein structure and function, crucial to determine the structural basis of protein misfolding. In this context, thermal unfolding studies represent a useful tool to get a molecular description of the conformational transitions governing the folding/unfolding equilibrium of a given protein. Here, we report the thermal folding/unfolding pathway of VEGFR1D2, a member of the immunoglobulin superfamily by means of a high‐resolution thermodynamic approach that combines differential scanning calorimetry with atomic‐level unfolding monitored by NMR. We show how VEGFR1D2 folding is driven by an oxidatively induced disulfide pairing: the key event in the achievement of its functional structure is the formation of a small hydrophobic core that surrounds a disulfide bridge. Such a ‘folding nucleus’ induces the cooperative transition to the properly folded conformation supporting the hypothesis that a disulfide bond can act as a folding nucleus that eases the folding process.

Rational Design of Protein-Specific Folding Modifiers.

Protein-folding can go wrong in vivo and in vitro, with significant consequences for the living organism and the pharmaceutical industry, respectively. Here we propose a design principle for small-peptide-based protein-specific folding modifiers. The principle is based on constructing a "xenonucleus", which is a prefolded peptide that mimics the folding nucleus of a protein. Using stopped-flow kinetics, NMR spectroscopy, Förster resonance energy transfer, single-molecule force measurements, and molecular dynamics simulations, we demonstrate that a xenonucleus can make the refolding of ubiquitin faster by 33 ± 5%, while variants of the same peptide have little or no effect. Our approach provides a novel method for constructing specific, genetically encodable folding catalysts for suitable proteins that have a well-defined contiguous folding nucleus.

The ubiquitous buried water in the beta-trefoil architecture contributes to the folding nucleus and ~20% of the folding enthalpy.

The beta-trefoil protein architecture is characterized by three repeating "trefoil" motifs related by rotational symmetry and postulated to have evolved via gene duplication and fusion events. Despite this apparent structural symmetry, the primary and secondary structural elements typically exhibit pronounced asymmetric features. A survey of this family of proteins has revealed that among the most conserved symmetric structural elements is a ubiquitous buried solvent which participates in a bridging H-bond with three different beta-strands in each of the trefoil motifs. A computational analysis reported that these waters are likely associated with a substantial enthalpic contribution to overall stability. In this report, a Pro mutation is used to disrupt one of the water H-bond interactions to a main chain amide, and the effects upon stability and folding kinetics are determined. Combined with Ala mutations, the separate effects upon side chain truncation and H-bond deletion are analyzed in terms of stability and folding kinetics. The results show that these buried waters act to assemble a central folding nucleus, and are responsible for ~20% of the overall favorable enthalpy of folding.

A Glimpse into the Structural Properties of the Intermediate and Transition State in the Folding of Bromodomain 2 Domain 2 by Φ Value Analysis.

Bromodomains (BRDs) are small protein interaction modules of about 110 amino acids that selectively recognize acetylated lysine in histones and other proteins. These domains have been identified in a variety of multi-domain proteins involved in transcriptional regulation or chromatin remodeling in eukaryotic cells. BRD inhibition is considered an attractive therapeutic approach in epigenetic disorders, particularly in oncology. Here, we present a Φ value analysis to investigate the folding pathway of the second domain of BRD2 (BRD2(2)). Using an extensive mutational analysis based on 25 site-directed mutants, we provide structural information on both the intermediate and late transition state of BRD2(2). The data reveal that the C-terminal region represents part of the initial folding nucleus, while the N-terminal region of the domain consolidates its structure only later in the folding process. Furthermore, only a small number of native-like interactions have been identified, suggesting the presence of a non-compact, partially folded state with scarce native-like characteristics. Taken together, these results indicate that, in BRD2(2), a hierarchical mechanism of protein folding can be described with non-native interactions that play a significant role in folding.

A conserved folding nucleus sculpts the free energy landscape of bacterial and archaeal orthologs from a divergent TIM barrel family

The amino acid sequences of proteins have evolved over billions of years, preserving their structures and functions while responding to evolutionary forces. Are there conserved sequence and structural elements that preserve the protein folding mechanisms? The functionally diverse and ancient (βα)1-8 TIM barrel motif may answer this question. We mapped the complex six-state folding free energy surface of a ∼3.6 billion y old, bacterial indole-3-glycerol phosphate synthase (IGPS) TIM barrel enzyme by equilibrium and kinetic hydrogen-deuterium exchange mass spectrometry (HDX-MS). HDX-MS on the intact protein reported exchange in the native basin and the presence of two thermodynamically distinct on- and off-pathway intermediates in slow but dynamic equilibrium with each other. Proteolysis revealed protection in a small (α1β2) and a large cluster (β5α5β6α6β7) and that these clusters form cores of stability in Ia and Ibp The strongest protection in both states resides in β4α4 with the highest density of branched aliphatic side chain contacts in the folded structure. Similar correlations were observed previously for an evolutionarily distinct archaeal IGPS, emphasizing a key role for hydrophobicity in stabilizing common high-energy folding intermediates. A bioinformatics analysis of IGPS sequences from the three superkingdoms revealed an exceedingly high hydrophobicity and surprising α-helix propensity for β4, preceded by a highly conserved βα-hairpin clamp that links β3 and β4. The conservation of the folding mechanisms for archaeal and bacterial IGPS proteins reflects the conservation of key elements of sequence and structure that first appeared in the last universal common ancestor of these ancient proteins.

Folding nucleus and unfolding dynamics of protein 2GB1* *Project supported by the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDA17010504) and the National Natural Science Foundation of China (Grant No. 11947302).

The folding of many small proteins is kinetically a two-state process with one major free-energy barrier to overcome, which can be roughly regarded as the inverse process of unfolding. In this work, we first use a Gaussian network model to predict the folding nucleus corresponding to the major free-energy barrier of protein 2GB1, and find that the folding nucleus is located in the β-sheet domain. High-temperature molecular dynamics simulations are then used to investigate the unfolding process of 2GB1. We draw free-energy surface from unfolding simulations, taking RMSD and contact number as reaction coordinates, which confirms that the folding of 2GB1 is kinetically a two-state process. The comparison of the contact maps before and after the free energy barrier indicates that the transition from native to non-native structure of the protein is kinetically caused by the destruction of the β-sheet domain, which manifests that the folding nucleus is indeed located in the β-sheet domain. Moreover, the constrained MD simulation further confirms that the destruction of the secondary structures does not alter the topology of the protein retained by the folding nucleus. These results provide vital information for upcoming researchers to further understand protein folding in similar systems.

Conserved buried water molecules enable the β-trefoil architecture.

Available high-resolution crystal structures for the family of β-trefoil proteins in the structural databank were queried for buried waters. Such waters were classified as either: (a) unique to a particular domain, family, or superfamily or (b) conserved among all β-trefoil folds. Three buried waters conserved among all β-trefoil folds were identified. These waters are related by the threefold rotational pseudosymmetry characteristic of this protein architecture (representing three instances of an identical structural environment within each repeating trefoil-fold motif). The structural properties of this buried water are remarkable and include: residing in a cavity space no larger than a single water molecule, exhibiting a positional uncertainty (i.e., normalized B-factor) substantially lower than the average Cα atom, providing essentially ideal H-bonding geometry with three solvent-inaccessible main chain groups, simultaneously serving as a bridging H-bond for three different β-strands at a point of secondary structure divergence, and orienting conserved hydrophobic side chains to form a nascent core-packing group. Other published work supports an interpretation that these interactions are key to the formation of an efficient folding nucleus and folded thermostability. The fundamental threefold symmetric structural element of the β-trefoil fold is therefore, surprisingly, a buried water molecule.

Oligomerization of a symmetric β-trefoil protein in response to folding nucleus perturbation.

Gene duplication and fusion events in protein evolution are postulated to be responsible for the common protein folds exhibiting internal rotational symmetry. Such evolutionary processes can also potentially yield regions of repetitive primary structure. Repetitive primary structure offers the potential for alternative definitions of critical regions, such as the folding nucleus (FN). In principle, more than one instance of the FN potentially enables an alternative folding pathway in the face of a subsequent deleterious mutation. We describe the targeted mutation of the carboxyl-terminal region of the (internally located) FN of the de novo designed purely-symmetric β-trefoil protein Symfoil-4P. This mutation involves wholesale replacement of a repeating trefoil-fold motif with a "blade" motif from a β-propeller protein, and postulated to trap that region of the Symfoil-4P FN in a nonproductive folding intermediate. The resulting protein (termed "Bladefoil") is shown to be cooperatively folding, but as a trimeric oligomer. The results illustrate how symmetric protein architectures have potentially diverse folding alternatives available to them, including oligomerization, when preferred pathways are perturbed.

Do folding elements trade‐off with function in the human mitochondrial metabolite transporter?

In humans, transmembrane β‐barrels are found exclusively in the outer mitochondrial membrane (OMM). They comprise metabolite transporters, translocases, and chaperones, which together control cellular homeostasis. The metabolite transporters — known as voltage‐dependent anion channels (VDACs) – are the most abundant OMM barrels. They are 19‐stranded β‐barrel structures that also form important pharmacological targets due to their roles in metabolite flux, steroidogenesis, gametogenesis, and ROS regulation. Despite long‐standing efforts, molecular factors that regulate VDAC function and folding at the atomistic level are not known. Here, we map the folding nucleus and the assembly pathway of human VDAC isoform 2 in phosphocholine vesicles. We carry out a comprehensive measurement of hVDAC2 folding kinetics, and equilibrium thermodynamic contribution, of each of the 270 non‐Ala residues of the 294‐residue barrel. Interestingly, we obtain a two‐state thermodynamic equilibrium that is attained through two detectable on‐pathway early folding intermediates, which are assembled in milliseconds – seconds. These intermediates undergo slow conformational restructuring to form the folded VDAC2 barrel. We find that specific charged residues that form a part of these early intermediates also positively regulate the gating characteristics of VDAC2. This observation is unprecedented, as directed evolution of protein sequences in favor of function usually place a stability burden on its structure that additionally impedes protein assembly. Our findings provide molecular insight on the assembly pathway of human VDAC2, and directly link these folding elements with VDAC2 function as the metabolite transporter of human mitochondria.Support or Funding InformationThis work was supported by the Wellcome Trust/DBT India Alliance Fellowship [grant number IA/I/14/1/501305] awarded to RM.

Linking Folding Landscape with Function in the Human Mitochondrial VDAC2

Transmembrane β-barrels in humans are exclusive to the mitochondrial outer membrane (OMM). They include translocases, porins, and chaperones, which together regulate cellular homeostasis. The most abundant OMM proteins are the voltage-dependent anion channels (VDACs), which are 19-stranded multifunctional metabolite flux porins that are also important pharmacological targets. In humans, VDAC2 is anti-apoptotic and functions additionally in metabolite flux, steroidogenesis, gametogenesis, and ROS regulation, making this isoform indispensable for the cell. Despite long-standing efforts, molecular elements regulating VDAC function and folding are not known. Here, we map the folding nucleus and the assembly pathway of human VDAC2 in phosphocholine vesicles. We carry out a comprehensive measurement of folding kinetics and equilibrium thermodynamic contribution of each of the 270 non-Ala residues in the 294-residue barrel. Interestingly, we obtain a two-state thermodynamic equilibrium that is achieved through two detectable on-pathway early folding intermediates, which drive barrel assembly. Detailed analysis of the transition state structure reveals that specific strands at the N-terminal and central region of the barrel assemble first, which is followed by structuring of the C-terminal region. Residues causing frustration in the folding (Φ-value < 0) are more populated towards the C-terminal strands. Furthermore, specific charged residues that regulate the gating characteristics of human VDAC2 are surprisingly a part of the folding nucleus. The only interesting exception is D169, which is functionally important in vivo for cell survival, but establishes non-native contacts during VDAC2 folding. Relating our findings with the proposed mechanism of chaperone–assisted folding in the OMM provides molecular insight that links the assembly pathway of human VDAC2 with function.

Shared unfolding pathways of unrelated immunoglobulin-like β-sandwich proteins.

The Dynameomics project contains native state and unfolding simulations of 807 protein domains, where each domain is representative of a different metafold; these metafolds encompass ~97% of protein fold space. There is a long-standing question in structural biology as to whether proteins in the same fold family share the same folding/unfolding characteristics. Using molecular dynamics simulations from the Dynameomics project, we conducted a detailed study of protein unfolding/folding pathways for 5 protein domains from the immunoglobulin (Ig)-like β-sandwich metafold (the highest ranked metafold in our database). The domains have sequence similarities ranging from 4 to 15% and are all from different SCOP superfamilies, yet they share the same overall Ig-like topology. Despite having very different amino acid sequences, the dominant unfolding pathway is very similar for the 5 proteins, and the secondary structures that are peripheral to the aligned, shared core domain add variability to the unfolding pathway. Aligned residues in the core domain display consensus structure in the transition state primarily through conservation of hydrophobic positions. Commonalities in the obligate folding nucleus indicate that insights into the major events in the folding/unfolding of other domains from this metafold may be obtainable from unfolding simulations of a few representative proteins.

Ubiquitin folds via a flip-twist-lock mechanism

To perform specific functional activities, the majority of proteins should fold into their distinct three-dimensional conformations. However, the biologically active conformation of a protein is generally found to be marginally stable than the other conformations that the chain can adopt. How a protein finds its native conformation from its post-synthesis unfolded structure in a complex conformational landscape is the unsolved question that still drives the protein folding community. Here, we report the folding mechanism of a globular protein, ubiquitin, from its chemically denatured state using all-atom molecular dynamics simulations. From the kinetic analysis of the simulated trajectories we show that the folding process can be described by the hydrophobic collapse mechanism, initiated by the “dewetting transition”, and subsequently assisted by the origination of an N-terminal folding nucleus, and finally supported by a native salt-bridge interaction between K11 and E34. We show that ubiquitin folds via an intermediate. Finally, we confirm the presence of “biological water” and explain its role to the folding process.