Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

Chiral molecular self-assembly of phospholipid tubules: A circular dichroism study

We report on spectroscopic studies of the chiral structure in phospholipid tubules formed in mixtures of alcohol and water. Synthetic phospholipids containing diacetylenic moieties in the acyl chains self-assemble into hollow, cylindrical tubules in appropriate conditions. Circular dichroism provides a direct measure of chirality of the molecular structure. We find that the CD spectra of tubules formed in mixtures of alcohol and water depends strongly on the alcohol used and the lipid concentration. The relative spectral intensity of different circular dichroism bands correlates with the number of bilayers observed using microscopy. The results provide experimental evidence that tubule formation is based on chiral packing of the lipid molecules and that interbilayer interactions are important to the tubule structure.

Site-specific conformational determination in thermal unfolding studies of helical peptides using vibrational circular dichroism with isotopic substitution

Understanding the detailed mechanism of protein folding requires dynamic, site-specific stereochemical information. The short time response of vibrational spectroscopies allows evaluation of the distribution of populations in rapid equilibrium as the peptide unfolds. Spectral shifts associated with isotopic labels along with local stereochemical sensitivity of vibrational circular dichroism (VCD) allow determination of the segment sequence of unfolding. For a series of alanine-rich peptides that form α-helices in aqueous solution, we used isotopic labeling and VCD to demonstrate that the α-helix noncooperatively unwinds from the ends with increasing temperature. For these blocked peptides, the C-terminal is frayed at 5°C. Ab initio level theoretical simulations of the IR and VCD band shapes are used to analyze the spectra and to confirm the conformation of the labeled components. The VCD signals associated with the labeled residues are amplified by coupling to the nonlabeled parts of the molecule. Thus small labeled segments are detectable and stereochemically defined in moderately large peptides in this report of site-specific peptide VCD conformational analysis.

Spectroscopic recognition of guanine dimeric hairpin quadruplexes by a carbocyanine dye.

Isolated guanine quadruplex structures have been described at high resolution both in solution and in the solid state. The existence of this unusual DNA structure in vivo and its biological significance remain to be determined. We describe the binding of 3,3'-diethyloxadicarbocyanine to dimeric hairpin guanine quadruplexes. This interaction results in a set of unique spectrophotometric signatures, none of which arises from binding to single strands or Watson-Crick duplexes. These unique signatures include a new absorbance peak (lambda max = 534 nm), an induced circular dichroism (lambda = 534-626 nm), a quenching of the dye fluorescence upon excitation with visible light, and strong energy transfer from DNA. This last effect provides the basis for detecting hairpin quadruplex structures in the presence of excess amounts of nonquadruplex DNA structures, such as single strands and Watson-Crick duplexes. The mechanism of quadruplex recognition by this dye is discussed, along with the possibility of using this dye as a probe for hairpin quadruplex structures in vitro and in vivo.

Highly conserved features of DNA binding between two divergent members of the Myb family of transcription factors

Bas1p, a divergent yeast member of the Myb family of transcription factors, shares with the proteins of this family a highly conserved cysteine residue proposed to play a role in redox regulation. Substitutions of this residue in Bas1p (C153) allowed us to establish that, despite its very high conservation, it is not strictly required for Bas1p function: its substitution with a small hydrophobic residue led to a fully functional protein in vitro and in vivo. C153 was accessible to an alkylating agent in the free protein but was protected by prior exposure to DNA. The reactivity of cysteines in the first and third repeats was much lower than in the second repeat, suggesting a more accessible conformation of repeat 2. Proteolysis protection, fluorescence quenching and circular dichroism experiments further indicated that DNA binding induces structural changes making Bas1p less accessible to modifying agents. Altogether, our results strongly suggest that the second repeat of the DNA-binding domain of Bas1p behaves similarly to its Myb counterpart, i.e. a DNA-induced conformational change in the second repeat leads to formation of a full helix–turn–helix-related motif with the cysteine packed in the hydrophobic core of the repeat.

Alpha-helical, but not beta-sheet, propensity of proline is determined by peptide environment.

