Nanopore unzipping of individual DNA hairpin molecules

We have used the nanometer scale alpha-Hemolysin pore to study the unzipping kinetics of individual DNA hairpins under constant force or constant loading rate. Using a dynamic voltage control method, the entry rate of polynucleotides into the pore and the voltage pattern applied to induce hairpin unzipping are independently set. Thus, hundreds of unzipping events can be tested in a short period of time (few minutes), independently of the unzipping voltage amplitude. Because our method does not entail the physical coupling of the molecules under test to a force transducer, very high throughput can be achieved. We used our method to study DNA unzipping kinetics at small forces, which have not been accessed before. We find that in this regime the static unzipping times decrease exponentially with voltage with a characteristic slope that is independent of the duplex region sequence, and that the intercept depends strongly on the duplex region energy. We also present the first nanopore dynamic force measurements (time varying force). Our results are in agreement with the approximately logV dependence at high V (where V is the loading rate) observed by other methods. The extension of these measurements to lower loading rates reveals a much weaker dependence on V.     

Stretching short sequences of DNA with constant force axial optical tweezers

Single-molecule techniques for stretching DNA of contour lengths less than a kilobase are fraught with experimental difficulties. However, many interesting biological events such as histone binding and protein-mediated looping of DNA, occur on this length scale. In recent years, the mechanical properties of DNA have been shown to play a significant role in fundamental cellular processes like the packaging of DNA into compact nucleosomes and chromatin fibers. Clearly, it is then important to understand the mechanical properties of short stretches of DNA. In this paper, we provide a practical guide to a single-molecule optical tweezing technique that we have developed to study the mechanical behavior of DNA with contour lengths as short as a few hundred basepairs. The major hurdle in stretching short segments of DNA is that conventional optical tweezers are generally designed to apply force in a direction lateral to the stage (see Fig. 1). In this geometry, the angle between the bead and the coverslip, to which the DNA is tethered, becomes very steep for submicron length DNA. The axial position must now be accounted for, which can be a challenge, and, since the extension drags the microsphere closer to the coverslip, steric effects are enhanced. Furthermore, as a result of the asymmetry of the microspheres, lateral extensions will generate varying levels of torque due to rotation of the microsphere within the optical trap since the direction of the reactive force changes during the extension. Alternate methods for stretching submicron DNA run up against their own unique hurdles. For instance, a dual-beam optical trap is limited to stretching DNA of around a wavelength, at which point interference effects between the two traps and from light scattering between the microspheres begin to pose a significant problem. Replacing one of the traps with a micropipette would most likely suffer from similar challenges. While one could directly use the axial potential to stretch the DNA, an active feedback scheme would be needed to apply a constant force and the bandwidth of this will be quite limited, especially at low forces. We circumvent these fundamental problems by directly pulling the DNA away from the coverslip by using a constant force axial optical tweezers. This is achieved by trapping the bead in a linear region of the optical potential, where the optical force is constant-the strength of which can be tuned by adjusting the laser power. Trapping within the linear region also serves as an all optical force-clamp on the DNA that extends for nearly 350 nm in the axial direction. We simultaneously compensate for thermal and mechanical drift by finely adjusting the position of the stage so that a reference microsphere stuck to the coverslip remains at the same position and focus, allowing for a virtually limitless observation period.     

Mechanical unfolding of RNA: from hairpins to structures with internal multiloops

Mechanical unfolding of RNA structures, ranging from hairpins to ribozymes, using laser optical tweezer experiments have begun to reveal the features of the energy landscape that cannot be easily explored using conventional experiments. Upon application of constant force (f), RNA hairpins undergo cooperative transitions from folded to unfolded states whereas subdomains of ribozymes unravel one at a time. Here, we use a self-organized polymer model and Brownian dynamics simulations to probe mechanical unfolding at constant force and constant-loading rate of four RNA structures of varying complexity. For simple hairpins, such as P5GA, application of constant force or constant loading rate results in bistable cooperative transitions between folded and unfolded states without populating any intermediates. The transition state location (DeltaxFTS) changes dramatically as the loading rate is varied. At loading rates comparable to those used in laser optical tweezer experiments, the hairpin is plastic, with DeltaxFTS being midway between folded and unfolded states; whereas at high loading rates, DeltaxFTS moves close to the folded state, i.e., RNA is brittle. For the 29-nucleotide TAR RNA with the three-nucleotide bulge, unfolding occurs in a nearly two-state manner with an occasional pause in a high free energy metastable state. Forced unfolding of the 55 nucleotides of the Hepatitis IRES domain IIa, which has a distorted L-shaped structure, results in well-populated stable intermediates. The most stable force-stabilized intermediate represents straightening of the L-shaped structure. For these structures, the unfolding pathways can be predicted using the contact map of the native structures. Unfolding of a RNA motif with internal multiloop, namely, the 109-nucleotide prohead RNA that is part of the 29 DNA packaging motor, at constant value of rf occurs with three distinct rips that represent unraveling of the paired helices. The rips represent kinetic barriers to unfolding. Our work shows 1), the response of RNA to force is largely determined by the native structure; and 2), only by probing mechanical unfolding over a wide range of forces can the underlying energy landscape be fully explored.     

