A comparison of the folding of two knotted proteins: YbeA and YibK

The extraordinary topology of proteins belonging to the alpha/beta-knot superfamily of proteins is unexpected, due to the apparent complexities involved in the formation of a deep trefoil knot in a polypeptide backbone. Despite this, an increasing number of knotted structures are being identified; how such proteins fold remains a mystery. Studies on the dimeric protein YibK from Haemophilus influenzae have led to the characterisation of its folding pathway in some detail. To complement research into the folding of YibK, and to address whether folding pathways are conserved for members of the alpha/beta-knot superfamily, the structurally similar knotted protein YbeA from Escherichia coli has been studied. A comprehensive thermodynamic and kinetic analysis of the folding of YbeA is presented here, and compared to that of YibK. Both fold via an intermediate state populated under equilibrium conditions that is monomeric and considerably structured. The unfolding/refolding kinetics of YbeA are simpler than those found for YibK and involve two phases attributed to the formation of a monomeric intermediate state and a dimerisation step. In contrast to YibK, a change in the rate-determining step on the unfolding pathway for YbeA is observed with a changing concentration of urea. Despite this difference, both proteins fold by a mechanism involving at least one sequential monomeric intermediate that has properties similar to that observed during the equilibrium unfolding. The rate of dimerisation observed for YbeA and YibK is very similar, as is the rate constant for formation of the kinetic monomeric intermediate that precedes dimerisation. The findings suggest that relatively slow folding and dimerisation may be common attributes of knotted proteins.     

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.     

Probing nature's knots: the folding pathway of a knotted homodimeric protein

The homodimeric protein YibK from Haemophilus influenzae belongs to a recently discovered superfamily of knotted proteins that has brought about a new protein-folding conundrum. Members of the alpha/beta-knot clan form deep trefoil knots in their native backbone structure, a topological feature that is currently unexplained in the protein-folding field. To help solve the puzzle of how a polypeptide chain can efficiently knot itself, the folding kinetics of YibK have been studied extensively and the results are reported here. Folding was monitored using probes for changes in both secondary and tertiary structure, and the monomer-dimer equilibrium was perturbed with a variety of solution conditions to allow characterisation of otherwise inaccessible states. Multiphasic kinetics were observed in the unfolding and refolding reactions of YibK, and under conditions where the dimer is favoured, dissociation and association were rate-limiting, respectively. A folding model consistent with all kinetic data is proposed: YibK appears to fold via two parallel pathways, partitioned by proline isomerisation events, to two distinct monomeric intermediates. These form a common third intermediate that is able to fold to native dimer. Kinetic simulations suggest that all intermediates are on-pathway. These results provide the valuable groundwork required to further understand how Nature codes for knot formation.     

The folding mechanics of a knotted protein

An increasing number of proteins are being discovered with a remarkable and somewhat surprising feature, a knot in their native structures. How the polypeptide chain is able to "knot" itself during the folding process to form these highly intricate protein topologies is not known. Here we perform a computational study on the 160-amino-acid homodimeric protein YibK, which, like other proteins in the SpoU family of MTases, contains a deep trefoil knot in its C-terminal region. In this study, we use a coarse-grained C(alpha)-chain representation and Langevin dynamics to study folding kinetics. We find that specific, attractive nonnative interactions are critical for knot formation. In the absence of these interactions, i.e., in an energetics driven entirely by native interactions, knot formation is exceedingly unlikely. Further, we find, in concert with recent experimental data on YibK, two parallel folding pathways that we attribute to an early and a late formation of the trefoil knot, respectively. For both pathways, knot formation occurs before dimerization. A bioinformatics analysis of the SpoU family of proteins reveals further that the critical nonnative interactions may originate from evolutionary conserved hydrophobic segments around the knotted region.     

The dimerization of an alpha/beta-knotted protein is essential for structure and function

alpha/beta-Knotted proteins are an extraordinary example of biological self-assembly; they contain a deep topological trefoil knot formed by the backbone polypeptide chain. Evidence suggests that all are dimeric and function as methyltransferases, and the deep knot forms part of the active site. We investigated the significance of the dimeric structure of the alpha/beta-knot protein, YibK, from Haemophilus influenzae by the design and engineering of monomeric versions of the protein, followed by examination of their structural, functional, stability, and kinetic folding properties. Monomeric forms of YibK display similar characteristics to an intermediate species populated during the formation of the wild-type dimer. However, a notable loss in structure involving disruption to the active site, rendering it incapable of cofactor binding, is observed in monomeric YibK. Thus, dimerization is vital for preservation of the native structure and, therefore, activity of the protein.     

