Downhill protein folding: evolution meets physics

Hydrogen bonding heterogeneity correlates with protein folding transition state passage time as revealed by data sonification

Protein-protein and protein-water hydrogen bonding interactions play essential roles in the way a protein passes through the transition state during folding or unfolding, but the large number of these interactions in molecular dynamics (MD) simulations makes them difficult to analyze. Here, we introduce a state space representation and associated "rarity" measure to identify and quantify transition state passage (transit) events. Applying this representation to a long MD simulation trajectory that captured multiple folding and unfolding events of the GTT WW domain, a small protein often used as a model for the folding process, we identified three transition categories: Highway (faster), Meander (slower), and Ambiguous (intermediate). We developed data sonification and visualization tools to analyze hydrogen bond dynamics before, during, and after these transition events. By means of these tools, we were able to identify characteristic hydrogen bonding patterns associated with "Highway" versus "Meander" versus "Ambiguous" transitions and to design algorithms that can identify these same folding pathways and critical protein-water interactions directly from the data. Highly cooperative hydrogen bonding can either slow down or speed up transit. Furthermore, an analysis of protein-water hydrogen bond dynamics at the surface of WW domain shows an increase in hydrogen bond lifetime from folded to unfolded conformations with Ambiguous transitions as an outlier. In summary, hydrogen bond dynamics provide a direct window into the heterogeneity of transits, which can vary widely in duration (by a factor of 10) due to a complex energy landscape.

Distinguishing between concerted, sequential and barrierless conformational changes: Folding versus allostery

Characterization of transition and intermediate states of reactions provides insights into their mechanisms and is often achieved through analysis of linear free energy relationships. Such an approach has been used extensively in protein folding studies but less so for analyzing allosteric transitions. Here, we point out analogies in ways to characterize pathways and intermediates in folding and allosteric transitions. Achieving an understanding of the mechanisms by which proteins undergo allosteric switching is important in many cases for obtaining insights into how they function.

Effects of sequence-dependent non-native interactions in equilibrium and kinetic folding properties of knotted proteins

Determining the role of non-native interactions in folding dynamics, kinetics, and mechanisms is a classic problem in protein folding. More recently, this question has witnessed a renewed interest in light of the hypothesis that knotted proteins require the assistance of non-native interactions to fold efficiently. Here, we conduct extensive equilibrium and kinetic Monte Carlo simulations of a simple off-lattice C-alpha model to explore the role of non-native interactions in the thermodynamics and kinetics of three proteins embedding a trefoil knot in their native structure. We find that equilibrium knotted conformations are stabilized by non-native interactions that are non-local, and proximal to native ones, thus enhancing them. Additionally, non-native interactions increase the knotting frequency at high temperatures, and in partially folded conformations below the transition temperatures. Although non-native interactions clearly enhance the efficiency of transition from an unfolded conformation to a partially folded knotted one, they are not required to efficiently fold a knotted protein. Indeed, a native-centric interaction potential drives the most efficient folding transition, provided that the simulation temperature is well below the transition temperature of the considered model system.

Emerging role of heat shock proteins in cardiovascular diseases

Heat Shock Proteins (HSPs) are evolutionarily conserved proteins from prokaryotes to eukaryotes. They are ubiquitous proteins involved in key physiological and cellular pathways (viz. inflammation, immunity and apoptosis). Indeed, the survivability of the cells under various stressful conditions depends on appropriate levels of HSP expression. There is a growing line of evidence for the role of HSPs in regulating cardiovascular diseases (CVDs) (viz. hypertension, atherosclerosis, atrial fibrillation, cardiomyopathy and heart failure). Furthermore, studies indicate that a higher concentration of circulatory HSP antibodies correlate to CVDs; some are even potential markers for CVDs. The multifaceted roles of HSPs in regulating cellular signaling necessitate unraveling their links to pathophysiology of CVDs. This review aims to consolidate our understanding of transcriptional (via multiple transcription factors including HSF-1, NF-κB, CREB and STAT3) and post-transcriptional (via microRNAs including miR-1, miR-21 and miR-24) regulation of HSPs. The cytoprotective nature of HSPs catapults them to the limelight as modulators of cell survival. Yet another attractive prospect is the development of new therapeutic strategies against cardiovascular diseases (from hypertension to heart failure) by targeting the regulation of HSPs. Moreover, this review provides insights into how genetic variation of HSPs can contribute to the manifestation of CVDs. It would also offer a bird's eye view of the evolving role of different HSPs in the modulation and manifestation of cardiovascular disease.