Efficient golden gate assembly of DNA constructs for single molecule force spectroscopy and imaging

Single-molecule techniques such as optical tweezers and fluorescence imaging are powerful tools for probing the biophysics of DNA and DNA-protein interactions. The application of these methods requires efficient approaches for creating designed DNA structures with labels for binding to a surface or microscopic beads. In this paper, we develop a simple and fast technique for making a diverse range of such DNA constructs by combining PCR amplicons and synthetic oligonucleotides using golden gate assembly rules. We demonstrate high yield fabrication of torsionally-constrained duplex DNA up to 10 kbp in length and a variety of DNA hairpin structures. We also show how tethering to a cross-linked antibody substrate significantly enhances measurement lifetime under high force. This rapid and adaptable fabrication method streamlines the assembly of DNA constructs for single molecule biophysics.     

Heterogeneous folding of the trpzip hairpin: full atom simulation and experiment

The beta-hairpin trpzip2 can be tuned continuously from a two-state folder to folding on a rough energy landscape without a dominant refolding barrier. At high denaturant concentration, this extremely stable peptide exhibits a single apparent "two-state" transition temperature when monitored by different spectroscopic techniques. However, under optimal folding conditions the hairpin undergoes an unusual folding process with three clusters of melting transitions ranging from 15 degrees C to 160 degrees C, as monitored by 12 different experimental and computational observables. We explain this behavior in terms of a rough free energy landscape of the unfolded peptide caused by multiple tryptophan interactions and alternative backbone conformations. The landscape is mapped out by potentials of mean force derived from replica-exchange molecular dynamics simulations. Implications for deducing cooperativity from denaturant titrations, for the origin of folding cooperativity, and for the folding of thermophilic proteins are pointed out. trpzip is an excellent small tunable model system for the glass-like folding transitions predicted by landscape theory.     

Non-specific binding of Na+ and Mg2+ to RNA determined by force spectroscopy methods

RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to both monovalent and multivalent ions. Despite recent advances, we do not have general methods to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, and the thermodynamic parameters for secondary structure prediction have only been derived at 1M [Na(+)]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed monovalent/Mg(2+) salt conditions. We measure the unfolding and folding rupture forces and apply Kramers theory to extract accurate information about the hairpin free energy landscape under tension at a wide range of ionic conditions. We obtain non-specific corrections for the free energy of formation of the RNA hairpin and measure how the distance of the transition state to the folded state changes with force and ionic strength. We experimentally validate the Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test the approximate rule by which the non-specific binding affinity of divalent cations at a given concentration is equivalent to that of monovalent cations taken at 100-fold concentration for small molecular constructs.     

A shorter peptide model from staphylococcal nuclease for the folding-unfolding equilibrium of a beta-hairpin shows that unfolded state has significant contribution from compact conformational states

It is important to understand the conformational features of the unfolded state in equilibrium with folded state under physiological conditions. In this paper, we consider a short peptide model LMYKGQPM from staphylococcal nuclease to model the conformational equilibrium between a hairpin conformation and its unfolded state using molecular dynamics simulation under NVT conditions at 300K using GROMOS96 force field. The free energy landscape has overall funnel-like shape with hairpin conformations sampling the minima. The "unfolded" state has a higher free energy of approximately 12kJ/mol with respect to native hairpin minimum and occupies a plateau region. We find that the unfolded state has significant contributions from compact conformations. Many of these conformations have hairpin-like topology. Further, these compact conformational forms are stabilized by hydrophobic interactions. Conversion between native and non-native hairpins occurs via unfolded states. Frequent conversions between folded and unfolded hairpins are observed with single exponential kinetics. We compare our results with the emerging picture of unfolded state from both experimental and theoretical studies.     

40 Probing single molecule RNA folding using temperature and force

