Nanosecond folding dynamics of a three-stranded beta-sheet

We studied conformational stability and folding kinetics of a three-stranded beta-sheet containing two rigid turns. Static infrared measurements indicate that this beta-sheet undergoes a broad but cooperative thermal unfolding transition with a midpoint at approximately 53 degrees C. Interestingly, time-resolved infrared experiments show that its relaxation kinetics in response to a temperature-jump (T-jump) occur on the nanosecond time scale (e.g., the relaxation time is approximately 140 ns at 35.0 degrees C), thereby suggesting that the conformational relaxation encounters only a small free energy barrier or even proceeds in a downhill manner. Further Langevin dynamics simulations suggest that the observed T-jump relaxation kinetics could be modeled by a conformational diffusion process along a single-well free energy profile, which allowed us to determine the effective diffusion constant and also the roughness of the folding energy landscape.     

Theoretical characterization of alpha-helix and beta-hairpin folding kinetics

By means of the conformational free energy surface and corresponding diffusion coefficients, as obtained by long time scale atomistic molecular dynamics simulations (mus time scale), we model the folding kinetics of alpha-helix and beta-hairpin peptides as a diffusive process over the free energy surface. The two model systems studied in this paper (the alpha-helical temporin L and the beta-hairpin prion protein H1 peptide) exhibit a funnel-like almost barrierless free energy profile, leading to nonexponential folding kinetics matching rather well the available experimental data. Moreover, using the free energy profile provided by Mu?oz et al. [Mu?oz et al. Nature 1997, 390: 196-199], this model was also applied to reproduce the two-state folding kinetics of the C-terminal beta-hairpin of protein GB1, yielding an exponential folding kinetics with a time constant (approximately 5 micros) in excellent agreement with the experimentally observed one (approximately 6 micros). Finally, the folding kinetics obtained by solving the diffusion equation, considering either a one-dimensional or a two-dimensional free energy surface, are also compared in order to understand the relevance of the possible kinetic coupling between conformational degrees of freedom in the folding process.     

Folding kinetics of a naturally occurring helical peptide: implication of the folding speed limit of helical proteins

The folding mechanism and dynamics of a helical protein may strongly depend on how quickly its constituent alpha-helices can fold independently. Thus, our understanding of the protein folding problem may be greatly enhanced by a systematic survey of the folding rates of individual alpha-helical segments derived from their parent proteins. As a first step, we have studied the relaxation kinetics of the central helix (L9:41-74) of the ribosomal protein L9 from the bacterium Bacillus stearothermophilus , in response to a temperature-jump ( T-jump) using infrared spectroscopy. L9:41-74 has been shown to exhibit unusually high helicity in aqueous solution due to a series of side chain-side chain interactions, most of which are electrostatic in nature, while still remaining monomeric over a wide concentration range. Thus, this peptide represents an excellent model system not only for examining how the folding rate of naturally occurring helices differs from that of the widely studied alanine-based peptides, but also for estimating the folding speed limit of (small) helical proteins. Our results show that the T-jump induced relaxation rate of L9:41-74 is significantly slower than that of alanine-based peptides. For example, at 11 degrees C its relaxation time constant is about 2 micros, roughly seven times slower than that of SPE(5), an alanine-rich peptide of similar chain length. In addition, our results show that the folding rate of a truncated version of L9:41-74 is even slower. Taken together, these results suggest that individual alpha-helical segments in proteins may fold on a time scale that is significantly slower than the folding time of alanine-based peptides. Furthermore, we argue that the relaxation rate of L9:41-74 measured between 8 and 45 degrees C provides a realistic estimate of the ultimate folding rate of (small) helical proteins over this temperature range.     

Helix–coil kinetics of two 14-residue peptides

The helix–coil transition kinetics of two 14-residue helical peptides of different stability were studied by time-resolved infrared (IR) spectroscopy coupled with laser-induced temperature-jump (T-jump) technique. The T-jump induced relaxation kinetics of both peptides show strong dependence on the final temperature, implying the existence of an enthalpic barrier for the nucleation process. In addition, the peptide with end-capping groups, which is more stable, folds faster. Together, these results suggest that the overall helical stability plays an important role in controlling the kinetics of the helix–coil transition, in agreement with results of early theoretical studies.     

