The three states of globular proteins: acid denaturation.

We describe statistical mechanical theory that aims to predict protein stabilities as a function of temperature, pH, and salt concentration, from the physical properties of the constituent amino acids: (1) the number of nonpolar groups, (2) the chain length, (3) the temperature-dependent free energy of transfer, (4) the pKa's (including those in the native state) and their temperature dependencies. We calculate here the phase diagrams for apomyoglobin and hypothetical variant proteins. The theory captures essential features of protein stability including myoglobin's Tm vs pH as measured by P. L. Privalov [(1979) Advances in Protein Chemistry, Vol. 33, pp. 167-241] and its ionic strength vs pH phase diagram as measured by Y. Goto and A. L. Fink [(1990) Journal of Molecular Biology, Vol. 214, pp. 803-805]. The main predictions here are the following: (1) There are three stable states, corresponding to native (N), compact denatured (C), and highly unfolded (U), with transitions between them. (2) In agreement with experiments, the compact denatured state is predicted to have enthalpy closer to U than N because even though there is considerable hydrophobic "clustering" in C, this nevertheless represents a major loss of hydrophobic contacts relative to configurations (N) that have a hydrophobic "core." (3) C becomes more prominent in the phase diagram with increasing nonpolar content or decreasing chain length, perhaps thus accounting for (a) why lysozyme and alpha-lactalbumin differ in their denatured states, and (b) why shortened Staph nuclease molecules are compact. (4) Of major importance for protein calorimetry is Privalov's observation that the enthalpy of folding, delta H (T, pH) is independent of pH. The theory accounts for this through the prediction that the main electrostatic contribution to stability is not enthalpic; the main contribution is the entropy, mainly due to the different distributions of protons and small ions in the native and denatured states.     

Alanine-scanning of the 50's loop in the Clostridium beijerinckii flavodoxin: evaluation of additivity and the importance of interactions provided by the main chain in the modulation of the oxidation-reduction potentials.

The four-residue reverse turn -Met56-Gly-Asp-Glu59- in the Clostridium beijerinckii flavodoxin provides the majority of the critical interactions with the isoalloxazine ring of the flavin mononucleotide (FMN) cofactor that contribute to the binding and the differential stabilization of its three redox states. Direct side chain contacts include the sulfur-ring interaction of Met56, which primarily influences the oxidized and hydroquinone states, and the hydrogen bond by Glu59 with the N3H, which directly (and indirectly through its "anchoring" function) influences all three states to various extents. Involving a novel redox-dependent conformational change, the hydrogen bond formed between the carbonyl group of Gly57 and the N5H of the reduced cofactor strongly influences the stability of the semiquinone state. In this study, the sequential elimination of all side chain interactions in various combinations through a systematic alanine-scanning mutagenesis approach was conducted to more completely understand the functional inter-relationships as well as any synergistic interactions that might occur within the loop. In general, additive effects for each side chain on the midpoint potentials for both couples were observed except for the hydroquinone state where some degree of nonadditivity was noted in multiple mutants involving Glu59. The study concluded with the generation of the triple mutant -Ala56-Gly-Ala-Ala59- in which all side chain interactions are removed. Gly57 was left unchanged because of its critical conformational contribution. Remarkably, this mutant retained the ability to bind the FMN and to thermodynamically stabilize the semiquinone state despite the absence of all side chain interactions. Collectively, these observations emphasize the overriding importance of the main chain interactions with the N5H of the FMN and the associated redox-dependent conformational change in this loop and leaves little doubt as to its role in the thermodynamic stabilization of the neutral semiquinone state of the FMN cofactor.     

Transition states for folding of circular-permuted proteins

