NOE data demonstrating a compact unfolded state for an SH3 domain under non-denaturing conditions

Helicity of a tardigrade disordered protein contributes to its protective function during desiccation

To survive extreme drying (anhydrobiosis), many organisms, spanning every kingdom of life, accumulate intrinsically disordered proteins (IDPs). For decades, the ability of anhydrobiosis-related IDPs to form transient amphipathic helices has been suggested to be important for promoting desiccation tolerance. However, evidence empirically supporting the necessity and/or sufficiency of helicity in mediating anhydrobiosis is lacking. Here, we demonstrate that the linker region of CAHS D, a desiccation-related IDP from the tardigrade Hypsibius exemplaris, that contains significant helical structure, is the protective portion of this protein. Perturbing the sequence composition and grammar of the linker region of CAHS D, through the insertion of helix-breaking prolines, modulating the identity of charged residues, or replacement of hydrophobic amino acids with serine or glycine residues results in variants with different degrees of helical structure. Importantly, correlation of protective capacity and helical content in variants generated through different helix perturbing modalities does not show as strong a trend, suggesting that while helicity is important, it is not the only property that makes a protein protective during desiccation. These results provide direct evidence for the decades-old theory that helicity of desiccation-related IDPs is linked to their anhydrobiotic capacity.

The Analytical Flory Random Coil Is a Simple-to-Use Reference Model for Unfolded and Disordered Proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory random coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google Colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

The analytical Flory random coil is a simple-to-use reference model for unfolded and disordered proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory Random Coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

Atomic view of cosolute-induced protein denaturation probed by NMR solvent paramagnetic relaxation enhancement

The cosolvent effect arises from the interaction of cosolute molecules with a protein and alters the equilibrium between native and unfolded states. Denaturants shift the equilibrium toward the latter, while osmolytes stabilize the former. The molecular mechanism whereby cosolutes perturb protein stability is still the subject of considerable debate. Probing the molecular details of the cosolvent effect is experimentally challenging as the interactions are very weak and transient, rendering them invisible to most conventional biophysical techniques. Here, we probe cosolute-protein interactions by means of NMR solvent paramagnetic relaxation enhancement together with a formalism we recently developed to quantitatively describe, at atomic resolution, the energetics and dynamics of cosolute-protein interactions in terms of a concentration normalized equilibrium average of the interspin distance, [Formula: see text], and an effective correlation time, τ_c The system studied is the metastable drkN SH3 domain, which exists in dynamic equilibrium between native and unfolded states, thereby permitting us to probe the interactions of cosolutes with both states simultaneously under the same conditions. Two paramagnetic cosolute denaturants were investigated, one neutral and the other negatively charged, differing in the presence of a carboxyamide group versus a carboxylate. Our results demonstrate that attractive cosolute-protein backbone interactions occur largely in the unfolded state and some loop regions in the native state, electrostatic interactions reduce the [Formula: see text] values, and temperature predominantly impacts interactions with the unfolded state. Thus, destabilization of the native state in this instance arises predominantly as a consequence of interactions of the cosolutes with the unfolded state.

The denatured state of HIV-1 protease under native conditions

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high‐resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV‐1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant‐free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV‐1 protease under native conditions displays rich patterns of transient native and non‐native structures, which could be of relevance to its guidance through a complex folding process.

The Cold-Unfolded State Is Expanded but Contains Long- and Medium-Range Contacts and Is Poorly Described by Homopolymer Models