Backbone assignments for the SPOUT methyltransferase MTT<Subscript><Emphasis Type="Italic">Tm</Emphasis></Subscript>, a knotted protein from <Emphasis Type="Italic">Thermotoga maritima</Emphasis>

The SPOUT family of methyltransferase proteins is noted for containing a deep trefoil knot in their defining backbone fold. This unique fold is of high interest for furthering the understanding of knots in proteins. Here, we report the ¹H, ¹³C, ¹⁵N assignments for MTT _Tm , a canonical member of the SPOUT family. This protein is unique, as it is one of the smallest members of the family, making it an ideal system for probing the unique properties of the knot. Our present work represents the foundation for further studies into the topology of MTT _Tm , and understanding how its structure affects both its folding and function.     

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.     

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering,Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.     

Diversity in the disulfide folding pathways of cystine knot peptides

Summary The plant cyclotides are a fascinating family of circular proteins that contain a cyclic cystine knot motif (CCK). This unique family was discovered only recently but contains over 50 known sequences to date. Various biological activities are associated with these peptides including antimicrobial and insecticidal activity. The knotted topology and cyclic nature of the cyclotides poses interesting questions about the folding mechanisms and how the knotted arrangement of disulfide bonds is formed. Some studies have been performed on related inhibitor cystine knot (ICK) containing peptides, but little is known about the folding mechanisms of CCK molecules. We have examined the oxidative refolding and reductive unfolding of the prototypic member of the cyclotide family, kalata B1. Analysis of the rates of formation of the intermediates along the reductive unfolding pathway highlights the stability conferred by the cystine knot motif. Significant differences are observed between the folding of kalata B1 and an acyclic cystine knot protein, EETI-II, suggesting that the circular backbone has a significant influence in directing the folding pathway.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Simulation and experiment conspire to reveal cryptic intermediates and a slide from the nucleation-condensation to framework mechanism of folding

There is a change from three-state to two-state kinetics of folding across the homeodomain superfamily of proteins as the mechanism slides from framework to nucleation-condensation. The tendency for framework folding in this family correlates with inherent helical propensity. The cellular myeloblastis protein (c-Myb) falls in the mechanistic transition region. An earlier, preliminary report of protein engineering experiments and molecular dynamics simulations (MD) showed that the folding mechanism for this protein has aspects of both the nucleation-condensation and framework models. In the more in-depth analysis of the MD trajectories presented here, we find that folding may be attributed to both of these mechanisms in different regions of the protein. The folding of the loop, middle helix, and turn is best described by nucleation-condensation, whereas folding of the N and C-terminal helices may be described by the framework model. Experimentally, c-Myb folds by apparent two-state kinetics, but the MD simulations predict that the kinetics hide a high-energy intermediate. We stabilized this hypothetical folding intermediate by deleting a residue (P174) in the loop between its second and third helices, and the mutant intermediate is long-lived in the simulations. Equilibrium and kinetic experiments demonstrate that folding of the DeltaP174 mutant is indeed three-state. The presence and shape of the intermediate observed in the simulations were confirmed by small angle X-ray scattering experiments.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

Phi-analysis at the experimental limits: mechanism of beta-hairpin formation

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.     

Characterization of the protein unfolding processes induced by urea and temperature

Correct folding is critical for the biological activities of proteins. As a contribution to a better understanding of the protein (un)folding problem, we studied the effect of temperature and of urea on peptostreptococcal Protein L destructuration. We performed standard molecular dynamics simulations at 300 K, 350 K, 400 K, and 480 K, both in 10 M urea and in water. Protein L followed at least two alternative unfolding pathways. Urea caused the loss of secondary structure acting preferentially on the β-sheets, while leaving the α-helices almost intact; on the contrary, high temperature preserved the β-sheets and led to a complete loss of the α-helices. These data suggest that urea and high temperature act through different unfolding mechanisms, and protein secondary motives reveal a differential sensitivity to various denaturant treatments. As further validation of our results, replica-exchange molecular dynamics simulations of the temperature-induced unfolding process in the presence of urea were performed. This set of simulations allowed us to compute the thermodynamical parameters of the process and confirmed that, in the configurational space of Protein L unfolding, both of the above pathways are accessible, although to a different relative extent.     

(Un)Folding Mechanisms of the FBP28 WW Domain in Explicit Solvent Revealed by Multiple Rare Event Simulation Methods

