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.     

Cytochrome c unfolding on an anionic surface

It is now well accepted that the adsorption of proteins to solid supports sometimes involves surface-mediated unfolding. A detailed understanding of the adsorption and surface-mediated unfolding process is lacking. We selected a well studied protein, horse heart cytochrome c, and a weakly ionic support to examine some of the characteristics of protein adsorption under near-physiological conditions. We used high-performance liquid chromatography (HPLC) to investigate the effect of temperature on surface-mediated unfolding. Samples of cytochrome c were introduced to an anionic support, and a NaCl gradient was used to desorb the protein at different times and temperatures. The profiles and retention times were monitored to examine the adhesive properties of cytochrome c to the anionic support. We found that protein retention increased with time at temperatures as low as 0 degrees C, and a significant loss of cytochrome c occurred between 55 degrees C and 70 degrees C. The loss of recovery of cytochrome c indicates irreversible surface-mediated unfolding. The changes in retention time may indicate more subtle transitions, including reversible surface-mediated unfolding of cytochrome c. These results suggest that perturbations in the structure as well as unfolding of cytochrome c can be detected at a lower temperature on an anionic surface than in solution thereby acting like a catalyst for protein unfolding.     

Molecular folding screen: folding and unfolding of 1,8-anthrylene-ethynylene oligomers by photochemical cycloaddition and thermal cycloreversion

Molecular folding screens consisting of anthracene plates and acetylene linkers stereoselectively fold into a zigzag form by [4 + 4]photocycloaddition, and unfold by thermal cycloreversion.     

Detailed unfolding and folding of gaseous ubiquitin ions characterized by electron capture dissociation

The unfolding enthalpy of the native state of ubiquitin in solution is 5 to 8 times that of its gaseous ions, as determined by electron capture dissociation (ECD) mass spectrometry. Although two-state folding occurs in solution, the three-state gaseous process proposed for this by Clemmer and co-workers based on ion mobility data is supported in general by ECD mass spectra, including relative product yields, distinct Delta H(unfolding) values between states, site-specific melting temperatures, and folding kinetics indicating a cooperative process. ECD also confirms that the 13+ ions represent separate conformers, possibly with side-chain solvated alpha-helical structures. However, the ECD data on the noncovalent bonding in the 5+ to 13+ ions, determined overall in 69 of the 75 interresidue sites, shows that thermal unfolding proceeds via a diversity of intermediates whose conformational characteristics also depend on charge site locations. As occurs with increased acidity in solution, adding 6 protons to the 5+ ions completely destroys their tertiary noncovalent bonding. However, solvation of the newly protonated sites to the backbone instead increases the stability of the secondary structure (possibly an alpha-helix) of these gaseous ions, while in solution these new sites aid denaturation by solvation in the aqueous medium. Extensive ion equilibration can lead to even more compact and diverse conformers. The three-state unfolding of gaseous ubiquitin appears to involve ensembles of individual chain conformations in a "folding funnel" of parallel reaction paths. This also provides a further caution for characterizing solution conformers from their gas-phase behavior.     

Modeling of folding and unfolding mechanisms in alanine-based alpha-helical polypeptides

alpha-Helix formation is known to be opposed by the entropy loss due to the folding and favored by the energy of molecular interactions. However, the underlying mechanism of these factors is still being discussed. Here we have used the experimental and calculation data for short alanine-based peptides embedded in water to model the mechanism of helix folding and unfolding and to calculate microscopically the free energy factors of alanine in the frame of helix coil conformational integrals. Classical helix-coil transition theories take into account the interactions in a peptide chain only if the i, i + 3 peptide bond participates in hydrogen bonding. But quantum mechanical calculations showed that interactions of the i, i + 2 peptide bond play an important role in helix folding too. We also included the short-range repulsive interactions due to molecular steric clashes and the end effects due to polar/hydrogen-bonding interactions at the N and C termini. The helix and coil regions of peptide conformational space were defined using an experimental steric criterion for hydrogen bonding. Arginine helix propensity was discussed and estimated. Monte Carlo numerical simulations of thermodynamics and kinetics for the 21 amino acid alpha-helical polypeptide Ac-A5(AAARA)3A-NMe were carried out and found to be in an agreement with the experimental results.     

Design,synthesis, and surface activity of amphiphilic perfluorinated oxazoline polymers