Insights into Fluctuations of Structure of Proteins: Significance of Intermediary States in Regulating Biological Functions

Proteins are indispensable to cellular communication and metabolism. The structure on which cells and tissues are developed is deciphered from proteins. To perform functions, proteins fold into a three-dimensional structural design, which is specific and fundamentally determined by their characteristic sequence of amino acids. Few of them have structural versatility, allowing them to adapt their shape to the task at hand. The intermediate states appear momentarily, while protein folds from denatured (D) ⇔ native (N), which plays significant roles in cellular functions. Prolific effort needs to be taken in characterizing these intermediate species if detected during the folding process. Protein folds into its native structure through definite pathways, which involve a limited number of transitory intermediates. Intermediates may be essential in protein folding pathways and assembly in some cases, as well as misfolding and aggregation folding pathways. These intermediate states help to understand the machinery of proper folding in proteins. In this review article, we highlight the various intermediate states observed and characterized so far under in vitro conditions. Moreover, the role and significance of intermediates in regulating the biological function of cells are discussed clearly.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

A Comparative Study to Decipher the Structural and Dynamics Determinants Underlying the Activity and Thermal Stability of GH-11 Xylanases

With the growing need for renewable sources of energy, the interest for enzymes capable of biomass degradation has been increasing. In this paper, we consider two different xylanases from the GH-11 family: the particularly active GH-11 xylanase from Neocallimastix patriciarum, NpXyn11A, and the hyper-thermostable mutant of the environmentally isolated GH-11 xylanase, EvXyn11^TS. Our aim is to identify the molecular determinants underlying the enhanced capacities of these two enzymes to ultimately graft the abilities of one on the other. Molecular dynamics simulations of the respective free-enzymes and enzyme–xylohexaose complexes were carried out at temperatures of 300, 340, and 500 K. An in-depth analysis of these MD simulations showed how differences in dynamics influence the activity and stability of these two enzymes and allowed us to study and understand in greater depth the molecular and structural basis of these two systems. In light of the results presented in this paper, the thumb region and the larger substrate binding cleft of NpXyn11A seem to play a major role on the activity of this enzyme. Its lower thermal stability may instead be caused by the higher flexibility of certain regions located further from the active site. Regions such as the N-ter, the loops located in the fingers region, the palm loop, and the helix loop seem to be less stable than in the hyper-thermostable EvXyn11^TS. By identifying molecular regions that are critical for the stability of these enzymes, this study allowed us to identify promising targets for engineering GH-11 xylanases. Eventually, we identify NpXyn11A as the ideal host for grafting the thermostabilizing traits of EvXyn11^TS.

Studying protein folding in health and disease using biophysical approaches

Protein folding is crucial for normal physiology including development and healthy aging, and failure of this process is related to the pathology of diseases including neurodegeneration and cancer. Early thermodynamic and kinetic studies based on the unfolding and refolding equilibrium of individual proteins in the test tube have provided insight into the fundamental principles of protein folding, although the problem of predicting how any given protein will fold remains unsolved. Protein folding within cells is a more complex issue than folding of purified protein in isolation, due to the complex interactions within the cellular environment, including post-translational modifications of proteins, the presence of macromolecular crowding in cells, and variations in the cellular environment, for example in cancer versus normal cells. Development of biophysical approaches including fluorescence resonance energy transfer (FRET) and nuclear magnetic resonance (NMR) techniques and cellular manipulations including microinjection and insertion of noncanonical amino acids has allowed the study of protein folding in living cells. Furthermore, biophysical techniques such as single-molecule fluorescence spectroscopy and optical tweezers allows studies of simplified systems at the single molecular level. Combining in-cell techniques with the powerful detail that can be achieved from single-molecule studies allows the effects of different cellular components including molecular chaperones to be monitored, providing us with comprehensive understanding of the protein folding process. The application of biophysical techniques to the study of protein folding is arming us with knowledge that is fundamental to the battle against cancer and other diseases related to protein conformation or protein–protein interactions.

