Is the continuity of the domains required for the correct folding of a two-domain protein?

Unraveling the Mechanical Unfolding Pathways of a Multidomain Protein: Phosphoglycerate Kinase

Phosphoglycerate kinase (PGK) is a highly conserved enzyme that is crucial for glycolysis. PGK is a monomeric protein composed of two similar domains and has been the focus of many studies for investigating interdomain interactions within the native state and during folding. Previous studies used traditional biophysical methods (such as circular dichroism, tryptophan fluorescence, and NMR) to measure signals over a large ensemble of molecules, which made it difficult to observe transient changes in stability or structure during unfolding and refolding of single molecules. Here, we unfold single molecules of PGK using atomic force spectroscopy and steered molecular dynamic computer simulations to examine the conformational dynamics of PGK during its unfolding process. Our results show that after the initial forced separation of its domains, yeast PGK (yPGK) does not follow a single mechanical unfolding pathway; instead, it stochastically follows two distinct pathways: unfolding from the N-terminal domain or unfolding from the C-terminal domain. The truncated yPGK N-terminal domain unfolds via a transient intermediate, whereas the structurally similar isolated C-terminal domain has no detectable intermediates throughout its mechanical unfolding process. The N-terminal domain in the full-length yPGK displays a strong unfolding intermediate 13% of the time, whereas the truncated domain (yPGKNT) transitions through the intermediate 81% of the time. This effect indicates that the mechanical properties of yPGK cannot be simply deduced from the mechanical properties of its constituents. We also find that Escherichia coli PGK is significantly less mechanically stable as compared to yPGK, contrary to bulk unfolding measurements. Our results support the growing body of observations that the folding behavior of multidomain proteins is difficult to predict based solely on the studies of isolated domains.

The relationship between chain connectivity and domain stability in the equilibrium and kinetic folding mechanisms of dihydrofolate reductase from E.coli

The role of domains in defining the equilibrium and kinetic folding properties of dihydrofolate reductase (DHFR) from Escherichia coli was probed by examining the thermodynamic and kinetic properties of a set of variants in which the chain connectivity in the discontinuous loop domain (DLD) and the adenosine-binding domain (ABD) was altered by permutation. To test the concept that chain cleavage can selectively destabilize the domain in which the N- and C-termini are resident, permutations were introduced at one position within the ABD, one within the DLD and one at a boundary between the domains. The results demonstrated that a continuous ABD is required for a stable thermal intermediate and a continuous DLD is required for a stable urea intermediate. The permutation at the domain interface had both a thermal and urea intermediate. Strikingly, the observable kinetic folding responses of all three permuted proteins were very similar to the wild-type protein. These results demonstrate a crucial role for stable domains in defining the energy surface for the equilibrium folding reaction of DHFR. If domain connectivity affects the kinetic mechanism, the effects must occur in the sub-millisecond time range.

Asymmetric effect of domain interactions on the kinetics of folding in yeast phosphoglycerate kinase

The aim of this work is to shed more light on the effect of domain-domain interactions on the kinetics and the pathway of protein folding. A model protein system consisting of several single-tryptophan variants of the two-domain yeast phosphoglycerate kinase (PGK) and its individual domains was studied. Refolding was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 msec to 1000 sec. Denaturant titrations of both individual domains showed apparent two-state unfolding transitions. Refolding kinetics of the individual domains from different denaturant concentrations, however, revealed the presence of intermediate structures during titration for both domains. Refolding of the same domains within the complete protein showed that domain-domain interactions direct the folding of both domains, but in an asymmetric way. Folding of the N domain was already altered within 1 msec, while detectable changes in the folding of the C domain occurred only 60-100 msec after initiating refolding. All mutants showed a hyperfluorescent kinetic intermediate. Both the disappearance of this intermediate and the completion of the folding were significantly faster in the individual N domain than in the complete protein. On the contrary, folding of the individual C domain was slower than in the complete protein. The presence of the C domain directs the refolding of the N domain along a completely different pathway than that of the individual N domain, while folding of the individual C domain follows the same path as within the complete protein.

Tuning the Heterogeneous Early Folding Dynamics of Phosphoglycerate Kinase

We recently reported stretched kinetics during the formation of a collapsed, long-lived intermediate state of the large two-domain enzyme phosphoglycerate kinase (PGK). It was postulated that intrinsic roughness of the energy landscape on the way downhill to the intermediate causes the lack of a single time-scale. Here, we investigate several alternative explanations for stretched refolding dynamics in more detail: tyrosine fluorescence, multiple tryptophan probes, and rate differences between independently folding domains. To this end, we systematically simplify PGK in several steps from the full protein with two tryptophan residues and all tyrosine residues probed, to a single domain with only one tryptophan residue and no tyrosine residue probed. The kinetics in the 10 micros to 10 ms range are revealed by laser-induced temperature-jump relaxation experiments. The isolated N-terminal domain forms an intermediate by nearly single-exponential kinetics, but the isolated C-terminal domain shows strongly non-exponential kinetics. Thus, domain interaction and a cis-proline residue between the two domains are ruled out as the sole contributors to heterogeneity during the earliest folding dynamics of the C-terminal domain. We apply two limiting models for the roughness of the energy landscape. A sequential three-state model lumps all the roughness into a single trap. The "strange kinetics" model with logarithmic oscillations developed by Klafter and co-workers distributes the roughness over a larger number of states. Both models explain our data about equally well, but the coincidental values of rate constants in all of our double-exponential fits, and the absence of a spectroscopic signature distinct from the endpoints of the folding process favors more roughness than can be explained by just a single trap.