We report a numerical study of the (un)folding routes of the truncated FBP28 WW domain at ambient conditions using a combination of four advanced rare event molecular simulation techniques. We explore the free energy landscape of the native state, the unfolded state, and possible intermediates, with replica exchange molecular dynamics. Subsequent application of bias-exchange metadynamics yields three tentative unfolding pathways at room temperature. Using these paths to initiate a transition path sampling simulation reveals the existence of two major folding routes, differing in the formation order of the two main hairpins, and in hydrophobic side-chain interactions. Having established that the hairpin strand separation distances can act as reasonable reaction coordinates, we employ metadynamics to compute the unfolding barriers and find that the barrier with the lowest free energy corresponds with the most likely pathway found by transition path sampling. The unfolding barrier at 300 K is ∼17 k_BT ≈ 42 kJ/mol, in agreement with the experimental unfolding rate constant. This work shows that combining several powerful simulation techniques provides a more complete understanding of the kinetic mechanism of protein folding.     

Thermal unfolding simulations of a multimeric protein--transition state and unfolding pathways

The folding of an oligomeric protein poses an extra challenge to the folding problem because the protein not only has to fold correctly; it has to avoid nonproductive aggregation. We have carried out over 100 molecular dynamics simulations using an implicit solvation model at different temperatures to study the unfolding of one of the smallest known tetramers, p53 tetramerization domain (p53tet). We found that unfolding started with disruption of the native tetrameric hydrophobic core. The transition state for the tetramer to dimer transition was characterized as a diverse ensemble of different structures using Phi value analysis in quantitative agreement with experimental data. Despite the diversity, the ensemble was still native-like with common features such as partially exposed tetramer hydrophobic core and shifts in the dimer-dimer arrangements. After passing the transition state, the secondary and tertiary structures continued to unfold until the primary dimers broke free. The free dimer had little secondary structure left and the final free monomers were random-coil like. Both the transition states and the unfolding pathways from these trajectories were very diverse, in agreement with the new view of protein folding. The multiple simulations showed that the folding of p53tet is a mixture of the framework and nucleation-condensation mechanisms and the folding is coupled to the complex formation. We have also calculated the entropy and effective energy for the different states along the unfolding pathway and found that the tetramerization is stabilized by hydrophobic interactions.     

The yeast ribosome synthesis factor Emg1 is a novel member of the superfamily of alpha/beta knot fold methyltransferases

Emg1 was previously shown to be required for maturation of the 18S rRNA and biogenesis of the 40S ribosomal subunit. Here we report the determination of the crystal structure of Emg1 at 2 Å resolution in complex with the methyl donor, S-adenosyl-methionine (SAM). This structure identifies Emg1 as a novel member of the alpha/beta knot fold methyltransferase (SPOUT) superfamily. In addition to the conserved SPOUT core, Emg1 has two unique domains that form an extended surface, which we predict to be involved in binding of RNA substrates. A point mutation within a basic patch on this surface almost completely abolished RNA binding in vitro. Three point mutations designed to disrupt the interaction of Emg1 with SAM each caused>100-fold reduction in SAM binding in vitro. Expression of only Emg1 with these mutations could support growth and apparently normal ribosome biogenesis in strains genetically depleted of Emg1. We conclude that the catalytic activity of Emg1 is not essential and that the presence of the protein is both necessary and sufficient for ribosome biogenesis.     

About the structural role of disulfide bridges in serum albumins: evidence from protein simulated unfolding

The role of the 17 disulfide (S-S) bridges in preserving the native conformation of human serum albumin (HSA) is investigated by performing classical molecular dynamics (MD) simulations on protein structures with intact and, respectively, reduced S-S bridges. The thermal unfolding simulations predict a clear destabilization of the protein secondary structure upon reduction of the S-S bridges as well as a significant distortion of the tertiary structure that is revealed by the changes in the protein native contacts fraction. The effect of the S-S bridges reduction on the protein compactness was tested by calculating Gibbs free energy profiles with respect to the protein gyration radius. The theoretical results obtained using the OPLS-AA and the AMBER ff03 force fields are in agreement with the available experimental data. Beyond the validation of the simulation method, the results here reported provide new insights into the mechanism of the protein reductive/oxidative unfolding/folding processes. It is predicted that in the native conformation of the protein, the thiol (-SH) groups belonging to the same reduced S-S bridge are located in potential wells that maintain them in contact. The -SH pairs can be dispatched by specific conformational transitions of the peptide chain located in the neighborhood of the cysteine residues.