To explore the role of entropy and chain connectivity in protein folding, a particularly interesting scheme, namely, the circular permutation, has been used. Recently, experimental observations showed that there are large differences in the folding mechanisms between the wild-type proteins and their circular permutants. These differences are strongly related to the change in the intrachain connectivity. Some results obtained by molecular dynamics simulations also showed a good agreement with the experimental findings. Here, we use a topology-based free-energy functional method to study the role of the chain connectivity in folding by comparing features of transition states of the wild-type proteins with those of their circular permutants. We concentrate our study on 3 small globular proteins, namely, the alpha-spectrin SH3 domain (SH3), the chymotrypsin inhibitor 2 (CI2), and the ribosomal protein S6, and obtain exciting results that are consistent with the available experimental and simulation results. A heterogeneity of the interaction energies between contacts for protein CI2 and for protein S6 is also introduced, which characterizes the strong interactions between contacts with long loops, as speculated from experiments for protein S6. The comparison between the folding nucleus of the wild-type proteins and those of their circular permutants indicates that chain connectivity affects remarkably the shapes of the energy profiles and thus the folding mechanism. Further comparisons between our theoretical calculated phi(th) values and the experimental observed phi(exp) values for the 3 proteins and their permutants show that our results are in good agreement with experimental ones and that correlations between them are high. These indicate that the free-energy functional method really provides a way to analyze the folding behavior of the circular-permuted proteins and therefore the folding mechanism of the wild-type proteins.     

Sequence determinants of a protein folding pathway

Local hydrophobic collapse of the polypeptide chain and transient long-range interactions in unfolded states of apomyoglobin appear to occur in regions of the amino acid sequence which, upon folding, bury an above-average area of hydrophobic surface. To explore the role of these interactions in protein folding, we prepared and characterized apomyoglobins with compensating point mutations designed to change the average buried surface area in local regions of the sequence, while conserving as much as possible the constitution of the hydrophobic core. The behavior of the mutants in quench-flow experiments to determine the folding pathway was exactly as predicted by the changes in the buried surface area parameter calculated from the amino acid sequence. In addition, spin label experiments with acid-unfolded mutant apomyoglobin showed that the transient long-range contacts that occur in the wild-type protein are abolished in the mutant, while new contacts are observed between areas that now have above-average buried surface area. We conclude that specific groupings of amino acid side-chains, which can be predicted from the sequence, are responsible for early hydrophobic interactions in the first phase of folding in apomyoglobin, and that these early interactions determine the subsequent course of the folding process.     

Interaction of cytochrome c with cytochrome c oxidase: an NMR study on two soluble fragments derived from Paracoccus denitrificans

The functional interactions between the various components of the respiratory chain are relatively short-lived, thus allowing high turnover numbers but at the same time complicating the structural analysis of the complexes. Chemical shift mapping by NMR spectroscopy is a useful tool to investigate such transient contacts, since it can monitor changes in the electron-shielding properties of a protein as the result of temporary contacts with a reaction partner. In this study, we investigated the molecular interaction between two components of the electron-transfer chain from Paracoccus denitrificans: the engineered, water-soluble fragment of cytochrome c(552) and the Cu(A) domain from the cytochrome c oxidase. Comparison of [(15)N,(1)H]-TROSY spectra of the [(15)N]-labeled cytochrome c(552) fragment in the absence and in the presence of the Cu(A) fragment showed chemical shift changes for the backbone amide groups of several, mostly uncharged residues located around the exposed heme edge in cytochrome c(552). The detected contact areas on the cytochrome c(552) surface were comparable under both fully reduced and fully oxidized conditions, suggesting that the respective chemical shift changes represent biologically relevant protein-protein interactions.     

15N NMR assignments and chemical shift analysis of uniformly labeled 15N calbindin D9k in the apo, (Cd2+)1 and (Ca2+)2 states.

15N has been uniformly incorporated into the EF-hand Ca(2+)-binding protein calbindin D9k so that heteronuclear experiments can be used to further characterize the structure and dynamics of the apo, (Cd2+)1 and (Ca2+)2 states of the protein. The 15N NMR resonances were assigned by 2D 15N-resolved 1H experiments, which also allowed the identification of a number of sequential and medium-range 1H-1H contacts that are obscured by chemical shift degeneracy in homonuclear experiments. The 15N chemical shifts are analyzed with respect to correlations with protein secondary structure. In addition, the changes in 15N chemical shift found for the apo----(Cd2+)1----(Ca2+)2 binding sequence confirm that the effects on the protein are mainly associated with chelation of the first ion.     