Cold unfolding of proteins is predicted by the Gibbs-Helmholtz equation and is thought to be driven by a strongly temperature-dependent interaction of protein nonpolar groups with water. Studies of the cold-unfolded state provide insight into protein energetics, partially structured states, and folding cooperativity and are of practical interest in biotechnology. However, structural characterization of the cold-unfolded state is much less extensive than studies of thermally or chemically denatured unfolded states, in large part because the midpoint of the cold unfolding transition is usually below freezing. We exploit a rationally designed point mutation (I98A) in the hydrophobic core of the C-terminal domain of the ribosomal protein L9 that allows the cold denatured state ensemble to be observed above 0 °C at near neutral pH and ambient pressure in the absence of added denaturants. A combined approach consisting of paramagnetic relaxation enhancement measurements, analysis of small-angle X-ray scattering data, all-atom simulations, and polymer theory provides a detailed description of the cold-unfolded state. Despite a globally expanded ensemble, as determined by small-angle X-ray scattering, sequence-specific medium- and long-range interactions in the cold-unfolded state give rise to deviations from homopolymer-like behavior. Our results reveal that the cold-denatured state is heterogeneous with local and long-range intramolecular interactions that may prime the folded state and also demonstrate that significant long-range interactions are compatible with expanded unfolded ensembles. The work also highlights the limitations of homopolymer-based descriptions of unfolded states of proteins.

Extended Experimental Inferential Structure Determination Method in Determining the Structural Ensembles of Disordered Protein States

Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. We introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, that calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R _h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (R_h) of biomolecules in solution. For unfolded proteins, the R_h provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the R_h of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by ¹³C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that R_h values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the R_h of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Binding-induced folding under unfolding conditions: Switching between induced fit and conformational selection mechanisms

The chemistry of protein-ligand binding is the basis of virtually every biological process. Ligand binding can be essential for a protein to function in the cell by stabilizing or altering the conformation of a protein, particularly for partially or completely unstructured proteins. However, the mechanisms by which ligand binding impacts disordered proteins or influences the role of disorder in protein folding is not clear. To gain insight into this question, the mechanism of folding induced by the binding of a Pro-rich peptide ligand to the SH3 domain of phosphatidylinositol 3-kinase unfolded in the presence of urea has been studied using kinetic methods. Under strongly denaturing conditions, folding was found to follow a conformational selection (CS) mechanism. However, under mildly denaturing conditions, a ligand concentration-dependent switch in the mechanism was observed. The folding mechanism switched from being predominantly a CS mechanism at low ligand concentrations to being predominantly an induced fit (IF) mechanism at high ligand concentrations. The switch in the mechanism manifests itself as an increase in the reaction flux along the IF pathway at high ligand concentrations. The results indicate that, in the case of intrinsically disordered proteins too, the folding mechanism is determined by the concentration of the ligand that induces structure formation.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Cosolutes, Crowding, and Protein Folding Kinetics

Long accepted as the most important interaction, recent work shows that steric repulsions alone cannot explain the effects of macromolecular cosolutes on the equilibrium thermodynamics of protein stability. Instead, chemical interactions have been shown to modulate, and even dominate, crowding-induced steric repulsions. Here, we use 19F NMR to examine the effects of small and large cosolutes on the kinetics of protein folding and unfolding using the metastable 7 kDa N-terminal SH3 domain of the Drosophila signaling protein drk (SH3), which folds by a two-state mechanism. The small cosolutes consist of trimethylamine N-oxide and sucrose, which increase equilibrium protein stability, and urea, which destabilizes proteins. The macromolecules comprise the stabilizing sucrose polymer, Ficoll, and the destabilizing globular protein, lysozyme. We assessed the effects of these cosolutes on the differences in free energy between the folded state and the transition state and between the unfolded ensemble and the transition state. We then examined the temperature dependence to assess changes in activation enthalpy and entropy. The enthalpically mediated effects are more complicated than suggested by equilibrium measurements. We also observed enthalpic effects with the supposedly inert sucrose polymer, Ficoll, that arise from its macromolecular nature. Assessment of activation entropies shows important contributions from solvent and cosolute, in addition to the configurational entropy of the protein that, again, cannot be gleaned from equilibrium data. Comparing the effects of Ficoll to those of the more physiologically relevant cosolute lysozyme reveals that synthetic polymers are not appropriate models for understanding the kinetics of protein folding in cells.

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.

Probing the Action of Chemical Denaturant on an Intrinsically Disordered Protein by Simulation and Experiment

Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.

Steric trapping reveals a cooperativity network in the intramembrane protease GlpG