UV Raman spatially resolved melting dynamics of isotopically labeled polyalanyl peptide: slow alpha-helix melting follows 3(10)-helices and pi-bulges premelting

We used UV resonance Raman (UVRR) to examine the spatial dependence of the T-jump secondary structure relaxation of an isotopically labeled 21-residue mainly Ala peptide, AdP. The AdP penultimate Ala residues were perdeuterated, leaving the central residues hydrogenated, to allow separate monitoring of melting of the middle versus the end peptide bonds. For 5 to 30 degrees C T-jumps, the central peptide bonds show a approximately 2-fold slower relaxation time (189 +/- 31 ns) than do the exterior peptide bonds (97 +/- 15 ns). In contrast, for a 20 to 40 degrees C T-jump, the central peptide bond relaxation appears to be faster (56 +/- 6 ns) than that of the penultimate peptide bonds (131 +/- 46 ns). We show that, if the data are modeled as a two-state transition, we find that only exterior peptide bonds show anti-Arrhenius folding behavior; the middle peptide bonds show both normal Arrhenius-like folding and unfolding. This anti-Arrhenius behavior results from the involvement of pi-bulges/helices and 3(10)-helix states in the melting. The unusual temperature dependence of the (un)folding rates of the interior and exterior peptide bonds is due to the different relative (un)folding rates of 3(10)-helices, alpha-helices, and pi-bulges/helices. Pure alpha-helix unfolding rates are approximately 12-fold slower (approximately 1 micros) than that of pi-bulges and 3(10)-helices. In addition, we also find that the alpha-helix is most stable at the AdP N-terminus where eight consecutive Ala occur, whereas the three hydrophilic Arg located in the middle and at the C-terminus destabilize the alpha-helix in these regions and induce defects such as pi-bulges and 3(10)-helices.     

Protein folding kinetics: barrier effects in chemical and thermal denaturation experiments

Recent experimental work on fast protein folding brings about an intriguing paradox. Microsecond-folding proteins are supposed to fold near or at the folding speed limit (downhill folding), but yet their folding behavior seems to comply with classical two-state analyses, which imply the crossing of high free energy barriers. However, close inspection of chemical and thermal denaturation kinetic experiments in fast-folding proteins reveals systematic deviations from two-state behavior. Using a simple one-dimensional free energy surface approach we find that such deviations are indeed diagnostic of marginal folding barriers. Furthermore, the quantitative analysis of available fast-kinetic data indicates that many microsecond-folding proteins fold downhill in native conditions. All of these proteins are then promising candidates for an atom-by-atom analysis of protein folding using nuclear magnetic resonance.1 We also find that the diffusion coefficient for protein folding is strongly temperature dependent, corresponding to an activation energy of approximately 1 kJ.mol-1 per protein residue. As a consequence, the folding speed limit at room temperature is about an order of magnitude slower than the approximately 1 micros estimates from high-temperature T-jump experiments. Our analysis is quantitatively consistent with the available thermodynamic and kinetic data on slow two-state folding proteins and provides a straightforward explanation for the apparent fast-folding paradox.     

A designed beta-hairpin peptide for molecular recognition of ATP in water

A designed 12-residue beta-hairpin peptide with a diagonal tryptophan (Trp) pair was shown to bind ATP in water through a combination of aromatic and electrostatic interactions. The affinity for ATP was 5800 M-1 (DeltaG approximately -5.0 kcal/mol), a remarkable affinity for a short, structured peptide in water, consisting of entirely natural amino acid residues. Proton NMR measurements indicate that the adenine ring of the nucleotide is intercalated between the diagonal tryptophans in the bound state. Delineation of the contributions to ATP binding to the hairpin suggest that aromatic interactions contribute approximately -1.8 kcal/mol, while individual electrostatic interactions involving the ATP phosphates and positively charged side chains of the hairpin contribute approximately -1 kcal/mol each. The designed beta-hairpin receptor presents a novel minimalist system to investigate the energetic contributions to protein-nucleic acid recognition through the surface of a beta-sheet.     

Beta-hairpin folding by a model amyloid peptide in solution and at an interface

