Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Denaturant-induced movement of the transition state of protein folding revealed by high-pressure stopped-flow measurements

The small all-β protein tendamistat folds and unfolds with two-state kinetics. We determined the volume changes associated with the folding process by performing kinetic and equilibrium measurements at variable pressure between 0.1 and 100 MPa (1 to 1,000 bar). GdmCl-induced equilibrium unfolding transitions reveal that the volume of the native state is increased by 41.4 ± 2.0 cm3/mol relative to the unfolded state. This value is virtually independent of denaturant concentration. The use of a high-pressure stopped-flow instrument enabled us to measure the activation volumes for the refolding (ΔVf0‡) and unfolding reaction (ΔVu0‡) over a broad range of GdmCl concentrations. The volume of the transition state is 60% native-like (ΔVf0‡ = 25.0 ± 1.2 cm3/mol) in the absence of denaturant, indicating partial solvent accessibility of the core residues. The volume of the transition state increases linearly with denaturant concentration and exceeds the volume of the native state above 6 M GdmCl. This result argues for a largely desolvated transition state with packing deficiencies at high denaturant concentrations and shows that the structure of the transition state depends strongly on the experimental conditions.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

Guanidine hydrochloride blocks a critical step in the propagation of the prion-like determinant [PSI+] of Saccharomyces cerevisiae

The cytoplasmic heritable determinant [PSI+] of the yeast Saccharomyces cerevisiae reflects the prion-like properties of the chromosome-encoded protein Sup35p. This protein is known to be an essential eukaryote polypeptide release factor, namely eRF3. In a [PSI+] background, the prion conformer of Sup35p forms large oligomers, which results in the intracellular depletion of functional release factor and hence inefficient translation termination. We have investigated the process by which the [PSI+] determinant can be efficiently eliminated from strains, by growth in the presence of the protein denaturant guanidine hydrochloride (GuHCl). Strains are “cured” of [PSI+] by millimolar concentrations of GuHCl, well below that normally required for protein denaturation. Here we provide evidence indicating that the elimination of the [PSI+] determinant is not derived from the direct dissolution of self-replicating [PSI+] seeds by GuHCl. Although GuHCl does elicit a moderate stress response, the elimination of [PSI+] is not enhanced by stress, and furthermore, exhibits an absolute requirement for continued cell division. We propose that GuHCl inhibits a critical event in the propagation of the prion conformer and demonstrate that the kinetics of curing by GuHCl fit a random segregation model whereby the heritable [PSI+] element is diluted from a culture, after the total inhibition of prion replication by GuHCl.

Role of metallothionein in nitric oxide signaling as revealed by a green fluorescent fusion protein

Although the function of metallothionein (MT), a 6- to 7-kDa cysteine-rich metal binding protein, remains unclear, it has been suggested from in vitro studies that MT is an important component of intracellular redox signaling, including being a target for nitric oxide (NO). To directly study the interaction between MT and NO in live cells, we generated a fusion protein consisting of MT sandwiched between two mutant green fluorescent proteins (GFPs). In vitro studies with this chimera (FRET-MT) demonstrate that fluorescent resonance energy transfer (FRET) can be used to follow conformational changes indicative of metal release from MT. Imaging experiments with live endothelial cells show that agents that increase cytoplasmic Ca2+ act via endogenously generated NO to rapidly and persistently release metal from MT. A role for this interaction in intact tissue is supported by the finding that the myogenic reflex of mesenteric arteries is absent in MT knockout mice (MT−/−) unless endogenous NO synthesis is blocked. These results are the first application of intramolecular green fluorescent protein (GFP)-based FRET in a native protein and demonstrate the utility of FRET-MT as an intracellular surrogate indicator of NO production. In addition, an important role of metal thiolate clusters of MT in NO signaling in vascular tissue is revealed.

Specificity of native-like interhelical hydrophobic contacts in the apomyoglobin intermediate

On exposure to mildly acidic conditions, apomyoglobin forms a partially folded intermediate, I. The A, B, G, and H helices are significantly structured in this equilibrium intermediate, whereas the remainder of the protein is largely unfolded. We report here the effects of mutations at helix pairing sites on the stability of I in three classes of mutants that: (i) truncate hydrophobic side chains in native helix packing sites, (ii) truncate hydrophobic side chains not involved in interhelical contacts, and (iii) extend hydrophobic side chains at residues not involved in interhelical contacts. Class I mutants significantly decrease the stability and cooperativity of folding of the intermediate. Class II and III mutants show smaller effects on stability and have little effect on cooperativity. Qualitatively similar results to those found in I were obtained for all three classes of mutants in native myoglobin (N), demonstrating that hydrophobic burial is fairly specific to native helix packing sites in I as well as in N. These results suggest that hydrophobic burial along native-like interhelical contacts is important for the formation of the cooperatively folded intermediate.