Proline is established as a potent breaker of both alpha-helical and beta-sheet structures in soluble (globular) proteins. Thus, the frequent occurrence of the Pro residue in the putative transmembrane helices of integral membrane proteins, particularly transport proteins, presents a structural dilemma. We propose that this phenomenon results from the fact that the structural propensity of a given amino acid may be altered to conform to changes imposed by molecular environment. To test this hypothesis on proline, we synthesized model peptides of generic sequence H2N-(Ser-LyS)2-Ala- Leu-Z-Ala-Leu-Z-Trp-Ala-Leu-Z-(Lys-Ser)3-OH (Z = Ala and/or Pro). Peptide conformations were analyzed by circular dichroism spectroscopy in aqueous buffer, SDS, lysophosphatidylglycerol micelles, and organic solvents (methanol, trifluoroethanol, and 2-propanol). The helical propensity of Pro was found to be greatly enhanced in the membrane-mimetic environments of both lipid micelles and organic solvents. Proline was found to stabilize the alpha-helical conformation relative to Ala at elevated temperatures in 2-propanol, an observation that argues against the doctrine that Pro is the most potent alpha-helix breaker as established in aqueous media. Parallel studies in deoxycholate micelles of the temperature-induced conformational transitions of the single-spanning membrane bacteriophage IKe major coat protein, in which the Pro-containing wild type was compared with Pro30 --> Ala mutant, Pro was found to protect the helix, but disrupt the beta-sheet structure as effectively as it does to model peptides in water. The intrinsic capacity of Pro to disrupt beta-sheets was further reflected in a survey of porins where Pro was found to be selectively excluded from the core of membrane-spanning beta-sheet barrels. The overall data provide a rationale for predicting and understanding the structural consequences when Pro occurs in the context of a membrane.

Irreversible inactivation of interleukin 2 in a pump-based delivery environment.

The physical stability of pharmaceutical proteins in delivery environments is a critical determinant of biological potency and treatment efficacy, and yet it is often taken for granted. We studied both the bioactivity and physical stability of interleukin 2 upon delivery via continuous infusion. We found that the biological activity of the delivered protein was dramatically reduced by approximately 90% after a 24-hr infusion program. Only a portion of these losses could be attributed to direct protein deposition on the delivery surfaces. Analysis of delivered protein by size exclusion chromatography gave no indication of insulin-like, surface-induced aggregation phenomena. Examination of the secondary and tertiary structure of both adsorbed and delivered protein via Fourier-transform infrared spectroscopy, circular dichroism, and fluorescence spectroscopy indicated that transient surface association of interleukin 2 with the catheter tubing resulted in profound, irreversible structural changes that were responsible for the majority of the biological activity losses.

Stopped-flow NMR spectroscopy: real-time unfolding studies of 6-19F-tryptophan-labeled Escherichia coli dihydrofolate reductase.

Escherichia coli dihydrofolate reductase (DHFR; EC 1.5.1.3) contains five tryptophan residues that have been replaced with 6-19F-tryptophan. The 19F NMR assignments are known in the native, unliganded form and the unfolded form. We have used these assignments with stopped-flow 19F NMR spectroscopy to investigate the behavior of specific regions of the protein in real time during urea-induced unfolding. The NMR data show that within 1.5 sec most of the intensities of the native 19F resonances of the protein are lost but only a fraction (approximately 20%) of the intensities of the unfolded resonances appears. We postulate that the early disappearance of the native resonances indicates that most of the protein rapidly forms an intermediate in which the side chains have considerable mobility. Stopped-flow far-UV circular dichroism measurements indicate that this intermediate retains native-like secondary structure. Eighty percent of the intensities of the NMR resonances assigned to the individual tryptophans in the unfolded state appear with similar rate constants (k approximately 0.14 sec-1), consistent with the major phase of unfolding observed by stopped-flow circular dichroism (representing 80% of total amplitude). These data imply that after formation of the intermediate, which appears to represent an expanded structural form, all regions of the protein unfold at the same rate. Stopped-flow measurements of the fluorescence and circular dichroism changes associated with the urea-induced unfolding show a fast phase (half-time of about 1 sec) representing 20% of the total amplitude in addition to the slow phase mentioned above. The NMR data show that approximately 20% of the total intensity for each of the unfolded tryptophan resonances is present at 1.5 sec, indicating that these two phases may represent the complete unfolding of the two different populations of the native protein.

Structure and stability of a second molten globule intermediate in the apomyoglobin folding pathway.