The change of conditions does not affect Ros87 downhill folding mechanism

Downhill folding has been defined as a unique thermodynamic process involving a conformations ensemble that progressively loses structure with the decrease of protein stability. Downhill folders are estimated to be rather rare in nature as they miss an energetically substantial folding barrier that can protect against aggregation and proteolysis. We have previously demonstrated that the prokaryotic zinc finger protein Ros87 shows a bipartite folding/unfolding process in which a metal binding intermediate converts to the native structure through a delicate barrier-less downhill transition. Significant variation in folding scenarios can be detected within protein families with high sequence identity and very similar folds and for the same sequence by varying conditions. For this reason, we here show, by means of DSC, CD and NMR, that also in different pH and ionic strength conditions Ros87 retains its partly downhill folding scenario demonstrating that, at least in metallo-proteins, the downhill mechanism can be found under a much wider range of conditions and coupled to other different transitions. We also show that mutations of Ros87 zinc coordination sphere produces a different folding scenario demonstrating that the organization of the metal ion core is determinant in the folding process of this family of proteins.

Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics

As theory and experiment have shown, protein dehydration is a major contributor to protein folding. Dehydration upon folding can be characterized directly by all-atom simulations of fast pressure drops, which create desolvated pockets inside the nascent hydrophobic core. Here, we study pressure-drop refolding of three λ-repressor fragment (λ_6-85) mutants computationally and experimentally. The three mutants report on tertiary structure formation via different fluorescent helix-helix contact pairs. All-atom simulations of pressure drops capture refolding and unfolding of all three mutants by a similar mechanism, thus validating the nonperturbative nature of the fluorescent contact probes. Analysis of simulated interprobe distances shows that the α-helix 1-3 pair distance displays a slower characteristic time scale than the 1-2 or 3-2 pair distance. To see whether slow packing of α-helices 1 and 3 is reflected in the rate-limiting folding step, fast pressure-drop relaxation experiments captured refolding on a millisecond time scale. These experiments reveal that refolding monitored by 1-3 contact formation indeed is much slower than when monitored by 1-2 or 3-2 contact formation. Unlike the case of the two-state folder [three-α-helix bundle (α₃D)], whose drying and core formation proceed in concert, λ_6-85 repeatedly dries and rewets different local tertiary contacts before finally forming a solvent-excluded core, explaining the non-two-state behavior observed during refolding in molecular dynamics simulations. This work demonstrates that proteins can explore desolvated pockets and dry globular states numerous times before reaching the native conformation.

Protein Folding Cooperativity and Thermodynamic Barriers of the Simplest β-Sheet Fold: A Survey of WW Domains

Theory and experiments have shown that microsecond folding proteins exhibit characteristic thermodynamic properties that reflect the limited cooperativity of folding over marginal barriers (downhill folding). Those studies have mostly focused on proteins with large α-helical contents and small size, which tend to be the fastest folders. A key open question is whether such properties are also present in the fastest all-β proteins. We address this issue by investigating the unfolding thermodynamics of a collection of WW domains as representatives of the simplest β-sheet fold. WW domains are small microsecond folders, although they do not fold as fast as their α-helical counterparts. In previous work on the NEDD4-WW4 domain, we reported deviations from two-state thermodynamics that were less apparent and thus suggestive of an incipient downhill scenario. Here we investigate the unfolding thermodynamics of four other WW domains (NEDD4-WW3, YAP65-WW1(L30K), FBP11-WW1, and FBP11-WW2) by performing all of the thermodynamic tests for downhill folding that have been previously developed on α-helical proteins. This set of five WW domains shares low sequence identity and include examples from two specificity classes, thus providing a comprehensive survey. Thermodynamic analysis of the four new WW domains consistently reveals all of the properties of downhill folding equilibria, which are in all cases more marked than what we found before in NEDD4-WW4. Our results show that fast-folding all-β proteins do share limited cooperativity and gradual unfolding thermodynamics with fast α-helical proteins and suggest that the free energy barrier to folding of natural proteins is mostly determined by size and fold topology and much less by the specific amino acid sequence.

