Hydrogen exchange in a large 29 kD protein and characterization of molten globule aggregation by NMR

The nature of denatured ensembles of the enzyme human carbonic anhydrase (HCA) has been extensively studied by various methods in the past. The protein constitutes an interesting model for folding studies that does not unfold by a simple two-state transition, instead a molten globule intermediate is highly populated at 1.5 M GuHCl. In this work, NMR and H/D exchange studies have been conducted on one of the isozymes, HCA I. The H/D exchange studies, which were enabled by the previously obtained resonance assignment of HCA I, have been used to identify unfolded forms that are accessible from the native state. In addition, the GuHCl-induced unfolded states of HCA I have also been characterized by NMR at GuHCl concentrations in the 0-5 M range. The most important findings in this work are as follows: (1) Amide protons located in the center of the beta-sheet require global unfolding events for efficient H/D exchange. (2) The molten globule and the native state give similar protection against H/D exchange for all of the observable amide protons (i.e., water seems not to efficiently penetrate the interior of the molten globule). (3) At high protein concentrations, the molten globule can form large aggregates, which are not detectable by solution-state NMR methods. (4) The unfolded state (U), present at GuHCl concentrations above 2 M, is composed of an ensemble of conformations having residual structures with different stabilities.     

Describing Partially Unfolded States of Proteins from Sparse NMR Ddata

Proteins involved in signal transduction can usually be present in two states: an inactive and an active (signaling) state. In the case of photoreceptors such as PYP, it has been shown that the signaling state has a large degree of structural and dynamic disorder. Conventional structural NMR approaches present difficulties in describing such partially unfolded states. Owing to the disordered dynamical and transient nature of such states classical NOE-based information, when present, is sparse. Chemical shift changes upon partial unfolding can, however, be easily monitored from HSQC spectra. We show here that such states can be modeled by defining native-like inter-residue contacts for those residues that do not shift significantly upon partial unfolding. The feasibility of this approach is demonstrated using lysozyme as a test case and applied to model the partially unfolded signaling state (pB) of a truncated form of the photoactive yellow protein for which a “classical” NOE-based structure is available for validation. This approach should be generally applicable to systems in which part of the structure remains in a well-defined native-like conformation.     

Test for cooperativity in the early kinetic intermediate in lysozyme folding

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.     

Predicting the equilibrium protein folding pathway: Structure-based analysis of staphylococcal nuclease

The equilibrium folding pathway of staphylococcal nucleas (SNase) has been approximated using a statistical thermodynamic formalism that utilizes the high-resolution structure of the native state as a template to generate a large ensemble of partially folded states. Close to 400,000 different states ranging from the native to the completely unfolded states were included in the analysis. The probability of each state was estimated using an empirical structural parametrization of the folding energetics. It is shown that this formalism predicts accurately the stability of the protein, the cooperativity of the folding/unfolding transition observed by differential scanning calorimetry (DSC) or urea denaturation and the thermodynamic parameters for unfolding. More importantly, this formalism provides a quantitative account of the experimental hydrogen exchange protection factors measured under native conditions for SNase. These results suggest that the computer-generated distribution of states approximates well the ensemble of conformations existing in solution. Furthermore, this formalism represents the first model capable of quantitatively predicting within a unified framework the probability distribution of states seen under native conditions and its change upon unfolding. © 1997 Wiley-Liss, Inc.     

Study of multiple unfolding trajectories and unfolded states of the protein GB1 under the physical property space