Purification and structural characterization of the central hydrophobic domain of oleosin

The oil bodies of rapeseeds contain a triacylglycerol matrix surrounded by a monolayer of phospholipids embedded with abundant structural alkaline proteins termed oleosins and some other minor proteins. Oleosins are unusual proteins because they contain a 70-80-residue uninterrupted nonpolar domain flanked by relatively polar C- and N-terminal domains. Although the hydrophilic N-terminal domain had been studied, the structural feature of the central hydrophobic domain remains unclear due to its high hydrophobicity. In the present study, we reported the generation, purification, and characterization of a 9-kDa central hydrophobic domain from rapeseed oleosin (19 kDa). The 9-kDa central hydrophobic domain was produced by selectively degrading the N and C termini with enzymes and then purifying the digest by SDS-PAGE and electroelution. We have also reconstituted the central domain into liposomes and synthetic oil bodies to determine the secondary structure of the domain using CD and Fourier transform infrared (FTIR) spectroscopy. The spectra obtained from CD and FTIR were analyzed with reference to structural information of the N-terminal domain and the full-length rapeseed oleosin. Both CD and FTIR analysis revealed that 50-63% of the domain was composed of beta-sheet structure. Detailed analysis of the FTIR spectra indicated that 80% of the beta-sheet structure, present in the central domain, was arranged in parallel to the intermolecular beta-sheet structure. Therefore, interactions between adjacent oleosin proteins would give rise to a stable beta-sheet structure that would extend around the surface of the seed oil bodies stabilizing them in emulsion systems. The strategies used in our present study are significant in that it could be generally used to study difficult proteins with different independent structural domains, especially with long hydrophobic domains.

Protein folding in the post-genomic era

Protein folding is a topic of fundamental interest since it concerns the mechanisms by which the genetic message is translated into the three-dimensional and functional structure of proteins. In these post-genomic times, the knowledge of the fundamental principles are required in the exploitation of the information contained in the increasing number of sequenced genomes. Protein folding also has practical applications in the understanding of different pathologies and the development of novel therapeutics to prevent diseases associated with protein misfolding and aggregation. Significant advances have been made ranging from the Anfinsen postulate to the "new view" which describes the folding process in terms of an energy landscape. These new insights arise from both theoretical and experimental studies. The problem of folding in the cellular environment is briefly discussed. The modern view of misfolding and aggregation processes that are involved in several pathologies such as prion and Alzheimer diseases. Several approaches of structure prediction, which is a very active field of research, are described.

Role of loops in the folding and stability of yeast phosphoglycerate kinase

Yeast phosphoglycerate kinase (yPGK) is a monomeric two domain protein used as folding model representative of large proteins. We inserted short unstructured sequences (four Gly or four Thr) into the connections between secondary structure elements and studied the consequences of these insertions on the folding process and stability of yPGK. All the mutated proteins can refold efficiently. The effect per residue on stability is larger for the first inserted residue. Insertion in two long betaalpha loops (at residue positions 71 and 129) is more destabilizing than an insertion in a short alphabeta loop (at residue position 89) located on the opposite side of the N-terminal domain. The effect on stability is mainly due to a large increase of the unfolding rate rather than a decrease of the folding rate. This suggests that these connections between secondary structure elements do not play an active role in directing the folding process. Insertion into the short alphabeta loop (position 89) has limited effects on stability and results in the detection of a kinetic phase not previously seen with the wild-type protein, suggesting that insertions in this particular loop do qualitatively affect the folding process without a large effect on folding efficiency. For the two long betaalpha loops (positions 71 and 129) located in the inner surface of the N-terminal domain, the effects on stability are possibly associated with decoupling of the two domains as observed by differential scanning calorimetry during thermal unfolding.

Testing the role of chain connectivity on the stability and structure of dihydrofolate reductase from E. coli: fragment complementation and circular permutation reveal stable, alternatively folded forms