Pressure- and heat-induced protein unfolding in bacterial cells: crowding vs. sticking

In-cell protein stability is increased by crowding, but can be reduced by destabilizing surface interactions. Will different denaturation techniques yield similar trends? Here, we apply pressure and thermal denaturation to green fluorescent protein/ReAsH-labeled yeast phosphoglycerate kinase (PGK) in Escherichia coli cells. Pressure denaturation is more two state-like in E. coli than in vitro, stabilizing the native state. Thermal denaturation destabilizes PGK in E. coli, unlike in mammalian cells. Results in wild-type MG1655 strain are corroborated in pressure-resistant J1 strain, where PGK is less prone to aggregation. Thus, destabilizing surface interactions overcome stabilizing crowding in the E. coli cytoplasm under thermal denaturation, but not under pressure denaturation.

Folding mechanisms steer the amyloid fibril formation propensity of highly homologous proteins

Understanding the molecular determinants of fibrillogenesis by studying the aggregation propensities of high homologous proteins with different folding pathways.

Significant advances in the understanding of the molecular determinants of fibrillogenesis can be expected from comparative studies of the aggregation propensities of proteins with highly homologous structures but different folding pathways. Here, we fully characterize, by means of stopped-flow, T-jump, CD and DSC experiments, the unfolding mechanisms of three highly homologous proteins, zinc binding Ros87 and Ml1_53–149 and zinc-lacking Ml4_52–151. The results indicate that the three proteins significantly differ in terms of stability and (un)folding mechanisms. Particularly, Ros87 and Ml1_53–149 appear to be much more stable to guanidine denaturation and are characterized by folding mechanisms including the presence of an intermediate. On the other hand, metal lacking Ml4_52–151 folds according to a classic two-state model. Successively, we have monitored the capabilities of Ros87, Ml4_52–151 and Ml1_53–149 to form amyloid fibrils under native conditions. Particularly, we show, by CD, fluorescence, DLS, TEM and SEM experiments, that after 168 hours, amyloid formation of Ros87 has started, while Ml1_53–149 has formed only amorphous aggregates and Ml4_52–151 is still monomeric in solution. This study shows how metal binding can influence protein folding pathways and thereby control conformational accessibility to aggregation-prone states, which in turn changes aggregation kinetics, shedding light on the role of metal ions in the development of protein deposition diseases.

Be positive: optimizing pramlintide from microcanonical analysis of amylin isoforms

Amylin, or human islet amyloid polypeptide (hIAPP), is a 37-residue hormone synergistic to insulin and co-secreted with it by β-cells in the pancreas. The deposition of its cytotoxic amyloid fibrils is strongly related to the progression of Type II diabetes (T2D) and islet graft failures. Notably, isoforms from some mammalian species, such as rats (rIAPP) and porcine (pIAPP), present a few key mutations preventing aggregation. This has lead to biotechnological development of drugs for adjunct therapies of T2D, such as pramlintide, a variant of hIAPP inspired by rIAPP whose proline substitutions have β-strand fibril-breaking properties. Ideally, such a drug should be formulated with insulin and co-administered, but this has been prevented by a poor solubility profile at the appropriate pH. Hopefully, this could be improved with appropriate point mutations, increasing the molecular net charge. Despite experimental progress, preliminary screening during rational drug design can greatly benefit from thermodynamic insight derived from molecular simulations. So we introduce microcanonical thermostatistics analysis of multicanonical (MUCA) simulations of wild-type amylin isoforms as a systematic assessment of protein thermostability. As a consequence of this comprehensive investigation, the most suitable single-point mutations able to optimize pramlintide are located among the wild-type amylin isoforms. In particular, we find that aggregation inhibition and increased solubility are inherited by pramlintide through further S20R substitution typical of pIAPP. Thus, we provide a consistent thermostatistical methodology to aid the design of improved adjunct therapies for T2D according to current clinical knowledge.

Slow domain reconfiguration causes power-law kinetics in a two-state enzyme

