Backbone NMR assignments of a topologically knotted protein in urea-denatured state

YbeA is a 3-methylpseudoridine methyltransferase from Escherichia coli that forms a stable homodimer in solution. It is one of the deeply trefoil 3₁ knotted proteins, of which the knot encompasses the C-terminal helix that threads through a long loop. Recent studies on the knotted protein folding pathways using YbeA have suggested that the protein knot remains present under chemically denaturing conditions. Here, we report ¹H, ¹³C and ¹⁵N chemical shift assignments for urea-denatured YbeA, which will serve as the basis for further structural characterisations using solution state NMR spectroscopy with paramagnetic spin labeled and partial alignment media.     

How Difficult Is It to Fold a Knotted Protein? In Silico Insights from Surface-Tethered Folding Experiments

We explore the effect of surface tethering on the folding process of a lattice protein that contains a trefoil knot in its native structure via Monte Carlo simulations. We show that the outcome of the tethering experiment depends critically on which terminus is used to link the protein to a chemically inert plane. In particular, if surface tethering occurs at the bead that is closer to the knotted core the folding rate becomes exceedingly slow and the protein is not able to find the native structure in all the attempted folding trajectories. Such low folding efficiency is also apparent from the analysis of the probability of knot formation, p_knot, as a function of nativeness. Indeed, p_knot increases abruptly from ∼0 to ∼1 only when the protein has more than 80% of its native contacts formed, showing that a highly compact conformation must undergo substantial structural re-arrangement in order to get effectively knotted. When the protein is surface tethered by the bead that is placed more far away from the knotted core p_knot is higher than in the other folding setups (including folding in the bulk), especially if conformations are highly native-like. These results show that the mobility of the terminus closest to the knotted core is critical for successful folding of trefoil proteins, which, in turn, highlights the importance of a knotting mechanism that is based on a threading movement of this terminus through a knotting loop. The results reported here predict that if this movement is blocked, knotting occurs via an alternative mechanism, the so-called spindle mechanism, which is prone to misfolding. Our simulations show that in the three considered folding setups the formation of the knot is typically a late event in the folding process. We discuss the implications of our findings for co-translational folding of knotted trefoils.     

Unique Optimal Foldings of Proteins on a Triangular Lattice

Background: A problem for unique protein folding was raised in 1998: are there proteins having unique optimal foldings for all lengths in the hydrophobic-hydrophilic (hydrophobic-polar; HP) model? To such a question, it was proved that on a square lattice there are (i) closed chains of monomers having unique optimal foldings for all even lengths and (ii) open monomer chains having unique optimal foldings for all lengths divisible by four. In this article, we aim to extend the previous work on a square lattice to the optimal foldings of proteins on a triangular lattice by examining the uniqueness property or stability of HP chain folding. Method: We consider this protein folding problem on a triangular lattice using graph theory. For an HP chain with length n > 13, generally it is very time-consuming to enumerate all of its possible folding conformations. Hence, one can hardly know whether or not it has a unique optimal folding. A natural problem is to determine for what value of n there is an n-node HP chain that has a unique optimal folding on a triangular lattice. Results and conclusion: Using graph theory, this article proves that there are both closed and open chains having unique optimal foldings for all lengths >19 in a triangular lattice. This result is not only general from the theoretical viewpoint, but also can be expected to apply to areas of protein structure prediction and protein design because of their close relationship with the concept of energy state and designability.     

The Unique Cysteine Knot Regulates the Pleotropic Hormone Leptin

Leptin plays a key role in regulating energy intake/expenditure, metabolism and hypertension. It folds into a four-helix bundle that binds to the extracellular receptor to initiate signaling. Our work on leptin revealed a hidden complexity in the formation of a previously un-described, cysteine-knotted topology in leptin. We hypothesized that this unique topology could offer new mechanisms in regulating the protein activity. A combination of in silico simulation and in vitro experiments was used to probe the role of the knotted topology introduced by the disulphide-bridge on leptin folding and function. Our results surprisingly show that the free energy landscape is conserved between knotted and unknotted protein, however the additional complexity added by the knot formation is structurally important. Native state analyses led to the discovery that the disulphide-bond plays an important role in receptor binding and thus mediate biological activity by local motions on distal receptor-binding sites, far removed from the disulphide-bridge. Thus, the disulphide-bridge appears to function as a point of tension that allows dissipation of stress at a distance in leptin.     