A specific transition state for S-peptide combining with folded S-protein and then refolding

We measured the folding and unfolding kinetics of mutants for a simple protein folding reaction to characterize the structure of the transition state. Fluorescently labeled S-peptide analogues combine with S-protein to form ribonuclease S analogues: initially, S-peptide is disordered whereas S-protein is folded. The fluorescent probe provides a convenient spectroscopic probe for the reaction. The association rate constant, kon, and the dissociation rate constant, koff, were both determined for two sets of mutants. The dissociation rate constant is measured by adding an excess of unlabeled S-peptide analogue to a labeled complex (RNaseS*). This strategy allows kon and koff to be measured under identical conditions so that microscopic reversibility applies and the transition state is the same for unfolding and refolding. The first set of mutants tests the role of the α-helix in the transition state. Solvent-exposed residues Ala-6 and Gln-11 in the α-helix of native RNaseS were replaced by the helix destabilizing residues glycine or proline. A plot of log kon vs. log Kd for this series of mutants is linear over a very wide range, with a slope of −0.3, indicating that almost all of the molecules fold via a transition state involving the helix. A second set of mutants tests the role of side chains in the transition state. Three side chains were investigated: Phe-8, His-12, and Met-13, which are known to be important for binding S-peptide to S-protein and which also contribute strongly to the stability of RNaseS*. Only the side chain of Phe-8 contributes significantly, however, to the stability of the transition state. The results provide a remarkably clear description of a folding transition state.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Mapping the parameters of prion-induced neuropathology

We present a theoretical framework that enables us to dissect out the parametric dependencies of the pathogenesis of prion diseases. We are able to determine the influence of both host-dependent factors (connectivity, cell density, protein synthesis rate, and cell death) and strain-dependent factors (cell tropism, virulence, and replication rate). We use a model based on a linked system of differential equations on a lattice to explore how the regional distribution of central nervous system pathology in Creutzfeldt-Jakob disease, Gerstmann-Sträussler-Scheinker syndrome, and fatal familial insomnia relates to each of these factors. The model then is used to make qualitative predictions about the pathology for two possible hypothetical triggers of neuronal loss in prion diseases. Pathological progression in overexpressing mouse models has been shown to depend on the site of initial infection. The model allows us to compare the pathologies resulting from different inoculation routes.

Statistical significance of protein structure prediction by threading

In this study, we estimate the statistical significance of structure prediction by threading. We introduce a single parameter ɛ that serves as a universal measure determining the probability that the best alignment is indeed a native-like analog. Parameter ɛ takes into account both length and composition of the query sequence and the number of decoys in threading simulation. It can be computed directly from the query sequence and potential of interactions, eliminating the need for sequence reshuffling and realignment. Although our theoretical analysis is general, here we compare its predictions with the results of gapless threading. Finally we estimate the number of decoys from which the native structure can be found by existing potentials of interactions. We discuss how this analysis can be extended to determine the optimal gap penalties for any sequence-structure alignment (threading) method, thus optimizing it to maximum possible performance.

Novel folded protein domains generated by combinatorial shuffling of polypeptide segments

It has been proposed that the architecture of protein domains has evolved by the combinatorial assembly and/or exchange of smaller polypeptide segments. To investigate this proposal, we fused DNA encoding the N-terminal half of a β-barrel domain (from cold shock protein CspA) with fragmented genomic Escherichia coli DNA and cloned the repertoire of chimeric polypeptides for display on filamentous bacteriophage. Phage displaying folded polypeptides were selected by proteolysis; in most cases the protease-resistant chimeric polypeptides comprised genomic segments in their natural reading frames. Although the genomic segments appeared to have no sequence homologies with CspA, one of the originating proteins had the same fold as CspA, but another had a different fold. Four of the chimeric proteins were expressed as soluble polypeptides; they formed monomers and exhibited cooperative unfolding. Indeed, one of the chimeric proteins contained a set of very slowly exchanging amides and proved more stable than CspA itself. These results indicate that native-like proteins can be generated directly by combinatorial segment assembly from nonhomologous proteins, with implications for theories of the evolution of new protein folds, as well as providing a means of creating novel domains and architectures in vitro.

The biochemical properties of the ATPase activity of a 70-kDa heat shock protein (Hsp70) are governed by the C-terminal domains

