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 and unfolding kinetics of DNA hairpins in flowing solution by multiparameter fluorescence correlation spectroscopy

Dynamic equilibrium between the folded and unfolded conformations of single stranded DNA hairpin molecules containing polythymine hairpin loops was investigated using simultaneous two-beam fluorescence cross-correlation spectroscopy and single beam autocorrelation spectroscopy. The hairpins were end-labeled with a fluorescent dye and a quencher, such that folding and unfolding of the DNA hairpin primary structure caused the dye fluorescence to fluctuate on the same characteristic time scale as the folding and unfolding reaction. These fluctuations were observed as the molecules flowed sequentially between two spatially offset, microscopic detection volumes. Cross-correlation analysis of fluorescence from the two detection volumes revealed the translational diffusion and flow properties of the hairpins, as well as the average molecular occupancy of the two volumes. Autocorrelation analysis of the fluorescence from the individual detection volumes revealed the kinetics of hairpin folding and unfolding, with the parameters relating to diffusion, flow, and molecular occupancy constrained to the values determined from the cross-correlation analysis. This allowed unambiguous characterization of the folding and unfolding kinetics, without the need to determine the hydrodynamic properties by analyzing a separate control sample. The analysis revealed nonexponential relaxation kinetics and DNA size-dependent folding times characteristic of dynamic heterogeneity in the DNA hairpin-forming mechanism.     

DNA as helical ruler: exciton-coupled circular dichroism in DNA conjugates

The structure and properties of oligonucleotide conjugates possessing stilbenedicarboxamide chromophores at both ends of a poly(dA):poly(dT) base-pair domain of variable length have been investigated using a combination of spectroscopic and computational methods. These conjugates form capped hairpin structures in which one stilbene serves as a hairpin linker and the other as a hydrophobic end-cap. The capping stilbene stabilizes the hairpin structures by ca. 2 kcal/mol, making possible the formation of a stable folded structure containing a single A:T base pair. Exciton coupling between the stilbene chromophores has little effect on the absorption bands of capped hairpins. However, exciton-coupled circular dichroism (EC-CD) can be observed for capped hairpins possessing as many as 11 base pairs. Both the sign and intensity of the EC-CD spectrum are sensitive to the number of base pairs separating the stilbene chromophores, as a consequence of the distance and angular dependence of exciton coupling. Calculated spectra obtained using a static vector model based on canonical B-DNA are in good agreement with the experimental spectra. Molecular dynamics simulations show that conformational fluctuations of the capped hairpins result in large deviations of the averaged spectra in both the positive and negative directions. These results demonstrate for the first time the ability of B-DNA to serve as a helical ruler for the study of electronic interactions between aligned chromophores. Furthermore, they provide important tests for atomistic theoretical models of DNA.     

A three-state mechanism for DNA hairpin folding characterized by multiparameter fluorescence fluctuation spectroscopy

The folding of a dye-quencher labeled DNA hairpin molecule was investigated using fluorescence autocorrelation and cross-correlation spectroscopy (FCS) and photon counting histogram analysis (PCH). The autocorrelation and cross-correlation measurements revealed the flow and diffusion times of the DNA molecules through two spatially offset detection volumes, the relaxation time of the folding reaction, and the total concentration of DNA molecules participating in the reaction. The PCH measurements revealed the equilibrium distribution of DNA molecules in folded and unfolded conformations and the specific brightnesses of the fluorophore in each conformational state. These measurements were carried out over a range of NaCl concentrations, from those that favored the open form of the DNA hairpin to those that favored the closed form. DNA melting curves obtained from each sample were also analyzed for comparison. It was found that the reactant concentrations were depleted as the reaction progressed and that the equilibrium distributions measured by FCS and PCH deviated from those obtained from the melting curve analyses. These observations suggest a three-state mechanism for the DNA hairpin folding reaction that involves a stable intermediate form of the DNA hairpin. The reaction being probed by FCS and PCH is suggested to be a rapid equilibrium between open and intermediate conformations. Formation of the fully closed DNA hairpin is suggested to occur on a much longer time scale than the FCS and PCH measurement time. The closed form of the hairpin thus serves as a sink into which the reactants are depleted as the reaction progresses.     

