Light-triggered beta-hairpin folding and unfolding.
Ontology highlight
ABSTRACT:
SUBMITTER: Schrader TE
PROVIDER: S-EPMC1993841 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:2-Aminopurine (2AP) is a fluorescent adenine analog that probes mainly base stacking in nucleic acids. We labeled the loop or the stem of the RNA hairpin gacUACGguc with 2AP to study folding thermodynamics and kinetics at both loci. Thermal melts and fast laser temperature jumps detected by 2AP fluorescence monitored the stability and folding/unfolding kinetics. The observed thermodynamic and kinetic traces of the stem and loop mutants, though strikingly different at a first glance, can be fitted to the same free-energy landscape. The differences between the two probe locations arise because base stacking decreases upon unfolding in the stem, whereas it increases in the loop. We conclude that 2AP is a conservative adenine substitution for mapping out the contributions of different RNA structural elements to the overall folding process. Molecular dynamics (MD) totaling 0.6 ?sec were performed to look at the conformations populated by the RNA at different temperatures. The combined experimental data, and MD simulations lead us to propose a minimal four-state free-energy landscape for the RNA hairpin. Analysis of this landscape shows that a sequential folding model is a good approximation for the full folding dynamics. The frayed state formed initially from the native state is a heterogeneous ensemble of structures whose stem is frayed either from the end or from the loop.
| S-EPMC2995403 | biostudies-literature
Project description:We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active "docked" conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.
| S-EPMC3508787 | biostudies-literature
Project description:We examine the dynamical folding pathways of the C-terminal beta-hairpin of protein G-B1 in explicit solvent at room temperature by means of a transition-path sampling algorithm. In agreement with previous free-energy calculations, the resulting path ensembles reveal a folding mechanism in which the hydrophobic residues collapse first followed by backbone hydrogen-bond formation, starting with the hydrogen bonds inside the hydrophobic core. In addition, the path ensembles contain information on the folding kinetics, including solvent motion. Using the recently developed transition interface sampling technique, we calculate the rate constant for unfolding of the protein fragment and find it to be in reasonable agreement with experiments. The results support the validation of using all-atom force fields to study protein folding.
| S-EPMC218724 | biostudies-literature
Project description:All-atom force fields are now routinely used for more detailed understanding of protein folding mechanisms. However, it has been pointed out that use of all-atom force fields does not guarantee more accurate representations of proteins; in fact, sometimes it even leads to biased structural distributions. Indeed, several issues remain to be solved in force field developments, such as accurate treatment of implicit solvation for efficient conformational sampling and proper treatment of backbone interactions for secondary structure propensities. In this study, we first investigate the quality of several recently improved backbone interaction schemes in AMBER for folding simulations of a beta-hairpin peptide, and further study their influences on the peptide's folding mechanism. Due to the significant number of simulations needed for a thorough analysis of tested force fields, the implicit Poisson-Boltzmann solvent was used in all simulations. The chosen implicit solvent was found to be reasonable for studies of secondary structures based on a set of simulations of both alpha-helical and beta-hairpin peptides with the TIP3P explicit solvent as benchmark. Replica exchange molecular dynamics was also utilized for further efficient conformational sampling. Among the tested AMBER force fields, ff03 and a revised ff99 force field were found to produce structural and thermodynamic data in comparably good agreement with the experiment. However, detailed folding pathways, such as the order of backbone hydrogen bond zipping and the existence of intermediate states, are different between the two force fields, leading to force field-dependent folding mechanisms.
| S-EPMC2242415 | biostudies-literature
Project description:Single-molecule mechanical unfolding experiments have the potential to provide insights into the details of protein folding pathways. To investigate the relationship between force-extension unfolding curves and microscopic events, we performed molecular dynamics simulations of the mechanical unfolding of the C-terminal hairpin of protein G. We have studied the dependence of the unfolding pathway on pulling speed, cantilever stiffness, and attachment points. Under conditions that generate low forces, the unfolding trajectory mimics the untethered, thermally accessible pathway previously proposed based on high-temperature studies. In this stepwise pathway, complete breakdown of backbone hydrogen bonds precedes dissociation of the hydrophobic cluster. Under more extreme conditions, the cluster and hydrogen bonds break simultaneously. Transitions between folding intermediates can be identified in our simulations as features of the calculated force-extension curves.
| S-EPMC1300663 | biostudies-other
Project description:Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces beta-hairpin folding of the glycoprotein Ibalpha beta-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-micros molecular dynamics simulations, we investigate the kinetics of flow-induced beta-hairpin folding. Simulations sampling different flow velocities reveal that under flow, beta-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded beta-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the beta-switch has two stable conformations imprinted on it, namely, loop and hairpin--with flow inducing a transition between the two.
| S-EPMC2920744 | biostudies-literature
Project description:We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.
| S-EPMC1301720 | biostudies-other
Project description: Not available
| S-EPMC8366679 | biostudies-literature
Project description:We examine the dynamical (un)folding pathways of the C-terminal beta-hairpin of protein G-B1 at room temperature in explicit solvent, by employing transition path sampling algorithms. The path ensembles contain information on the folding kinetics, including solvent motion. We determine the transition state ensembles for the two main transitions: 1), the hydrophobic collapse; and 2), the backbone hydrogen bond formation. In both cases the transition state ensembles are characterized by a layer (1) or a strip (2) of water molecules in between the two hairpin strands, supporting the hypothesis of the solvent as lubricant in the folding process. The transition state ensembles do not correspond with saddle points in the equilibrium free-energy landscapes. The kinetic pathways are thus not completely determined by the free-energy landscape. This phenomenon can occur if the order parameters obey different timescales. Using the transition interface sampling technique, we calculate the rate constants for (un)folding and find them in reasonable agreement with experiments, thus supporting the validation of using all-atom force fields to study protein folding.
| S-EPMC1305027 | biostudies-literature
Project description:We have performed the first unbiased folding simulations of the GB1 hairpin in explicit solvent, using hundreds of microsecond-long molecular dynamics simulations (total time: 0.7 ms). Our simulations are initiated from two sets of structures. Starting from an equilibrium unfolded state, we obtain single-exponential folding kinetics with rate coefficients in good agreement (T=350 K) or within an order of magnitude (T=300 K) of the experimental values. However, simulations initiated from unfolded configurations lacking secondary structure result in biexponential kinetics with an additional fast nanosecond kinetic mode. This mode can strongly bias the folding rate estimated from the mean first passage time, when the trials are much shorter than the folding time. We find that the mechanism of the hairpin folding is insensitive to the details of the initial unfolded ensemble and is initiated by correct formation of the turn of the hairpin, followed by the formation of the native hydrogen bonds and hydrophobic contacts, consistent with experimental -value analysis. Subsequent native interactions can be formed either from the turn or from the hairpin termini, helping to explain an apparent discrepancy in experimental results. From our simulations, we also obtain the transition path durations, a critical parameter for single molecule experiments aiming to resolve events along folding pathways. The lengths of transition paths span a wide range, from 50 ps to 140 ns, at 300 K.
| S-EPMC3131341 | biostudies-literature