Membrane proteins are assembled through balanced interactions among protein, lipids and water. Studying their folding while maintaining the native lipid environment is necessary but challenging. Here we present methods for analyzing key elements in membrane protein folding including thermodynamic stability, compactness of the unfolded state and folding cooperativity under native conditions. The methods are based on steric trapping which couples unfolding of a doubly-biotinylated protein to binding of monovalent streptavidin (mSA). We further advanced this technology for general application by developing versatile biotin probes possessing spectroscopic reporters that are sensitized by mSA binding or protein unfolding. By applying these methods to an intramembrane protease GlpG of Escherichia coli, we elucidated a widely unraveled unfolded state, subglobal unfolding of the region encompassing the active site, and a network of cooperative and localized interactions to maintain the stability. These findings provide crucial insights into the folding energy landscape of membrane proteins.

Crystallographic studies on protein misfolding: Domain swapping and amyloid formation in the SH3 domain

Oligomerization by 3D domain swapping is found in a variety of proteins of diverse size, fold and function. In the early 1960s this phenomenon was postulated for the oligomers of ribonuclease A, but it was not until the 1990s that X-ray diffraction provided the first experimental evidence of this special manner of oligomerization. Nowadays, structural information has allowed the identification of these swapped oligomers in over one hundred proteins. Although the functional relevance of this phenomenon is not clear, this alternative folding of protomers into intertwined oligomers has been related to amyloid formation. Studies on proteins that develop 3D domain swapping might provide some clues on the early stages of amyloid formation. The SH3 domain is a small modular domain that has been used as a model to study the basis of protein folding. Among SH3 domains, the c-Src-SH3 domain emerges as a helpful model to study 3D domain swapping and amyloid formation.

Determination of pseudocontact shifts of low-populated excited states by NMR chemical exchange saturation transfer

Angular and distance restraints for low populated excited conformations are studied using PCS–CEST NMR spectroscopy.

Quantitative assessments of the distinct contributions of polypeptide backbone amides versus side chain groups to chain expansion via chemical denaturation

In aqueous solutions with high concentrations of chemical denaturants such as urea and guanidinium chloride (GdmCl) proteins expand to populate heterogeneous conformational ensembles. These denaturing environments are thought to be good solvents for generic protein sequences because properties of conformational distributions align with those of canonical random coils. Previous studies showed that water is a poor solvent for polypeptide backbones, and therefore, backbones form collapsed globular structures in aqueous solvents. Here, we ask if polypeptide backbones can intrinsically undergo the requisite chain expansion in aqueous solutions with high concentrations of urea and GdmCl. We answer this question using a combination of molecular dynamics simulations and fluorescence correlation spectroscopy. We find that the degree of backbone expansion is minimal in aqueous solutions with high concentrations of denaturants. Instead, polypeptide backbones sample conformations that are denaturant-specific mixtures of coils and globules, with a persistent preference for globules. Therefore, typical denaturing environments cannot be classified as good solvents for polypeptide backbones. How then do generic protein sequences expand in denaturing environments? To answer this question, we investigated the effects of side chains using simulations of two archetypal sequences with amino acid compositions that are mixtures of charged, hydrophobic, and polar groups. We find that side chains lower the effective concentration of backbone amides in water leading to an intrinsic expansion of polypeptide backbones in the absence of denaturants. Additional dilution of the effective concentration of backbone amides is achieved through preferential interactions with denaturants. These effects lead to conformational statistics in denaturing environments that are congruent with those of canonical random coils. Our results highlight the role of side chain-mediated interactions as determinants of the conformational properties of unfolded states in water and in influencing chain expansion upon denaturation.

Inactivation and unfolding of protein tyrosine phosphatase from Thermus thermophilus HB27 during urea and guanidine hydrochloride denaturation