Forty nine molecular dynamics simulations of unfolding trajectories of the segment B1 of streptococcal protein G (GB1) provide a direct demonstration of the diversity of unfolding pathway and give a statistically utmost unfolding pathway under the physical property space. Twelve physical properties of the protein were chosen to construct a 12-dimensional property space. Then the 12-dimensional property space was reduced to a 3-dimensional principle component property space. Under the property space, the multiple unfolding trajectories look like "trees", which have some common characters. The "root of the tree" corresponds to the native state, the "bole" homologizes the partially unfolded conformations, and the "crown" is in correspondence to the unfolded state. These unfolding trajectories can be divided into three types. The first one has the characters of straight "bole" and "crown" corresponding to a fast two-state unfolding pathway of GB1. The second one has the character of "the standstill in the middle tree bole", which may correspond to a three-state unfolding pathway. The third one has the character of "the circuitous bole" corresponding to a slow two-state unfolding pathway. The fast two-state unfolding pathway is a statistically utmost unfolding pathway or preferred pathway of GB1, which occupies 53% of 49 unfolding trajectories. In the property space all the unfolding trajectories construct a thermal unfolding pathway ensemble of GB1. The unfolding pathway ensemble resembles a funnel that is gradually emanative from the native state ensemble to the unfolded state ensemble. In the property space, the thermal unfolded state distribution looks like electronic cloud in quantum mechanics. The unfolded states of the independent unfolding simulation trajectories have substantial overlaps, indicating that the thermal unfolded states are confined by the physical property values, and the number of protein unfolded state are much less than that was believed before.     

Abstracts: Albany 2007, The 15th Conversation

Abstract

Forty nine molecular dynamics simulations of unfolding trajectories of the segment B1 of streptococcal protein G (GB1) provide a direct demonstration of the diversity of unfolding pathway and give a statistically utmost unfolding pathway under the physical property space. Twelve physical properties of the protein were chosen to construct a 12-dimensional property space. Then the 12-dimentional property space was reduced to a 3-dimentional principle component property space. Under the property space, the multiple unfolding trajectories look like “trees”, which have some common characters. The “root of the tree” corresponds to the native state, the “bole” homologizes the partially unfolded conformations, and the “crown” is in correspondence to the unfolded state. These unfolding trajectories can be divided into three types. The first one has the characters of straight “bole” and “crown” corresponding to a fast two-state unfolding pathway of GB1. The second one has the character of “the standstill in the middle tree bole”, which may correspond to a three-state unfolding pathway. The third one has the character of “the circuitous bole” corresponding to a slow two-state unfolding pathway. The fast two-state unfolding pathway is a statistically utmost unfolding pathway or preferred pathway of GB1, which occupies 53% of 49 unfolding trajectories. In the property space all the unfolding trajectories construct a thermal unfolding pathway ensemble of GB1. The unfolding pathway ensemble resembles a funnel that is gradually emanative from the native state ensemble to the unfolded state ensemble. In the property space, the thermal unfolded state distribution looks like electronic cloud in quantum mechanics. The unfolded states of the independent unfolding simulation trajectories have substantial overlaps, indicating that the thermal unfolded states are confined by the physical property values, and the number of protein unfolded state are much less than that was believed before.     

Insights into conformation and dynamics of protein GB1 during folding and unfolding by NMR

Understanding protein stability requires characterization of structural determinants of the folded and unfolded states. Many proteins are capable of populating partially folded states under specific solution conditions. Occasionally, coexistence of the folded and an unfolded state under non- or mildly denaturing conditions can be observed by NMR, allowing us to structurally probe these states under identical conditions. Here we report on a destabilized mutant of the B1 domain of protein G (GB1) whose equilibrium unfolding was systematically investigated. Backbone amide residual dipolar couplings (RDCs), the tryptophan Nepsilon-H resonance and the amide nitrogen transverse relaxation rates (R2s) for varying pH values and different temperatures were measured. The backbone amide RDCs indicate that prior to complete unfolding, two melting hot spots are formed at the turn around T11, L12 and K13 and the N terminus of the helix at A24 and T25. The RDCs for the low pH, thermally unfolded state of GB1 are very small and do not indicate the presence of any native-like structure. Amide nitrogen transverse relaxation rates for GB1 in the folded state at different temperatures exhibit large contributions from exchange processes and the associated dynamics display considerable heterogeneity. Our data provide clear evidence for intermediate conformations and multi-state equilibrium un/folding for this GB1 variant.     

Native-like partially folded conformations and folding process revealed in the N-terminal large fragments of staphylococcal nuclease: a study by NMR spectroscopy