Nucleic acids can be unfolded either by temperature, such as in UV melting, or by mechanical force using optical tweezers. In UV melting experiments, the folding free energy of nucleic acids at mesophilic temperatures are extrapolated from unfolding occurring at elevated temperatures. Additionally, single molecule unfolding experiments are typically performed only at room temperature, preventing calculation of changes in enthalpy and entropy. Here, we present temperature-controlled optical tweezers suitable for studying folding of single RNA molecules at physiological temperatures. Constant temperatures between 22 and 37?°C are maintained with an accuracy of 0.1?°C, whereas the optical tweezers display a spatial resolution of ～1?nm over the temperature range. Using this instrument, we measured the folding thermodynamics and kinetics of a 20-base-pair RNA hairpin by force-ramp and constant force experiments. Between 22 and 37?°C, the hairpin unfolds and refolds in a single step. Increasing temperature decreases the stability of the hairpin and thus decreases the force required to unfold it. The equilibrium force, at which unfolding and refolding rates are equal, drops ～1?pN as temperature increases every 5?°C. At each temperature, the folding energy can be quantified by reversible work done to unfold the RNA and from the equilibrium constant at constant forces. Over the experimental temperature range, the folding free energy of the hairpin depends linearly on temperature, indicating that ΔH is constant. The measured folding thermodynamics are further compared with the nearest neighbor calculations using Turner’s parameters of nucleic acid folding energetics.     

Elastic light scattering from single cells: orientational dynamics in optical trap

Light-scattering diagrams (phase functions) from single living cells and beads suspended in an optical trap were recorded with 30-ms time resolution. The intensity of the scattered light was recorded over an angular range of 0.5-179.5 degrees using an optical setup based on an elliptical mirror and rotating aperture. Experiments revealed that light-scattering diagrams from biological cells exhibit significant and complex time dependence. We have attributed this dependence to the cell's orientational dynamics within the trap. We have also used experimentally measured phase function information to calculate the time dependence of the optical radiation pressure force on the trapped particle and show how it changes depending on the orientation of the particle. Relevance of these experiments to potential improvement in the sensitivity of label-free flow cytometry is discussed.     

Single-molecule manipulation of macromolecules on GUV or SUV membranes using optical tweezers

Despite their wide applications in soluble macromolecules, optical tweezers have rarely been used to characterize the dynamics of membrane proteins, mainly due to the lack of model membranes compatible with optical trapping. Here, we examined optical trapping and mechanical properties of two potential model membranes, giant and small unilamellar vesicles (GUVs and SUVs, respectively) for studies of membrane protein dynamics. We found that optical tweezers can stably trap GUVs containing iodixanol with controlled membrane tension. The trapped GUVs with high membrane tension can serve as a force sensor to accurately detect reversible folding of a DNA hairpin or membrane binding of synaptotagmin-1 C2AB domain attached to the GUV. We also observed that SUVs are rigid enough to resist large pulling forces and are suitable for detecting protein conformational changes induced by force. Our methodologies may facilitate single-molecule manipulation studies of membrane proteins using optical tweezers.     

A coarse-grained model for force-induced protein deformation and kinetics

Force-induced changes in protein conformation are thought to be responsible for certain cellular responses to mechanical force. Changes in conformation subsequently initiate a biochemical response by alterations in, for example, binding affinity to another protein or enzymatic activity. Here, a model of protein extension under external forcing is created inspired by Kramers' theory for reaction rate kinetics in liquids. The protein is assumed to have two distinct conformational states: a relaxed state, C(1), preferred in the absence of external force, and an extended state, C(2), favored under force application. In the context of mechanotransduction, the extended state is a conformation from which the protein can initiate signaling. Appearance and persistence of C(2) are assumed to lead to transduction of the mechanical signal into a chemical one. The protein energy landscape is represented by two harmonic wells of stiffness kappa(1) and kappa(2), whose minima correspond to conformations C(1) and C(2). First passage time t(f) from C(1) to C(2) is determined from the Fokker-Plank equation employing several different approaches found in the literature. These various approaches exhibit significant differences in behavior as force increases. Although the level of applied force and the energy difference between states largely determine equilibrium, the dominant influence on t(f) is the height of the transition state. Distortions in the energy landscape due to force can also have a significant influence, however, exhibiting a weaker force dependence than exponential as previously reported, approaching a nearly constant value at a level of force that depends on the ratio kappa(1)/kappa(2). Two model systems are used to demonstrate the utility of this approach: a short alpha-helix undergoing a transition between two well-defined states and a simple molecular motor.     

Measuring Molecular Motor Forces In Vivo: Implications for Tug-of-War Models of Bidirectional Transport