The development of specific agents against amyloidoses requires an understanding of the conformational distribution of fibrillogenic peptides at a microscopic level. Here, I present molecular dynamics simulations of the model amyloid peptide LSFD with sequence LSFDNSGAITIG-NH2 in explicit water and at a water/vapor interface for a total time scale of approximately 1.8 micros. An extended structure was used as initial peptide configuration. At approximately 290 K, solvated LSFD was kinetically trapped in diverse misfolded beta-sheet/coil conformations. At 350 K, in contrast, the same type II' beta-hairpin in equilibrium with less ordered but also U-shaped conformations was observed for the core residues DNSGAITI in solution and at the interface in multiple independent simulations. The most stable structural unit of the beta-hairpin was the two residue turn (GA). The core residues exhibited a well-defined folded state in which the beta-hairpin was stabilized by a hydrogen bond between the side chain of Asn-385 and the main chain carbonyl group of Gly-387. My results suggest that beta-sheet conformations indicated from previous Fourier-transform infrared spectroscopy measurements immediately after preparation of the peptide solution may not arise from protofilaments as speculated by others but are a property of LSFD monomers. In addition, combined with previous results from X-ray scattering, my findings suggest that interfacial aggregation of LSFD implies a transition from U-shaped to extended peptide conformations. This work including the first simulations of reversible beta-hairpin folding at an interface is an essential step toward a microscopic understanding of interfacial peptide folding and self-assembly. Knowledge of the main conformation of the peptide core may facilitate the design of possible inhibitors of LSFD aggregation as a test ground for future computational therapeutic strategies against amyloid diseases.     

Thermodynamic analysis of beta-hairpin-forming peptides from the thermal dependence of (1)H NMR chemical shifts

The temperature dependence of the (1)H chemical shifts of six designed peptides previously shown to adopt beta-hairpin structures in aqueous solution has been analyzed in terms of two-state (beta-hairpin left arrow over right arrow coil) equilibrium. The stability of the beta-hairpins formed by these peptides, as derived from their T(m) (midpoint transition temperature) values, parallels in general their ability to adopt those structures as deduced from independent NMR parameters: NOEs, Deltadelta(C)(alpha)(H), Deltadelta(C)(alpha), and Deltadelta(C)(beta) values. The observed T(m) values are dependent on the particular position within the beta-hairpin that is probed, indicating that their folding to a beta-hairpin conformation deviates from a "true" two-state transition. To obtain individual T(m) values for each hairpin region in each peptide, a simplified model of a successive uncoupled two-state equilibrium covering the entire process has been applied. The distribution of T(m) values obtained for the different beta-hairpin regions (turn, strands, backbone, side chains) in the six analyzed peptides reveals a similar pattern. A model for beta-hairpin folding is proposed on the basis of this pattern and the reasonable assumption that regions showing higher T(m) values are the last ones to unfold and, presumably, the first to form. With this assumption, the analysis suggests that turn formation is the first event in beta-hairpin folding. This is consistent with previous results on the essential role of the turn sequence in beta-hairpin folding.     

Single-molecule electrophoresis of beta-hairpin peptides by electrical recordings and Langevin dynamics simulations

We used single-channel electrical recordings and Langevin molecular dynamics simulations to explore the electrophoretic translocation of various beta-hairpin peptides across the staphylococcal alpha-hemolysin (alphaHL) protein pore at single-molecule resolution. The beta-hairpin peptides, which varied in their folding properties, corresponded to the C terminal residues of the B1 domain of protein G. The translocation time was strongly dependent on the electric force and was correlated with the folding features of the beta-hairpin peptides. Highly unfolded peptides entered the pore in an extended conformation, resulting in fast single-file translocation events. In contrast, the translocation of the folded beta-hairpin peptides occurred more slowly. In this case, the beta-hairpin peptides traversed the alphaHL pore in a misfolded or fully folded conformation. This study demonstrates that the interaction between a polypeptide and a beta-barrel protein pore is dependent on the folding features of the polypeptide.     

Characterization of the stereochemical selectivity of beta-hairpin formation by molecular dynamics simulations