Protein dynamics are typically captured well by rate equations that predict exponential decays for two-state reactions. Here, we describe a remarkable exception. The electron-transfer enzyme quiescin sulfhydryl oxidase (QSOX), a natural fusion of two functionally distinct domains, switches between open- and closed-domain arrangements with apparent power-law kinetics. Using single-molecule FRET experiments on time scales from nanoseconds to milliseconds, we show that the unusual open-close kinetics results from slow sampling of an ensemble of disordered domain orientations. While substrate accelerates the kinetics, thus suggesting a substrate-induced switch to an alternative free energy landscape of the enzyme, the power-law behavior is also preserved upon electron load. Our results show that the slow sampling of open conformers is caused by a variety of interdomain interactions that imply a rugged free energy landscape, thus providing a generic mechanism for dynamic disorder in multidomain enzymes.

Globular Protein Folding In Vitro and In Vivo

In vitro, computational, and theoretical studies of protein folding have converged to paint a rich and complex energy landscape. This landscape is sensitively modulated by environmental conditions and subject to evolutionary pressure on protein function. Of these environments, none is more complex than the cell itself, where proteins function in the cytosol, in membranes, and in different compartments. A wide variety of kinetic and thermodynamics experiments, ranging from single-molecule studies to jump kinetics and from nuclear magnetic resonance to imaging on the microscope, have elucidated how protein energy landscapes facilitate folding and how they are subject to evolutionary constraints and environmental perturbation. Here we review some recent developments in the field and refer the reader to some original work and additional reviews that cover this broad topic in protein science.

Chaperone-client interactions: Non-specificity engenders multifunctionality

Here, we provide an overview of the different mechanisms whereby three different chaperones, Spy, Hsp70, and Hsp60, interact with folding proteins, and we discuss how these chaperones may guide the folding process. Available evidence suggests that even a single chaperone can use many mechanisms to aid in protein folding, most likely due to the need for most chaperones to bind clients promiscuously. Chaperone mechanism may be better understood by always considering it in the context of the client's folding pathway and biological function.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.

Probing the origins of two-state folding

Many protein systems fold in a two-state manner. Random models, however, rarely display two-state kinetics and thus such behavior should not be accepted as a default. While theories for the prevalence of two-state kinetics have been presented, none sufficiently explain the breadth of experimental observations. A model, making minimal assumptions, is introduced that suggests two-state behavior is likely for any system with an overwhelmingly populated native state. We show two-state folding is a natural consequence of such two-state thermodynamics, and is strengthened by increasing the population of the native state. Further, the model exhibits hub-like behavior, with slow interconversions between unfolded states. Despite this, the unfolded state equilibrates quickly relative to the folding time. This apparent paradox is readily understood through this model. Finally, our results compare favorable with measurements of folding rates as a function of chain length and Keq, providing new insight into these relations.

Filling up the heme pocket stabilizes apomyoglobin and speeds up its folding

Wild type apomyoglobin folds in at least two steps: the ABGH core rapidly, followed much later by the heme-binding CDEF core. We hypothesize that the evolved heme-binding function of the CDEF core frustrates its folding: it has a smaller contact order and is no more complex topologically than ABGH, and thus, it should be able to fold faster. Therefore, filling up the empty heme cavity of apomyoglobin with larger, hydrophobic side chains should significantly stabilize the protein and increase its folding rate. Molecular dynamics simulations allowed us to design four different mutants with bulkier side chains that increase the native bias of the CDEF region. In vitro thermal denaturation shows that the mutations increase folding stability and bring the protein closer to two-state behavior, as judged by the difference of fluorescence- and circular dichroism-detected protein stability. Millisecond stopped flow measurements of the mutants exhibit refolding kinetics that are over 4 times faster than the wild type’s. We propose that myoglobin-like proteins not evolved to bind heme are equally stable, and find an example. Our results illustrate how evolution for function can force proteins to adapt frustrated folding mechanisms, despite having simple topologies.

Evolutionary biochemistry: revealing the historical and physical causes of protein properties