Electrostatics of nucleic acid folding under conformational constraint

RNA folding is enabled by interactions between the nucleic acid and its ion atmosphere, the mobile sheath of aqueous ions that surrounds and stabilizes it. Understanding the ion atmosphere requires the interplay of experiment and theory. However, even an apparently simple experiment to probe the ion atmosphere, measuring the dependence of DNA duplex stability upon ion concentration and identity, suffers from substantial complexity, because the unfolded ensemble contains many conformational states that are difficult to treat accurately with theory. To minimize this limitation, we measured the unfolding equilibrium of a DNA hairpin using a single-molecule optical trapping assay, in which the unfolded state is constrained to a limited set of elongated conformations. The unfolding free energy increased linearly with the logarithm of monovalent cation concentration for several cations, such that smaller cations tended to favor the folded state. Mg(2+) stabilized the hairpin much more effectively at low concentrations than did any of the monovalent cations. Poisson-Boltzmann theory captured trends in hairpin stability measured for the monovalent cation titrations with reasonable accuracy, but failed to do so for the Mg(2+) titrations. This finding is consistent with previous work, suggesting that Poisson-Boltzmann and other mean-field theories fail for higher valency cations where ion-ion correlation effects may become significant. The high-resolution data herein, because of the straightforward nature of both the folded and the unfolded states, should serve as benchmarks for the development of more accurate electrostatic theories that will be needed for a more quantitative and predictive understanding of nucleic acid folding.     

Temperature jump and fast photochemical oxidation probe submillisecond protein folding

We report a new mass-spectrometry-based approach for studying protein-folding dynamics on the submillisecond time scale. The strategy couples a temperature jump with fast photochemical oxidation of proteins (FPOP), whereby folding/unfolding is followed by changes in oxidative modifications by OH radical reactions. Using a flow system containing the protein barstar as a model, we altered the protein's equilibrium conformation by applying the temperature jump and demonstrated that its reactivity with OH free radicals serves as a reporter of the conformational change. Furthermore, we found that the time-dependent increase in mass resulting from free-radical oxidation is a measure of the rate constant for the transition from the unfolded to the first intermediate state. This advance offers the promise that, when extended with mass-spectrometry-based proteomic analysis, the sites and kinetics of folding/unfolding can also be followed on the submillisecond time scale.     

Effects of turn stability on the kinetics of refolding of a hairpin in a beta-sheet

As part of our continuing study of the effects of the turn sequence on the conformational stability as well as the mechanism of folding of a beta-sheet structure, we have undertaken a parallel investigation of the solution structure, conformational stability, and kinetics of refolding of the beta-sheet VFIVDGOTYTEV(D)PGOKILQ. The latter peptide is an analogue of the original Gellman beta-sheet VFITS(D)PGKTYTEV(D)PGOKILQ, wherein the TS(D)PGK turn sequence in the first hairpin has been replaced by VDGO. Thermodynamics studies revealed comparable conformational stability of the two peptides. However, unlike the Gellman peptide, which showed extremely rapid refolding of the first hairpin, early kinetic events associated with the refolding of the corresponding hairpin in the VDGO mutant were found to be significantly slower. A detailed study of the conformation of the modified peptide suggested that hydrophobic interactions might be contributing to its stability. Accordingly, we surmise that the early kinetic events are sensitive to whether the formation of the hairpin is nucleated at the turn or by sequestering of the hydrophobic residues across the strand, before structural rearrangements to produce the nativelike topology. Nucleation of the hairpin at the turn is expected to be intrinsically rapid for a strong turn. However, if the process must involve collapse of hydrophobic side chains, the nucleation should be slower as solvent molecules must be displaced to sequester the hydrophobic residues. These findings reflect the contribution of different forces toward nucleation of hairpins in the mechanism of folding of beta-sheets.     

