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.

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.

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.

The pKa of His-24 in the folding transition state of apomyoglobin

In native apomyoglobin, His-24 cannot be protonated, although at pH 4 the native protein forms a molten globule folding intermediate in which the histidine residues are readily protonated. The inability to protonate His-24 in the native protein dramatically affects the unfolding/refolding kinetics, as demonstrated by simulations for a simple model. Kinetic data for wild type and for a mutant lacking His-24 are analyzed. The pKa values of histidine residues in native apomyoglobin are known from earlier studies, and the average histidine pKa in the molten globule is determined from the pH dependence of the equilibrium between the native and molten globule forms. Analysis of the pH-dependent unfolding/refolding kinetics reveals that the average pKa of the histidine residues, including His-24, is closely similar in the folding transition state to the value found in the molten globule intermediate. Consequently, protonation of His-24 is not a barrier to refolding of the molten globule to the native protein. Instead, the normal pKa of His-24 in the transition state, coupled with its inaccessibility in the native state, promotes fast unfolding at low pH. The analysis of the wild-type results is confirmed and extended by using the wild-type parameters to fit the unfolding kinetics of a mutant lacking His-24.

Simple mechanochemistry describes the dynamics of kinesin molecules

Recently, Block and coworkers [Visscher, K., Schnitzer, M. J., & Block, S. M. (1999) Nature (London) 400, 184–189 and Schnitzer, M. J., Visscher, K. & Block, S. M. (2000) Nat. Cell Biol. 2, 718–723] have reported extensive observations of individual kinesin molecules moving along microtubules in vitro under controlled loads, F = 1 to 8 pN, with [ATP] = 1 μM to 2 mM. Their measurements of velocity, V, randomness, r, stalling force, and mean run length, L, reveal a need for improved theoretical understanding. We show, presenting explicit formulae that provide a quantitative basis for comparing distinct molecular motors, that their data are satisfactorily described by simple, discrete-state, sequential stochastic models. The simplest (N = 2)-state model with fixed load-distribution factors and kinetic rate constants concordant with stopped-flow experiments, accounts for the global (V, F, L, [ATP]) interdependence and, further, matches relative acceleration observed under assisting loads. The randomness, r(F,[ATP]), is accounted for by a waiting-time distribution, ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{1}^{+}}}\end{equation*}\end{document}(t), [for the transition(s) following ATP binding] with a width parameter ν ≡ 〈t〉2/〈(Δt)2〉≃2.5, indicative of a dispersive stroke of mechanicity ≃0.6 or of a few (≳ν − 1) further, kinetically coupled states: indeed, N = 4 (but not N = 3) models do well. The analysis reveals: (i) a substep of d0 = 1.8–2.1 nm on ATP binding (consistent with structurally based suggestions); (ii) comparable load dependence for ATP binding and unbinding; (iii) a strong load dependence for reverse hydrolysis and subsequent reverse rates; and (iv) a large (≳50-fold) increase in detachment rate, with a marked load dependence, following ATP binding.

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.