Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant “internal friction” arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

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.

Solid-state synthesis and mechanical unfolding of polymers of T4 lysozyme

Recent advances in single molecule manipulation methods offer a novel approach to investigating the protein folding problem. These studies usually are done on molecules that are naturally organized as linear arrays of globular domains. To extend these techniques to study proteins that normally exist as monomers, we have developed a method of synthesizing polymers of protein molecules in the solid state. By introducing cysteines at locations where bacteriophage T4 lysozyme molecules contact each other in a crystal and taking advantage of the alignment provided by the lattice, we have obtained polymers of defined polarity up to 25 molecules long that retain enzymatic activity. These polymers then were manipulated mechanically by using a modified scanning force microscope to characterize the force-induced reversible unfolding of the individual lysozyme molecules. This approach should be general and adaptable to many other proteins with known crystal structures. For T4 lysozyme, the force required to unfold the monomers was 64 ± 16 pN at the pulling speed used. Refolding occurred within 1 sec of relaxation with an efficiency close to 100%. Analysis of the force versus extension curves suggests that the mechanical unfolding transition follows a two-state model. The unfolding forces determined in 1 M guanidine hydrochloride indicate that in these conditions the activation barrier for unfolding is reduced by 2 kcal/mol.

A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

The "allosteric three-site model" of elongation cannot be confirmed in a well-defined ribosome system from Escherichia coli.

For the functional role of the ribosomal tRNA exit (E) site, two different models have been proposed. It has been suggested that transient E-site binding of the tRNA leaving the peptidyl (P) site promotes elongation factor G (EF-G)-dependent translocation by lowering the energetic barrier of tRNA release [Lill, R., Robertson, J. M. & Wintermeyer, W. (1989) EMBO J. 8, 3933-3938]. The alternative "allosteric three-site model" [Nierhaus, K.H. (1990) Biochemistry 29, 4997-5008] features stable, codon-dependent tRNA binding to the E site and postulates a coupling between E and aminoacyl (A) sites that regulates the tRNA binding affinity of the two sites in an anticooperative manner. Extending our testing of the two conflicting models, we have performed translocation experiments with fully active ribosomes programmed with heteropolymeric mRNA. The results confirm that the deacylated tRNA released from the P site is bound to the E site in a kinetically labile fashion, and that the affinity of binding, i.e., the occupancy of the E site, is increased by Mg2+ or polyamines. At conditions of high E-site occupancy in the posttranslocation complex, filling the A site with aminoacyl-tRNA had no influence on the E site, i.e., there was no detectable anticooperative coupling between the two sites, provided that second-round translocation was avoided by removing EF-G. On the basis of these results, which are entirely consistent with our previous results, we consider the allosteric three-site model of elongation untenable. Rather, as proposed earlier, the E site-bound state of the leaving tRNA is a transient intermediate and, as such, is a mechanistic feature of the classic two-state model of the elongating ribosome.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

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.

NMR investigation of ferricytochrome c unfolding: Detection of an equilibrium unfolding intermediate and residual structure in the denatured state

Horse ferricytochrome c (cyt c) undergoes exchange of one of its axial heme ligands (Met-80) for one or more non-native ligands under denaturing conditions. We have used 1H NMR spectroscopy to detect two conformations of paramagnetic cyt c with non-native heme ligation through a range of urea concentrations. One non-native form is an equilibrium unfolding intermediate observed under partially denaturing conditions and is attributed to replacement of Met-80 with one or more Lys side chains. The second non-native form, in which the native Met ligand is replaced by a His, is observed under strongly denaturing conditions. Thermodynamic analysis of these data indicates a relatively small ΔG (17 kJ/mol) for the transition from native to the Lys-ligated intermediate and a significantly larger ΔG (47 kJ/mol) for the transition from native to the His-ligated species. Although CD and fluorescence data indicate that the equilibrium unfolding of cyt c is a two-state process, these NMR results implicate an intermediate with His-Lys ligation.

A model for amplification of hair-bundle motion by cyclical binding of Ca2+ to mechanoelectrical-transduction channels

Amplification of auditory stimuli by hair cells augments the sensitivity of the vertebrate inner ear. Cell-body contractions of outer hair cells are thought to mediate amplification in the mammalian cochlea. In vertebrates that lack these cells, and perhaps in mammals as well, active movements of hair bundles may underlie amplification. We have evaluated a mathematical model in which amplification stems from the activity of mechanoelectrical-transduction channels. The intracellular binding of Ca2+ to channels is posited to promote their closure, which increases the tension in gating springs and exerts a negative force on the hair bundle. By enhancing bundle motion, this force partially compensates for viscous damping by cochlear fluids. Linear stability analysis of a six-state kinetic model reveals Hopf bifurcations for parameter values in the physiological range. These bifurcations signal conditions under which the system’s behavior changes from a damped oscillatory response to spontaneous limit-cycle oscillation. By varying the number of stereocilia in a bundle and the rate constant for Ca2+ binding, we calculate bifurcation frequencies spanning the observed range of auditory sensitivity for a representative receptor organ, the chicken’s cochlea. Simulations using prebifurcation parameter values demonstrate frequency-selective amplification with a striking compressive nonlinearity. Because transduction channels occur universally in hair cells, this active-channel model describes a mechanism of auditory amplification potentially applicable across species and hair-cell types.