The N-terminal large fragments of staphylococcal nuclease (SNase), SNase110 (1-110 residues), SNase121 (1-121 residues), and SNase135 (1-135 residues), and the fragment mutants G88W110, G88W121, V66W110 and V66W121 were studied by heteronuclear multidimensional NMR spectroscopy. Ensembles of co-existent native-like partially folded and unfolded states were observed for fragments. The persistent native-like tertiary interaction drives fragments to be in partially folded states, which reveal native-like beta-barrel conformations. G88W and V66W mutations modulate the extent of inherent native-like tertiary interaction in fragment molecules, and in consequence, fragment mutants fold into native-like beta-subdomain conformations. In cooperation with the inherent tertiary interaction, 2 M TMAO (trimethylamine N-oxide) can promote the folding reaction of fragments through the changes of unfolding free energy, and a native-like beta-subdomain conformation is observed when the chain length contains 135 residues. Heterogeneous partially folded conformations of 1-121 and 1-135 fragments due to cis and trans X-prolyl bond of Lys116-Pro117 make a non-unique folding pathway of fragments. The folding reaction of fragments can be characterized as a hierarchical process.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

The Ac-RGD-NH2 peptide as a probe of slow conformational exchange of short linear peptides in DMSO

According to general belief, the conformational information on short linear peptides in solution derived at ambient temperature from NMR spectrometry represents a population-weighted average over all members of an ensemble of rapidly interconverting conformations. Usually the search for discrete conformations is concentrated at low temperatures especially when sharp NMR resonances are detected at room temperature. Using the peptide Ac-RGD-NH(2) (Ac-Arg-Gly-Asp-NH(2), Ac: acetyl) as a model system and following a new approach, we have been able to demonstrate that short linear peptides can adopt discrete conformational states in DMSO-d(6) (DMSO: dimethylsulfoxide) which vary in a way critically dependent on the reconstitution conditions used before their dissolution in DMSO-d(6). The conformers are stabilized by intramolecular hydrogen bonds, which persist at high temperatures and undergo a very slow exchange with their extended structures in the NMR chemical shift time scale. The reported findings provide clear evidence for the occurrence of solvent-induced conformational exchange and point to DMSO as a valuable medium for folding studies of short linear peptides.     

Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.     

The effect of the polyproline II (PPII) conformation on the denatured state entropy

Polyproline II (PPII) is reported to be a dominant conformation in the unfolded state of peptides, even when no prolines are present in the sequence. Here we use isothermal titration calorimetry (ITC) to investigate the PPII bias in the unfolded state by studying the binding of the SH3 domain of SEM-5 to variants of its putative PPII peptide ligand, Sos. The experimental system is unique in that it provides direct access to the conformational entropy change of the substituted amino acids. Results indicate that the denatured ensemble can be characterized by at least two thermodynamically distinct states, the PPII conformation and an unfolded state conforming to the previously held idea of the denatured state as a random collection of conformations determined largely by hard-sphere collision. The probability of the PPII conformation in the denatured states for Ala and Gly were found to be significant, approximately 30% and approximately 10%, respectively, resulting in a dramatic reduction in the conformational entropy of folding.     

The effect of a DeltaK280 mutation on the unfolded state of a microtubule-binding repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer's disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a DeltaK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and DeltaK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of beta-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.     

Dynamics in the unfolded state of beta2-microglobulin studied by NMR

Many proteins form amyloid-like fibrils in vitro under conditions that favour the population of partially folded conformations or denatured state ensembles. Characterising the structural and dynamic properties of these states is crucial towards understanding the mechanisms of self-assembly in amyloidosis. The aggregation of beta2-microglobulin (beta2m) into amyloid fibrils in vivo occurs in the condition known as dialysis-related amyloidosis (DRA) and the protein has been shown to form amyloid-like fibrils under acidic conditions in vitro. We have used a number of 1H-15N nuclear magnetic resonance (NMR) experiments in conjunction with site-directed mutagenesis to study the acid-unfolded state of beta2m. 15N NMR transverse relaxation experiments reveal that the acid-denatured ensemble, although predominantly unfolded at the N and C termini, contains substantial non-native structure in the central region of the polypeptide chain, stabilised by long-range interactions between aromatic residues and by the single disulphide bond. Relaxation dispersion studies indicate that the acid-unfolded ensemble involves two or more distinct species in conformational equilibrium on the micro- to millisecond time-scale. One of these species appears to be hydrophobically collapsed, as mutations in an aromatic-rich region of the protein, including residues that are solvent-exposed in the native protein, disrupt this structure and cause a consequent decrease in the population of this conformer. Thus, acid-unfolded beta2m consists of a heterogeneous ensemble of rapidly fluctuating species, some of which contain stable, non-native hydrophobic clusters. Given that amyloid assembly of beta2m proceeds with lag kinetics under the conditions of this study, a rarely populated species such as a conformer with non-native aromatic clustering could be key to the initiation of amyloidosis.     