The stability of secondary structure motifs found in proteins is influenced by the choice of the configuration of the chiral centers present in the amino acid residues (i.e., D vs L). Experimental studies showed that the structural properties of the tetrapeptide (L)V(L)P(L)A(L)L (all-L) are drastically altered upon mutating the L-proline and the L-alanine by their d-enantiomers [J. Am. Chem. Soc. 1996, 118, 6975]. The all-L diastereomer is unstructured, experiencing little or no beta-hairpin formation, while the (L)V(D)P(D)A(L)L peptide exhibits a substantial population of beta-hairpin conformation. In this study, we perform molecular dynamics simulations to investigate the folding propensity of these two model peptides. The results confirm the experimental findings, namely, that the presence of d-amino acids in the loop region strongly induces beta-hairpin formation (a population increase from about 1.5% to 50% is observed). The major factor determining the different behavior is found to be the large difference in energy between the two diastereomers, approximately 22 kJ/mol, when they adopt a beta-hairpin structure. The higher energy observed for the all-L peptide is a consequence of none-ideal hydrogen bond formation and of steric repulsions. The results suggest that selective incorporation of D-amino acids in proteins can be used to enhance certain secondary structure elements. The kinetic behavior of the folding process observed in the simulations is also investigated. We find that the decay rate of the folded structure fits to a biexponential function, suggesting that the folding/unfolding process of a beta-hairpin is governed by two different mechanisms.     

Characterizing the kinetics of volume recovery in glasses by instantaneous temperature-jump experiments

This article proposes a temperature-jump (T-jump) approach for characterizing the kinetics of volume recovery in glassy materials. The kinetic characterization is based on the Kovacs-Aklonis model. This incorporates a retardation-time spectrum which shifts according to both the temperature and the instantaneous volume. The proposed experiments involve measuring the change in recovery rate caused by an abrupt temperature jump. Although an analogous procedure has been used to determine the activation energy for linear viscoelastic creep, the analysis for volume recovery is complicated by its inherent nonlinearity. Nevertheless, accounting for the nonlinearity by a reduction of the time scale permits the T-jump results to be analyzed. In particular, the T-jump approach can be used to: (i) test a particular functional form for the shift factor and (ii) determine the previously unmeasurable parameter x, which defines the relative importance of the temperature dependence and the volume dependence in this function. In addition, numerical simulations indicate that the proposed method can be implemented in the laboratory.     

Experimental resolution of early steps in protein folding: testing molecular dynamics simulations

Time-resolved Tyr fluorescence spectroscopy coupled with a laser-induced temperature-jump (T-jump) was employed to follow the folding relaxation dynamics of the B-domain of Staphylococcal protein A. The single Tyr is located in helix 1 (H1) and is a sensitive probe of the structure of this helix and the overall helical bundle structure. The results from this study were compared to those from a complementary infrared T-jump study on this protein [Vu, D. M.; Myers, J. K.; Oas, T. G.; Dyer, R. B. Biochemistry 2004, 43, 3582]. Both methods detect a microsecond process that follows the cooperative relaxation of the helical bundle core. However, a fast process (10-7 s) that follows the relaxation of the individual helices was observed only with the infrared probe. Thus, fast formation of H1 is not observed, but rather H1 forms in the microsecond phase, concomitantly with the docking to (and stabilization by) the other two helices to form the helical bundle structure. This observation validates the results of several previous molecular dynamics simulations that predict H1 formation only in the final assembly of the helix bundle.     

Sequential merocyanine product isomerization following femtosecond UV excitation of a spiropyran

The ring-opening dynamics of the photochromic switch 1',3'-dihydro-1',3',3'-trimethyl-6-nitrospiro[2H-1-benzopyran-2,2'-(2H)-indole] in tetrachloroethene is studied with both femtosecond time-resolved ultraviolet (UV)/visible and UV/mid-infrared (IR) pump-probe spectroscopy. During the first picosecond we identify two new transient features in the UV/vis experiments, the first of which we assign to spiropyran S1 --> S(n) absorption (lifetime < or = 0.2 ps). The second feature (lifetime 0.5 +/- 0.2 ps) we tentatively assign to the merocyanine T2 state. After 1 ps both probing methods show biexponential merocyanine formation kinetics, with average time constants of 17 +/- 3 and 350 +/- 20 ps. In the UV/IR experiments, the initial dynamics show more dispersion in formation times than in the UV/vis measurements, whereas the slower time constant is the same in both. A weak transient IR signal at approximately 1360 cm(-1) demonstrates that this biexponentiality is caused by a sequential isomerization between two merocyanine species. Lifetimes provide evidence that the merocyanine S1 state is not involved in the photochemical reaction.     