Folding of a stable DNA motif involves a highly cooperative network of interactions

Hairpins are structural elements that play important roles in the folding and function of RNA and DNA. The extent of cooperativity in folding is an important aspect of the RNA folding problem. We reasoned that an investigation into the origin of cooperativity might be best carried out on a stable nucleic acid system with a limited number of interactions, such as a stable DNA hairpin loop. The stable d(cGNAg) hairpin loop motif (closing base pair in lower case; loop in upper case; N = A, C, G, or T) is stabilized through only three interactions: two loop-loop hydrogen bonds in a sheared GA base pair and a loop-closing base pair interaction. Herein, we investigate this network of interactions and test whether the loop-loop and loop-closing base pair interactions communicate. Thermodynamic measurements of nucleotide analogue substituted oligonucleotides were used to probe the additivity of the interactions. On the basis of double mutant cycles, all interactions were found to be nonadditive and interdependent, suggesting that loop-loop and loop-closing base pair interactions form in a highly cooperative manner. When double mutant cycles were repeated in the absence of the other interaction, nonadditivity was significantly reduced suggesting that coupling is indirect and requires all three interactions in order to be optimal. A cooperative network of interactions helps explain the structural and energetic bases of stability in certain DNA hairpins and paves the way for similar studies in more complex nucleic acid systems.     

Direct Measurement of Single‐Molecule DNA Hybridization Dynamics with Single‐Base Resolution

Herein, we report label‐free detection of single‐molecule DNA hybridization dynamics with single‐base resolution. By using an electronic circuit based on point‐decorated silicon nanowires as electrical probes, we directly record the folding/unfolding process of individual hairpin DNAs with sufficiently high signal‐to‐noise ratio and bandwidth. These measurements reveal two‐level current oscillations with strong temperature dependence, enabling us to determine the thermodynamic and kinetic properties of hairpin DNA hybridization. More importantly, successive, stepwise increases and decreases in device conductance at low temperature on a microsecond timescale are successfully observed, indicating a base‐by‐base unfolding/folding process. The process demonstrates a kinetic zipper model for DNA hybridization/dehybridization at the single base‐pair level. This measurement capability promises a label‐free single‐molecule approach to probe biomolecular interactions with fast dynamics.     

Helical Folding Competing with Unfolded Aggregation in Phenylene Ethynylene Foldamers

The folding and aggregation behavior of a pair of oligo(phenylene ethynylene) (OPE) foldamers are investigated by means of UV/Vis absorption and circular dichroism spectroscopy. With identical OPE backbones, two foldamers, 1 with alkyl side groups and 2 with triethylene glycol side chains, manifest similar helical conformations in solutions in n‐hexane and methanol, respectively. However, disparate and competing folding and aggregation processes are observed in alternative solvents. In cyclohexane, oligomer 1 initially adopts the helical conformation, but the self‐aggregation of unfolded chains, as a minor component, gradually drives the folding–unfolding transition eventually to the unfolded aggregate state completely. In contrast, in aqueous solution (CH₃OH/H₂O) both folded and unfolded oligomer 2 appear to undergo self‐association; aggregates of the folded chains are thermodynamically more stable. In solutions with a high H₂O content, self‐aggregation among unfolded oligomers is kinetically favored; these oligomers very slowly transform into aggregates of helical structures with greater thermodynamic stability. The folded–unfolded conformational switch thus takes place with the free (nonaggregated) molecules, and the very slow folding transition is due to the low concentration of molecularly dispersed oligomers.     

beta-Hairpins,alpha-helices,and the intermediates among the secondary structures in the energy landscape of a peptide from a distal beta-hairpin of SH3 domain

