Protein unfolding rates correlate as strongly as folding rates with native structure

Although the folding rates of proteins have been studied extensively, both experimentally and theoretically, and many native state topological parameters have been proposed to correlate with or predict these rates, unfolding rates have received much less attention. Moreover, unfolding rates have generally been thought either to not relate to native topology in the same manner as folding rates, perhaps depending on different topological parameters, or to be more difficult to predict. Using a dataset of 108 proteins including two-state and multistate folders, we find that both unfolding and folding rates correlate strongly, and comparably well, with well-established measures of native topology, the absolute contact order and the long range order, with correlation coefficient values of 0.75 or higher. In addition, compared to folding rates, the absolute values of unfolding rates vary more strongly with native topology, have a larger range of values, and correlate better with thermodynamic stability. Similar trends are observed for subsets of different protein structural classes. Taken together, these results suggest that choosing a scaffold for protein engineering may require a compromise between a simple topology that will fold sufficiently quickly but also unfold quickly, and a complex topology that will unfold slowly and hence have kinetic stability, but fold slowly. These observations, together with the established role of kinetic stability in determining resistance to thermal and chemical denaturation as well as proteases, have important implications for understanding fundamental aspects of protein unfolding and folding and for protein engineering and design.     

Plakophilin3 loss leads to an increase in PRL3 levels promoting K8 dephosphorylation, which is required for transformation and metastasis

The desmosome anchors keratin filaments in epithelial cells leading to the formation of a tissue wide IF network. Loss of the desmosomal plaque protein plakophilin3 (PKP3) in HCT116 cells, leads to an increase in neoplastic progression and metastasis, which was accompanied by an increase in K8 levels. The increase in levels was due to an increase in the protein levels of the Phosphatase of Regenerating Liver 3 (PRL3), which results in a decrease in phosphorylation on K8. The increase in PRL3 and K8 protein levels could be reversed by introduction of an shRNA resistant PKP3 cDNA. Inhibition of K8 expression in the PKP3 knockdown clone S10, led to a decrease in cell migration and lamellipodia formation. Further, the K8 PKP3 double knockdown clones showed a decrease in colony formation in soft agar and decreased tumorigenesis and metastasis in nude mice. These results suggest that a stabilisation of K8 filaments leading to an increase in migration and transformation may be one mechanism by which PKP3 loss leads to tumor progression and metastasis.     

Designing cooperativity into the designed protein Top7

The topology of the designed protein Top7 is not found in natural proteins. Top7 shows signatures of non‐cooperative folding in both experimental studies and computer simulations. In particular, molecular dynamics of coarse‐grained structure‐based models of Top7 show a well‐populated C‐terminal folding‐intermediate. Since most similarly sized globular proteins are cooperative folders, the non‐natural topology of Top7 has been suggested as a reason for its non‐cooperative folding. Here, we computationally examine the folding of Top7 with the intent of making it cooperative. We find that its folding cooperativity can be increased in two ways: (a) Optimization of packing interactions in the N‐terminal half of the protein enables further folding of the C‐terminal intermediate. (b) Reduction in the packing density of the C‐terminal region destabilizes the intermediate. In practice, these strategies are implemented in our Top7 model through modifications to the contact‐map. These modifications do not alter the topology of Top7 but result in cooperative folding. Amino‐acid mutations that mimic these modifications also lead to a significant increase in folding cooperativity. Finally, we devise a method to randomize the sizes of amino‐acids within the same topology, and confirm that the structure of Top7 makes its folding sensitive to packing interactions. In contrast, the ribosomal protein S6, which has secondary structure similar to Top7, has folding which is much less sensitive to packing perturbations. Thus, it should be possible to make a sequence fold cooperatively to the structure of Top7, but to do so its side‐chain packing needs to be carefully designed. Proteins 2014; 82:364–374. © 2013 Wiley Periodicals, Inc.     

Conformational transitions in adenylate kinase. Allosteric communication reduces misligation