The low-pH unfolded state of the C-terminal domain of the ribosomal protein L9 contains significant secondary structure in the absence of denaturant but is no more compact than the low-pH urea unfolded state

There is considerable interest in the properties of the unfolded states of proteins, particularly unfolded states which can be populated in the absence of high concentrations of denaturants. Interest in the unfolded state ensemble reflects the fact that it is the starting point for protein folding as well as the reference state for protein stability studies and can be the starting state for pathological aggregation. The unfolded state of the C-terminal domain (residues 58-149) of the ribosomal protein L9 (CTL9) can be populated in the absence of denaturant at low pH. CTL9 is a 92-residue globular alpha, beta protein. The low-pH unfolded state contains more secondary structure than the low-pH urea unfolded state, but it is not a molten globule. Backbone ( (1)H, (13)C, and (15)N) NMR assignments as well as side chain (13)C beta and (1)H beta assignments and (15)N R 2 values were obtained for the pH 2.0 unfolded form of CTL9 and for the urea unfolded state at pH 2.5. Analysis of the deviations of the chemical shifts from random coil values indicates that residues that comprise the two helices in the native state show a clear preference for adopting helical phi and psi angles in the pH 2.0 unfolded state. There is a less pronounced but nevertheless clear tendency for residues 107-124 to preferentially populate helical phi and psi values in the unfolded state. The urea unfolded state has no detectable tendency to populate any type of secondary structure even though it is as compact as the pH 2.0 unfolded state. Comparison of the two unfolded forms of CTL9 provides direct experimental evidence that states which differ significantly in their secondary structure can have identical hydrodynamic properties. This in turn demonstrates that global parameters such as R h or R g are very poor indicators of "random coil" behavior.     

Close correspondence between the motions from principal component analysis of multiple HIV-1 protease structures and elastic network modes

The large number of available HIV-1 protease structures provides a remarkable sampling of conformations of the different conformational states, which can be viewed as direct structural information about the dynamics of the HIV-1 protease. After structure matching, we apply principal component analysis (PCA) to obtain the important apparent motions for both bound and unbound structures. There are significant similarities between the first few key motions and the first few low-frequency normal modes calculated from a static representative structure with an elastic network model (ENM), strongly suggesting that the variations among the observed structures and the corresponding conformational changes are facilitated by the low-frequency, global motions intrinsic to the structure. Similarities are also found when the approach is applied to an NMR ensemble, as well as to molecular dynamics (MD) trajectories. Thus, a sufficiently large number of experimental structures can directly provide important information about protein dynamics, but ENM can also provide similar sampling of conformations.     

Design and Structure of an Equilibrium Protein Folding Intermediate: A Hint into Dynamical Regions of Proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.     

Single molecule force spectroscopy reveals a weakly populated microstate of the FnIII domains of tenascin

The native states of proteins exist as an ensemble of conformationally similar microstates. The fluctuations among different microstates are of great importance for the functions and structural stability of proteins. Here, we demonstrate that single molecule atomic force microscopy (AFM) can be used to directly probe the existence of multiple folded microstates. We used the AFM to repeatedly stretch and relax a recombinant tenascin fragment TNfnALL to allow the fibronectin type III (FnIII) domains to undergo repeated unfolding/refolding cycles. In addition to the native state, we discovered that some FnIII domains can refold from the unfolded state into a previously unrecognized microstate, N* state. This novel state is conformationally similar to the native state, but mechanically less stable. The native state unfolds at approximately 120 pN, while the N* state unfolds at approximately 50 pN. These two distinct populations of microstates constitute the ensemble of the folded states for some FnIII domains. An unfolded FnIII domain can fold into either one of the two microstates via two distinct folding routes. These results reveal the dynamic and heterogeneous picture of the folded ensemble for some FnIII domains of tenascin, which may carry important implications for the mechanical functions of tenascins in vivo.