Apomyoglobin folding proceeds through a molten globule intermediate (low-salt form; I1) that has been characterized by equilibrium (pH 4) and kinetic (pH 6) folding experiments. Of the eight alpha-helices in myoglobin, three (A, G, and H) are structured in I1, while the rest appear to be unfolded. Here we report on the structure and stability of a second intermediate, the trichloroacetate form of the molten globule intermediate (I2), which is induced either from the acid-unfolded protein or from I1 by > or = 5 mM sodium trichloroacetate. Circular dichroism measurements monitoring urea- and acid-induced unfolding indicate that I2 is more highly structured and more stable than I1. Although I2 exhibits properties closer to those of the native protein, one-dimensional NMR spectra show that it maintains the lack of fixed side-chain structure that is the hallmark of a molten globule. Amide proton exchange and 1H-15N two-dimensional NMR experiments are used to identify the source of the extra helicity observed in I2. The results reveal that the existing A, G, and H helices present in I1 have become more stable in I2 and that a fourth helix--the B helix--has been incorporated into the molten globule. Available evidence is consistent with I2 being an on-pathway intermediate. The data support the view that apomyoglobin folds in a sequential fashion through a single pathway populated by intermediates of increasing structure and stability.

Amyloid fibril formation by an SH3 domain

The SH3 domain is a well characterized small protein module with a simple fold found in many proteins. At acid pH, the SH3 domain (PI3-SH3) of the p85α subunit of bovine phosphatidylinositol 3-kinase slowly forms a gel that consists of typical amyloid fibrils as assessed by electron microscopy, a Congo red binding assay, and x-ray fiber diffraction. The soluble form of PI3-SH3 at acid pH (the A state by a variety of techniques) from which fibrils are generated has been characterized. Circular dichroism in the far- and near-UV regions and 1H NMR indicate that the A state is substantially unfolded relative to the native protein at neutral pH. NMR diffusion measurements indicate, however, that the effective hydrodynamic radius of the A state is only 23% higher than that of the native protein and is 20% lower than that of the protein denatured in 3.5 M guanidinium chloride. In addition, the A state binds the hydrophobic dye 1-anilinonaphthalene-8-sulfonic acid, which suggests that SH3 in this state has a partially formed hydrophobic core. These results indicate that the A state is partially folded and support the hypothesis that partially folded states formed in solution are precursors of amyloid deposition. Moreover, that this domain aggregates into amyloid fibrils suggests that the potential for amyloid deposition may be a common property of proteins, and not only of a few proteins associated with disease.

Structural basis of an embryonically lethal single Ala → Thr mutation in the vnd/NK-2 homeodomain

The structural and DNA binding behavior is described for an analog of the vnd/NK-2 homeodomain, which contains a single amino acid residue alanine to threonine replacement in position 35 of the homeodomain. Multidimensional nuclear magnetic resonance, circular dichroism, and electrophoretic gel retardation assays were carried out on recombinant 80-aa residue proteins that encompass the wild-type and mutant homeodomains. The mutant A35T vnd/NK-2 homeodomain is unable to adopt a folded conformation free in solution at temperatures down to −5°C in contrast to the behavior of the corresponding wild-type vnd/NK-2 homeodomain, which is folded into a functional three-dimensional structure below 25°C. The A35T vnd/NK-2 binds specifically to the vnd/NK-2 target DNA sequence, but with an affinity that is 50-fold lower than that of the wild-type homeodomain. Although the three-dimensional structure of the mutant A35T vnd/NK-2 in the DNA bound state shows characteristic helix–turn–helix behavior similar to that of the wild-type homeodomain, a notable structural deviation in the mutant A35T analog is observed for the amide proton of leucine-40. The wild-type homeodomain forms an unusual i,i-5 hydrogen bond with the backbone amide oxygen of residue 35. In the A35T mutant this amide proton resonance is shifted upfield by 1.27 ppm relative to the resonance frequency for the wild-type analog, thereby indicating a significant alteration of this i,i-5 hydrogen bond.

Ultrafast signals in protein folding and the polypeptide contracted state

To test the significance of ultrafast protein folding signals (≪1 msec), we studied cytochrome c (Cyt c) and two Cyt c fragments with major C-terminal segments deleted. The fragments remain unfolded under all conditions and so could be used to define the unfolded baselines for protein fluorescence and circular dichroism (CD) as a function of denaturant concentration. When diluted from high to low denaturant in kinetic folding experiments, the fragments readjust to their new baseline values in a “burst phase” within the mixing dead time. The fragment burst phase reflects a contraction of the polypeptide from a more extended unfolded condition at high denaturant to a more contracted unfolded condition in the poorer, low denaturant solvent. Holo Cyt c exhibits fluorescence and CD burst phase signals that are essentially identical to the fragment signals over the whole range of final denaturant concentrations, evidently reflecting the same solvent-dependent, relatively nonspecific contraction and not the formation of a specific folding intermediate. The significance of fast folding signals in Cyt c and other proteins is discussed in relation to the hypothesis of an initial rate-limiting search-nucleation-collapse step in protein folding [Sosnick, T. R., Mayne, L. & Englander, S. W. (1996) Proteins Struct. Funct. Genet. 24, 413–426].