Partially formed native tertiary interactions in the A-state of cytochrome c.

Considerable insight into protein structure, stability, and folding has been obtained from studies of non-native states. We have studied the extent of native tertiary contacts in one such molecule, the A-state of yeast iso-1-ferricytochrome c. Previously, we showed that the interface between the N and C-terminal helices is completely formed in the A-state. Here, we focus on interactions essential for forming the heme pocket of eukaryotic cytochromes c. To determine the extent of these interactions, we used saturation mutagenesis at the evolutionarily invariant residue leucine 68, and measured the free energy of denaturation for the native states and the A-states of functional variants. We show that, unlike the interaction between the terminal helices, the native interactions between the 60s helix and the rest of the protein are not completely formed in the A-state.     

Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques

Pulse field gradient NMR methods have been used to determine the effective hydrodynamic radii of a range of native and nonnative protein conformations. From these experimental data, empirical relationships between the measured hydrodynamic radius (R(h)) and the number of residues in the polypeptide chain (N) have been established; for native folded proteins R(h) = 4.75N (0.29)A and for highly denatured states R(h) = 2.21N (0.57)A. Predictions from these equations agree well with experimental data from dynamic light scattering and small-angle X-ray or neutron scattering studies reported in the literature for proteins ranging in size from 58 to 760 amino acid residues. The predicted values of the hydrodynamic radii provide a framework that can be used to analyze the conformational properties of a range of nonnative states of proteins. Several examples are given here to illustrate this approach including data for partially structured molten globule states and for proteins that are unfolded but biologically active under physiological conditions. These reveal evidence for significant coupling between local and global features of the conformational ensembles adopted in such states. In particular, the effective dimensions of the polypeptide chain are found to depend significantly on the level of persistence of regions of secondary structure or features such as hydrophobic clusters within a conformational ensemble.     

Particle Polymorphism Caused by Deletion of a Peptide Molecular Switch in a Quasiequivalent Icosahedral Virus

The capsid of flock house virus is composed of 180 copies of a single type of coat protein which forms a T=3 icosahedral shell. High-resolution structural analysis has shown that the protein subunits, although chemically identical, form different contacts across the twofold axes of the virus particle. Subunits that are related by icosahedral twofold symmetry form flat contacts, whereas subunits that are related by quasi-twofold symmetry form bent contacts. The flat contacts are due to the presence of ordered genomic RNA and an ordered peptide arm which is inserted in the groove between the subunits and prevents them from forming the dihedral angle observed at the bent quasi-twofold contacts. We hypothesized that by deleting the residues that constitute the ordered peptide arm, formation of flat contacts should be impossible and therefore result in assembly of particles with only bent contacts. Such particles would have T=1 symmetry. To test this hypothesis we generated two deletion mutants in which either 50 or 31 residues were eliminated from the N terminus of the coat protein. We found that in the absence of residues 1 to 50, assembly was completely inhibited, presumably because the mutation removed a cluster of positively charged amino acids required for neutralization of encapsidated RNA. When the deletion was restricted to residues 1 to 31, assembly occurred, but the products were highly heterogeneous. Small bacilliform-like structures and irregular structures as well as wild-type-like T=3 particles were detected. The anticipated T=1 particles, on the other hand, were not observed. We conclude that residues 20 to 30 are not critical for formation of flat protein contacts and formation of T=3 particles. However, the N terminus of the coat protein appears to play an essential role in regulating assembly such that only one product, T=3 particles, is synthesized.     

Equilibrium unfolding of DLC8 monomer by urea and guanidine hydrochloride: Distinctive global and residue level features