The effects of urea and guanidine hydrochloride (GdnHCl) on the activity, conformation and unfolding process of protein tyrosine phosphatase (PTPase), a thermostable low molecular weight protein from Thermus thermophilus HB27, have been studied. Enzymatic activity assays showed both urea and GdnHCl resulted in the inactivation of PTPase in a concentration and time-dependent manner. Inactivation kinetics analysis suggested that the inactivation of PTPase induced by urea and GdnHCl were both monophasic and reversible processes, and the effects of urea and GdnHCl on PTPase were similar to that of mixed-type reversible inhibitors. Far-ultraviolet (UV) circular dichroism (CD), Tryptophan and 1-anilinonaphthalene -8-sulfonic acid (ANS) fluorescence spectral analyses indicated the existence of a partially active and an inactive molten globule-like intermediate during the unfolding processes induced by urea and GdnHCl, respectively. Based on the sequence alignment and the homolog Tt1001 protein structure, we discussed the possible conformational transitions of PTPase induced by urea and GdnHCl and compared the conformations of these unfolding intermediates with the transient states in bovine PTPase and its complex structures in detail. Our results may be able to provide some valuable clues to reveal the relationship between the structure and enzymatic activity, and the unfolding pathway and mechanism of PTPase.

Energetically significant networks of coupled interactions within an unfolded protein

Unfolded and partially unfolded proteins participate in a wide range of biological processes from pathological aggregation to the regulation of normal cellular activity. Unfolded states can be populated under strongly denaturing conditions, but the ensemble which is relevant for folding, stability, and aggregation is that populated under physiological conditions. Characterization of nonnative states is critical for the understanding of these processes, yet comparatively little is known about their energetics and their structural propensities under native conditions. The standard view is that energetically significant coupled interactions involving multiple residues are generally not present in the denatured state ensemble (DSE) or in intrinsically disordered proteins. Using the N-terminal domain of the ribosomal protein L9, a small α-β protein, as an experimental model system, we demonstrate that networks of energetically significant, coupled interactions can form in the DSE of globular proteins, and can involve residues that are distant in sequence and spatially well separated in the native structure. X-ray crystallography, NMR, dynamics studies, native state pKa measurements, and thermodynamic analysis of more than 25 mutants demonstrate that residues are energetically coupled in the DSE. Altering these interactions by mutation affects the stability of the domain. Mutations that alter the energetics of the DSE can impact the analysis of cooperativity and folding, and may play a role in determining the propensity to aggregate.

Towards a true protein movie: a perspective on the potential impact of the ensemble-based structure determination using exact NOEs

Confined by the Boltzmann distribution of the energies of the states, a multitude of structural states are inherent to biomolecules. For a detailed understanding of a protein's function, its entire structural landscape at atomic resolution and insight into the interconversion between all the structural states (i.e. dynamics) are required. Whereas dedicated trickery with NMR relaxation provides aspects of local dynamics, and 3D structure determination by NMR is well established, only recently have several attempts been made to formulate a more comprehensive description of the dynamics and the structural landscape of a protein. Here, a perspective is given on the use of exact NOEs (eNOEs) for the elucidation of structural ensembles of a protein describing the covered conformational space.

Denatured state ensembles with the same radii of gyration can form significantly different long-range contacts

Defining the structural, dynamic, and energetic properties of the unfolded state of proteins is critical for an in-depth understanding of protein folding, protein thermodynamics, and protein aggregation. Here we analyze long-range contacts and compactness in two apparently fully unfolded ensembles of the same protein: the acid unfolded state of the C-terminal domain of ribosomal protein L9 in the absence of high concentrations of urea as well as the urea unfolded state at low pH. Small angle X-ray scattering reveals that the two states are expanded with values of Rg differing by <7%. Paramagnetic relaxation enhancement (PRE) nuclear magnetic resonance studies, however, reveal that the acid unfolded state samples conformations that facilitate contacts between residues that are distant in sequence while the urea unfolded state ensemble does not. The experimental PRE profiles for the acid unfolded state differ significantly from these predicted using an excluded volume limit ensemble, but these long-range contacts are largely eliminated by the addition of 8 M urea. The work shows that expanded unfolded states can sample very different distributions of long-range contacts yet still have similar radii of gyration. The implications for protein folding and for the characterization of unfolded states are discussed.

Mining residue contacts in proteins using local structure predictions

In this paper we develop data mining techniques to predict 3D contact potentials among protein residues (or amino acids) based on the hierarchical nucleation-propagation model of protein folding. We apply a hybrid approach, using a hidden Markov model to extract folding initiation sites, and then apply association mining to discover contact potentials. The new hybrid approach achieves accuracy results better than those reported previously.