Large conformational changes in the LID and NMP domains of adenylate kinase (AKE) are known to be key to ligand binding and catalysis, yet the order of binding events and domain motion is not well understood. Combining the multiple available structures for AKE with the energy landscape theory for protein folding, a theoretical model was developed for allostery, order of binding events, and efficient catalysis. Coarse-grained models and nonlinear normal mode analysis were used to infer that intrinsic structural fluctuations dominate LID motion, whereas ligand-protein interactions and cracking (local unfolding) are more important during NMP motion. In addition, LID-NMP domain interactions are indispensable for efficient catalysis. LID domain motion precedes NMP domain motion, during both opening and closing. These findings provide a mechanistic explanation for the observed 1:1:1 correspondence between LID domain closure, NMP domain closure, and substrate turnover. This catalytic cycle has likely evolved to reduce misligation, and thus inhibition, of AKE. The separation of allosteric motion into intrinsic structural fluctuations and ligand-induced contributions can be generalized to further our understanding of allosteric transitions in other proteins.     

Molecular phylogeny of Glyphochloa (Poaceae,Panicoideae), an endemic grass genus from the Western Ghats,India

The genus Glyphochloa (Poaceae: Panicoideae: Andropogoneae: Rottboellinae) is endemic to peninsular India and is distributed on lateritic plateaus of low and high altitude in and around Western Ghats and the Malabar Coast. The genus presumably originated and diversified in the Western Ghats. Species relationships in the genus Glyphochloa were deduced here based on molecular phylogenies inferred using nuclear ribosomal ITS sequences and plastid intergenic spacer regions (atpB-rbcL, trnT-trnL, trnL-trnF), and new observations were made of spikelet morphology, caryopsis morphology and meiotic chromosome counts. We observed two distinct clades of Glyphochloa s.l. One of these (‘group I’) includes Ophiuros bombaiensis, and is characterized by a single-awned lower glume and a base chromosome number of 6; it grows in low elevation coastal areas. The other clade (‘group II’) has a double-awned lower glume, a base chromosome number of 7, and is restricted to higher elevation lateritic plateaus; G. ratnagirica may belong to the group II clade, or may be a third distinct lineage in the genus. A sister-group relationship between group I and II taxa (with or without G. ratnagirica) is not well supported, although the genus is recovered as monophyletic in shortest trees inferred using ITS or concatenated plastid data. We present a key to species of Glyphochloa and make a new combination for O. bombaiensis.     

Topological frustration and the folding of interleukin-1 beta

The cytokine, interleukin-1beta (IL-1beta), adopts a beta-trefoil fold. It is known to be much slower folding than similarly sized proteins, despite having a low contact order. Proteins are sufficiently well designed that their folding is not dominated by local energetic traps. Therefore, protein models that encode only the folded structure and are energetically unfrustrated (Gō-type), can capture the essentials of the folding routes. We investigate the folding thermodynamics of IL-1beta using such a model and molecular dynamics (MD) simulations. We develop an enhanced sampling technique (a modified multicanonical method) to overcome the sampling problem caused by the slow folding. We find that IL-1beta has a broad and high free energy barrier. In addition, the protein fold causes intermediate unfolding and refolding of some native contacts within the protein along the folding trajectory. This "backtracking" occurs around the barrier region. Complex folds like the beta-trefoil fold and functional loops like the beta-bulge of IL-1beta can make some of the configuration space unavailable to the protein and cause topological frustration.     

Folding circular permutants of IL-1β: route selection driven by functional frustration