We present circular dichroism (CD), steady state fluorescence and multidimensional NMR investigations on the equilibrium unfolding of monomeric dynein light chain protein (DLC8) by urea and guanidine hydrochloride (GdnHCl). Quantitative analysis of the CD and fluorescence denaturation curves reveals that urea unfolding is a two-state process, whereas guanidine unfolding is more complex. NMR investigations in the native state and in the near native states created by low denaturant concentrations enabled residue level characterization of the early structural and dynamic perturbations by the two denaturants. Firstly, (15)N transverse relaxation rates in the native state indicate that the regions around N10, Q27, the loop between beta2 and beta4 strands, and K87 at the C-terminal are potential unfolding initiation sites in the protein. Amide and (15)N chemical shift perturbations indicate different accessibilities of the residues along the chain and help identify locations of the early perturbations by the two denaturants. Guanidine and urea are seen to interact at several sites some of which are different in the two cases. Notable among the common interaction site is that around K87 which is in close proximity to W54 on the protein structure, but the interaction modes of the two denaturants are different. The secondary chemical shifts indicate that the structural perturbation by 1M urea is small, compared to that by guanidine which is more encompassing over the length of the chain. The probable (phi, psi) changes at the individual residues have been calculated using the TALOS algorithm. It appears that the helices in the protein are significantly perturbed by guanidine. Further, comparison of the spectral density functions of the native and the two near native states in the two denaturants implicate greater loosening of the structure by guanidine as compared to that by urea, even though the structures are still in the native state ensemble. These differences in the early perturbations of the native state structure and dynamics by the two denaturants might direct the protein along different pathways, as the unfolding progresses on further increasing the denaturant concentration.     

Enzymatic transition states: thermodynamics, dynamics and analogue design

Kinetic isotope effects and computational chemistry have defined the transition state structures for several members of the N-ribosyltransferase family. Transition state analogues designed to mimic their cognate transition state structures are among the most powerful enzyme inhibitors. In complexes of N-ribosyltransferases with their transition state analogues, the dynamic nature of the transition state is converted to an ordered, thermodynamic structure closely related to the transition state. This phenomenon is documented by peptide bond H/D exchange, crystallography and computational chemistry. Complexes with substrate, transition state and product analogues reveal reaction coordinate motion and catalytic interactions. Isotope-edited spectroscopic analysis and binding specificity of these complexes provides information about specific enzyme-transition state contacts. In combination with protein dynamic QM/MM models, it is proposed that the transition state is reached by stochastic dynamic excursions of the protein groups near the substrates in the closed conformation. Examples from fully dissociated (D(N) *A(N)), hybrid (D(N)A(N)) and symmetric nucleophilic displacement (A(N)D(N)) transition states are found in the N-ribosyltransferases. The success of transition state analogue inhibitor design based on kinetic isotope effects validates this approach to understanding enzymatic transition states.     

Van der Waals locks: loop-n-lock structure of globular proteins

In a globular protein the polypeptide chain returns to itself many times, making numerous chain-to-chain contacts. The stability of these contacts is maintained primarily by van der Waals interactions. In this work we isolated and analysed van der Waals contacts that stabilise spatial structures of nine major folds. We suggest a specific way to identify the tightest contacts of prime importance for the stability of a given crystallized protein and introduce the notion of the van der Waals lock. The loops closed by the van der Waals interactions provide a basically novel view of protein globule organization: the loop-n-lock structure. This opens a new perspective in understanding protein folding as well: the consecutive looping of the polypeptide chain and the locking of the loop ends by tight van der Waals interactions.     

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.     

Alternatively folded states of an immunoglobulin

Well-defined, non-native protein structures of low stability have been increasingly observed as intermediates in protein folding or as equilibrium structures populated under specific solvent conditions. These intermediate structures, frequently referred to as molten globule states, are characterized by the presence of secondary structure, a lack of significant tertiary contacts, increased hydrophobicity and partial specific volume as compared to native structures, and low cooperativity in thermal unfolding. The present study demonstrates that under acidic conditions (pH less than 3) the antibody MAK33 can assume a folded stable conformation. This A-state is characterized by a high degree of secondary structure, increased hydrophobicity, a native-like maximum wavelength of fluorescence emission, and a tendency toward slow aggregation. A prominent feature of this low-pH conformation is the stability against denaturant and thermal unfolding that is manifested in highly cooperative reversible phase transitions indicative of the existence of well-defined tertiary contacts. These thermodynamic results are corroborated by the kinetics of folding from the completely unfolded chain to the alternatively folded state at pH 2. The given data suggest that MAK33 at pH 2 adopts a cooperative structure that differs from the native immunoglobulin fold at pH 7. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its extraordinary structural stability that is characteristic for native protein structures.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