Ensemble modeling of protein disordered states: experimental restraint contributions and validation

Disordered states of proteins include the biologically functional intrinsically disordered proteins and the unfolded states of normally folded proteins. In recent years, ensemble-modeling strategies using various experimental measurements as restraints have emerged as powerful means for structurally characterizing disordered states. However, these methods are still in their infancy compared with the structural determination of folded proteins. Here, we have addressed several issues important to ensemble modeling using our ENSEMBLE methodology. First, we assessed how calculating ensembles containing different numbers of conformers affects their structural properties. We find that larger ensembles have very similar properties to smaller ensembles fit to the same experimental restraints, thus allowing a considerable speed improvement in our calculations. In addition, we analyzed the contributions of different experimental restraints to the structural properties of calculated ensembles, enabling us to make recommendations about the experimental measurements that should be made for optimal ensemble modeling. The effects of different restraints, most significantly from chemical shifts, paramagnetic relaxation enhancements and small-angle X-ray scattering, but also from other data, underscore the importance of utilizing multiple sources of experimental data. Finally, we validate our ENSEMBLE methodology using both cross-validation and synthetic experimental restraints calculated from simulated ensembles. Our results suggest that secondary structure and molecular size distribution can generally be modeled very accurately, whereas the accuracy of calculated tertiary structure is dependent on the number of distance restraints used.

Analysis of electrostatic interactions in the denatured state ensemble of the N-terminal domain of L9 under native conditions

The pH dependence of protein stability is defined by the difference in the number of protons bound to the folded state and to the denatured state ensemble (DSE) as a function of pH. In many cases, the protonation behavior can be described as the sum of a set of independently titrating residues; in this case, the pH dependence of stability reflects differences in folded and DSE pK(a)'s. pH dependent stability studies have shown that there are energetically important interactions involving charged residues in the DSE of the N-terminal domain of L9 (NTL9), which affect significantly the stability of the protein. The DSE of wild type NTL9 cannot be directly characterized under native conditions because of its high stability. A destabilized double mutant of NTL9, V3AI4A, significantly populates the folded state and the DSE in the absence of denaturant. The two states are in slow exchange on the nuclear magnetic resonance time scale, and diffusion measurements indicate that the DSE is compact. The DSE pK(a)'s of all of the acidic residues were directly determined. The DSE pK(a) of Asp8 and Asp23 are depressed relative to model compounds values. Use of the mutant DSE pK(a)'s together with known native state pK(a)'s leads to a significantly improved agreement between the measured pH dependent stability and that predicted by the Tanford-Wyman linkage relationship. An analysis of the literature suggests that DSE interactions involving charged residues are relatively common and should be considered in discussions of protein stability.

Visualization of early events in acetic acid denaturation of HIV-1 protease: a molecular dynamics study

Protein denaturation plays a crucial role in cellular processes. In this study, denaturation of HIV-1 Protease (PR) was investigated by all-atom MD simulations in explicit solvent. The PR dimer and monomer were simulated separately in 9 M acetic acid (9 M AcOH) solution and water to study the denaturation process of PR in acetic acid environment. Direct visualization of the denaturation dynamics that is readily available from such simulations has been presented. Our simulations in 9 M AcOH reveal that the PR denaturation begins by separation of dimer into intact monomers and it is only after this separation that the monomer units start denaturing. The denaturation of the monomers is flagged off by the loss of crucial interactions between the α-helix at C-terminal and surrounding β-strands. This causes the structure to transit from the equilibrium dynamics to random non-equilibrating dynamics. Residence time calculations indicate that denaturation occurs via direct interaction of the acetic acid molecules with certain regions of the protein in 9 M AcOH. All these observations have helped to decipher a picture of the early events in acetic acid denaturation of PR and have illustrated that the α-helix and the β-sheet at the C-terminus of a native and functional PR dimer should maintain both the stability and the function of the enzyme and thus present newer targets for blocking PR function.

Infrared study of the folding mechanism of a helical hairpin: porcine PYY