Dissection of the ion-induced folding of the hammerhead ribozyme using 19F NMR

We have used 19F NMR to analyze the metal ion-induced folding of the hammerhead ribozyme by selective incorporation of 5fluorouridine. We have studied the chemical shift and linewidths of 19F resonances of 5-fluorouridine at the 4 and 7 positions in the ribozyme core as a function of added Mg2+. The data fit well to a simple two-state model whereby the formation of domain 1 is induced by the noncooperative binding of Mg2+ with an association constant in the range of 100 to 500 M−1, depending on the concentration of monovalent ions present. The results are in excellent agreement with data reporting on changes in the global shape of the ribozyme. However, the NMR experiments exploit reporters located in the center of the RNA sections undergoing the folding transitions, thereby allowing the assignment of specific nucleotides to the separate stages. The results define the folding pathway at high resolution and provide a time scale for the first transition in the millisecond range.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Thermal unfolding of human high-density apolipoprotein A-1: implications for a lipid-free molten globular state.

Apolipoprotein A-1 (apoA-1) in complex with high-density lipoprotein is critically involved in the transport and metabolism of cholesterol and in the pathogenesis of atherosclerosis. We reexamined the thermal unfolding of lipid-free apoA-1 in low-salt solution at pH approximately 7, by using differential scanning calorimetry and circular dichroism. At protein concentrations <5 mg/ml, thermal unfolding of apoA-1 is resolved as an extended peak (25 degrees C-90 degrees C) that can be largely accounted for by a single reversible non-two-state transition with midpoint Tm 57 +/- 1 degree C, calorimetric enthalpy deltaH(Tm)= 200 +/- 20 kcal/mol (1 kcal = 4.18 kJ), van't Hoff enthalpy deltaHv(Tm) approximately 32.5 kcal/mol, and cooperativity deltaHv(Tm)/deltaH(Tm) approximately 0.16. The enthalpy deltaH(Tm) can be accounted for by melting of the alpha-helical structure that is inferred by CD to constitute approximately 60% of apoA-1 amino acids. Farand near-UV CD spectra reveal noncoincident melting of the secondary and tertiary structural elements and indicate a well-defined secondary structure but a largely melted tertiary structure for apoA-1 at approximately 37 degrees C and pH 7. This suggests a molten globular-like state for lipid-free apoA-1 under near-physiological conditions. Our results suggest that in vivo lipid binding by apoA-1 may be mediated via the molten globular apolipoprotein state in plasma.

Simple model of protein folding kinetics.

A simple model of the kinetics of protein folding is presented. The reaction coordinate is the "correctness" of a configuration compared with the native state. The model has a gap in the energy spectrum, a large configurational entropy, a free energy barrier between folded and partially folded states, and a good thermodynamic folding transition. Folding kinetics is described by a master equation. The folding time is estimated by means of a local thermodynamic equilibrium assumption and then is calculated both numerically and analytically by solving the master equation. The folding time has a maximum near the folding transition temperature and can have a minimum at a lower temperature.

A CDKN2-like polymorphism in Xiphophorus LG V is associated with UV-B-induced melanoma formation in platyfish–swordtail hybrids

The genetic basis of spontaneous melanoma formation in spotted dorsal (Sd) Xiphophorus platyfish–swordtail hybrids has been studied for decades, and is adequately explained by a two-gene inheritance model involving a sex-linked oncogene, Xmrk, and an autosomal tumor suppressor, DIFF. The Xmrk oncogene encodes a receptor tyrosine kinase related to EGFR; the nature of the DIFF tumor suppressor gene is unknown. We analyzed the genetic basis of UV-B-induced melanoma formation in closely related, spotted side platyfish–swordtail hybrids, which carry a different sex-linked pigment pattern locus, Sp. We UV-irradiated spotted side Xiphophorus platyfish–swordtail backcross hybrids to induce melanomas at frequencies 6-fold higher than occur spontaneously in unirradiated control animals. To identify genetic determinants of melanoma susceptibility in this UV-inducible Xiphophorus model, we genotyped individual animals from control and UV-irradiated experimental regimes using allozyme and DNA restriction fragment length polymorphisms and tested for joint segregation of genetic markers with pigmentation phenotype and UV-induced melanoma formation. Joint segregation results show linkage of a CDKN2-like DNA polymorphism with UV-B-induced melanoma formation in these hybrids. The CDKN2-like polymorphism maps to Xiphophorus linkage group V and exhibits recombination fractions with ES1 and MDH2 allozyme markers consistent with previous localization of the DIFF tumor suppressor locus. Our results indicate that the CDKN2-like sequence we have cloned and mapped is a candidate for the DIFF tumor suppressor gene.