The cytosolic 70-kDa heat shock proteins (Hsp70s), Ssa and Ssb, of Saccharomyces cerevisiae are functionally distinct. Here we report that the ATPase activities of these two classes of Hsp70s exhibit different kinetic properties. The Ssa ATPase has properties similar to those of other Hsp70s studied, such as DnaK and Hsc70. Ssb, however, has an unusually low steady-state affinity for ATP but a higher maximal velocity. In addition, the ATPase activity of Hsp70s, like that of Ssa1, depends on the addition of K+ whereas Ssb activity does not. Suprisingly, the isolated 44-kDa ATPase domain of Ssb has a Km and Vmax for ATP hydrolysis similar to those of Ssa, rather than those of full length Ssb. Analysis of Ssa/Ssb fusion proteins demonstrates that the Ssb peptide-binding domain fused to the Ssa ATPase domain generates an ATPase of relatively high activity and low steady-state affinity for ATP similar to that of native Ssb. Therefore, at least some of the biochemical differences between the ATPases of these two classes of Hsp70s are not intrinsic to the ATPase domain itself. The differential influence of the peptide-binding domain on the ATPase domain may, in part, explain the functional uniqueness of these two classes of Hsp70s.

The barrier-to-autointegration protein is a host factor for HIV type 1 integration

In vivo, retroviral integration is mediated by a large nucleoprotein complex, termed the preintegration complex (PIC). PICs isolated from infected cells display in vitro integration activity. Here, we analyze the roles of different host cell factors in the structure and function of HIV type 1 (HIV-1) PICs. PICs purified by size exclusion after treatment with high salt lost their integration activity, and adding back an extract from uninfected cells restored this activity. In parallel, the native protein–DNA intasome structure detected at the ends of HIV-1 by Mu-mediated PCR footprinting was abolished by high salt and restored by the crude cell extract. Various purified proteins previously implicated in retroviral PIC function then were analyzed for their effects on the structure and function of salt-treated HIV-1 PICs. Whereas relatively low amounts (5–20 nM) of human barrier-to-autointegration factor (BAF) protein restored integration activity, substantially more (5–10 μM) human host factor HMG I(Y) was required. Similarly high levels (3–8 μM) of bovine RNase A, a DNA-binding protein used as a nonspecific control, also restored activity. Mu-mediated PCR footprinting revealed that of these three purified proteins, only BAF restored the native structure of the HIV-1 protein–DNA intasome. We suggest that BAF is a natural host cofactor for HIV-1 integration.

Production, specificity, and functionality of monoclonal antibodies to specific peptide–major histocompatibility complex class II complexes formed by processing of exogenous protein

Several unanswered questions in T cell immunobiology relating to intracellular processing or in vivo antigen presentation could be approached if convenient, specific, and sensitive reagents were available for detecting the peptide–major histocompatibility complex (MHC) class I or class II ligands recognized by αβ T cell receptors. For this reason, we have developed a method using homogeneously loaded peptide–MHC class II complexes to generate and select specific mAb reactive with these structures using hen egg lysozyme (HEL) and I-Ak as a model system. mAbs specific for either HEL-(46–61)–Ak or HEL-(116–129)–Ak have been isolated. They cross-react with a small subset of I-Ak molecules loaded with self peptides but can nonetheless be used for flow cytometry, immunoprecipitation, Western blotting, and intracellular immunofluorescence to detect specific HEL peptide–MHC class II complexes formed by either peptide exposure or natural processing of native HEL. An example of the utility of these reagents is provided herein by using one of the anti-HEL-(46–61)–Ak specific mAbs to visualize intracellular compartments where I-Ak is loaded with HEL-derived peptides early after antigen administration. Other uses, especially for in vivo tracking of specific ligand-bearing antigen-presenting cells, are discussed.

Identification of human DNA helicase V with the far upstream element-binding protein

The properties of human DNA helicase V (HDH V) were studied in greater detail following an improved purification procedure. From 450 g of cultured cells, <0.1 mg of pure protein was isolated. HDH V unwinds DNA unidirectionally by moving in the 3′ to 5′ direction along the bound strand in an ATP- and Mg2+-dependent fashion. The enzyme is not processive and can also unwind partial RNA–RNA duplexes such as HDH IV and HDH VIII. The Mr determined by SDS–PAGE (66 kDa) corresponds to that measured under native conditions, suggesting that HDH V exists as a monomer in the nucleus. Microsequencing of the purified HDH V shows that this enzyme is identical to the far upstream element-binding protein (FBP), a protein that stimulates the activity of the c-myc gene by binding specifically to the ‘FUSE’ DNA region localized upstream of its promoter. The sequence of HDH V/FBP contains RGG motifs like HDH IV/nucleolin, HDH VIII/G3BP as well as other human RNA and DNA helicases identified by other laboratories.