From snapshot to movie: φ analysis of protein folding transition states taken one step further

Kinetic anomalies in protein folding can result from changes of the kinetic ground states (D, I, and N), changes of the protein folding transition state, or both. The 102-residue protein U1A has a symmetrically curved chevron plot which seems to result mainly from changes of the transition state. At low concentrations of denaturant the transition state occurs early in the folding reaction, whereas at high denaturant concentration it moves close to the native structure. In this study we use this movement to follow continuously the formation and growth of U1A's folding nucleus by φ analysis. Although U1A's transition state structure is generally delocalized and displays a typical nucleation–condensation pattern, we can still resolve a sequence of folding events. However, these events are sufficiently coupled to start almost simultaneously throughout the transition state structure.

Chemically distinct transition states govern rapid dissociation of single L-selectin bonds under force

Carbohydrate–protein bonds interrupt the rapid flow of leukocytes in the circulation by initiation of rolling and tethering at vessel walls. The cell surface carbohydrate ligands are glycosylated proteins like the mucin P-selectin glycoprotein ligand-1 (PSGL-1), which bind ubiquitously to the family of E-, P-, and L-selectin proteins in membranes of leukocytes and endothelium. The current view is that carbohydrate–selectin bonds dissociate a few times per second, and the unbinding rate increases weakly with force. However, such studies have provided little insight into how numerous hydrogen bonds, a Ca2+ metal ion bond, and other interactions contribute to the mechanical strength of these attachments. Decorating a force probe with very dilute ligands and controlling touch to achieve rare single-bond events, we have varied the unbinding rates of carbohydrate–selectin bonds by detachment with ramps of force/time from 10 to 100,000 pN/sec. Testing PSGL-1, its outer 19 aa (19FT), and sialyl LewisX (sLeX) against L-selectin in vitro on glass microspheres and in situ on neutrophils, we found that the unbinding rates followed the same dependence on force and increased by nearly 1,000-fold as rupture forces rose from a few to ≈200 pN. Plotted on a logarithmic scale of loading rate, the rupture forces reveal two prominent energy barriers along the unbinding pathway. Strengths above 75 pN arise from rapid detachment (<0.01 sec) impeded by an inner barrier that requires a Ca2+ bond between a single sLeX and the lectin domain. Strengths below 75 pN occur under slow detachment (>0.01 sec) impeded by the outer barrier, which appears to involve an array of weak (putatively hydrogen) bonds.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

Non-enzymatic and enzymatic hydrolysis of alkyl halides: A haloalkane dehalogenation enzyme evolved to stabilize the gas-phase transition state of an SN2 displacement reaction

The semiempirical PM3 method, calibrated against ab initio HF/6–31+G(d) theory, has been used to elucidate the reaction of 1,2-dichloroethane (DCE) with the carboxylate of Asp-124 at the active site of haloalkane dehalogenase of Xanthobacter autothropicus. Asp-124 and 13 other amino acid side chains that make up the active site cavity (Glu-56, Trp-125, Phe-128, Phe-172, Trp-175, Leu-179, Val-219, Phe-222, Pro-223, Val-226, Leu-262, Leu-263, and His-289) were included in the calculations. The three most significant observations of the present study are that: (i) the DCE substrate and Asp-124 carboxylate, in the reactive ES complex, are present as an ion-molecule complex with a structure similar to that seen in the gas-phase reaction of AcO− with DCE; (ii) the structures of the transition states in the gas-phase and enzymatic reaction are much the same where the structure formed at the active site is somewhat exploded; and (iii) the enthalpies in going from ground states to transition states in the enzymatic and gas-phase reactions differ by only a couple kcal/mol. The dehalogenase derives its catalytic power from: (i) bringing the electrophile and nucleophile together in a low-dielectric environment in an orientation that allows the reaction to occur without much structural reorganization; (ii) desolvation; and (iii) stabilizing the leaving chloride anion by Trp-125 and Trp-175 through hydrogen bonding.

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.

Visualization of intermediate and transition-state structures in protein-tyrosine phosphatase catalysis.

