Electrostatic effects on funneled landscapes and structural diversity in denatured protein ensembles

Dramatic Shape Changes Occur as Cytochrome c Folds

Extensive experimental studies on the folding of cytochrome c (Cyt c) make this small protein an ideal target for atomic detailed simulations for the purposes of quantitatively characterizing the structural transitions and the associated time scales for folding to the native state from an ensemble of unfolded states. We use previously generated atomically detailed folding trajectories by the stochastic difference equation in length to calculate the time-dependent changes in the small-angle X-ray scattering (SAXS) profiles. Excellent agreement is obtained between experiments and simulations for the time-dependent SAXS spectra, allowing us to identify the structures of the folding intermediates, which shows that Cyt c reaches the native state by a sequential folding mechanism. Using the ensembles of structures along the folding pathways, we show that compaction and the sphericity of Cyt c change dramatically from the prolate ellipsoid shape in the unfolded state to the spherical native state. Our data, which are in unprecedented quantitative agreement with all aspects of time-resolved SAXS experiments, show that hydrophobic collapse and amide group protection coincide on the 100 microseconds time scale, which is in accordance with ultrafast hydrogen/deuterium exchange studies. Based on these results, we propose that compaction of polypeptide chains, accompanied by dramatic shape changes, is a universal characteristic of globular proteins, regardless of the underlying folding mechanism.

Unfolding cytochromes c-b562 and Rd apo b562

We have analyzed the early stages of unfolding of cytochromes c-b₅₆₂ (PDB ID: 2BC5) and Rd apo b₅₆₂ (PDB ID: 1YYJ). Our geometrical approach proceeds from an analysis of the crystal structure reported for each protein. We quantify, residue-by-residue and region-by-region, the spatial and angular changes in the structure as the protein denatures, and quantify differences that result from the seven residues that differ in the two proteins. Using two independent analyses, one based on spatial metrics and the second on angular metrics, we establish the order of unfolding of the five helices in cyt c-b₅₆₂ and the four helices in the apo protein. For the two helices nearest the N-terminal end of both proteins, the ones in the apo protein unfold first. For the two helices nearest the C-terminal end, the interior helix of the apo protein unfolds first, whereas the terminal helix of the holo protein unfolds first. Excluded-volume effects (repulsive interactions) are minimized in turning regions; the overall range in Δ values is Δ = 36.3 Å³ for cyt c-b₅₆₂ and Δ = 36.6 Å³ for the apo protein, whereas the span for all 20 amino acids is Δ = 167.7 Å³. As our work indicates that the interior helix of cytochrome c-b₅₆₂ is the first to fold, we suggest that this helix protects the heme from misligation, consistent with ultrafast folding over a minimally frustrated funneled landscape.

Funneled angle landscapes for helical proteins