Polymers derived from 2-[2-(perfluoroalkyl)ethyl]-2-oxazolines have been prepared in very good yields with molecular weights up to 50000. The polymers are crystalline with melting temperatures (Tm's) varying from 180 to 232°C depending on the length and type of the perfluorinated side chain, have good thermal stability and are soluble in fluorinated solvents. They are effective as modifiers of surfaces of hydrophilic polymers, such as nylon-6,6, and give acceptable complement activation parameters. Factor C3a (μg/ml) released is 5.3 in contrast to positive control value of 45.2 (activated sepharose) and negative control = 4.8. A 1:1 mixture of 2-octyl-2-oxazoline and 2-[2-(perfluorohexyl)ethyl]-2-oxazoline copolymerize to give a block copolymer that exhibits two separate Tm's. The partially perfluorinated polymers also form stable monolayers at air-water interface.     

Bayesian analysis of folding and unfolding time series of single-forced RNAs

On the basis of a coarse-grain physical model of the folding and unfolding of single-forced RNAs conducted in light tweezer experiments, we theoretically investigate the feasibility of inferring the RNA's intrinsic kinetic parameters from the noisy time series of the molecule's extension. A Bayesian approach using Monte Carlo Markov Chain is proposed. We prove that this statistical approach is efficient and accurate in inferring the molecule's physical parameters, even if the experimental data are yielded under a narrow range of forces.     

Direct UV Raman monitoring of 3(10)-helix and pi-bulge premelting during alpha-helix unfolding

We used UV resonance Raman (UVRR) spectroscopy exciting at approximately 200 nm within the peptide bond pi --> pi* transitions to selectively study the amide vibrations of peptide bonds during alpha-helix melting. The dependence of the amide frequencies on their Psi Ramachandran angles and hydrogen bonding enables us, for the first time, to experimentally determine the temperature dependence of the peptide bond Psi Ramachandran angle population distribution of a 21-residue mainly alanine peptide. These Psi distributions allow us to easily discriminate between alpha-helix, 3(10)-helix and pi-helix/bulge conformations, obtain their individual melting curves, and estimate the corresponding Zimm and Bragg parameters. A striking finding is that alpha-helix melting is more cooperative and shows a higher melting temperature than previously erroneously observed. These Psi distributions also enable the experimental determination of the Gibbs free energy landscape along the Psi reaction coordinate, which further allows us to estimate the free energy barriers along the AP melting pathway. These results will serve as a benchmark for the numerous untested theoretical studies of protein and peptide folding.     

Millisecond time-scale folding and unfolding of DNA hairpins using rapid-mixing stopped-flow kinetics

We report stopped-flow kinetics experiments to study the folding and unfolding of 5 base-pair stem and 21 nucleotide polythymidine loop DNA hairpins over various concentrations of NaCl. The reactions occurred on a time scale of milliseconds, considerably longer than the microsecond time scale suggested by previous kinetics studies of similar-sized hairpins. In comparison to a recent fluorescence correlation spectroscopy study (J. Am. Chem. Soc. 2006, 128, 1240-1249), we suggest the microsecond time-scale reactions are due to intermediate states and the millisecond time-scale reactions reported here are due to the formation of the fully folded DNA hairpin. These results support our view that DNA hairpin folding occurs via a minimum three-state mechanism.     

Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations

Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.     

Determining the folding and unfolding rate constants of nucleic acids by biosensor. Application to telomere G-quadruplex

Nucleic acid molecules may fold into secondary structures, and the formation of such structures is involved in many biological processes and technical applications. The folding and unfolding rate constants define the kinetics of conformation interconversion and the stability of these structures and is important in realizing their functions. We developed a method to determine these kinetic parameters using an optical biosensor based on surface plasmon resonance. The folding and unfolding of a nucleic acid is coupled with a hybridization reaction by immobilization of the target nucleic acid on a sensor chip surface and injection of a complementary probe nucleic acid over the sensor chip surface. By monitoring the time course of duplex formation, both the folding and unfolding rate constants for the target nucleic acid and the association and dissociation rate constants for the target-probe duplex can all be derived from the same measurement. We applied this method to determine the folding and unfolding rate constants of the G-quadruplex of human telomere sequence (TTAGGG)(4) and its association and dissociation rate constants with the complementary strand (CCCTAA)(4). The results show that both the folding and unfolding occur on the time scale of minutes at physiological concentration of K(+). We speculate that this property might be important for telomere elongation. A complete set of the kinetic parameters for both of the structures allows us to study the competition between the formation of the quadruplex and the duplex. Calculations indicate that the formation of both the quadruplex and the duplex is strand concentration-dependent, and the quadruplex can be efficiently formed at low strand concentration. This property may provide the basis for the formation of the quadruplex in vivo in the presence of a complementary strand.     

Biomolecular and amphiphilic films probed by surface sensitive X-ray and neutron scattering

In this review article we discuss the thin film analytical techniques of interface sensitive X-ray and neutron scattering applied to aligned stacks of amphiphilic bilayers, in particular phospholipid membranes in the fluid L phase. We briefly discuss how the structure, composition, fluctuations and interactions in lipid or synthetic membranes can be studied by modern surface sensitive scattering techniques, using X-rays or neutrons as a probe. These techniques offer an in-situ approach to study lipid bilayer systems in different environments over length scales extending from micrometer to nanometer, both with and without additional membrane-active molecules such as amphiphilic peptides or membrane proteins.     

Simulation studies of protein folding/unfolding equilibrium under polar and nonpolar confinement

We study the equilibrium folding/unfolding thermodynamics of a small globular miniprotein, the Trp cage, that is confined to the interior of a 2 nm radius fullerene ball. The interactions of the fullerene surface are changed from nonpolar to polar to mimic the interior of the GroEL/ES chaperonin that assists proteins to fold in vivo. We find that nonpolar confinement stabilizes the folded state of the protein due to the effects of volume reduction that destabilize the unfolded state and also due to interactions with the fullerene surface. For the Trp cage, polar confinement has a net destabilizing effect that results from the stabilizing confinement and the competitive exclusion effect that keeps the protein away from the surface hydration shell and stronger interactions between charged side chains in the protein and the polar surface that compete against the formation of an ion pair that stabilizes the protein folded state. We show that confinement effects due to volume reduction can be overcome by sequence-specific interactions of the protein side chains with the encapsulating surface. This study shows that there is a complex balance among many competing effects that determine the mechanism of GroEL chaperonin in enhancing the folding rate of polypeptide inside its cavity.     

New insight into surface properties of LB film of an amphiphilic terpolymer

Amphiphilic terpolymer(TPR) exhibits good film-forming behavior on pure water observed by means ofπ-A isotherms.To gain insight into physical properties of TPR,the films have been deposited on silicon substrates at different surface pressure by Langmuir-Blodgett(LB) technique.It was found that the increase in peak intensities of stretching mode was due to orderly packing of the films.The contact angles increased with increasing surface pressure,indicating an increase in hydrophobicity due to dense packing of chains of TPR.The cyclic voltammetric measurements indicated that TPR showed good current shielding effect for electron-transfer.In a word,LB films of TPR can produce a variety of structures with varied topography,enabling us to control not only the functionality of the surface,but also the interfacial transport characteristics.     

Reversible folding/unfolding of small a-helix in explicit solvent investigated by ABEEMσπ/MM

Science China Chemistry - We have performed molecular dynamics simulations on the reversible folding/unfolding of small α-helix (short Ala based peptide Ala5) in explicit water solvent in...     

Integrated prediction of protein folding and unfolding rates from only size and structural class

Protein stability, folding and unfolding rates are all determined by the multidimensional folding free energy surface, which in turn is dictated by factors such as size, structure, and amino-acid sequence. Work over the last 15 years has highlighted the role of size and 3D structure in determining folding rates, resulting in many procedures for their prediction. In contrast, unfolding rates are thought to depend on sequence specifics and be much more difficult to predict. Here we introduce a minimalist physics-based model that computes one-dimensional folding free energy surfaces using the number of aminoacids (N) and the structural class (α-helical, all-β, or α-β) as only protein-specific input. In this model N sets the overall cost in conformational entropy and the net stabilization energy, whereas the structural class defines the partitioning of the stabilization energy between local and non-local interactions. To test its predictive power, we calibrated the model empirically and implemented it into an algorithm for the PREdiction of Folding and Unfolding Rates (PREFUR). We found that PREFUR predicts the absolute folding and unfolding rates of an experimental database of 52 proteins with accuracies of ±0.7 and ±1.4 orders of magnitude, respectively (relative to experimental spans of 6 and 8 orders of magnitude). Such prediction uncertainty for proteins vastly varying in size and structure is only two-fold larger than the differences in folding (±0.34) and unfolding rates (±0.7) caused by single-point mutations. Moreover, PREFUR predicts protein stability with an accuracy of ±6.3 kJ mol(-1), relative to the 5 kJ mol(-1) average perturbation induced by single-point mutations. The remarkable performance of our simplistic model demonstrates that size and structural class are the major determinants of the folding landscapes of natural proteins, whereas sequence variability only provides the final 10-20% tuning. PREFUR is thus a powerful bioinformatic tool for the prediction of folding properties and analysis of experimental data.