The helical hairpin motif plays a key role as a receptor site in DNA binding and protein-protein interactions. Thus, various helical hairpins have recently been developed to assess the factors that control the DNA and/or protein binding affinities of this structural motif and to form synthetic templates for protein and drug design. In addition, several lines of evidence suggest that rapid acquisition of a helical hairpin structure from the unfolded ensemble may guide the rapid formation of helical proteins. Despite its importance as a crucial structural element in protein folding and binding, the folding mechanism of the helical hairpin motif has not been thoroughly studied. Herein, we investigate the structural determinants of the folding kinetics of a naturally occurring helical hairpin (porcine PYY) that is free of disulfide bonds and metal ion-induced cross-links using an infrared temperature-jump technique. It is found that mutations in the turn region predominantly increase the barrier of folding irrespective of the temperature, whereas the effect of mutations that perturb the hydrophobic interactions between the two helices is temperature-dependent. At low temperatures, deletion of hydrophobic side chains is found to predominantly affect the unfolding rate, while the opposite is observed at high temperatures. These results are interpreted in terms of a folding mechanism in which the turn is formed in the transition state and also based on the assumption that cross-strand hydrophobic contacts exist in the thermally unfolded state of PYY.

Protein folding, protein collapse, and tanford's transfer model: lessons from single-molecule FRET

The essential and nontrivial role of the denatured state of proteins in their folding reaction is being increasingly scrutinized in recent years. Single molecule FRET (smFRET) experiments show that the denatured state undergoes a continuous collapse (or coil-to-globule) transition as the concentration of a chemical denaturant is decreased, suggesting that conformational entropy of the denatured state is an important part of the free energy of folding. Such observations question the validity of the classical Tanford transfer model, which suggests that the folding free energy can be understood solely based on the difference in amino acid solvation between the folded state and a fixed unfolded state. An alternative to the transfer model is obtained here from a polymer theoretical analysis of a series of published smFRET data. The analysis shows that the free energy of denatured-state collapse has a linear dependence on denaturant concentration, an outcome of the interplay between enthalpic and entropic contributions. Surprisingly, the slope of the free energy of collapse agrees very well with the respective slope of the free energy of folding. This conformity of values obtained from two very different measurements shows that it is the collapse transition in the denatured state which mediates the effect of denaturants on folding. The energetics of folding are thus governed by the competition of solvation and conformational entropy in the denatured state.

The cold denatured state of the C-terminal domain of protein L9 is compact and contains both native and non-native structure

Cold denaturation is a general property of globular proteins, and the process provides insight into the origins of the cooperativity of protein folding and the nature of partially folded states. Unfortunately, studies of protein cold denaturation have been hindered by the fact that the cold denatured state is normally difficult to access experimentally. Special conditions such as addition of high concentrations of denaturant, encapsulation into reverse micelles, the formation of emulsified solutions, high pressure, or extremes of pH have been applied, but these can perturb the unfolded state of proteins. The cold denatured state of the C-terminal domain of the ribosomal protein L9 can be populated under native-like conditions by taking advantage of a destabilizing point mutation which leads to cold denaturation at temperatures above 0 degrees C. This state is in slow exchange with the native state on the NMR time scale. Virtually complete backbone (15)N, (13)C, and (1)H as well as side-chain (13)C(beta) and (1)H(beta) chemical shift assignments were obtained for the cold denatured state at pH 5.7, 12 degrees C. Chemical shift analysis, backbone N-H residual dipolar couplings, amide proton NOEs, and R(2) relaxation rates all indicate that the cold denatured state of CTL9 (the C-terminal domain of the ribosomal protein L9) not only contains significant native-like secondary structure but also non-native structure. The regions corresponding to the two native alpha-helices show a strong tendency to populate helical Phi and Psi angles. The segment which connects alpha-helix 2 and beta-strand 2 (residues 107-124) in the native state exhibits a significant preference to form non-native helical structure in the cold denatured state. The structure observed in the cold denatured state of the I98A mutant is similar to that observed in the pH 3.8 unfolded state of wild type CTL9 at 25 degrees C, suggesting that it is a robust feature of the denatured state ensemble of this protein. The implications for protein folding and for studies of cold denatured states are discussed.