We use crystallographic data for four helical iron proteins (cytochrome c-b₅₆₂, cytochrome c', sperm whale myoglobin, human cytoglobin) to calculate radial and angular signatures as each unfolds from the native state stepwise though four unfolded states. From these data we construct an angle phase diagram to display the evolution of each protein from its native state; and, in turn, the phase diagram is used to construct a funneled angle landscape for comparison with the topography of its folding energy landscape. We quantify the departure of individual helical and turning regions from the areal, angular profile of corresponding regions of the native state. This procedure allows us to identify the similarities and differences among individual helical and turning regions in the early stages of unfolding of the four helical heme proteins.

Are Parallel Proliferation Pathways Redundant?

Are the receptor tyrosine kinase (RTK) and JAK-STAT-driven proliferation pathways 'parallel' or 'redundant'? And what about those of K-Ras4B versus N-Ras? 'Parallel' proliferation pathways accomplish a similar drug resistance outcome. Thus, are they 'redundant'? In this paper, it is argued that there is a fundamental distinction between 'parallel' and 'redundant'. Cellular proliferation pathways are influenced by the genome sequence, 3D organization and chromatin accessibility, and determined by protein availability prior to cancer emergence. In the opinion presented, if they operate the same downstream protein families, they are redundant; if evolutionary-independent, they are parallel. Thus, RTK and JAK-STAT-driven proliferation pathways are parallel; those of Ras isoforms are redundant. Our Precision Medicine Call to map cancer proliferation pathways is vastly important since it can expedite effective therapeutics.

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.

Synergistic approaches unraveling regulation and aggregation of intrinsically disordered β-amyloids implicated in Alzheimer's disease

Alzheimer's disease is a severe brain illness that causes vast numbers of nerve cells in the brain to die, driven by the production and deposition of amyloid beta (Aβ) peptides. Intrinsically disordered proteins (IDPs) generally lack stable structures and are abundant in nature. Aβ peptide is a well-known IDP with a wide range of oligomeric forms. Dysfunctions in Aβ lead to oligomerization, formation of fibrils, and neurodegenerative disorders or other forms of dementia. In this study, we used replica exchange molecular dynamics (REMD) to elucidate the roles of different osmolytes, particularly urea and trimethylamine N-oxide (TMAO), to study shifts in IDP populations. REMD samples the conformational space efficiently and at physiologically relevant temperatures, compared to conventional molecular dynamics that sample at a constant temperature. Urea is known to minimize the aggregation process, while TMAO is beneficial for its stabilizing action. The two osmolytes displayed characteristic effects on Aβ peptides and resulted in progressive modulation of conformations. The present study underlines the hypothesis of "modulation of conformational ensembles" to explain the regulation and aggregation of IDPs.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Conformational Ensembles of α-Synuclein Derived Peptide with Different Osmolytes from Temperature Replica Exchange Sampling

Intrinsically disordered proteins (IDP) are a class of proteins that do not have a stable three-dimensional structure and can adopt a range of conformations playing various vital functional role. Alpha-synuclein is one such IDP which can aggregate into toxic protofibrils and has been associated largely with Parkinson's disease (PD) along with other neurodegenerative diseases. Osmolytes are small organic compounds that can alter the environment around the proteins by acting as denaturants or protectants for the proteins. In the present study, we have conducted a series of replica exchange molecular dynamics simulations to explore the role of osmolytes, urea which is a denaturant and TMAO (trimethylamine N-oxide), a protecting osmolyte, in aggregation and conformations of the synuclein peptide. We observed that both the osmolytes have significantly distinct impacts on the peptide and led to transitions of the conformations of the peptide from one state to other. Our findings highlighted that urea attenuated peptide aggregation and resulted in the formation of extended peptide structures whereas TMAO led to compact and folded forms of the peptide.

Opposing Intermolecular Tuning of Ca2+ Affinity for Calmodulin by Neurogranin and CaMKII Peptides

We investigated the impact of bound calmodulin (CaM)-target compound structure on the affinity of calcium (Ca²⁺) by integrating coarse-grained models and all-atomistic simulations with nonequilibrium physics. We focused on binding between CaM and two specific targets, Ca²⁺/CaM-dependent protein kinase II (CaMKII) and neurogranin (Ng), as they both regulate CaM-dependent Ca²⁺ signaling pathways in neurons. It was shown experimentally that Ca²⁺/CaM (holoCaM) binds to the CaMKII peptide with overwhelmingly higher affinity than Ca²⁺-free CaM (apoCaM); the binding of CaMKII peptide to CaM in return increases the Ca²⁺ affinity for CaM. However, this reciprocal relation was not observed in the Ng peptide (Ng_13–49), which binds to apoCaM or holoCaM with binding affinities of the same order of magnitude. Unlike the holoCaM-CaMKII peptide, whose structure can be determined by crystallography, the structural description of the apoCaM-Ng_13–49 is unknown due to low binding affinity, therefore we computationally generated an ensemble of apoCaM-Ng_13–49 structures by matching the changes in the chemical shifts of CaM upon Ng_13–49 binding from nuclear magnetic resonance experiments. Next, we computed the changes in Ca²⁺ affinity for CaM with and without binding targets in atomistic models using Jarzynski’s equality. We discovered the molecular underpinnings of lowered affinity of Ca²⁺ for CaM in the presence of Ng_13–49 by showing that the N-terminal acidic region of Ng peptide pries open the β-sheet structure between the Ca²⁺ binding loops particularly at C-domain of CaM, enabling Ca²⁺ release. In contrast, CaMKII peptide increases Ca²⁺ affinity for the C-domain of CaM by stabilizing the two Ca²⁺ binding loops. We speculate that the distinctive structural difference in the bound complexes of apoCaM-Ng_13–49 and holoCaM-CaMKII delineates the importance of CaM’s progressive mechanism of target binding on its Ca²⁺ binding affinities.

Genomic Energy Landscapes

Energy landscape theory, developed in the context of protein folding, provides, to our knowledge, a new perspective on chromosome architecture. We review what has been learned concerning the topology and structure of both the interphase and mitotic chromosomes from effective energy landscapes constructed using Hi-C data. Energy landscape thinking raises new questions about the nonequilibrium dynamics of the chromosome and gene regulation.

How EF-Tu can contribute to efficient proofreading of aa-tRNA by the ribosome

It has long been recognized that the thermodynamics of mRNA–tRNA base pairing is insufficient to explain the high fidelity and efficiency of aminoacyl-tRNA (aa-tRNA) selection by the ribosome. To rationalize this apparent inconsistency, Hopfield proposed that the ribosome may improve accuracy by utilizing a multi-step kinetic proofreading mechanism. While biochemical, structural and single-molecule studies have provided a detailed characterization of aa-tRNA selection, there is a limited understanding of how the physical–chemical properties of the ribosome enable proofreading. To this end, we probe the role of EF-Tu during aa-tRNA accommodation (the proofreading step) through the use of energy landscape principles, molecular dynamics simulations and kinetic models. We find that the steric composition of EF-Tu can reduce the free-energy barrier associated with the first step of accommodation: elbow accommodation. We interpret this effect within an extended kinetic model of accommodation and show how EF-Tu can contribute to efficient and accurate proofreading.

The translation of mRNA by the ribosome is governed by a series of large-scale conformational transitions. Here the authors use MD simulations to demonstrate how the rate of dissociation of elongation factor Tu affects the dynamics of tRNA accommodation and proofreading.

Cytochrome unfolding pathways from computational analysis of crystal structures

We have developed a model to study the role of geometrical factors in influencing the early stages of unfolding in three cytochromes: cyt c', cyt c-b562 and cyt c. Each stage in unfolding is quantified by the spatial extension λ̂i of n-residue segments, and by their angular extension 〈βn〉. Similarities and differences between and among the three cytochromes in the unfolding of helical and non-helical regions can be determined by analyzing the data for each signature separately. Definite conclusions can be drawn when spatial and angular changes are considered in tandem. To facilitate comparisons, we present graphical portraits of the three cytochromes at the same stage of unfolding, and in relation to their native state structures. We also display specific segments at different stages of unfolding to illustrate differences in stability of defined domains thereby allowing us to make specific predictions on the unfolding of corresponding internal and terminal helices in cyt c' and cyt c-b562. Our work accords with an earlier experimental report on the presence and persistence of a hydrophobic core in cyt c.

Mimicking Tyrosine Phosphorylation in Human Cytochrome c by the Evolved tRNA Synthetase Technique

Phosphorylation of tyrosine 48 of cytochrome c is related to a wide range of human diseases due to the pleiotropic role of the heme-protein in cell life and death. However, the structural conformation and physicochemical properties of phosphorylated cytochrome c are difficult to study as its yield from cell extracts is very low and its kinase remains unknown. Herein, we report a high-yielding synthesis of a close mimic of phosphorylated cytochrome c, developed by optimization of the synthesis of the non-canonical amino acid p-carboxymethyl-L-phenylalanine (pCMF) and its efficient site-specific incorporation at position 48. It is noteworthy that the Y48pCMF mutation significantly destabilizes the Fe-Met bond in the ferric form of cytochrome c, thereby lowering the pKa value for the alkaline transition of the heme-protein. This finding reveals the differential ability of the phosphomimic protein to drive certain events. This modified cytochrome c might be an important tool to investigate the role of the natural protein following phosphorylation.

Electrostatics, structure prediction, and the energy landscapes for protein folding and binding

While being long in range and therefore weakly specific, electrostatic interactions are able to modulate the stability and folding landscapes of some proteins. The relevance of electrostatic forces for steering the docking of proteins to each other is widely acknowledged, however, the role of electrostatics in establishing specifically funneled landscapes and their relevance for protein structure prediction are still not clear. By introducing Debye-Hückel potentials that mimic long-range electrostatic forces into the Associative memory, Water mediated, Structure, and Energy Model (AWSEM), a transferable protein model capable of predicting tertiary structures, we assess the effects of electrostatics on the landscapes of thirteen monomeric proteins and four dimers. For the monomers, we find that adding electrostatic interactions does not improve structure prediction. Simulations of ribosomal protein S6 show, however, that folding stability depends monotonically on electrostatic strength. The trend in predicted melting temperatures of the S6 variants agrees with experimental observations. Electrostatic effects can play a range of roles in binding. The binding of the protein complex KIX-pKID is largely assisted by electrostatic interactions, which provide direct charge-charge stabilization of the native state and contribute to the funneling of the binding landscape. In contrast, for several other proteins, including the DNA-binding protein FIS, electrostatics causes frustration in the DNA-binding region, which favors its binding with DNA but not with its protein partner. This study highlights the importance of long-range electrostatics in functional responses to problems where proteins interact with their charged partners, such as DNA, RNA, as well as membranes.

Multiscale enhanced sampling of intrinsically disordered protein conformations

In a recently developed multiscale enhanced sampling (MSES) technique, topology-based coarse-grained (CG) models are coupled to atomistic force fields to enhance the sampling of atomistic protein conformations. Here, the MSES protocol is refined by designing more sophisticated Hamiltonian/temperature replica exchange schemes that involve additional parameters in the MSES coupling restraint potential, to more carefully control how conformations are coupled between the atomistic and CG models. A specific focus is to derive an optimal MSES protocol for simulating conformational ensembles of intrinsically disordered proteins (IDPs). The efficacy of the refined protocols, referred to as MSES-soft asymptote (SA), was evaluated using two model peptides with various levels of residual helicities. The results show that MSES-SA generates more reversible helix-coil transitions and leads to improved convergence on various ensemble conformational properties. This study further suggests that more detailed CG models are likely necessary for more effective sampling of local conformational transition of IDPs. © 2015 Wiley Periodicals, Inc.

Topology, structures, and energy landscapes of human chromosomes

Chromosome conformation capture experiments provide a rich set of data concerning the spatial organization of the genome. We use these data along with a maximum entropy approach to derive a least-biased effective energy landscape for the chromosome. Simulations of the ensemble of chromosome conformations based on the resulting information theoretic landscape not only accurately reproduce experimental contact probabilities, but also provide a picture of chromosome dynamics and topology. The topology of the simulated chromosomes is probed by computing the distribution of their knot invariants. The simulated chromosome structures are largely free of knots. Topologically associating domains are shown to be crucial for establishing these knotless structures. The simulated chromosome conformations exhibit a tendency to form fibril-like structures like those observed via light microscopy. The topologically associating domains of the interphase chromosome exhibit multistability with varying liquid crystalline ordering that may allow discrete unfolding events and the landscape is locally funneled toward "ideal" chromosome structures that represent hierarchical fibrils of fibrils.

Regulation and aggregation of intrinsically disordered peptides

Intrinsically disordered proteins (IDPs) are a unique class of proteins that have no stable native structure, a feature that allows them to adopt a wide variety of extended and compact conformations that facilitate a large number of vital physiological functions. One of the most well-known IDPs is the microtubule-associated tau protein, which regulates microtubule growth in the nervous system. However, dysfunctions in tau can lead to tau oligomerization, fibril formation, and neurodegenerative disease, including Alzheimer's disease. Using a combination of simulations and experiments, we explore the role of osmolytes in regulating the conformation and aggregation propensities of the R2/wt peptide, a fragment of tau containing the aggregating paired helical filament (PHF6*). We show that the osmolytes urea and trimethylamine N-oxide (TMAO) shift the population of IDP monomer structures, but that no new conformational ensembles emerge. Although urea halts aggregation, TMAO promotes the formation of compact oligomers (including helical oligomers) through a newly proposed mechanism of redistribution of water around the perimeter of the peptide. We put forth a "superposition of ensembles" hypothesis to rationalize the mechanism by which IDP structure and aggregation is regulated in the cell.

Modulation of folding energy landscape by charge-charge interactions: linking experiments with computational modeling

The kinetics of folding-unfolding of a structurally diverse set of four proteins optimized for thermodynamic stability by rational redesign of surface charge-charge interactions is characterized experimentally. The folding rates are faster for designed variants compared with their wild-type proteins, whereas the unfolding rates are largely unaffected. A simple structure-based computational model, which incorporates the Debye-Hückel formalism for the electrostatics, was used and found to qualitatively recapitulate the experimental results. Analysis of the energy landscapes of the designed versus wild-type proteins indicates the differences in refolding rates may be correlated with the degree of frustration of their respective energy landscapes. Our simulations indicate that naturally occurring wild-type proteins have frustrated folding landscapes due to the surface electrostatics. Optimization of the surface electrostatics seems to remove some of that frustration, leading to enhanced formation of native-like contacts in the transition-state ensembles (TSE) and providing a less frustrated energy landscape between the unfolded and TS ensembles. Macroscopically, this results in faster folding rates. Furthermore, analyses of pairwise distances and radii of gyration suggest that the less frustrated energy landscapes for optimized variants are a result of more compact unfolded and TS ensembles. These findings from our modeling demonstrates that this simple model may be used to: (i) gain a detailed understanding of charge-charge interactions and their effects on modulating the energy landscape of protein folding and (ii) qualitatively predict the kinetic behavior of protein surface electrostatic interactions.

Mechanisms of cellular proteostasis: insights from single-molecule approaches

Cells employ a variety of strategies to maintain proteome homeostasis. Beginning during protein biogenesis, the translation machinery and a number of molecular chaperones promote correct de novo folding of nascent proteins even before synthesis is complete. Another set of molecular chaperones helps to maintain proteins in their functional, native state. Polypeptides that are no longer needed or pose a threat to the cell, such as misfolded proteins and aggregates, are removed in an efficient and timely fashion by ATP-dependent proteases. In this review, we describe how applications of single-molecule manipulation methods, in particular optical tweezers, are shedding new light on the molecular mechanisms of quality control during the life cycles of proteins.

Learning To Fold Proteins Using Energy Landscape Theory

This review is a tutorial for scientists interested in the problem of protein structure prediction, particularly those interested in using coarse-grained molecular dynamics models that are optimized using lessons learned from the energy landscape theory of protein folding. We also present a review of the results of the AMH/AMC/AMW/AWSEM family of coarse-grained molecular dynamics protein folding models to illustrate the points covered in the first part of the article. Accurate coarse-grained structure prediction models can be used to investigate a wide range of conceptual and mechanistic issues outside of protein structure prediction; specifically, the paper concludes by reviewing how AWSEM has in recent years been able to elucidate questions related to the unusual kinetic behavior of artificially designed proteins, multidomain protein misfolding, and the initial stages of protein aggregation.

Cytochrome c: A Multifunctional Protein Combining Conformational Rigidity with Flexibility

Cytochrome has served as a model system for studying redox reactions, protein folding, and more recently peroxidase activity induced by partial unfolding on membranes. This review illuminates some important aspects of the research on this biomolecule. The first part summarizes the results of structural analyses of its active site. Owing to heme-protein interactions the heme group is subject to both in-plane and out-of-plane deformations. The unfolding of the protein as discussed in detail in the second part of this review can be induced by changes of pH and temperature and most prominently by the addition of denaturing agents. Both the kinetic and thermodynamic folding and unfolding involve intermediate states with regard to all unfolding conditions. If allowed to sit at alkaline pH (11.5) for a week, the protein does not return to its folding state when the solvent is switched back to neutral pH. It rather adopts a misfolded state that is prone to aggregation via domain swapping. On the surface of cardiolipin containing liposomes, the protein can adopt a variety of partially unfolded states. Apparently, ferricytochrome c can perform biological functions even if it is only partially folded.

Reduced model captures Mg(2+)-RNA interaction free energy of riboswitches

The stability of RNA tertiary structures depends heavily on Mg(2+). The Mg(2+)-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ2+, the number of excess Mg(2+) associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg(2+) ions interact closely and strongly with the RNA, unlike monovalent ions such as K(+), suggesting that an explicit treatment of Mg(2+) is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg(2+)-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg(2+) through the explicit representation of Mg(2+) ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K(+) and models its competition with condensed Mg(2+). The model contains one fitted parameter, the number of condensed K(+) ions in the absence of Mg(2+). Values of Γ2+ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ2+ to the Mg(2+)-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg(2+) and RNA dynamics.

Preferential water exclusion in protein unfolding

Association of water with protein plays a central role in the latter's folding, structure acquisition, ligand binding, catalytic reactivity, oligomerization, and crystallization. Because these phenomena are also influenced by the net charge content on the protein, the present study examines the association of water with cytochrome c held at different pH values so as to allow its side chains to ionize to variable extents. Equilibrium unfolding of differently charged cytochrome c molecules in water-methanol binary mixtures, where the alcohol acts as the cosolvent denaturant, was used to quantify the preferential exclusion of water during the unfolding transition. The extent of exclusion was found to be related to the net-charge-dependent molecular expansion of the protein in an alcohol-free aqueous medium. The degree of water exclusion was also found to be linearly related to the observed rate of protein unfolding, where the net charge contents of the initial and final states are the same. The results suggest that side-chain ionization, molecular expansion due to charge repulsion, and hence the loss of tertiary contacts lead to additional water-protein association. Protein unfolding rates appear to be linearly correlated with the effective number of water molecules excluded across the end states of unfolding equilibria.

Quantifying the topography of the intrinsic energy landscape of flexible biomolecular recognition

Biomolecular functions are determined by their interactions with other molecules. Biomolecular recognition is often flexible and associated with large conformational changes involving both binding and folding. However, the global and physical understanding for the process is still challenging. Here, we quantified the intrinsic energy landscapes of flexible biomolecular recognition in terms of binding-folding dynamics for 15 homodimers by exploring the underlying density of states, using a structure-based model both with and without considering energetic roughness. By quantifying three individual effective intrinsic energy landscapes (one for interfacial binding, two for monomeric folding), the association mechanisms for flexible recognition of 15 homodimers can be classified into two-state cooperative "coupled binding-folding" and three-state noncooperative "folding prior to binding" scenarios. We found that the association mechanism of flexible biomolecular recognition relies on the interplay between the underlying effective intrinsic binding and folding energy landscapes. By quantifying the whole global intrinsic binding-folding energy landscapes, we found strong correlations between the landscape topography measure Λ (dimensionless ratio of energy gap versus roughness modulated by the configurational entropy) and the ratio of the thermodynamic stable temperature versus trapping temperature, as well as between Λ and binding kinetics. Therefore, the global energy landscape topography determines the binding-folding thermodynamics and kinetics, crucial for the feasibility and efficiency of realizing biomolecular function. We also found "U-shape" temperature-dependent kinetic behavior and a dynamical cross-over temperature for dividing exponential and nonexponential kinetics for two-state homodimers. Our study provides a unique way to bridge the gap between theory and experiments.

Snapshots of a protein folding intermediate

We have investigated the folding dynamics of Thermus thermophilus cytochrome c(552) by time-resolved fluorescence energy transfer between the heme and each of seven site-specific fluorescent probes. We have found both an equilibrium unfolding intermediate and a distinct refolding intermediate from kinetics studies. Depending on the protein region monitored, we observed either two-state or three-state denaturation transitions. The unfolding intermediate associated with three-state folding exhibited native contacts in β-sheet and C-terminal helix regions. We probed the formation of a refolding intermediate by time-resolved fluorescence energy transfer between residue 110 and the heme using a continuous flow mixer. The intermediate ensemble, a heterogeneous mixture of compact and extended polypeptides, forms in a millisecond, substantially slower than the ∼100-μs formation of a burst-phase intermediate in cytochrome c. The surprising finding is that, unlike for cytochrome c, there is an observable folding intermediate, but no microsecond burst phase in the folding kinetics of the structurally related thermostable protein.

Origin of the conformational heterogeneity of cardiolipin-bound cytochrome C

Interactions of cytochrome c (cyt c) with cardiolipin (CL) partially unfold the protein, activating its peroxidase function, a critical event in the execution of apoptosis. However, structural features of the altered protein species in the heterogeneous ensemble are difficult to probe with ensemble averaging. Analyses of the dye-to-heme distance distributions P(r) from time-resolved FRET (TR-FRET) have uncovered two distinct types of CL-bound cyt c conformations, extended and compact. We have combined TR-FRET, fluorescence correlation spectroscopy (FCS), and biolayer interferometry to develop a systematic understanding of the functional partitioning between the two conformations. The two subpopulations are in equilibrium with each other, with a submillisecond rate of conformational exchange reflecting the protein folding into a compact non-native state, as well as protein interactions with the lipid surface. Electrostatic interactions with the negatively charged lipid surface that correlate with physiologically relevant changes in CL concentrations strongly affect the kinetics of cyt c binding and conformational exchange. A predominantly peripheral binding mechanism, rather than deep protein insertion into the membrane, provides a rationale for the general denaturing effect of the CL surface and the large-scale protein unfolding. These findings closely relate to cyt c folding dynamics and suggest a general strategy for extending the time window in monitoring the kinetics of folding.

Nonnative interactions regulate folding and switching of myristoylated protein

We present an integrated experimental and computational study of the molecular mechanisms by which myristoylation affects protein folding and function, which has been little characterized to date. Myristoylation, the covalent linkage of a hydrophobic C14 fatty acyl chain to the N-terminal glycine in a protein, is a common modification that plays a critical role in vital regulated cellular processes by undergoing reversible energetic and conformational switching. Coarse-grained folding simulations for the model pH-dependent actin- and membrane-binding protein hisactophilin reveal that nonnative hydrophobic interactions of the myristoyl with the protein as well as nonnative electrostatic interactions have a pronounced effect on folding rates and thermodynamic stability. Folding measurements for hydrophobic residue mutations of hisactophilin and atomistic simulations indicate that the nonnative interactions of the myristoyl group in the folding transition state are nonspecific and robust, and so smooth the energy landscape for folding. In contrast, myristoyl interactions in the native state are highly specific and tuned for sensitive control of switching functionality. Simulations and amide hydrogen exchange measurements provide evidence for increases as well as decreases in stability localized on one side of the myristoyl binding pocket in the protein, implicating strain and altered dynamics in switching. The effects of folding and function arising from myristoylation are profoundly different from the effects of other post-translational modifications.

Importance of electrostatic interactions in the association of intrinsically disordered histone chaperone Chz1 and histone H2A.Z-H2B

Histone chaperones facilitate assembly and disassembly of nucleosomes. Understanding the process of how histone chaperones associate and dissociate from the histones can help clarify their roles in chromosome metabolism. Some histone chaperones are intrinsically disordered proteins (IDPs). Recent studies of IDPs revealed that the recognition of the biomolecules is realized by the flexibility and dynamics, challenging the century-old structure-function paradigm. Here we investigate the binding between intrinsically disordered chaperone Chz1 and histone variant H2A.Z-H2B by developing a structure-based coarse-grained model, in which Debye-Hückel model is implemented for describing electrostatic interactions due to highly charged characteristic of Chz1 and H2A.Z-H2B. We find that major structural changes of Chz1 only occur after the rate-limiting electrostatic dominant transition state and Chz1 undergoes folding coupled binding through two parallel pathways. Interestingly, although the electrostatic interactions stabilize bound complex and facilitate the recognition at first stage, the rate for formation of the complex is not always accelerated due to slow escape of conformations with non-native electrostatic interactions at low salt concentrations. Our studies provide an ionic-strength-controlled binding/folding mechanism, leading to a cooperative mechanism of “local collapse or trapping” and “fly-casting” together and a new understanding of the roles of electrostatic interactions in IDPs' binding.

Histone chaperones facilitate the assembly and disassembly of nucleosome by interacting with the corresponding histones or histone variants. As the biomolecules in nucleosome are highly charged, electrostatic interactions are particularly important in these processes. The experiments have explored that the histone chaperon Chz1 as an intrinsically disordered protein (IDP) can fold by binding to its histone variants H2A.Z-H2B. Here, we developed a molecular simulation program that treated electrostatic interactions with Debye-Hückel model to study the mechanism of the association. We found that the inter-chain electrostatic interactions facilitate the coupled folding and binding transitions, consistent with the kinetic experiments and microscopic structural perspectives. Furthermore, we show that the intra-chain electrostatic interactions collapse Chz1 and slow the binding rate. The collapsed structure in IDPs caused by intra-chain electrostatic interactions has been widely found in experiments and the effect in binding is well studied in our work. Our theoretical approach shed new light on the role of electrostatics on inter-chain and intra-chain interactions for IDPs' binding and is applicable to the binding/folding of other IDPs to their targets.

Folding without charges

Surface charges of proteins have in several cases been found to function as "structural gatekeepers," which avoid unwanted interactions by negative design, for example, in the control of protein aggregation and binding. The question is then if side-chain charges, due to their desolvation penalties, play a corresponding role in protein folding by avoiding competing, misfolded traps? To find out, we removed all 32 side-chain charges from the 101-residue protein S6 from Thermus thermophilus. The results show that the charge-depleted S6 variant not only retains its native structure and cooperative folding transition, but folds also faster than the wild-type protein. In addition, charge removal unleashes pronounced aggregation on longer timescales. S6 provides thus an example where the bias toward native contacts of a naturally evolved protein sequence is independent of charges, and point at a fundamental difference in the codes for folding and intermolecular interaction: specificity in folding is governed primarily by hydrophobic packing and hydrogen bonding, whereas solubility and binding relies critically on the interplay of side-chain charges.

The ribosome modulates nascent protein folding

Proteins are synthesized by the ribosome and generally must fold to become functionally active. Although it is commonly assumed that the ribosome affects the folding process, this idea has been extremely difficult to demonstrate. We have developed an experimental system to investigate the folding of single ribosome-bound stalled nascent polypeptides with optical tweezers. In T4 lysozyme, synthesized in a reconstituted in vitro translation system, the ribosome slows the formation of stable tertiary interactions and the attainment of the native state relative to the free protein. Incomplete T4 lysozyme polypeptides misfold and aggregate when free in solution, but they remain folding-competent near the ribosomal surface. Altogether, our results suggest that the ribosome not only decodes the genetic information and synthesizes polypeptides, but also promotes efficient de novo attainment of the native state.

Conformational properties of cardiolipin-bound cytochrome c

Interactions of cytochrome c (cyt c) with cardiolipin (CL) are important for both electron transfer and apoptotic functions of this protein. A sluggish peroxidase in its native state, when bound to CL, cyt c catalyzes CL peroxidation, which contributes to the protein apoptotic release. The heterogeneous CL-bound cyt c ensemble is difficult to characterize with traditional structural methods and ensemble-averaged probes. We have employed time-resolved FRET measurements to evaluate structural properties of the CL-bound protein in four dansyl (Dns)-labeled variants of horse heart cyt c. The Dns decay curves and extracted Dns-to-heme distance distributions P(r) reveal a conformational diversity of the CL-bound cyt c ensemble with distinct populations of the polypeptide structures that vary in their degree of protein unfolding. A fraction of the ensemble is substantially unfolded, with Dns-to-heme distances resembling those in the guanidine hydrochloride-denatured state. These largely open cyt c structures likely dominate the peroxidase activity of the CL-bound cyt c ensemble. Site variations in P(r) distributions uncover structural features of the CL-bound cyt c, rationalize previous findings, and implicate the prime role of electrostatic interactions, particularly with the protein C terminus, in the CL-induced unfolding.

Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.

Synergistic folding of two intrinsically disordered proteins: searching for conformational selection

Intrinsically disordered proteins (IDPs) lack stable structures under physiological conditions but often fold into stable structures upon specific binding. These coupled binding and folding processes underlie the organization of cellular regulatory networks, and a mechanistic understanding is thus of fundamental importance. Here, we investigated the synergistic folding of two IDPs, namely, the NCBD domain of transcription coactivator CBP and the p160 steroid receptor coactivator ACTR, using a topology-based model that was carefully calibrated to balance intrinsic folding propensities and intermolecular interactions. As one of the most structured IDPs, NCBD is a plausible candidate that interacts through conformational selection-like mechanisms, where binding is mainly initiated by pre-existing folded-like conformations. Indeed, the simulations demonstrate that, even though binding and folding of both NCBD and ACTR is highly cooperative on the baseline level, the tertiary folding of NCBD is best described by the "extended conformational selection" model that involves multiple stages of selection and induced folding. The simulations further predict that the NCBD/ACTR recognition is mainly initiated by forming a mini folded core that includes the second and third helices of NCBD and ACTR. These predictions are fully consistent with independent physics-based atomistic simulations as well as a recent experimental mapping of the H/D exchange protection factors. The current work thus adds to the limited number of existing mechanistic studies of coupled binding and folding of IDPs, and provides a first direct demonstration of how conformational selection might contribute to efficient recognition of IDPs. Interestingly, even for highly structured IDPs like NCBD, the recognition is initiated by the more disordered C-terminal segment and with substantial contribution from induced folding. Together with existing studies of IDP interaction mechanisms, this argues that induced folding is likely prevalent in IDP-protein interaction, and emphasizes the importance of understanding how IDPs manage to fold efficiently upon (nonspecific) binding. Success of the current study also further supports the notion that, with careful calibration, topology-based models can be effective tools for mechanistic study of IDP interaction and regulation, especially when combined with physics-based atomistic simulations and experiments.

Multi-scaled explorations of binding-induced folding of intrinsically disordered protein inhibitor IA3 to its target enzyme

Biomolecular function is realized by recognition, and increasing evidence shows that recognition is determined not only by structure but also by flexibility and dynamics. We explored a biomolecular recognition process that involves a major conformational change – protein folding. In particular, we explore the binding-induced folding of IA3, an intrinsically disordered protein that blocks the active site cleft of the yeast aspartic proteinase saccharopepsin (YPrA) by folding its own N-terminal residues into an amphipathic alpha helix. We developed a multi-scaled approach that explores the underlying mechanism by combining structure-based molecular dynamics simulations at the residue level with a stochastic path method at the atomic level. Both the free energy profile and the associated kinetic paths reveal a common scheme whereby IA3 binds to its target enzyme prior to folding itself into a helix. This theoretical result is consistent with recent time-resolved experiments. Furthermore, exploration of the detailed trajectories reveals the important roles of non-native interactions in the initial binding that occurs prior to IA3 folding. In contrast to the common view that non-native interactions contribute only to the roughness of landscapes and impede binding, the non-native interactions here facilitate binding by reducing significantly the entropic search space in the landscape. The information gained from multi-scaled simulations of the folding of this intrinsically disordered protein in the presence of its binding target may prove useful in the design of novel inhibitors of aspartic proteinases.

The intrinsically disordered peptide IA3 is the endogenous inhibitor for the enzyme named yeast aspartic proteinase saccharopepsin (YPrA). In the presence of YPrA, IA3 folds itself into an amphipathic helix that blocks the active site cleft of the enzyme. We developed a multi-scaled approach to explore the underlying mechanism of this binding-induced ordering transition. Our approach combines a structure-based molecular dynamics model at the residue level with a stochastic path method at the atomic level. Our simulations suggest that IA3 inhibits YPrA through an induced-fit mechanism where the enzyme (YPrA) induces conformational change of its inhibitor (IA3). This expands the definition of an induced-fit model from its original meaning that the binding of substrate (IA3) drives conformational change in the protein (YPrA). Our result is consistent with recent kinetic experiments and provides a microscopic explanation for the underlying mechanism. We also discuss the important roles of non-native interactions and backtracking. These results enrich our understanding of the enzyme-inhibition mechanism and may have value in the design of drugs.

The origin of nonmonotonic complex behavior and the effects of nonnative interactions on the diffusive properties of protein folding

We present a method for calculating the configurational-dependent diffusion coefficient of a globular protein as a function of the global folding process. Using a coarse-grained structure-based model, we determined the diffusion coefficient, in reaction coordinate space, as a function of the fraction of native contacts formed Q for the cold shock protein (TmCSP). We find nonmonotonic behavior for the diffusion coefficient, with high values for the folded and unfolded ensembles and a lower range of values in the transition state ensemble. We also characterized the folding landscape associated with an energetically frustrated variant of the model. We find that a low-level of frustration can actually stabilize the native ensemble and increase the associated diffusion coefficient. These findings can be understood from a mechanistic standpoint, in that the transition state ensemble has a more homogeneous structural content when frustration is present. Additionally, these findings are consistent with earlier calculations based on lattice models of protein folding and more recent single-molecule fluorescence measurements.

From the Cover: Charge interactions can dominate the dimensions of intrinsically disordered proteins

Many eukaryotic proteins are disordered under physiological conditions, and fold into ordered structures only on binding to their cellular targets. Such intrinsically disordered proteins (IDPs) often contain a large fraction of charged amino acids. Here, we use single-molecule Förster resonance energy transfer to investigate the influence of charged residues on the dimensions of unfolded and intrinsically disordered proteins. We find that, in contrast to the compact unfolded conformations that have been observed for many proteins at low denaturant concentration, IDPs can exhibit a prominent expansion at low ionic strength that correlates with their net charge. Charge-balanced polypeptides, however, can exhibit an additional collapse at low ionic strength, as predicted by polyampholyte theory from the attraction between opposite charges in the chain. The pronounced effect of charges on the dimensions of unfolded proteins has important implications for the cellular functions of IDPs.

Denatured states of low-complexity polypeptide sequences differ dramatically from those of foldable sequences

How the primary sequence of a protein encodes conformational preferences that operate early in folding to promote efficient formation of the correct native topology is still poorly understood. To address this issue, we have prepared a set of yeast iso-1-cytochrome c variants that contain polyalanine inserts ranging from 6 to 30 residues in length near the N terminus of the protein. We study the thermodynamics and kinetics of His-heme loop formation in the denatured state at 3 and 6 M guanidine-HCl concentration. We find that polyalanine closely approximates a random coil with excluded volume giving scaling exponents, nu(3), for equilibrium loop formation of 2.26 +/- 0.13 and 1.97 +/- 0.04 in 3 and 6 M guanidine-HCl, respectively. The rate of loop breakage initially decreases and then becomes independent of loop size as would be expected for a random coil. Comparison with previously reported data for denatured state His-heme loop formation for iso-1-cytochrome c and Rhodopseudomonas palustris cytochrome c', shows that foldable sequences deviate significantly from random coil behavior and that the deviation is fold-dependent.

The folding energy landscape and free energy excitations of cytochrome c

The covalently bound heme cofactor plays a dominant role in the folding of cytochrome c. Because of the complicated inorganic chemistry of the heme, some might consider the folding of cytochrome c to be a special case, following principles different from those used to describe the folding of proteins without cofactors. Recent investigations, however, demonstrate that common models describing folding for many proteins work well for cytochrome c when heme is explicitly introduced, generally providing results that agree with experimental observations. In this Account, we first discuss results from simple native structure-based models. These models include attractive interactions between nonadjacent residues only if they are present in the crystal structure at pH 7. Because attractive nonnative contacts are not included in native structure-based models, their energy landscapes can be described as "perfectly funneled". In other words, native structure-based models are energetically guided towards the native state and contain no energetic traps that would hinder folding. Energetic traps are denoted sources of "frustration", which cause specific transient intermediates to be populated. Native structure-based models do, however, include repulsion between residues due to excluded volume. Nonenergetic traps can therefore exist if the chain, which cannot cross over itself, must partially unfold so that folding can proceed. The ability of native structure-based models to capture this kind of motion is partly responsible for their successful predictions of folding pathways for many types of proteins. Models without frustration describe the sequence of folding events for cytochrome c well (as inferred from hydrogen-exchange experiments), thereby justifying their use as a starting point. At low pH, the experimentally observed folding sequence of cytochrome c deviates from that at pH 7 and from models with perfectly funneled energy landscapes. Here, alternate folding pathways are a result of "chemical frustration". This frustration arises because some regions of the protein are destabilized more than others due to the heterogeneous distribution of titratable residues that are protonated at low pH. Beginning with native structure-based terms, we construct more complex models by adding chemical frustration. These more complex models only modestly perturb the energy landscape, which remains, overall, well funneled. These perturbed models can accurately describe how alternative folding pathways are used at low pH. At alkaline pH, cytochrome c populates distinctly different structural ensembles. For instance, lysine residues are deprotonated and compete for the heme ligation site. The same models that can describe folding at low pH also predict well the structures and relative stabilities of intermediates populated at alkaline pH. The success of models based on funneled energy landscapes suggest that cytochrome c folding is driven primarily by native contacts. The presence of heme appears to add chemical complexity to the folding process, but it does not require fundamental modification of the general principles used to describe folding. Moreover, its added complexity provides a valuable means of probing the folding energy landscape in greater detail than is possible with simpler systems.

Chemical frustration in the protein folding landscape: grand canonical ensemble simulations of cytochrome c

A grand canonical formalism is developed to combine discrete simulations for chemically distinct species in equilibrium. Each simulation is based on a perturbed funneled landscape. The formalism is illustrated using the alkaline-induced transitions of cytochrome c as observed by FTIR spectroscopy and with various other experimental approaches. The grand canonical simulation method accounts for the acid/base chemistry of deprotonation, the inorganic chemistry of heme ligation and misligation, and the minimally frustrated folding energy landscape, thus elucidating the physics of protein folding involved with an acid/base titration of a protein. The formalism combines simulations for each of the relevant chemical species, varying by protonation and ligation states. In contrast to models based on perfectly funneled energy landscapes that contain only contacts found in the native structure, this study introduces "chemical frustration" from deprotonation and misligation that gives rise to many intermediates at alkaline pH. While the nature of these intermediates cannot be easily inferred from available experimental data, this study provides specific structural details of these intermediates, thus extending our understanding of how cytochrome c changes with an increase in pH. The results demonstrate the importance of chemical frustration for understanding biomolecular energy landscapes.

Modulating native-like residual structure in the fully denatured state of photoactive yellow protein affects its refolding

Residual structure in the fully unfolded state is a key element for understanding protein folding. We show that the residual structure in fully denatured photoactive yellow protein (PYP) is affected by isomerization of its p-coumaric acid (pCA) chromophore. The exposure of total surface area and hydrophobic surface area upon unfolding was quantified by denaturant m values and heat capacity changes (DeltaC(p)), respectively. The exposure of the buried surface area upon the unfolding of the acid-denatured state of PYP containing trans-pCA is approximately 20% smaller than that of the native state. In contrast, for the partially unfolded pB photocycle intermediate containing cis-pCA, unfolding-induced exposure of the surface area is not decreased. These results show that pCA photoisomerization reduces residual structure in the fully unfolded state. Thus, residual structure in the fully unfolded state of PYP is under direct experimental control by photoexcitation. The sensitivity of the unfolded state to pCA isomerization provides a novel criterion that residual structure in the unfolded state of PYP is native-like, involving native-like protein-chromophore interactions. A largely untested prediction is that native-like residual structure facilitates the conformational search during folding. In the case of PYP, refolding from the less disordered fully unfolded state containing trans-pCA indeed is substantially accelerated. The burial of hydrophobic surface area in the fully unfolded state suggests that a significant part of the hydrophobic collapse process already has occurred in the denatured state.

Urea-mediated protein denaturation: a consensus view

We have performed all-atom molecular dynamics simulations of three structurally similar small globular proteins in 8 M urea and compared the results with pure aqueous simulations. Protein denaturation is preceded by an initial loss of water from the first solvation shell and consequent in-flow of urea toward the protein. Urea reaches the first solvation shell of the protein mainly due to electrostatic interaction with a considerable contribution coming from the dispersion interaction. Urea shifts the equilibrium from the native to denatured ensemble by making the protein-protein contact less stable than protein-urea contact, which is just the reverse of the condition in pure water, where protein-protein contact is more stable than protein-water contact. We have also seen that water follows urea and reaches the protein interior at later stages of denaturation, while urea preferentially and efficiently solvates different parts of the protein. Solvation of the protein backbone via hydrogen bonding, favorable electrostatic interaction with hydrophilic residues, and dispersion interaction with hydrophobic residues are the key steps through which urea intrudes the core of the protein and denatures it. Why urea is preferred over water for binding to the protein backbone and how urea orients itself toward the protein backbone have been identified comprehensively. All the key components of intermolecular forces are found to play a significant part in urea-induced protein denaturation and also toward the stability of the denatured state ensemble. Changes in water network/structure and dynamical properties and higher degree of solvation of the hydrophobic residues validate the presence of "indirect mechanism" along with the "direct mechanism" and reinforce the effect of urea on protein.

Single-molecule spectroscopy of the temperature-induced collapse of unfolded proteins

We used single-molecule FRET in combination with other biophysical methods and molecular simulations to investigate the effect of temperature on the dimensions of unfolded proteins. With single-molecule FRET, this question can be addressed even under near-native conditions, where most molecules are folded, allowing us to probe a wide range of denaturant concentrations and temperatures. We find a compaction of the unfolded state of a small cold shock protein with increasing temperature in both the presence and the absence of denaturant, with good agreement between the results from single-molecule FRET and dynamic light scattering. Although dissociation of denaturant from the polypeptide chain with increasing temperature accounts for part of the compaction, the results indicate an important role for additional temperature-dependent interactions within the unfolded chain. The observation of a collapse of a similar extent in the extremely hydrophilic, intrinsically disordered protein prothymosin alpha suggests that the hydrophobic effect is not the sole source of the underlying interactions. Circular dichroism spectroscopy and replica exchange molecular dynamics simulations in explicit water show changes in secondary structure content with increasing temperature and suggest a contribution of intramolecular hydrogen bonding to unfolded state collapse.