Engineering site-specific amino acid substitutions into the protein-tyrosine phosphatase (PTPase) PTP1 and the dual-specific vaccinia H1-related phosphatase (VHR), has kinetically isolated the two chemical steps of the reaction and provided a rare opportunity for examining transition states and directly observing the phosphoenzyme intermediate. Changing serine to alanine in the active-site sequence motif HCXXGXXRS shifted the rate-limiting step from intermediate formation to intermediate hydrolysis. Using phosphorus 31P NMR, the covalent thiol-phosphate intermediate was directly observed during catalytic turnover. The importance of the conserved aspartic acid (D92 in VHR and D181 in PTP1) in both chemical steps was established. Kinetic analysis of D92N and D181N mutants indicated that aspartic acid acts as a general acid by protonating the leaving-group phenolic oxygen. Structure-reactivity experiments with native and aspartate mutant enzymes established that proton transfer is concomitant with P-O cleavage, such that no charge develops on the phenolic oxygen. Steady- and presteady-state kinetics, as well as NMR analysis of the double mutant D92N/S131A (VHR), suggested that the conserved aspartic acid functions as a general base during intermediate hydrolysis. As a general base, aspartate would activate a water molecule to facilitate nucleophilic attack. The amino acids involved in transition-state stabilization for cysteinylphosphate hydrolysis were confirmed by the x-ray structure of the Yersinia PTPase complexed with vanadate, a transition-state mimic that binds covalently to the active-site cysteine. Consistent with the NMR, x-ray, biochemical, and kinetic data, a unifying mechanism for catalysis is proposed.

Design of potent and selective human cathepsin K inhibitors that span the active site

Potent and selective active-site-spanning inhibitors have been designed for cathepsin K, a cysteine protease unique to osteoclasts. They act by mechanisms that involve tight binding intermediates, potentially on a hydrolytic pathway. X-ray crystallographic, MS, NMR spectroscopic, and kinetic studies of the mechanisms of inhibition indicate that different intermediates or transition states are being represented that are dependent on the conditions of measurement and the specific groups flanking the carbonyl in the inhibitor. The species observed crystallographically are most consistent with tetrahedral intermediates that may be close approximations of those that occur during substrate hydrolysis. Initial kinetic studies suggest the possibility of irreversible and reversible active-site modification. Representative inhibitors have demonstrated antiresorptive activity both in vitro and in vivo and therefore are promising leads for therapeutic agents for the treatment of osteoporosis. Expansion of these inhibitor concepts can be envisioned for the many other cysteine proteases implicated for therapeutic intervention.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.

Protein folding from a highly disordered denatured state: The folding pathway of chymotrypsin inhibitor 2 at atomic resolution

Previous experimental and theoretical studies have produced high-resolution descriptions of the native and folding transition states of chymotrypsin inhibitor 2 (CI2). In similar fashion, here we use a combination of NMR experiments and molecular dynamics simulations to examine the conformations populated by CI2 in the denatured state. The denatured state is highly unfolded, but there is some residual native helical structure along with hydrophobic clustering in the center of the chain. The lack of persistent nonnative structure in the denatured state reduces barriers that must be overcome, leading to fast folding through a nucleation–condensation mechanism. With the characterization of the denatured state, we have now completed our description of the folding/unfolding pathway of CI2 at atomic resolution.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

Extension of the thermodynamics of small systems to open metastable states: An example

By using a simplified model of small open liquid-like clusters with surface effects, in the gas phase, it is shown how the statistical thermodynamics of small systems can be extended to include metastable supersaturated gaseous states not too far from the gas–liquid equilibrium transition point. To accomplish this, one has to distinguish between mathematical divergence and physical convergence of the open-system partition function.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].

A signature of the T → R transition in human hemoglobin

Allosteric effects in hemoglobin arise from the equilibrium between at least two energetic states of the molecule: a tense state, T, and a relaxed state, R. The two states differ from each other in the number and energy of the interactions between hemoglobin subunits. In the T state, constraints between subunits oppose the structural changes resulting from ligand binding. In the R state, these constraints are released, thus enhancing ligand-binding affinity. In the present work, we report the presence of four sites in hemoglobin that are structurally stabilized in the R relative to the T state. These sites are Hisα103(G10) and Hisα122(H5) in each α subunit of hemoglobin. They are located at the α1β1 and α2β2 interfaces of the hemoglobin tetramer, where the histidine side chains form hydrogen bonds with specific residues from the β chains. We have measured the solvent exchange rates of side chain protons of Hisα103(G10) and Hisα122(H5) in both deoxygenated and ligated hemoglobin by NMR spectroscopy. The exchange rates were found to be higher in the deoxygenated-T than in ligated-R state. Analysis of exchange rates in terms of the local unfolding model revealed that the structural stabilization free energy at each of these two histidines is larger by ≈1.5 kcal/(mol tetramer) in the R relative to the T state. The location of these histidines at the intradimeric α1β1 and α2β2 interfaces also suggests a role for these interfaces in the allosteric equilibrium of hemoglobin.