Protein knots provide mechano-resilience to an AAA+ protease-mediated proteolysis with profound ATP energy expenses

Knotted proteins are some of the most fascinating examples of how linear polypeptide chains can achieve intricate topological arrangements efficiently and spontaneously. The entanglements of polypeptide chains could potentially enhance their folding stabilities. We recently reported the unprecedented mechanostability of the Gordian (5₂) knotted family of human ubiquitin C-terminal hydrolases (UCHs) in the context of withstanding the mechanical unfolding of the bacterial AAA+ proteasome, ClpXP; a green fluorescence protein (GFP) was fused to the N-terminus of various UCHs as a reporter of the unfolding and degradation of these topologically knotted substrates, but it also limited the ability to examine the effect of untying the knotted topology via N-terminal truncation. In this study, we directly monitored the ClpXP-mediated degradation of UCH variants by electrophoresis and quantitative imaging analyses. We demonstrated that untying of the 5₂ knot in UCHL1 via N-terminal truncation (UCHL1_Δ11) significantly reduces its mechanostability. We further quantified the ATP expenditures of degrading different UCH variants by ClpXP. The unknotted UCHL1_Δ11 underwent accelerated ClpXP-dependent proteolysis, with a 30-fold reduction in ATP consumption compared to the knotted wild type. Unlike all other known ClpXP substrates, UCHL5, which is the most resilient substrate known to date, significantly slowed down the ATP turnover rate by ClpXP. Furthermore, UCHL5 required 1000-fold more ATP to be fully degraded by ClpXP compared to GFP. Our results underscored how the complex, knotted folding topology in UCHs may interfere with the mechano-unfolding processes of the AAA+ unfoldase, ClpX.     

Fragmentation Follows Structure: Top-Down Mass Spectrometry Elucidates the Topology of Engineered Cystine-Knot Miniproteins

Over the last decades the field of pharmaceutically relevant peptides has enormously expanded. Among them, several peptide families exist that contain three or more disulfide bonds. In this context, elucidation of the disulfide patterns is extremely important as these motifs are often prerequisites for folding, stability, and activity. An example of this structure-determining pattern is a cystine knot which comprises three constrained disulfide bonds and represents a core element in a vast number of mechanically interlocked peptidic structures possessing different biological activities. Herein, we present our studies on disulfide pattern determination and structure elucidation of cystine-knot miniproteins derived from Momordica cochinchinensis peptide MCoTI-II, which act as potent inhibitors of human matriptase-1. A top-down mass spectrometric analysis of the oxidised and bioactive peptides is described. Following the detailed sequencing of the peptide backbone, interpretation of the MS³ spectra allowed for the verification of the knotted topology of the examined miniproteins. Moreover, we found that the fragmentation pattern depends on the knottin’s folding state, hence, tertiary structure, which to our knowledge has not been described for a top-down MS approach before.     

Characterizing a partially ordered miniprotein through folding molecular dynamics simulations: Comparison with the experimental data

The villin headpiece helical subdomain (HP36) is one of the best known model systems for computational studies of fast‐folding all‐α miniproteins. HP21 is a peptide fragment—derived from HP36—comprising only the first and second helices of the full domain. Experimental studies showed that although HP21 is mostly unfolded in solution, it does maintain some persistent native‐like structure as indicated by the analysis of NMR‐derived chemical shifts. Here we compare the experimental data for HP21 with the results obtained from a 15‐μs long folding molecular dynamics simulation performed in explicit water and with full electrostatics. We find that the simulation is in good agreement with the experiment and faithfully reproduces the major experimental findings, namely that (a) HP21 is disordered in solution with <10% of the trajectory corresponding to transiently stable structures, (b) the most highly populated conformer is a native‐like structure with an RMSD from the corresponding portion of the HP36 crystal structure of <1 Å, (c) the simulation‐derived chemical shifts—over the whole length of the trajectory—are in reasonable agreement with the experiment giving reduced χ² values of 1.6, 1.4, and 0.8 for the Δδ¹³C^α, Δδ¹³CO, and Δδ¹³C^β secondary shifts, respectively (becoming 0.8, 0.7, and 0.3 when only the major peptide conformer is considered), and finally, (d) the secondary structure propensity scores are in very good agreement with the experiment and clearly indicate the higher stability of the first helix. We conclude that folding molecular dynamics simulations can be a useful tool for the structural characterization of even marginally stable peptides.     