Atomic resolution structure of Moloney murine leukemia virus matrix protein and its relationship to other retroviral matrix proteins

Matrix proteins associated with the viral membrane are important in the formation of the viral particle and in virus maturation. The 1.0 A crystal structure of the ecotropic Gammaretrovirus Moloney murine leukemia virus (M-MuLV) matrix protein reveals the conserved topology of other retroviral matrix proteins, despite undetectable sequence similarity. The N terminus (normally myristylated) is exposed and adjacent to a basic surface patch, features likely to contribute to membrane binding. The four proteins in the asymmetric unit make varied contacts. The M-MuLV matrix structure is intermediate, between those of the lentiviruses and other retroviruses. The protein fold appears to be maintained, in part, by the conservation of side chain packing, which may provide a useful tool for searching for weak distant similarities in proteins.     

Progressive stabilization of intermediate and transition states in protein folding reactions by introducing surface hydrophobic residues

It can be argued from the principle of solvent exclusion that the introduction of hydrophobic residues onto the surface of a protein will not destabilize the folded state because the nonpolar side chain will be at least as exposed in the unfolded state as it is when the protein chain is folded. A comparison of the folding pathway of wild type and 11 site-directed mutants of CD2.d1 shows this to be true. In fact, owing to partial burial of nonpolar groups as folding proceeds, we find that the rapidly formed intermediate state and, to a greater extent, the transition state are generally stabilized by hydrophobic surface mutations. This effect is slightly moderated in the folded state presumably by the perturbation of van der Waals' contacts and/or local electrostatic interactions that have a greater influence in this fully compact structure. The fact that in all but one case we find that stabilization of the rapidly collapsed intermediate is accompanied by a faster acquisition of the folded state refutes the argument that I states are generally "off pathway" conformations or ensembles that lead to the inhibition of otherwise more rapid folding trajectories.     

On the optimal folding of protein molecules

The process of formation of a globular structure by a long molecular chain has been examined. In this process, various regions of the chain interact with one another. We classify the contacts thus formed as “correct” and “erroneous” ones. The correct contacts are those characteristic of the final native globular structure. All other contacts can be treated as erroneous. It is demonstrated that globule formation may proceed actually without formation and subsequent decay of erroneous contacts. Our model permits avoiding examination of numerous erroneous variants inasmuch as the regions of the chain that form correct contacts enter “long-range” interactions that at the same time can be highly selective. The existence of interactions of this kind facilitates the mutual approach and interaction of just those regions of the chain that yield correct contacts. Based on database analysis, it is shown that the model is valid not only for abstract structures but also for real polypeptide chains capable of forming protein globules and helical fibrils.     

The behaviour of polyamino acids reveals an inverse side chain effect in amyloid structure formation

Amyloid fibrils and prions are proteinaceous aggregates that are based on a unique form of polypeptide configuration, termed cross-beta structure. Using a group of chemically distinct polyamino acids, we show here that the existence of such a structure does not require the presence of specific side chain interactions or sequence patterns. These observations firmly establish that amyloid formation and protein folding represent two fundamentally different ways of organizing polypeptides into ordered conformations. Protein folding depends critically on the presence of distinctive side chain sequences and produces a unique globular fold. By contrast, the properties of different polyamino acids suggest that amyloid formation arises primarily from main chain interactions that are, in some environments, overruled by specific side chain contacts. This side chain effect can be thought of as the inverse of the one that characterizes protein folding. Conditions including Alzheimer's and Creutzfeldt-Jakob diseases represent, on this basis, pathological cases in which a natural polypeptide chain has aberrantly adopted the conformation that is primarily defined by main chain interactions and not the structure that is determined by specific side chain contacts that depend on the polypeptide sequence.