Calorimetric evidence for a two-state unfolding of the beta-hairpin peptide trpzip4

beta-Sheets are a common secondary structural element found in proteins. The difficulty in studying beta-sheet folding and stability is that their formation is often dependent on the tertiary structural environment within the protein. However, the discovery of water-soluble beta-hairpin peptides has allowed them to be used as model systems because they represent the smallest units of beta-sheet structure independent of tertiary structural context. Trpzip4 has been used as a model beta-hairpin peptide to study beta-hairpin folding and stability because it is highly soluble in aqueous solutions, maintains its monomeric state, and shows reversible cooperative thermal unfolding. The previously determined thermodynamic parameters for trpzip4 thermal unfolding vary depending on the spectroscopic probe used, which questions the assumption that trpzip4 unfolds in a two-state manner. Here we provide direct calorimetric evidence that the unfolding of trpzip4 follows a two-state unfolding mode. Furthermore, the thermal unfolding of trpzip4 monitored using near- and far-UV-CD yielded thermodynamic parameters similar to those determined calorimetrically, providing additional evidence for a two-state unfolding mode.     

Interstrand side chain--side chain interactions in a designed beta-hairpin: significance of both lateral and diagonal pairings

The contributions of interstrand side chain-side chain contacts to beta-sheet stability have been examined with an autonomously folding beta-hairpin model system. RYVEV(D)PGOKILQ-NH2 ((D)P = D-proline, O = ornithine) has previously been shown to adopt a beta-hairpin conformation in aqueous solution, with a two-residue loop at D-Pro-Gly. In the present study, side chains that display interstrand NOEs (Tyr-2, Lys-9, and Leu-11) are mutated to alanine or serine, and the conformational impact of the mutations is assessed. In the beta-hairpin conformation Tyr-2 and Leu-11 are directly across from one another (non-hydrogen bonded pair). This "lateral" juxtaposition of two hydrophobic side chains appears to contribute to beta-hairpin conformational stability, which is consistent with results from other beta-sheet model studies and with statistical analyses of interstrand residue contacts in protein crystal structures. Interaction between the side chains of Tyr-2 and Lys-9 also stabilizes the beta-hairpin conformation. Tyr-2/Lys-9 is a "diagonal" interstrand juxtaposition because these residues are not directly across from one another in terms of the hydrogen bonding registry between the strands. This diagonal interaction arises from the right-handed twist that is commonly observed among beta-sheets. Evidence of diagonal side chain-side chain contacts has been observed in other autonomously folding beta-sheet model systems, but we are not aware of other efforts to determine whether a diagonal interaction contributes to beta-sheet stability.     

Folding mechanisms of individual beta-hairpins in a Go model of Pin1 WW domain by all-atom molecular dynamics simulations

This paper examines the folding mechanism of an individual beta-hairpin in the presence of other hairpins by using an off-lattice model of a small triple-stranded antiparallel beta-sheet protein, Pin1 WW domain. The turn zipper model and the hydrophobic collapse model originally developed for a single beta-hairpin in literature is confirmed to be useful in describing beta-hairpins in model Pin1 WW domain. We find that the mechanism for folding a specific hairpin is independent of whether it folds first or second, but the formation process are significantly dependent on temperature. More specifically, beta1-beta2 hairpin folds via the turn zipper model at a low temperature and the hydrophobic collapse model at a high temperature, while the folding of beta2-beta3 hairpin follows the turn zipper model at both temperatures. The change in folding mechanisms is interpreted by the interplay between contact stability (enthalpy) and loop lengths (entropy), the effect of which is temperature dependent.     

A photocontrolled beta-hairpin peptide