The repertoire of proteins and nucleic acids in the living world is determined by evolution; their properties are determined by the laws of physics and chemistry. Explanations of these two kinds of causality - the purviews of evolutionary biology and biochemistry, respectively - are typically pursued in isolation, but many fundamental questions fall squarely at the interface of fields. Here we articulate the paradigm of evolutionary biochemistry, which aims to dissect the physical mechanisms and evolutionary processes by which biological molecules diversified and to reveal how their physical architecture facilitates and constrains their evolution. We show how an integration of evolution with biochemistry moves us towards a more complete understanding of why biological molecules have the properties that they do.

Stabilization of a protein conferred by an increase in folded state entropy

Entropic stabilization of native protein structures typically relies on strategies that serve to decrease the entropy of the unfolded state. Here we report, using a combination of experimental and computational approaches, on enhanced thermodynamic stability conferred by an increase in the configurational entropy of the folded state. The enhanced stability is observed upon modifications of a loop region in the enzyme acylphosphatase and is achieved despite significant enthalpy losses. The modifications that lead to increased stability, as well as those that result in destabilization, however, strongly compromise enzymatic activity, rationalizing the preservation of the native loop structure even though it does not provide the protein with maximal stability or kinetic foldability.

Sub-millisecond chain collapse of the Escherichia coli globin ApoHmpH

Myoglobins are ubiquitous proteins that play a seminal role in oxygen storage, transport, and NO metabolism. The folding mechanism of apomyoglobins from different species has been studied to a fair extent over the last two decades. However, integrated investigations of the entire process, including both the early (sub-ms) and late (ms-s) folding stages, have been missing. Here, we study the folding kinetics of the single-Trp Escherichia coli globin apoHmpH via a combination of continuous-flow microfluidic and stopped-flow approaches. A rich series of molecular events emerges, spanning a very wide temporal range covering more than 7 orders of magnitude, from sub-microseconds to tens of seconds. Variations in fluorescence intensity and spectral shifts reveal that the protein region around Trp120 undergoes a fast collapse within the 8 μs mixing time and gradually reaches a native-like conformation with a half-life of 144 μs from refolding initiation. There are no further fluorescence changes beyond ca. 800 μs, and folding proceeds much more slowly, up to 20 s, with acquisition of the missing helicity (ca. 30%), long after consolidation of core compaction. The picture that emerges is a gradual acquisition of native structure on a free-energy landscape with few large barriers. Interestingly, the single tryptophan, which lies within the main folding core of globins, senses some local structural consolidation events after establishment of native-like core polarity (i.e., likely after core dedydration). In all, this work highlights how the main core of the globin fold is capable of becoming fully native efficiently, on the sub-millisecond time scale.

Understanding the folding-function tradeoff in proteins

When an amino-acid sequence cannot be optimized for both folding and function, folding can get compromised in favor of function. To understand this tradeoff better, we devise a novel method for extracting the "function-less" folding-motif of a protein fold from a set of structurally similar but functionally diverse proteins. We then obtain the β-trefoil folding-motif, and study its folding using structure-based models and molecular dynamics simulations. CompariA protein sequence serves two purpson with the folding of wild-type β-trefoil proteins shows that function affects folding in two ways: In the slower folding interleukin-1β, binding sites make the fold more complex, increase contact order and slow folding. In the faster folding hisactophilin, residues which could have been part of the folding-motif are used for function. This reduces the density of native contacts in functional regions and increases folding rate. The folding-motif helps identify subtle structural deviations which perturb folding. These may then be used for functional annotation. Further, the folding-motif could potentially be used as a first step in the sequence design of function-less scaffold proteins. Desired function can then be engineered into these scaffolds.

Microsecond folding experiments and simulations: a match is made

For the past two decades, protein folding experiments have been speeding up from the second or millisecond time scale to the microsecond time scale, and full-atom simulations have been extended from the nanosecond to the microsecond and even millisecond time scale. Where the two meet, it is now possible to compare results directly, allowing force fields to be validated and refined, and allowing experimental data to be interpreted in atomistic detail. In this perspective we compare recent experiments and simulations on the microsecond time scale, pointing out the progress that has been made in determining native structures from physics-based simulations, refining experiments and simulations to provide more quantitative underlying mechanisms, and tackling the problems of multiple reaction coordinates, downhill folding, and complex underlying structure of unfolded or misfolded states.

Hysteresis as a Marker for Complex, Overlapping Landscapes in Proteins