Molecular motor proteins use the energy released from ATP hydrolysis to generate force and haul cargoes along cytoskeletal filaments. Thus, measuring the force motors generate amounts to directly probing their function. We report on optical trapping methodology capable of making precise in vivo stall-force measurements of individual cargoes hauled by molecular motors in their native environment. Despite routine measurement of motor forces in vitro, performing and calibrating such measurements in vivo has been challenging. We describe the methodology recently developed to overcome these difficulties, and used to measure stall forces of both kinesin-1 and cytoplasmic dynein-driven lipid droplets in Drosophila embryos. Critically, by measuring the cargo dynamics in the optical trap, we find that there is memory: it is more likely for a cargo to resume motion in the same direction—rather than reverse direction—after the motors transporting it detach from the microtubule under the force of the optical trap. This suggests that only motors of one polarity are active on the cargo at any instant in time and is not consistent with the tug-of-war models of bidirectional transport where both polarity motors can bind the microtubules at all times. We further use the optical trap to measure in vivo the detachment rates from microtubules of kinesin-1 and dynein-driven lipid droplets. Unlike what is commonly assumed, we find that dynein’s but not kinesin’s detachment time in vivo increases with opposing load. This suggests that dynein’s interaction with microtubules behaves like a catch bond.     

Quantitative measurements of force and displacement using an optical trap.

We combined a single-beam gradient optical trap with a high-resolution photodiode position detector to show that an optical trap can be used to make quantitative measurements of nanometer displacements and piconewton forces with millisecond resolution. When an external force is applied to a micron-sized bead held by an optical trap, the bead is displaced from the center of the trap by an amount proportional to the applied force. When the applied force is changed rapidly, the rise time of the displacement is on the millisecond time scale, and thus a trapped bead can be used as a force transducer. The performance can be enhanced by a feedback circuit so that the position of the trap moves by means of acousto-optic modulators to exert a force equal and opposite to the external force applied to the bead. In this case the position of the trap can be used to measure the applied force. We consider parameters of the trapped bead such as stiffness and response time as a function of bead diameter and laser beam power and compare the results with recent ray-optic calculations.     

A semiflexible polymer model applied to loop formation in DNA hairpins

A statistical mechanical "zipper" model is applied to describe the equilibrium melting of short DNA hairpins with poly(dT) loops ranging from 4 to 12 bases in the loop. The free energy of loop formation is expressed in terms of the persistence length of the chain. This method provides a new measurement of the persistence length of single-stranded DNA, which is found to be approximately 1.4 nm for poly(dT) strands in 100 mM NaCl. The free energy of the hairpin relative to the random coil state is found to scale with the loop size with an apparent exponent of > or = 7, much larger than the exponent of approximately 1.5-1.8 expected from considerations of loop entropy alone. This result indicates a strong dependence of the excess stability of the hairpins, from stacking interactions of the bases within the loop, on the size of the loop. We interpret this excess stability as arising from favorable hydrophobic interactions among the bases in tight loops and which diminish as the loops get larger. Free energy profiles along a generalized reaction coordinate are calculated from the equilibrium zipper model. The transition state for hairpin formation is identified as an ensemble of looped conformations with one basepair closing the loop, and with a lower enthalpy than the random coil state. The equilibrium model predicts apparent activation energy of approximately -11 kcal/mol for the hairpin closing step, in remarkable agreement with the value obtained from kinetics measurements.     

Analysis of force generation during flagellar assembly through optical trapping of free-swimming Chlamydomonas reinhardtii

Many studies have used velocity measurements, waveform analyses, and theoretical flagella models to investigate the establishment, maintenance, and function of flagella of the biflagellate green algae Chlamydomonas reinhardtii. We report the first direct measurement of Chlamydomonas flagellar swimming force. Using an optical trap ("optical tweezers") we detect a 75% decrease in swimming force between wild type (CC124) cells and mutants lacking outer flagellar dynein arms (oda1). This difference is consistent with previous estimates and validates the force measurement approach. To examine mechanisms underlying flagella organization and function, we deflagellated cells and examined force generation during flagellar regeneration. As expected, fully regenerated flagella are functionally equivalent to flagella of untreated wild type cells. However, analysis of swimming force vs. flagella length and the increase in force over regeneration time reveals intriguing patterns where increases in force do not always correspond with increases in length. These investigations of flagellar force, therefore, contribute to the understanding of Chlamydomonas motility, describe phenomena surrounding flagella regeneration, and demonstrate the advantages of the optical trapping technique in studies of cell motility.     

Single-molecule manipulation of double-stranded DNA using optical tweezers: interaction studies of DNA with RecA and YOYO-1.