Beta-hairpins constitute the smallest beta-type structures in peptides and proteins. The development of highly stable, yet monomeric beta-hairpins based on the tryptophan zipper motif was therefore a remarkable success [A. G. Cochran, N. J. Skelton, M. A. Starovasnik, Proc. Natl. Acad. Sci USA 2001, 98, 5578-5583]. We have been able to design, synthesize and characterize a hairpin based on this motif which incorporates an azobenzene-based photoswitch, that allows for time-resolved folding studies of beta-structures with unprecedented time resolution. At room temperature the trans-azo isomer exhibits a mostly disordered structure; however, light-induced isomerization to the cis-azo form leads to a predominantly extended and parallel conformation of the two peptide parts, which are linked by the novel photoswitch, [3-(3-aminomethyl)phenylazo]phenylacetic acid (AMPP). While in the original sequence the dipeptide Asn-Gly forms a type I' beta-turn which connects the two strands of the hairpin, this role is adopted by the AMPP chromophore in our photoresponsive beta-hairpin that can apparently act as a beta I'-turn mimetic. The beta-hairpin structure was determined and confirmed by NMR spectroscopy, but the folding process can be monitored by pronounced changes in the CD, IR and fluorescence spectra. Finally, incorporation of the structurally and functionally important beta-hairpin motif into proteins by chemical ligation might allow for the photocontrol of protein structures and/or functions.     

probing the gas-phase folding kinetics of peptide ions by IR activated DR-ECD

The effect of infrared (IR) irradiation on the electron capture dissociation (ECD) fragmentation pattern of peptide ions was investigated. IR heating increases the internal energy of the precursor ion, which often amplifies secondary fragmentation, resulting in the formation of w-type ions as well as other secondary fragments. Improved sequence coverage was observed with IR irradiation before ECD, likely due to the increased conformational heterogeneity upon IR heating, rather than faster breakdown of the initially formed product ion complex, as IR heating after ECD did not have similar effect. Although the ECD fragment ion yield of peptide ions does not typically increase with IR heating, in double resonance (DR) ECD experiments, fragment ion yield may be reduced by fast resonant ejection of the charge reduced molecular species, and becomes dependent on the folding state of the precursor ion. In this work, the fragment ion yield was monitored as a function of the delay between IR irradiation and the DR-ECD event to study the gas-phase folding kinetics of the peptide ions. Furthermore, the degree of intracomplex hydrogen transfer of the ECD fragment ion pair was used to probe the folding state of the precursor ion. Both methods gave similar refolding time constants of approximately 1.5 s(-1), revealing that gaseous peptide ions often refold in less than a second, much faster than their protein counterparts. It was also found from the IR-DR-ECD study that the intramolecular H. transfer rate can be an order of magnitude higher than that of the separation of the long-lived c/z product ion complexes, explaining the common observation of c. and z type ions in ECD experiments.     

Well-defined carboxyl-terminated alkyl monolayers grafted onto H-Si(111): packing density from a combined AFM and quantitative IR study

This work demonstrates that well-defined mixed carboxyl-terminated/methyl-terminated alkyl monolayers can be prepared in one step on H-terminated Si(111) via direct photochemical hydrosilylation of undecylenic acid and 1-decene mixtures. As evidenced by AFM imaging and IR spectroscopy, a final rinse in hot acetic acid leaves the functionalized surface atomically smooth and perfectly free of physisorbed contaminants while unwanted material remains atop the monolayer with most other common solvents. The compositional surface chemistry was determined from a truly quantitative IR (ATR geometry) study in the range of 900-4000 cm(-)(1). Results prove that neither surface oxidation nor grafting through the carboxyl end groups occurs. Monolayers are fairly dense for such bulky end groups, with a total molecular surface density of approximately 2.7 10(14) cm(-)(2) corresponding to a surface coverage of 0.35 (maximum theoretical value approximately 0.5). Careful analysis of the CH- and COOH-related IR bands reveals that the composition of the grafted layers is richer in acid chains than the starting grafting mixture. A simple model is presented that shows that the grafting kinetics is about twice as fast for undecylenic acid as for 1-decene. Complementary electrochemical impedance measurements indicate the excellent electronic properties of the interface with a very low density of gap states. They also show that the acid terminal groups promote the penetration of water in the outer part of the organic film.