Topologically complex proteins fold by multiple routes as a result of hard-to-fold regions of the proteins. Oftentimes these regions are introduced into the protein scaffold for function and increase frustration in the otherwise smooth-funneled landscape. Interestingly, while functional regions add complexity to folding landscapes, they may also contribute to a unique behavior referred to as hysteresis. While hysteresis is predicted to be rare, it is observed in various proteins, including proteins containing a unique peptide cyclization to form a fluorescent chromophore as well as proteins containing a knotted topology in their native fold. Here, hysteresis is demonstrated to be a consequence of the decoupling of unfolding events from the isomerization or hula-twist of a chromophore in one protein and the untying of the knot in a second protein system. The question now is- can hysteresis be a marker for the interplay of landscapes where complex folding and functional regions overlap?

Configuration-dependent diffusion dynamics of downhill and two-state protein folding

Configuration-dependent diffusion (CDD) is important for protein folding kinetics with small thermodynamic barriers. CDD can be even more crucial in downhill folding without thermodynamic barriers. We explored the CDD of a downhill protein (BBL), and a two-state protein (CI2). The hidden kinetic barriers due to CDD were revealed. The increased ~1 k(B)T kinetic barrier is in line with experimental value based on other fast folding proteins. Compared to that of CI2, the effective free-energy profile of BBL is found to be significantly influenced by CDD, and the kinetics are totally determined by diffusion. These findings are consistent with both earlier bulk and single-molecule fluorescence measurements. In addition, we found the temperature dependence of CDD. We also found that the ratio of folding transition temperature against optimal kinetic folding temperature can provide both a quantitative measure for the underlying landscape topography and an indicator for the possible appearance of downhill folding. Our study can help for a better understanding of the role of diffusion in protein folding dynamics.

Strand swapping regulates the iron-sulfur cluster in the diabetes drug target mitoNEET

MitoNEET is a recently identified diabetes drug target that coordinates a transferable 2Fe-2S cluster, and additionally contains an unusual strand swap. In this manuscript, we use a dual basin structure-based model to predict and characterize the folding and functionality of strand swapping in mitoNEET. We demonstrate that a strand unswapped conformation is kinetically accessible and that multiple levels of control are employed to regulate the conformational dynamics of the system. Environmental factors such as temperature can shift route preference toward the unswapped pathway. Additionally we see that a region recently identified as contributing to frustration in folding acts as a regulatory hinge loop that modulates conformational balance. Interestingly, strand unswapping transfers strain specifically to cluster-coordinating residues, opening the cluster-coordinating pocket. Strengthening contacts within the cluster-coordinating pocket opens a new pathway between the swapped and unswapped conformation that utilizes cracking to bypass the unfolded basin. These results suggest that local control within distinct regions affect motions important in regulating mitoNEET's 2Fe-2S clusters.

Extending Tlusty's rate distortion index theorem method to the glycome: do even 'low level' biochemical phenomena require sophisticated cognitive paradigms?

Unlike the universal genetic code and ordered protein folding, direct application of Tlusty's method to the glycome produces a reducto ad absurdum: From the beginning a complicated system of chemical cognition is needed so that external information constrains and tunes what would otherwise be a monstrously large 'glycan code error network'. Further, the glycan manufacture machinery itself must be regulated by yet other levels of chemical cognition to ensure that what is produced matches what was chosen for production. Application of a rate distortion index theorem/operator method at this second stage appears possible, permitting analytic characterization of the complicated 'glycan spectra' associated with cellular interactions and their dynamics. The regulation of 'low level' biochemical processes, ranging from gene expression and protein folding through the production of flexible glycan surface signalling fronds, appears to require systems of chemical cognition whose sophistication may rival that of high order neural process.

Unifying model for two-state and downhill protein folding

A protein-folding model is proposed at the amino acid level, in which the folding process is divided into two successive stages: the rate-determining step, dominated by the "stochastic interactions"of solvent molecules, and the rapid phase, dominated by the "order interactions"among atoms in polypeptide. The master equation approach is used to investigate the folding kinetics, and an analytical treatment of the master equation yields a simple three-parameter expression for folding time. It is found that both two-state and downhill protein-folding kinetics can be described by a unifying model.