Paramagnetic relaxation enhancements in unfolded proteins: theory and application to drkN SH3 domain

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as (1)H(N) is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the (1)H(N) residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between (1)H(N) spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.

The unfolded state of the C-terminal domain of the ribosomal protein L9 contains both native and non-native structure

Interest in the structural and dynamic properties of unfolded proteins has increased in recent years owing to continued interest in protein folding and misfolding. Knowledge of the unfolded state under native conditions is particularly important for obtaining a complete picture of the protein folding process. The C-terminal domain of protein L9 is a globular alpha, beta protein with an unusual mixed parallel and antiparallel beta-strand structure. The folding kinetics and equilibrium unfolding of CTL9 strongly depend on pH, and follow a simple two state model. Both the native and the unfolded state can be significantly populated at pH 3.8 in the absence of denaturant, allowing the native state and the unfolded state to be characterized under identical conditions. Backbone (15)N, (13)C, (1)H and side-chain (13)C(beta), (1)H(beta) chemical shifts, amide proton NOEs, and (15)N R(2) relaxation rates were obtained for the two conformational states at pH 3.8. All the data indicate that the pH 3.8 native state is well folded and is similar to the native state at neutral pH. There is significant residual structure in the pH 3.8 unfolded state. The regions corresponding to the two native state alpha-helices show strong preference to populate helical phi and psi angles. The segment that connects alpha-helix 2 and beta-strand 2 has a significant tendency to form non-native alpha-helical structure. Comparison with the pH 2.0 unfolded state and the urea unfolded state indicates that the tendency to adopt both native and non-native helical structure is stronger at pH 3.8, demonstrating that the unfolded state of CTL9 under native-like conditions is more structured. The implications for the folding of CTL9 are discussed.

An unfolded CH1 domain controls the assembly and secretion of IgG antibodies

A prerequisite for antibody secretion and function is their assembly into a defined quaternary structure, composed of two heavy and two light chains for IgG. Unassembled heavy chains are actively retained in the endoplasmic reticulum (ER). Here, we show that the C(H)1 domain of the heavy chain is intrinsically disordered in vitro, which sets it apart from other antibody domains. It folds only upon interaction with the light-chain C(L) domain. Structure formation proceeds via a trapped intermediate and can be accelerated by the ER-specific peptidyl-prolyl isomerase cyclophilin B. The molecular chaperone BiP recognizes incompletely folded states of the C(H)1 domain and competes for binding to the C(L) domain. In vivo experiments demonstrate that requirements identified for folding the C(H)1 domain in vitro, including association with a folded C(L) domain and isomerization of a conserved proline residue, are essential for antibody assembly and secretion in the cell.

Structure and disorder in an unfolded state under nondenaturing conditions from ensemble models consistent with a large number of experimental restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include (1)H-(15)N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O(2)-induced (13)C paramagnetic shifts, hydrogen-exchange protection factors, and (15)N R(2) data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu(2+)-Ni(2+) binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative alpha-helical structure and both native-like and nonnative tertiary structure.

Identification of the folding inhibitors of hen-egg lysozyme: gathering the right tools

The unfolded state of proteins displays a surprisingly rich amount of local native structure, which appears to be critical for driving the protein to its native state. Peptides with the same sequence of the corresponding structured segments can be used to interfere with the correct folding of the protein. Using model simulations, we investigate the folding of hen-egg lysozyme, identifying its key segments. Activity assays, NMR and circular dichroism experiments are used to screen the peptides which are able to inhibit the folding of lysozyme. Few peptides, corresponding to the segments of the protein which are structured in the unfolded state, are identified to have significant inhibitory effects.

Experimental characterization of the denatured state ensemble of proteins