Interleukin-1β (IL-1β) is the cytokine crucial to inflammatory and immune response. Two dominant routes are populated in the folding to native structure. These distinct routes are a result of the competition between early packing of the functional loops versus closure of the β-barrel to achieve efficient folding and have been observed both experimentally and computationally. Kinetic experiments on the WT protein established that the dominant route is characterized by early packing of geometrically frustrated functional loops. However, deletion of one of the functional loops, the β-bulge, switches the dominant route to an alternative, yet, as accessible, route, where the termini necessary for barrel closure form first. Here, we explore the effect of circular permutation of the WT sequence on the observed folding landscape with a combination of kinetic and thermodynamic experiments. Our experiments show that while the rate of formation of permutant protein is always slower than that observed for the WT sequence, the region of initial nucleation for all permutants is similar to that observed for the WT protein and occurs within a similar timescale. That is, even permutants with significant sequence rearrangement in which the functional-nucleus is placed at opposing ends of the polypeptide chain, fold by the dominant WT "functional loop-packing route", despite the entropic cost of having to fold the N- and C- termini early. Taken together, our results indicate that the early packing of the functional loops dominates the folding landscape in active proteins, and, despite the entropic penalty of coalescing the termini early, these proteins will populate an entropically unfavorable route in order to conserve function. More generally, circular permutation can elucidate the influence of local energetic stabilization of functional regions within a protein, where topological complexity creates a mismatch between energetics and topology in active proteins.     

An all-atom structure-based potential for proteins: bridging minimal models with all-atom empirical forcefields

Protein dynamics take place on many time and length scales. Coarse-grained structure-based (Go) models utilize the funneled energy landscape theory of protein folding to provide an understanding of both long time and long length scale dynamics. All-atom empirical forcefields with explicit solvent can elucidate our understanding of short time dynamics with high energetic and structural resolution. Thus, structure-based models with atomic details included can be used to bridge our understanding between these two approaches. We report on the robustness of folding mechanisms in one such all-atom model. Results for the B domain of Protein A, the SH3 domain of C-Src Kinase, and Chymotrypsin Inhibitor 2 are reported. The interplay between side chain packing and backbone folding is explored. We also compare this model to a C(alpha) structure-based model and an all-atom empirical forcefield. Key findings include: (1) backbone collapse is accompanied by partial side chain packing in a cooperative transition and residual side chain packing occurs gradually with decreasing temperature, (2) folding mechanisms are robust to variations of the energetic parameters, (3) protein folding free-energy barriers can be manipulated through parametric modifications, (4) the global folding mechanisms in a C(alpha) model and the all-atom model agree, although differences can be attributed to energetic heterogeneity in the all-atom model, and (5) proline residues have significant effects on folding mechanisms, independent of isomerization effects. Because this structure-based model has atomic resolution, this work lays the foundation for future studies to probe the contributions of specific energetic factors on protein folding and function.     

In the Multi-domain Protein Adenylate Kinase,Domain Insertion Facilitates Cooperative Folding while Accommodating Function at Domain Interfaces

Having multiple domains in proteins can lead to partial folding and increased aggregation. Folding cooperativity, the all or nothing folding of a protein, can reduce this aggregation propensity. In agreement with bulk experiments, a coarse-grained structure-based model of the three-domain protein, E. coli Adenylate kinase (AKE), folds cooperatively. Domain interfaces have previously been implicated in the cooperative folding of multi-domain proteins. To understand their role in AKE folding, we computationally create mutants with deleted inter-domain interfaces and simulate their folding. We find that inter-domain interfaces play a minor role in the folding cooperativity of AKE. On further analysis, we find that unlike other multi-domain proteins whose folding has been studied, the domains of AKE are not singly-linked. Two of its domains have two linkers to the third one, i.e., they are inserted into the third one. We use circular permutation to modify AKE chain-connectivity and convert inserted-domains into singly-linked domains. We find that domain insertion in AKE achieves the following: (1) It facilitates folding cooperativity even when domains have different stabilities. Insertion constrains the N- and C-termini of inserted domains and stabilizes their folded states. Therefore, domains that perform conformational transitions can be smaller with fewer stabilizing interactions. (2) Inter-domain interactions are not needed to promote folding cooperativity and can be tuned for function. In AKE, these interactions help promote conformational dynamics limited catalysis. Finally, using structural bioinformatics, we suggest that domain insertion may also facilitate the cooperative folding of other multi-domain proteins.