Distinct unfolding and refolding pathways of ribonuclease a revealed by heating and cooling temperature jumps

Heating and cooling temperature jumps (T-jumps) were performed using a newly developed technique to trigger unfolding and refolding of wild-type ribonuclease A and a tryptophan-containing variant (Y115W). From the linear Arrhenius plots of the microscopic folding and unfolding rate constants, activation enthalpy (ΔH^#), and activation entropy (ΔS^#) were determined to characterize the kinetic transition states (TS) for the unfolding and refolding reactions. The single TS of the wild-type protein was split into three for the Y115W variant. Two of these transition states, TS1 and TS2, characterize a slow kinetic phase, and one, TS3, a fast phase. Heating T-jumps induced protein unfolding via TS2 and TS3; cooling T-jumps induced refolding via TS1 and TS3. The observed speed of the fast phase increased at lower temperature, due to a strongly negative ΔH^# of the folding-rate constant. The results are consistent with a path-dependent protein folding/unfolding mechanism. TS1 and TS2 are likely to reflect X-Pro¹¹⁴ isomerization in the folded and unfolded protein, respectively, and TS3 the local conformational change of the β-hairpin comprising Trp¹¹⁵. A very fast protein folding/unfolding phase appears to precede both processes. The path dependence of the observed kinetics is suggestive of a rugged energy protein folding funnel.     

Folding studies on a knotted protein

YibK is a 160 residue homodimeric protein belonging to the SPOUT class of methyltransferases. Proteins in this group all display a unique topological feature; the backbone polypeptide chain folds to form a deep trefoil knot. Such knotted structures were completely unpredicted, it being thought impossible for a protein to fold efficiently in this way. However, they are becoming more common and there are now a growing number of examples in the Protein Data Bank. These intriguing knotted structures represent a new and significant challenge in the field of protein folding. Here, we present an initial characterisation of the folding of YibK, one of the smallest knotted proteins to be identified. This is the first detailed folding study on a knotted protein to be reported. We have established conditions under which the protein can be denatured reversibly in vitro using urea, thereby showing that molecular chaperones are not required for the efficient folding of this protein. A series of equilibrium unfolding experiments were performed over a 400-fold range of protein concentration. Both secondary and tertiary structural probes show a single, protein concentration-dependent unfolding transition, and data are most consistent with a three-state equilibrium denaturation model involving a monomeric intermediate. Thermodynamic parameters obtained from the fit of the data to this model indicate that the intermediate is a stable species with appreciable secondary and tertiary structure; whether the topological knot remains in the intermediate state is still to be shown. Together, these results demonstrate that, despite its complex knotted structure, YibK is able to fold efficiently and behaves remarkably similarly to other dimeric proteins under equilibrium conditions.     

Tackling Force-Field Bias in Protein Folding Simulations: Folding of Villin HP35 and Pin WW Domains in Explicit Water

The ability to fold proteins on a computer has highlighted the fact that existing force fields tend to be biased toward a particular type of secondary structure. Consequently, force fields for folding simulations are often chosen according to the native structure, implying that they are not truly “transferable.” Here we show that, while the AMBER ff03 potential is known to favor helical structures, a simple correction to the backbone potential (ff03^∗) results in an unbiased energy function. We take as examples the 35-residue α-helical Villin HP35 and 37 residue β-sheet Pin WW domains, which had not previously been folded with the same force field. Starting from unfolded configurations, simulations of both proteins in Amber ff03^∗ in explicit solvent fold to within 2.0 Å RMSD of the experimental structures. This demonstrates that a simple backbone correction results in a more transferable force field, an important requirement if simulations are to be used to interpret folding mechanism.     

Equilibrium folding of pro-HlyA from Escherichia coli reveals a stable calcium ion dependent folding intermediate

HlyA from Escherichia coli is a member of the repeats in toxin (RTX) protein family, produced by a wide range of Gram-negative bacteria and secreted by a dedicated Type 1 Secretion System (T1SS). RTX proteins are thought to be secreted in an unfolded conformation and to fold upon secretion by Ca^2 + binding. However, the exact mechanism of secretion, ion binding and folding to the correct native state remains largely unknown. In this study we provide an easy protocol for high-level pro-HlyA purification from E. coli. Equilibrium folding studies, using intrinsic tryptophan fluorescence, revealed the well-known fact that Ca^2 + is essential for stability as well as correct folding of the whole protein. In the absence of Ca^2 +, pro-HlyA adopts a non-native conformation. Such molecules could however be rescued by Ca^2 + addition, indicating that these are not dead-end species and that Ca^2 + drives pro-HlyA folding. More importantly, pro-HlyA unfolded via a two-state mechanism, whereas folding was a three-state process. The latter is indicative of the presence of a stable folding intermediate. Analysis of deletion and Trp mutants revealed that the first folding transition, at 6–7 M urea, relates to Ca^2 + dependent structural changes at the extreme C-terminus of pro-HlyA, sensed exclusively by Trp914. Since all Trp residues of HlyA are located outside the RTX domain, our results demonstrate that Ca^2 + induced folding is not restricted to the RTX domain. Taken together, Ca^2 + binding to the pro-HlyA RTX domain is required to drive the folding of the entire protein to its native conformation.     

HP1α mediates defective heterochromatin repair and accelerates senescence in Zmpste24-deficient cells

Heterochromatin protein 1 (HP1) interacts with various proteins, including lamins, to play versatile functions within nuclei, such as chromatin remodeling and DNA repair. Accumulation of prelamin A leads to misshapen nuclei, heterochromatin disorganization, genomic instability, and premature aging in Zmpste24-null mice. Here, we investigated the effects of prelamin A on HP1α homeostasis, subcellular distribution, phosphorylation, and their contribution to accelerated senescence in mouse embryonic fibroblasts (MEFs) derived from Zmpste24^−/− mice. The results showed that the level of HP1α was significantly increased in Zmpste24^−/− cells. Although prelamin A interacted with HP1α in a manner similar to lamin A, HP1α associated with the nuclease-resistant nuclear matrix fraction was remarkably increased in Zmpste24^−/− MEFs compared with that in wild-type littermate controls. In wild-type cells, HP1α was phosphorylated at Thr50, and the phosphorylation was maximized around 30 min, gradually dispersed 2 h after DNA damage induced by camptothecin. However, the peak of HP1α phosphorylation was significantly compromised and appeared until 2 h, which is correlated with the delayed maximal formation of γ-H2AX foci in Zmpste24^−/− MEFs. Furthermore, knocking down HP1α by siRNA alleviated the delayed DNA damage response and accelerated senescence in Zmpste24^−/− MEFs, evidenced by the rescue of the delayed γ-H2AX foci formation, downregulation of p16, and reduction of senescence-associated β-galactosidase activity. Taken together, these findings establish a functional link between prelamin A, HP1α, chromatin remodeling, DNA repair, and early senescence in Zmpste24-deficient mice, suggesting a potential therapeutic strategy for laminopathy-based premature aging via the intervention of HP1α.     

Folding network of villin headpiece subdomain

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 μs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.     

Biophysical characterisation of the small ankyrin repeat protein myotrophin

The 118 residue protein myotrophin is composed of four ankyrin repeats that stack linearly to form an elongated, predominantly α-helical structure. The protein folds via a two-state mechanism at equilibrium. The free energy change of unfolding in water (ΔG_U-N^H₂O) is 5.8 kcal.mol⁻¹. The chevron plot reveals that the folding reaction has a broad energy barrier and that it conforms to a two-state mechanism. The rate of folding in water (k_f^H₂O) of 95 s⁻¹ is several orders of magnitude slower than the value predicted by topological calculations. Proline mutants were used to show that the minor kinetic phases observed for myotrophin arise from heterogeneity of the ground states due to cis-trans isomerisation of prolyl as well as non-prolyl peptide bonds. Myotrophin is the first example of a naturally occurring ankyrin repeat protein that conforms to an apparent two-state mechanism at equilibrium and under kinetic conditions, making it highly suitable for high resolution protein folding studies.     

Protein folding rates and thermodynamic stability are key determinants for interaction with the Hsp70 chaperone system

The Hsp70 family of molecular chaperones participates in vital cellular processes including the heat shock response and protein homeostasis. E. coli''s Hsp70, known as DnaK, works in concert with the DnaJ and GrpE co-chaperones (K/J/E chaperone system), and mediates cotranslational and post-translational protein folding in the cytoplasm. While the role of the K/J/E chaperones is well understood in the presence of large substrates unable to fold independently, it is not known if and how K/J/E modulates the folding of smaller proteins able to fold even in the absence of chaperones. Here, we combine experiments and computation to evaluate the significance of kinetic partitioning as a model to describe the interplay between protein folding and binding to the K/J/E chaperone system. First, we target three nonobligatory substrates, that is, proteins that do not require chaperones to fold. The experimentally observed chaperone association of these client proteins during folding is entirely consistent with predictions from kinetic partitioning. Next, we develop and validate a computational model (CHAMP70) that assumes kinetic partitioning of substrates between folding and interaction with K/J/E. CHAMP70 quantitatively predicts the experimentally measured interaction of RNase H^D as it refolds in the presence of various chaperones. CHAMP70 shows that substrates are posed to interact with K/J/E only if they are slow-folding proteins with a folding rate constant k_f <50 s⁻¹, and/or thermodynamically unstable proteins with a folding free energy ΔG⁰_UN ≥−2 kcal mol⁻¹. Hence, the K/J/E system is tuned to use specific protein folding rates and thermodynamic stabilities as substrate selection criteria.     

Using SecM Arrest Sequence as a Tool to Isolate Ribosome Bound Polypeptides

Extensive research has provided ample evidences suggesting that protein folding in the cell is a co-translational process^1-5. However, the exact pathway that polypeptide chain follows during co-translational folding to achieve its functional form is still an enigma. In order to understand this process and to determine the exact conformation of the co-translational folding intermediates, it is essential to develop techniques that allow the isolation of RNCs carrying nascent chains of predetermined sizes to allow their further structural analysis.SecM (secretion monitor) is a 170 amino acid E. coli protein that regulates expression of the downstream SecA (secretion driving) ATPase in the secM-secA operon⁶. Nakatogawa and Ito originally found that a 17 amino acid long sequence (150-FSTPVWISQAQGIRAGP-166) in the C-terminal region of the SecM protein is sufficient and necessary to cause stalling of SecM elongation at Gly165, thereby producing peptidyl-glycyl-tRNA stably bound to the ribosomal P-site^7-9. More importantly, it was found that this 17 amino acid long sequence can be fused to the C-terminus of virtually any full-length and/or truncated protein thus allowing the production of RNCs carrying nascent chains of predetermined sizes⁷. Thus, when fused or inserted into the target protein, SecM stalling sequence produces arrest of the polypeptide chain elongation and generates stable RNCs both in vivo in E. coli cells and in vitro in a cell-free system. Sucrose gradient centrifugation is further utilized to isolate RNCs.The isolated RNCs can be used to analyze structural and functional features of the co-translational folding intermediates. Recently, this technique has been successfully used to gain insights into the structure of several ribosome bound nascent chains^10,11. Here we describe the isolation of bovine Gamma-B Crystallin RNCs fused to SecM and generated in an in vitro translation system.