By using optical tweezers and a specially designed flow cell with an integrated glass micropipette, we constructed a setup similar to that of Smith et al. (Science 271:795-799, 1996) in which an individual double-stranded DNA (dsDNA) molecule can be captured between two polystyrene beads. The first bead is immobilized by the optical tweezers and the second by the micropipette. Movement of the micropipette allows manipulation and stretching of the DNA molecule, and the force exerted on it can be monitored simultaneously with the optical tweezers. We used this setup to study elongation of dsDNA by RecA protein and YOYO-1 dye molecules. We found that the stability of the different DNA-ligand complexes and their binding kinetics were quite different. The length of the DNA molecule was extended by 45% when RecA protein was added. Interestingly, the speed of elongation was dependent on the external force applied to the DNA molecule. In experiments in which YOYO-1 was added, a 10-20% extension of the DNA molecule length was observed. Moreover, these experiments showed that a change in the applied external force results in a time-dependent structural change of the DNA-YOYO-1 complex, with a time constant of approximately 35 s (1/e2). Because the setup provides an oriented DNA molecule, we determined the orientation of the transition dipole moment of YOYO-1 within DNA by using fluorescence polarization. The angle of the transition dipole moment with respect to the helical axis of the DNA molecule was 69 degrees +/- 3.     

Preliminary dielectric measurement and analysis protocol for determining the melting temperature and binding energy of short sequences of DNA in solution

Measurement of the real dielectric constant of bulk buffer solutions containing short sequences of DNA as a function of temperature through the DNA melting or denaturiztion transition can be used to determine melting temperatures, T(m), and to estimate the binding energy of the complimentary strands. We describe a preliminary dielectric measurement and analysis protocol to determine these parameters and its application to two known short sequences. The relative real dielectric constant for the bulk solutions was determined over the frequency range of 50 Hz-20 kHz and temperature range of <40-65 degrees C. The measurements were performed on dilute solutions and utilized low electric field strengths. Based on fits to the data by modified sigmoid functions, the melting temperatures, width of transition, and binding energy for the two sequences in solution were estimated. It was observed that the order of the transition appeared to be second order. The results were then compared against predictions of a number of models from the literature that provide theoretical estimates for the melting temperatures of known short sequences of DNA.     

Exploring the energy landscape of a beta hairpin in explicit solvent

We studied the energy landscape of the peptide Ace-GEWTYDDATKTFTVTE-Nme, taken from the C-terminal fragment (41-56) of protein G, in explicit aqueous solution by a highly parallel replica-exchange approach that combines molecular dynamics trajectories with a temperature exchange Monte Carlo process. The combined trajectories in T and configurational space allow a replica to overcome a free energy barrier present at one temperature by increasing T, changing configurations, and cooling in a self-regulated manner, thus allowing sampling of broad regions of configurational space in short (nanoseconds) time scales. The free energy landscape of this system over a wide range of temperatures shows that the system preferentially adopts a beta hairpin structure. However, the peptide also samples other stable ensembles where the peptide adopts helices and helix-turn-helix states, among others. The helical states become increasingly stable at low temperatures, but are slightly less stable than the beta turn ensemble. The energy landscape is rugged at low T, where substates are separated by large energy barriers. These barriers disappear at higher T (approximately 330 K), where the system preferentially adopts a "molten globule" state with structures similar to the beta hairpin.     

Probing Position-Dependent Diffusion in Folding Reactions Using Single-Molecule Force Spectroscopy

Folding of proteins and nucleic acids involves a diffusive search over a multidimensional conformational energy landscape for the minimal-energy structure. When examining the projection of conformational motions onto a one-dimensional reaction coordinate, as done in most experiments, the diffusion coefficient D is generally position dependent. However, it has proven challenging to measure such position-dependence experimentally. We investigated the position-dependence of D in the folding of DNA hairpins as a simple model system in two ways: first, by analyzing the round-trip time to return to a given extension in constant-force extension trajectories measured by force spectroscopy, and second, by analyzing the fall time required to reach a given extension in force jump measurements. These methods yielded conflicting results: the fall time implied a fairly constant D, but the round-trip time implied variations of over an order of magnitude. Comparison of experiments with computational simulations revealed that both methods were strongly affected by experimental artifacts inherent to force spectroscopy measurements, which obscured the intrinsic position-dependence of D. Lastly, we applied Kramers’s theory to the kinetics of hairpins with energy barriers located at different positions along the hairpin stem, as a crude probe of D at different stem positions, and we found that D did not vary much as the barrier was moved along the reaction coordinate. This work underlines the difficulties faced when trying to deduce position-dependent diffusion coefficients from experimental folding trajectories.