In vitro evolution of horse heart myoglobin to increase peroxidase activity

Random mutagenesis and screening for enzymatic activity has been used to engineer horse heart myoglobin to enhance its intrinsic peroxidase activity. A chemically synthesized gene encoding horse heart myoglobin was subjected to successive cycles of PCR random mutagenesis. The mutated myoglobin gene was expressed in Escherichia coli LE392, and the variants were screened for peroxidase activity with a plate assay. Four cycles of mutagenesis and screening produced a series of single, double, triple, and quadruple variants with enhanced peroxidase activity. Steady-state kinetics analysis demonstrated that the quadruple variant T39I/K45D/F46L/I107F exhibits peroxidase activity significantly greater than that of the wild-type protein with k1 (for H2O2 oxidation of metmyoglobin) of 1.34 × 104 M−1 s−1 (≈25-fold that of wild-type myoglobin) and k3 [for reducing the substrate (2, 2′-azino-di-(3-ethyl)benzthiazoline-6-sulfonic acid] of 1.4 × 106 M−1 s−1 (1.6-fold that of wild-type myoglobin). Thermal stability of these variants as measured with circular dichroism spectroscopy demonstrated that the Tm of the quadruple variant is decreased only slightly compared with wild-type (74.1°C vs. 76.5°C). The rate constants for binding of dioxygen exhibited by the quadruple variant are identical to the those observed for wild-type myoglobin (kon, 22.2 × 10−6 M−1 s−1 vs. 22.3 × 10−6 M−1 s−1; koff, 24.3 s−1 vs. 24.2 s−1; KO2, 0.91 × 10−6 M−1 vs. 0.92 × 10−6 M−1). The affinity of the quadruple variant for CO is increased slightly (kon, 0.90 × 10−6 M−1s−1 vs. 0.51 × 10−6 M−1s−1; koff, 5.08 s−1 vs. 3.51 s−1; KCO, 1.77 × 10−7 M−1 vs. 1.45 × 10−7 M−1). All four substitutions are in the heme pocket and within 5 Å of the heme group.

Conversion of a c type cytochrome to a b type that spontaneously forms in vitro from apo protein and heme: Implications for c type cytochrome biogenesis and folding

Cytochrome c552 from Hydrogenobacter thermophilus, a thermophilic bacterium, has been converted into a b type cytochrome, after mutagenesis of both heme-binding cysteines to alanine and expression in the cytoplasm of Escherichia coli. The b type variant is less stable, with the guanidine hydrochloride unfolding midpoint occurring at a concentration 2 M lower than for the wild-type protein. The reduction potential is 75 mV lower than that of the recombinant wild-type protein. The heme can be removed from the b type variant, thus generating an apo protein that has, according to circular dichroism spectroscopy, an α-helical content different from that of the holo b type protein. The latter is readily reformed in vitro by addition of heme to the apo protein. This reforming suggests that previously observed assembly of cytochrome c552, which has the typical class I cytochrome c fold, in the E. coli cytoplasm is a consequence of spontaneous thioether bond formation after binding of heme to a prefolded polypeptide. These observations have implications for the general problem of c type cytochrome biogenesis.

Self-assembly of large, ordered lamellae from non-bilayer lipids and integral membrane proteins in vitro

In many biological membranes, the major lipids are “non-bilayer lipids,” which in purified form cannot be arranged in a lamellar structure. The structural and functional roles of these lipids are poorly understood. This work demonstrates that the in vitro association of the two main components of a membrane, the non-bilayer lipid monogalactosyldiacylglycerol (MGDG) and the chlorophyll-a/b light-harvesting antenna protein of photosystem II (LHCII) of pea thylakoids, leads to the formation of large, ordered lamellar structures: (i) thin-section electron microscopy and circular dichroism spectroscopy reveal that the addition of MGDG induces the transformation of isolated, disordered macroaggregates of LHCII into stacked lamellar aggregates with a long-range chiral order of the complexes; (ii) small-angle x-ray scattering discloses that LHCII perturbs the structure of the pure lipid and destroys the inverted hexagonal phase; and (iii) an analysis of electron micrographs of negatively stained 2D crystals indicates that in MGDG-LHCII the complexes are found in an ordered macroarray. It is proposed that, by limiting the space available for MGDG in the macroaggregate, LHCII inhibits formation of the inverted hexagonal phase of lipids; in thylakoids, a spatial limitation is likely to be imposed by the high concentration of membrane-associated proteins.