Energy landscape of a peptide, extracted from a distal beta-hairpin of src SH3 domain, in explicit water was obtained with the multicanonical molecular dynamics. A variety of beta-hairpins with various strand-strand hydrogen bonds were found in the energy landscape at 300 K. There was no energy barrier between random-coil and hairpins. Thus, the peptide conformation can easily change from the random-coil to the hairpins in the thermal fluctuations at 300 K. The landscape also included two clusters of alpha-helices, among which an energy barrier existed, and besides, these helix clusters were separated from the other conformations. Thus, the free-energy barrier exists among the helices and the other conformations. Intermediate clusters were found between the helix and the hairpin clusters. The current study showed that the isolated state of this peptide in water fluctuates among random-coil, beta-hairpin, and alpha-helix. In SH3 domain, which has a topology of mainly beta-protein, the whole-protein folding may proceed when the segment is folded in the beta-hairpin and the other parts of the protein are coupled with the beta-hairpin in an energetically or kinetically favorite way.     

Isomerization‐ and Temperature‐Jump‐Induced Dynamics of a Photoswitchable β‐Hairpin

Conformational changes in proteins and peptides can be initiated by diverse processes. This raises the question how the variation of initiation mechanisms is connected to differences in folding or unfolding processes. In this work structural dynamics of a photoswitchable β‐hairpin model peptide were initiated by two different mechanisms: temperature jump (T‐jump) and isomerization of a backbone element. In both experiments the structural changes were followed by time‐resolved IR spectroscopy in the nanosecond to microsecond range. When the photoisomerization of the azobenzene backbone switch initiated the folding reaction, pronounced absorption changes related to folding into the hairpin structure were found with a time constant of about 16 μs. In the T‐jump experiment kinetics with the same time constant were observed. For both initiation processes the reaction dynamics revealed the same strong dependence of the reaction time on temperature. The highly similar transients in the microsecond range show that the peptide dynamics induced by T‐jump and isomerization are both determined by the same mechanism and exclude a downhill‐folding process. Furthermore, the combination of the two techniques allows a detailed model for folding and unfolding to be presented: The isomerization‐induced folding process ends in a transition‐state reaction scheme, in which a high energetic barrier of 48 kJ mol^?1 separates unfolded and folded structures.     

Molecular dynamics simulation of the structure, dynamics, and thermostability of the RNA hairpins uCACGg and cUUCGg

Classical replica-exchange molecular dynamics simulations are performed to study structure, dynamics and thermostability of the 14-mer RNA hairpins uCACGg and cUUCGg. Despite of the different sequence and closing base pair of the two systems, recent NMR studies have shown that the tetraloop CACG is strikingly similar in overall geometry and hydrogen bonding to the canonical UUCG tetraloop. On the other hand, the two systems differ significantly in their functionality and thermostability. The simulations confirm the structural similarities of the two RNA hairpins at room temperature but also reveal that the UUCG loop is more flexible than the CACG loop. Concerning the functionality, the CACG loop shows a stronger attitude to donate hydrogens than the UUCG loop, although their global solvent accessible surface is quite similar. The simulations qualitatively reproduce the experimentally found difference in melting temperatures (20 K). In the case of the uCACGg hairpin, the thermal unfolding occurs cooperatively in an all-or-none fashion, while the cUUCGg hairpin shows less cooperativity but exhibits intermediate states during the unfolding process.     

Face-to-face and edge-to-face pi-pi interactions in a synthetic DNA hairpin with a stilbenediether linker

Synthetic conjugates possessing bis(2-hydroxyethyl)stilbene-4,4'-diether linkers (Sd2) form the most stable DNA hairpins reported to date. Factors that affect stability are length and flexibility of the linkers and pi-stacking of the stilbene moiety on the adjacent base pair. The crystal structure of the hairpin d(GT(4)G)-Sd2-d(CA(4)C) was determined at 1.5 A resolution. The conformations of the two molecules in the asymmetric unit differ both in the linker and the stem portions. One of them shows a planar stilbene that is stacked on the adjacent G:C base pair. The other displays considerable rotation between the phenyl rings and an unprecedented edge-to-face orientation of stilbene and base pair. The observation of considerable variations in the conformation of the Sd moiety in the crystal structure allows us to exclude restriction of motion as the reason for the absence of Sd photoisomerization in the hairpins. Conformational differences in the stem portion of the two hairpin molecules go along with different Mg(2+) binding modes. Most remarkable among them is the sequence-specific coordination of a metal ion in the narrow A-tract minor groove. The crystal structure provides unequivocal evidence that a fully hydrated Mg(2+) ion can penetrate the narrow A-tract minor groove, causing the groove to further contract. Overall, the structural data provide a better understanding of the origins of hairpin stability and their photochemical behavior in solution.     