The traditional view of the denatured state ensemble of proteins is that it behaves as a classic random coil. This model has important implications for the analysis of protein stability, protein folding, and cooperativity; namely that the effects of mutations on the free energy of the denatured state ensemble can be ignored. This assumption, which is still routinely made, at least at the implicit level, greatly simplifies the analysis of such experiments. However it has long been recognized that the denatured state ensemble (DSE) of real proteins is often quite different from a random coil and can exhibit significant structural preferences. In some cases parts of the chain can even adopt relatively well-defined conformations, particularly under native conditions. Well-studied examples of DSE interactions include elements of hydrogen-bonded secondary structure, particularly helices or turns, as well hydrophobic clusters, hydrophobic aromatic clusters, and more recently interactions involving charged residues. Deviations from random-coil behavior are of practical importance if they influence protein folding, stability, or function, or if they compromise our analysis and interpretation of experiments. The existence of residual structure in the DSE naturally leads to the question of its role in protein folding and stability, and raises the possibility that some mutations could exert a significant part of their effect by altering the DSE. Much of our understanding of the interactions governing protein stability and the folding process have been generated by mutational studies; thus, a detailed understanding of the denatured state ensemble is critical.

Identification and characterization of folding inhibitors of hen egg lysozyme: an example of a new paradigm of drug design

Studies of protein folding indicate the presence of native contacts in the denatured state, giving rise to folding elements which contribute to the accomplishment of the native state. The possibility of finding molecules which can interact with specific folding elements of a target protein preventing it from reaching its native state, and hence from becoming biologically active, is particularly attractive. The notion that folding elements not only provide molecular recognition directing the folding process, but also have conserved sequence, implies that targeting such elements will make protein folding inhibitors less susceptible to mutations which, in many cases, abrogate drug effects. The folding-inhibition strategy can lead to a truly novel and rational approach to drug design, aside from providing new insight into folding. This is illustrated in the case of hen egg lysozyme.

The denatured state of N-PGK is compact and predominantly disordered

The organisation of the structure present in the chemically denatured N-terminal domain of phosphoglycerate kinase (N-PGK) has been determined by paramagnetic relaxation enhancements (PREs) to define the conformational landscape accessible to the domain. Below 2.0 M guanidine hydrochloride (GuHCl), a species of N-PGK (denoted I(b)) is detected, distinct from those previously characterised by kinetic experiments [folded (F), kinetic intermediate (I(k)) and denatured (D)]. The transition to I(b) is never completed at equilibrium, because F predominates below 1.0 M GuHCl. Therefore, the ability of PREs to report on transient or low population species has been exploited to characterise I(b). Five single cysteine variants of N-PGK were labelled with the nitroxide electron spin-label MTSL [(1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-methyl)methanesulfonate] and the denaturant dependences of the relaxation properties of the amide NMR signals between 1.2 and 3.6 M GuHCl were determined. Significant PREs for I(b) were obtained, but these were distributed almost uniformly throughout the sequence. Furthermore, the PREs indicate that no specific short tertiary contacts persist. The data indicate a collapsed state with no coherent three-dimensional structure, but with a restricted radius beyond which the protein chain rarely reaches. The NMR characteristics of I(b) indicate that it forms from the fully denatured state within 100 micros, and therefore a rapid collapse is the initial stage of folding of N-PGK from its chemically denatured state. By extrapolation, I(b) is the predominant form of the denatured state under native conditions, and the non-specifically collapsed structure implies that many non-native contacts and chain reversals form early in protein folding and must be broken prior to attaining the native state topology.

A shorter peptide model from staphylococcal nuclease for the folding-unfolding equilibrium of a beta-hairpin shows that unfolded state has significant contribution from compact conformational states

It is important to understand the conformational features of the unfolded state in equilibrium with folded state under physiological conditions. In this paper, we consider a short peptide model LMYKGQPM from staphylococcal nuclease to model the conformational equilibrium between a hairpin conformation and its unfolded state using molecular dynamics simulation under NVT conditions at 300K using GROMOS96 force field. The free energy landscape has overall funnel-like shape with hairpin conformations sampling the minima. The "unfolded" state has a higher free energy of approximately 12kJ/mol with respect to native hairpin minimum and occupies a plateau region. We find that the unfolded state has significant contributions from compact conformations. Many of these conformations have hairpin-like topology. Further, these compact conformational forms are stabilized by hydrophobic interactions. Conversion between native and non-native hairpins occurs via unfolded states. Frequent conversions between folded and unfolded hairpins are observed with single exponential kinetics. We compare our results with the emerging picture of unfolded state from both experimental and theoretical studies.