The effects of chain cleavage and circular permutation on the structure, stability, and activity of dihydrofolate reductase (DHFR) from Escherichia coli were investigated by various spectroscopic and biochemical methods. Cleavage of the backbone after position 86 resulted in two fragments, (1--86) and (87--159) each of which are poorly structured and enzymatically inactive. When combined in a 1 : 1 molar ratio, however, the fragments formed a high-affinity (K(a) = 2.6 x 10(7) M(-1)) complex that displays a weakly cooperative urea-induced unfolding transition at micromolar concentrations. The retention of about 15% of the enzymatic activity of full-length DHFR is surprising, considering that the secondary structure in the complex is substantially reduced from its wild-type counterpart. In contrast, a circularly permuted form with its N-terminus at position 86 has similar overall stability to full-length DHFR, about 50% of its activity, substantial secondary structure, altered side-chain packing in the adenosine binding domain, and unfolds via an equilibrium intermediate not observed in the wild-type protein. After addition of ligand or the tight-binding inhibitor methotrexate, both the fragment complex and the circular permutant adopt more native-like secondary and tertiary structures. These results show that changes in the backbone connectivity can produce alternatively folded forms and highlight the importance of protein-ligand interactions in stabilizing the active site architecture of DHFR.

Functionally accepted insertions of proteins within protein domains

Experiments were designed to explore the tolerance of protein structure and folding to very large insertions of folded protein within a structural domain. Dihydrofolate reductase and beta-lactamase have been inserted in four different positions of phosphoglycerate kinase. The resultant chimeric proteins are all overexpressed, and the host as well as the inserted partners are functional. Although not explicitly designed, functional coupling between the two fused partners was observed in some of the chimeras. These results show that the tolerance of protein structures to very large structured insertions is more general than previously expected and supports the idea that the natural sequence continuity of a structural domain is not required for the folding process. These results directly suggest a new experimental approach to screen, for example, for folded protein in randomized polypeptide sequences.

Eye lens betaB2-crystallin: circular permutation does not influence the oligomerization state but enhances the conformational stability

The related vertebrate eye lens polypeptides, betaB2- and gammaB-crystallin, each fold into two similar beta-sheet domains. The main difference is the state of oligomerization resulting from intermolecular domain interactions in the oligomeric beta-crystallins and intramolecular contacts in the monomeric gamma-crystallins. The question arises whether it is possible to create a monomeric gammaB-like betaB2-molecule by protein engineering methods. We wanted to produce such a molecule by circularly permuting the domains of betaB2-crystallin. The new termini were created from the original connecting peptide, and the new linker from stumps of the original extensions, while the rest of the flexible extensions were deleted. As judged by circular dichroism and fluorescence, the permutation causes little change in the structure of the protein. The circularly permuted protein forms dimers as wild-type betaB2-crystallin. On the other hand, cpbetaB2 shows a slightly enhanced stability against urea with a midpoint of transition of 2.1 M urea versus 1.9 M for the wild-type protein lacking N and C-terminal arms, thus indicating stronger domain interactions. To our knowledge this is the first circularly permuted protein which exhibits a higher stability than the corresponding wild-type protein.

Structure of pvu II DNA-(cytosine N4) methyltransferase, an example of domain permutation and protein fold assignment

We have determined the structure of Pvu II methyltransferase (M. Pvu II) complexed with S -adenosyl-L-methionine (AdoMet) by multiwavelength anomalous diffraction, using a crystal of the selenomethionine-substituted protein. M. Pvu II catalyzes transfer of the methyl group from AdoMet to the exocyclic amino (N4) nitrogen of the central cytosine in its recognition sequence 5'-CAGCTG-3'. The protein is dominated by an open alpha/beta-sheet structure with a prominent V-shaped cleft: AdoMet and catalytic amino acids are located at the bottom of this cleft. The size and the basic nature of the cleft are consistent with duplex DNA binding. The target (methylatable) cytosine, if flipped out of the double helical DNA as seen for DNA methyltransferases that generate 5-methylcytosine, would fit into the concave active site next to the AdoMet. This M. Pvu IIalpha/beta-sheet structure is very similar to those of M. Hha I (a cytosine C5 methyltransferase) and M. Taq I (an adenine N6 methyltransferase), consistent with a model predicting that DNA methyltransferases share a common structural fold while having the major functional regions permuted into three distinct linear orders. The main feature of the common fold is a seven-stranded beta-sheet (6 7 5 4 1 2 3) formed by five parallel beta-strands and an antiparallel beta-hairpin. The beta-sheet is flanked by six parallel alpha-helices, three on each side. The AdoMet binding site is located at the C-terminal ends of strands beta1 and beta2 and the active site is at the C-terminal ends of strands beta4 and beta5 and the N-terminal end of strand beta7. The AdoMet-protein interactions are almost identical among M. Pvu II, M. Hha I and M. Taq I, as well as in an RNA methyltransferase and at least one small molecule methyltransferase. The structural similarity among the active sites of M. Pvu II, M. Taq I and M. Hha I reveals that catalytic amino acids essential for cytosine N4 and adenine N6 methylation coincide spatially with those for cytosine C5 methylation, suggesting a mechanism for amino methylation.