Effect of loop distortion on the stability and structural dynamics of DNA hairpin and dumbbell conjugates

The thermal stability and conformational dynamics of DNA hairpin and dumbbell conjugates having short A-tract base pair domains connected by tri- or hexa(ethylene glycol) linkers is reported. The formation of stable base-paired A-tract hairpins having oligo(ethylene glycol) linkers requires a minimum of four or five A-T base pairs. The formation of base-paired dumbbells having oligo(ethylene glycol) linkers by means of chemical ligation of nicked dumbbells requires a minimum of two A-T base pairs on either side of the nick. Molecular modeling indicates that the hexa(ethylene glycol) linker is sufficiently long to permit formation of strain-free loop regions and B-DNA base pair domains. In contrast, the tri(ethylene glycol) is too short to permit Watson-Crick base pairing between the bases attached to the linker. The shorter linker distorts the duplex, resulting in fluxional behavior in which the base pairs adjacent to the linker and at the open end of the hairpin dissociate on the nanosecond time scale. The loss of interstrand binding energy caused by these fluctuations leads to a difference of approximately 5 degrees C in melting temperature between EG3 and EG6 hairpins. An analysis of the fluxional behavior of the EG3 adjacent base-pair has been used to study the pathways for base flipping and base stacking, including the identification of rotated base (partially flipped) intermediates that have not been described previously for A-T base pairs.     

Kinetic Mechanism of Conformational Switch between Bistable RNA Hairpins

Transitions between the different conformational states play a critical role in many RNA catalytic and regulatory functions. In this study, we use the Kinetic Monte Carlo method to investigate the kinetic mechanism for the conformational switches between bistable RNA hairpins. We find three types of conformational switch pathways for RNA hairpins: refolding after complete unfolding, folding through basepair-exchange pathways and through pseudoknot-assisted pathways, respectively. The result of the competition between the three types of pathways depends mainly on the location of the rate-limiting base stacks (such as the GC base stacks) in the structures. Depending on the structural relationships between the two bistable hairpins, the conformational switch can follow single or multiple dominant pathways. The predicted folding pathways are supported by the activation energy results derived from the Arrhenius plot as well as the NMR spectroscopy data.     

温控分子动力学研究微管蛋白活性肽链的折叠机制

运用温控和常温分子动力学方法, 研究了微管蛋白活性部位Pep1-28肽链的折叠机制, 总模拟时间为380.0 ns. 对于温控分子动力学, 逐渐降温可以清晰显示Pep1-28肽链的折叠途径, 发生明显折叠的温度约为550 K, 其折叠和去折叠可逆机制为U(>1200 K)←→I1(1200-1000 K)←→I2(800 K)←→I3(600 K)←→I4(450 K)←→F1(400 K)←→F2(300 K), 其中U为去折叠态构象, I1、I2、I3和I4是折叠过程中的四个重要的中间态构象, F1和F2是两个结构相近的折叠态构象. 对于常温(300 K)分子动力学, 其构象转变和折叠过程相当迅速, 很难观察到有效、稳定的中间态构象. 尤其引人注意的是, 其折叠态结构陷入了能量局域极小点, 与温控(300 K)的有较大差异, 两者能量差高达297.53 kJ·mol-1. 可见, 温控分子动力学方法不仅清晰地显示多肽和蛋白质折叠过程的重要中间态构象, 为折叠和去折叠机制提供直接、可靠的依据, 而且还有助于跨越较高的构象能垒, 促使多肽和蛋白质折叠以形成全局能量最低的稳定结构.