Thermodynamics of DNA interactions from single molecule stretching experiments

Thymine DNA glycosylase combines sliding, hopping, and nucleosome interactions to efficiently search for 5-formylcytosine

Base excision repair is the main pathway involved in active DNA demethylation. 5-formylcytosine and 5-carboxylcytosine, two oxidized moieties of methylated cytosine, are recognized and removed by thymine DNA glycosylase (TDG) to generate an abasic site. Using single molecule fluorescence experiments, we study TDG in the presence and absence of 5-formylcytosine. TDG exhibits multiple modes of linear diffusion, including hopping and sliding, in search of base modifications. TDG active site variants and truncated N-terminus, reveals these variants alter base modification search and recognition mechanism of TDG. On DNA containing an undamaged nucleosome, TDG is found to either bypass, colocalize with, or encounter but not bypass the nucleosome. Truncating the N-terminus reduces the number of interactions with the nucleosome. Our findings provide mechanistic insights into how TDG searches for modified DNA bases in chromatin.

Plasmonic-Driven Regulation of Biomolecular Activity In Situ

Selective and remote manipulation of activity for biomolecules, including protein, DNA, and lipids, is crucial to elucidate their molecular function and to develop biomedical applications. While advances in tool development, such as optogenetics, have significantly impacted these directions, the requirement for genetic modification significantly limits their therapeutic applications. Plasmonic nanoparticle heating has brought new opportunities to the field, as hot nanoparticles are unique point heat sources at the nanoscale. In this review, we summarize fundamental engineering problems such as plasmonic heating and the resulting biomolecular responses. We highlight the biological responses and applications of manipulating biomolecules and provide perspectives for future directions in the field.

A predictive model for the thermomechanical melting transition of double stranded DNA

By extending the classical Peyrard-Bishop model, we are able to obtain a fully analytical description for the mechanical response of DNA under stretching at variable values of temperature, number of base pairs and intrachains and interchains bonds stiffness. In order to compare elasticity and temperature effects, we first analyze the system in the zero temperature mechanical limit, important to describe several experimental effects including possible hysteresis. We then analyze temperature effects in the framework of equilibrium Statistical Mechanics. In particular, we obtain an analytical expression for the temperature-dependent melting force and unzipping assigned displacement in the thermodynamical limit, also depending on the relative stability of intra vs. inter molecular bonds. Such results coincide with the purely mechanical model in the limit of zero temperature and with the denaturation temperature that we obtain with the classical transfer integral method. Based on our analytical results, we obtain explicitly phase diagrams and cooperativity parameters, where also discreteness effect can be accounted for. The obtained results are successfully applied in reproducing the thermomechanical experimental melting of DNA and the response of DNA hairpins. Due to the generality of the model, exemplified in the proposed analysis of both overstretching and unzipping experiments, we argue that the proposed approach can be extended to other thermomechanically induced molecular melting phenomena. STATEMENT OF SIGNIFICANCE: We obtain a fully analytical description of the complex wiggly energy landscape of two stranded macromolecules under unzipping loading. Based on Equilibrium Statistical Mechanics, we describe the combined thermomechanical effects and the melting transition of double stranded molecules such as nucleic acids. This is proved by quantitatively predicting the experimental behavior of both melting of DNA and DNA hairpins opening. While analytical results have been previously attained under special conditions on the relative stiffness of the covalent vs. non-covalent bonds of the base pairs, our model is completely general in this respect, thus representing a tool in the perspective of the design at the molecular scale. We show that the obtained model can be fully inscribed in the theory of phase transitions giving a new interpretation of the thermomechanical behavior of double stranded molecules.

Statistical physics of DNA hybridization

Deoxyribonucleic acid (DNA) hybridization is at the heart of countless biological and biotechnological processes. Its theoretical modeling played a crucial role, since it has enabled extracting the relevant thermodynamic parameters from systematic measurements of DNA melting curves. In this article, we propose a framework based on statistical physics to describe DNA hybridization and melting in an arbitrary mixture of DNA strands. In particular, we are able to analytically derive closed expressions of the system partition functions for any number N of strings and explicitly calculate them in two paradigmatic situations: (i) a system made of self-complementary sequences and (ii) a system comprising two mutually complementary sequences. We derive the melting curve in the thermodynamic limit (N→∞) of our description, which provides a full justification for the extra entropic contribution that in classic hybridization modeling was required to correctly describe within the same framework the melting of sequences either self-complementary or not. We thus provide a thorough study comprising limit cases and alternative approaches showing how our framework can give a comprehensive view of hybridization and melting phenomena.

A review on nonlinear DNA physics

The study and the investigation of structural and dynamical properties of complex systems have attracted considerable interest among scientists in general and physicists and biologists in particular. The present review paper represents a broad overview of the research performed over the nonlinear dynamics of DNA, devoted to some different aspects of DNA physics and including analytical, quantum and computational tools to understand nonlinear DNA physics. We review in detail the semi-discrete approximation within helicoidal Peyrard-Bishop model and show that localized modulated solitary waves, usually called breathers, can emerge and move along the DNA. Since living processes occur at submolecular level, we then discuss a quantum treatment to address the problem of how charge and energy are transported on DNA and how they may play an important role for the functioning of living cells. While this problem has attracted the attention of researchers for a long time, it is still poorly understood how charge and energy transport can occur at distances comparable to the size of macromolecules. Here, we review a theory based on the mechanism of 'self-trapping' of electrons due to their interaction with mechanical (thermal) oscillation of the DNA structure. We also describe recent computational models that have been developed to capture nonlinear mechanics of DNA in vitro and in vivo, possibly under topological constraints. Finally, we provide some conjectures on potential future directions for this field.

Environment-dependent single-chain mechanics of synthetic polymers and biomacromolecules by atomic force microscopy-based single-molecule force spectroscopy and the implications for advanced polymer materials

"The Tao begets the One. One begets all things of the world." This quote from Tao Te Ching is still inspiring for scientists in chemistry and materials science: The "One" can refer to a single molecule. A macroscopic material is composed of numerous molecules. Although the relationship between the properties of the single molecule and macroscopic material is not well understood yet, it is expected that a deeper understanding of the single-chain mechanics of macromolecules will certainly facilitate the development of materials science. Atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) has been exploited extensively as a powerful tool to study the single-chain behaviors of macromolecules. In this review, we summarize the recent advances in the emerging field of environment-dependent single-chain mechanics of synthetic polymers and biomacromolecules by means of AFM-SMFS. First, the single-chain inherent elasticities of several typical linear macromolecules are introduced, which are also confirmed by one of three polymer models with theoretical elasticities of the corresponding macromolecules obtained from quantum mechanical (QM) calculations. Then, the effects of the external environments on the single-chain mechanics of synthetic polymers and biomacromolecules are reviewed. Finally, the impacts of single-chain mechanics of macromolecules on the development of polymer science especially polymer materials are illustrated.

A Tour de Force on the Double Helix: Exploiting DNA Mechanics To Study DNA-Based Molecular Machines

DNA is both a fundamental building block of life and a fascinating natural polymer. The advent of single-molecule manipulation tools made it possible to exert controlled force on individual DNA molecules and measure their mechanical response. Such investigations elucidated the elastic properties of DNA and revealed its distinctive structural configurations across force regimes. In the meantime, a detailed understanding of DNA mechanics laid the groundwork for single-molecule studies of DNA-binding proteins and DNA-processing enzymes that bend, stretch, and twist DNA. These studies shed new light on the metabolism and transactions of nucleic acids, which constitute a major part of the cell's operating system. Furthermore, the marriage of single-molecule fluorescence visualization and force manipulation has enabled researchers to directly correlate the applied tension to changes in the DNA structure and the behavior of DNA-templated complexes. Overall, experimental exploitation of DNA mechanics has been and will continue to be a unique and powerful strategy for understanding how molecular machineries recognize and modify the physical state of DNA to accomplish their biological functions.

Demodulated standing solitary wave and DNA-RNA transcription

Nonlinear dynamics of DNA molecule at segments where DNA-RNA transcription occurs is studied. Our basic idea is that the solitary wave, moving along the chain, transforms into a demodulated one at these segments. The second idea is that the wave becomes a standing one due to interaction with DNA surrounding, e.g., RNA polymerase molecules. We explain why this is biologically convenient and show that our results match the experimental ones. In addition, we suggest how to experimentally determine crucial constant describing covalent bonds within DNA.

Environmentally Sensitive Fluorescent Nucleoside Analogues for Surveying Dynamic Interconversions of Nucleic Acid Structures

Nucleic acids are characterized by a variety of dynamically interconverting structures that play a major role in transcriptional and translational regulation as well as recombination and repair. To monitor these interconversions, Förster resonance energy transfer (FRET)-based techniques can be used, but require two fluorophores that are typically large and can alter the DNA/RNA structure and protein binding. Additionally, events that do not alter the donor/acceptor distance and/or angular relationship are frequently left undetected. A more benign approach relies on fluorescent nucleobases that can substitute their native counterparts with minimal perturbation, such as the recently developed 2-thienyl-3-hydroxychromone (3HCnt) and thienoguanosine (^th G). To demonstrate the potency of 3HCnt and ^th G in deciphering interconversion mechanisms, we used the conversion of the (-)DNA copy of the HIV-1 primer binding site (-)PBS stem-loop into (+)/(-)PBS duplex, as a model system. When incorporated into the (-)PBS loop, the two probes were found to be highly sensitive to the individual steps both in the absence and the presence of a nucleic acid chaperone, providing the first complete mechanistic description of this critical process in HIV-1 replication. The combination of the two distinct probes appears to be instrumental for characterizing structural transitions of nucleic acids under various stimuli.

A stretched conformation of DNA with a biological role?

We have discovered a well-defined extended conformation of double-stranded DNA, which we call Σ-DNA, using laser-tweezers force-spectroscopy experiments. At a transition force corresponding to free energy change ΔG = 1·57 ± 0·12 kcal (mol base pair)-1 60 or 122 base-pair long synthetic GC-rich sequences, when pulled by the 3'-3' strands, undergo a sharp transition to the 1·52 ± 0·04 times longer Σ-DNA. Intriguingly, the same degree of extension is also found in DNA complexes with recombinase proteins, such as bacterial RecA and eukaryotic Rad51. Despite vital importance to all biological organisms for survival, genome maintenance and evolution, the recombination reaction is not yet understood at atomic level. We here propose that the structural distortion represented by Σ-DNA, which is thus physically inherent to the nucleic acid, is related to how recombination proteins mediate recognition of sequence homology and execute strand exchange. Our hypothesis is that a homogeneously stretched DNA undergoes a 'disproportionation' into an inhomogeneous Σ-form consisting of triplets of locally B-like perpendicularly stacked bases. This structure may ensure improved fidelity of base-pair recognition and promote rejection in case of mismatch during homologous recombination reaction. Because a triplet is the length of a gene codon, we speculate that the structural physics of nucleic acids may have biased the evolution of recombinase proteins to exploit triplet base stacks and also the genetic code.

DNA partitions into triplets under tension in the presence of organic cations, with sequence evolutionary age predicting the stability of the triplet phase

Using atomistic simulations, we show the formation of stable triplet structure when particular GC-rich DNA duplexes are extended in solution over a timescale of hundreds of nanoseconds, in the presence of organic salt. We present planar-stacked triplet disproportionated DNA (Σ DNA) as a possible solution phase of the double helix under tension, subject to sequence and the presence of stabilising co-factors. Considering the partitioning of the duplexes into triplets of base pairs as the first step of operation of recombinase enzymes like RecA, we emphasise the structure-function relationship in Σ DNA. We supplement atomistic calculations with thermodynamic arguments to show that codons for 'phase 1' amino acids (those appearing early in evolution) are more likely than a lower entropy GC-rich sequence to form triplets under tension. We further observe that the four amino acids supposed (in the 'GADV world' hypothesis) to constitute the minimal set to produce functional globular proteins have the strongest triplet-forming propensity within the phase 1 set, showing a series of decreasing triplet propensity with evolutionary newness. The weak form of our observation provides a physical mechanism to minimise read frame and recombination alignment errors in the early evolution of the genetic code.

Protein-nucleic acid interactions of LINE-1 ORF1p

Long interspersed nuclear element 1 (LINE-1 or L1) is the dominant retrotransposon in mammalian genomes. L1 encodes two proteins ORF1p and ORF2p that are required for retrotransposition. ORF2p functions as the replicase. ORF1p is a coiled coil-mediated trimeric, high affinity RNA binding protein that packages its full- length coding transcript into an ORF2p-containing ribonucleoprotein (RNP) complex, the retrotransposition intermediate. ORF1p also is a nucleic acid chaperone that presumably facilitates the proposed nucleic acid remodeling steps involved in retrotransposition. Although detailed mechanistic understanding of ORF1p function in this process is lacking, recent studies showed that the rate at which ORF1p can form stable nucleic acid-bound oligomers in vitro is positively correlated with formation of an active L1 RNP as assayed in vivo using a cell culture-based retrotransposition assay. This rate was sensitive to minor amino acid changes in the coiled coil domain, which had no effect on nucleic acid chaperone activity. Additional studies linking the complex nucleic acid binding properties to the conformational changes of the protein are needed to understand how ORF1p facilitates retrotransposition.

Dissecting the Dynamic Pathways of Stereoselective DNA Threading Intercalation

DNA intercalators that have high affinity and slow kinetics are developed for potential DNA-targeted therapeutics. Although many natural intercalators contain multiple chiral subunits, only intercalators with a single chiral unit have been quantitatively probed. Dumbbell-shaped DNA threading intercalators represent the next order of structural complexity relative to simple intercalators, and can provide significant insights into the stereoselectivity of DNA-ligand intercalation. We investigated DNA threading intercalation by binuclear ruthenium complex [μ-dppzip(phen)4Ru2](4+) (Piz). Four Piz stereoisomers are defined by the chirality of the intercalating subunit (Ru(phen)2dppz) and the distal subunit (Ru(phen)2ip), respectively, each of which can be either right-handed (Δ) or left-handed (Λ). We used optical tweezers to measure single DNA molecule elongation due to threading intercalation, revealing force-dependent DNA intercalation rates and equilibrium dissociation constants. The force spectroscopy analysis provided the zero-force DNA binding affinity, the equilibrium DNA-ligand elongation Δxeq, and the dynamic DNA structural deformations during ligand association xon and dissociation xoff. We found that Piz stereoisomers exhibit over 20-fold differences in DNA binding affinity, from a Kd of 27 ± 3 nM for (Δ,Λ)-Piz to a Kd of 622 ± 55 nM for (Λ,Δ)-Piz. The striking affinity decrease is correlated with increasing Δxeq from 0.30 ± 0.02 to 0.48 ± 0.02 nm and xon from 0.25 ± 0.01 to 0.46 ± 0.02 nm, but limited xoff changes. Notably, the affinity and threading kinetics is 10-fold enhanced for right-handed intercalating subunits, and 2- to 5-fold enhanced for left-handed distal subunits. These findings demonstrate sterically dispersed transition pathways and robust DNA structural recognition of chiral intercalators, which are critical for optimizing DNA binding affinity and kinetics.

Single-molecule studies of high-mobility group B architectural DNA bending proteins

Protein–DNA interactions can be characterized and quantified using single molecule methods such as optical tweezers, magnetic tweezers, atomic force microscopy, and fluorescence imaging. In this review, we discuss studies that characterize the binding of high-mobility group B (HMGB) architectural proteins to single DNA molecules. We show how these studies are able to extract quantitative information regarding equilibrium binding as well as non-equilibrium binding kinetics. HMGB proteins play critical but poorly understood roles in cellular function. These roles vary from the maintenance of chromatin structure and facilitation of ribosomal RNA transcription (yeast high-mobility group 1 protein) to regulatory and packaging roles (human mitochondrial transcription factor A). We describe how these HMGB proteins bind, bend, bridge, loop and compact DNA to perform these functions. We also describe how single molecule experiments observe multiple rates for dissociation of HMGB proteins from DNA, while only one rate is observed in bulk experiments. The measured single-molecule kinetics reveals a local, microscopic mechanism by which HMGB proteins alter DNA flexibility, along with a second, much slower macroscopic rate that describes the complete dissociation of the protein from DNA.

Mechanical properties of DNA origami nanoassemblies are determined by Holliday junction mechanophores

DNA nanoassemblies have demonstrated wide applications in various fields including nanomaterials, drug delivery and biosensing. In DNA origami, single-stranded DNA template is shaped into desired nanostructure by DNA staples that form Holliday junctions with the template. Limited by current methodologies, however, mechanical properties of DNA origami structures have not been adequately characterized, which hinders further applications of these materials. Using laser tweezers, here, we have described two mechanical properties of DNA nanoassemblies represented by DNA nanotubes, DNA nanopyramids and DNA nanotiles. First, mechanical stability of DNA origami structures is determined by the effective density of Holliday junctions along a particular stress direction. Second, mechanical isomerization observed between two conformations of DNA nanotubes at 10–35 pN has been ascribed to the collective actions of individual Holliday junctions, which are only possible in DNA origami with rotational symmetric arrangements of Holliday junctions, such as those in DNA nanotubes. Our results indicate that Holliday junctions control mechanical behaviors of DNA nanoassemblies. Therefore, they can be considered as ‘mechanophores’ that sustain mechanical properties of origami nanoassemblies. The mechanical properties observed here provide insights for designing better DNA nanostructures. In addition, the unprecedented mechanical isomerization process brings new strategies for the development of nano-sensors and actuators.

Force measurements reveal how small binders perturb the dissociation mechanisms of DNA duplex sequences

The force-driven separation of double-stranded DNA is crucial to the accomplishment of cellular processes like genome transactions. Ligands binding to short DNA sequences can have a local stabilizing or destabilizing effect and thus severely affect these processes. Although the design of ligands that bind to specific sequences is a field of intense research with promising biomedical applications, so far, their effect on the force-induced strand separation has remained elusive. Here, by means of AFM-based single molecule force spectroscopy, we show the co-existence of two different mechanisms for the separation of a short DNA duplex and demonstrate how they are perturbed by small binders. With the support of Molecular Dynamics simulations, we evidence that above a critical pulling rate one of the dissociation pathways becomes dominant, with a dramatic effect on the rupture forces. Around the critical threshold, we observe a drop of the most probable rupture forces for ligand-stabilized duplexes. Our results offer a deep understanding of how a stable DNA-ligand complex behaves under force-driven strand separation.

Elastic Properties of Nucleic Acids by Single-Molecule Force Spectroscopy

We review the current knowledge on the use of single-molecule force spectroscopy techniques to extrapolate the elastic properties of nucleic acids. We emphasize the lesser-known elastic properties of single-stranded DNA. We discuss the importance of accurately determining the elastic response in pulling experiments, and we review the simplest models used to rationalize the experimental data as well as the experimental approaches used to pull single-stranded DNA. Applications used to investigate DNA conformational transitions and secondary structure formation are also highlighted. Finally, we provide an overview of the effects of salt and temperature and briefly discuss the effects of contour length and sequence dependence.

A ruthenium dimer complex with a flexible linker slowly threads between DNA bases in two distinct steps

Several multi-component DNA intercalating small molecules have been designed around ruthenium-based intercalating monomers to optimize DNA binding properties for therapeutic use. Here we probe the DNA binding ligand [μ-C4(cpdppz)₂(phen)₄Ru₂]⁴⁺, which consists of two Ru(phen)₂dppz²⁺ moieties joined by a flexible linker. To quantify ligand binding, double-stranded DNA is stretched with optical tweezers and exposed to ligand under constant applied force. In contrast to other bis-intercalators, we find that ligand association is described by a two-step process, which consists of fast bimolecular intercalation of the first dppz moiety followed by ∼10-fold slower intercalation of the second dppz moiety. The second step is rate-limited by the requirement for a DNA-ligand conformational change that allows the flexible linker to pass through the DNA duplex. Based on our measured force-dependent binding rates and ligand-induced DNA elongation measurements, we are able to map out the energy landscape and structural dynamics for both ligand binding steps. In addition, we find that at zero force the overall binding process involves fast association (∼10 s), slow dissociation (∼300 s), and very high affinity (K_d ∼10 nM). The methodology developed in this work will be useful for studying the mechanism of DNA binding by other multi-step intercalating ligands and proteins.

Rigorous study of molecular dynamics of a single dsDNA confined in a nanochannel: Introduction of a critical mobility behaviour

The essential and effective characteristics of a double-stranded DNA (dsDNA) confined in a nanochannel is revisited by employing the rigorous full numerical approach of Molecular Dynamics (MD). The deformation of dsDNA and wall-biomolecule interaction which is critical in highly confined regime has been precisely imposed in numerical simulations. The numerical approach has been justified against available theoretical outcomes. A new and general expression for DNA electrophoretic mobility versus DNA length is extracted from numerical simulation which is out of reach of experimental methods due to practical shortcomings. The newly derived expression suggests an essential correction in the previously proposed expression for the critical case of small DNA molecules and reveals an astonishingly unbeknown trend of small DNA's mobility. Sub-molecular phenomenon of dsDNA melting under the condition of large external force is also studied. Assuming strong electric field exertion, the MD approach aptly demonstrates the elaborate melting phenomenon for dsDNA in sub-molecular scale.

Extracting physical chemistry from mechanics: a new approach to investigate DNA interactions with drugs and proteins in single molecule experiments

In this review we focus on the idea of establishing connections between the mechanical properties of DNA–ligand complexes and the physical chemistry of DNA–ligand interactions.

Equilibrium and kinetics of DNA overstretching modeled with a quartic energy landscape

It is well known that the dsDNA molecule undergoes a phase transition from B-DNA into an overstretched state at high forces. For some time, the structure of the overstretched state remained unknown and highly debated, but recent advances in experimental techniques have presented evidence of more than one possible phase (or even a mixed phase) depending on ionic conditions, temperature, and basepair sequence. Here, we present a theoretical model to study the overstretching transition with the possibility that the overstretched state is a mixture of two phases: a structure with portions of inner strand separation (melted or M-DNA), and an extended phase that retains the basepair structure (S-DNA). We model the double-stranded DNA as a chain composed of n segments of length l, where the transition is studied by means of a Landau quartic potential with statistical fluctuations. The length l is a measure of cooperativity of the transition and is key to characterizing the overstretched phase. By analyzing the different values of l corresponding to a wide spectrum of experiments, we find that for a range of temperatures and ionic conditions, the overstretched form is likely to be a mix of M-DNA and S-DNA. For a transition close to a pure S-DNA state, where the change in extension is close to 1.7 times the original B-DNA length, we find l ? 25 basepairs regardless of temperature and ionic concentration. Our model is fully analytical, yet it accurately reproduces the force-extension curves, as well as the transient kinetic behavior, seen in DNA overstretching experiments.

Cavity approach for modeling and fitting polymer stretching

The mechanical properties of molecules are today captured by single molecule manipulation experiments, so that polymer features are tested at a nanometric scale. Yet devising mathematical models to get further insight beyond the commonly studied force-elongation relation is typically hard. Here we draw from techniques developed in the context of disordered systems to solve models for single and double-stranded DNA stretching in the limit of a long polymeric chain. Since we directly derive the marginals for the molecule local orientation, our approach allows us to readily calculate the experimental elongation as well as other observables at wish. As an example, we evaluate the correlation length as a function of the stretching force. Furthermore, we are able to fit successfully our solution to real experimental data. Although the model is admittedly phenomenological, our findings are very sound. For single-stranded DNA our solution yields the correct (monomer) scale and yet, more importantly, the right persistence length of the molecule. In the double-stranded case, our model reproduces the well-known overstretching transition and correctly captures the ratio between native DNA and overstretched DNA. Also in this case the model yields a persistence length in good agreement with consensus, and it gives interesting insights into the bending stiffness of the native and overstretched molecule, respectively.

Strong DNA deformation required for extremely slow DNA threading intercalation by a binuclear ruthenium complex

DNA intercalation by threading is expected to yield high affinity and slow dissociation, properties desirable for DNA-targeted therapeutics. To measure these properties, we utilize single molecule DNA stretching to quantify both the binding affinity and the force-dependent threading intercalation kinetics of the binuclear ruthenium complex Δ,Δ-[μ‐bidppz‐(phen)₄Ru₂]⁴⁺ (Δ,Δ-P). We measure the DNA elongation at a range of constant stretching forces using optical tweezers, allowing direct characterization of the intercalation kinetics as well as the amount intercalated at equilibrium. Higher forces exponentially facilitate the intercalative binding, leading to a profound decrease in the binding site size that results in one ligand intercalated at almost every DNA base stack. The zero force Δ,Δ-P intercalation K_d is 44 nM, 25-fold stronger than the analogous mono-nuclear ligand (Δ-P). The force-dependent kinetics analysis reveals a mechanism that requires DNA elongation of 0.33 nm for association, relaxation to an equilibrium elongation of 0.19 nm, and an additional elongation of 0.14 nm from the equilibrium state for dissociation. In cells, a molecule with binding properties similar to Δ,Δ-P may rapidly bind DNA destabilized by enzymes during replication or transcription, but upon enzyme dissociation it is predicted to remain intercalated for several hours, thereby interfering with essential biological processes.

Differential contribution of basic residues to HIV-1 nucleocapsid protein's nucleic acid chaperone function and retroviral replication

The human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein contains 15 basic residues located throughout its 55-amino acid sequence, as well as one aromatic residue in each of its two CCHC-type zinc finger motifs. NC facilitates nucleic acid (NA) rearrangements via its chaperone activity, but the structural basis for this activity and its consequences in vivo are not completely understood. Here, we investigate the role played by basic residues in the N-terminal domain, the N-terminal zinc finger and the linker region between the two zinc fingers. We use in vitro ensemble and single-molecule DNA stretching experiments to measure the characteristics of wild-type and mutant HIV-1 NC proteins, and correlate these results with cell-based HIV-1 replication assays. All of the cationic residue mutations lead to NA interaction defects, as well as reduced HIV-1 infectivity, and these effects are most pronounced on neutralizing all five N-terminal cationic residues. HIV-1 infectivity in cells is correlated most strongly with NC’s NA annealing capabilities as well as its ability to intercalate the DNA duplex. Although NC’s aromatic residues participate directly in DNA intercalation, our findings suggest that specific basic residues enhance these interactions, resulting in optimal NA chaperone activity.

Coarse-graining DNA for simulations of DNA nanotechnology

To simulate long time and length scale processes involving DNA it is necessary to use a coarse-grained description. Here we provide an overview of different approaches to such coarse-graining, focussing on those at the nucleotide level that allow the self-assembly processes associated with DNA nanotechnology to be studied. OxDNA, our recently-developed coarse-grained DNA model, is particularly suited to this task, and has opened up this field to systematic study by simulations. We illustrate some of the range of DNA nanotechnology systems to which the model is being applied, as well as the insights it can provide into fundamental biophysical properties of DNA.

Coarse-grained simulations of DNA overstretching

We use a recently developed coarse-grained model to simulate the overstretching of duplex DNA. Overstretching at 23 °C occurs at 74 pN in the model, about 6-7 pN higher than the experimental value at equivalent salt conditions. Furthermore, the model reproduces the temperature dependence of the overstretching force well. The mechanism of overstretching is always force-induced melting by unpeeling from the free ends. That we never see S-DNA (overstretched duplex DNA), even though there is clear experimental evidence for this mode of overstretching under certain conditions, suggests that S-DNA is not simply an unstacked but hydrogen-bonded duplex, but instead probably has a more exotic structure.

Aromatic residue mutations reveal direct correlation between HIV-1 nucleocapsid protein's nucleic acid chaperone activity and retroviral replication

The human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein plays an essential role in several stages of HIV-1 replication. One important function of HIV-1 NC is to act as a nucleic acid chaperone, in which the protein facilitates nucleic acid rearrangements important for reverse transcription and recombination. NC contains only 55 amino acids, with 15 basic residues and two zinc fingers, each having a single aromatic residue (Phe16 and Trp37). Despite its simple structure, HIV-1 NC appears to have optimal chaperone activity, including the ability to strongly aggregate nucleic acids, destabilize nucleic acid secondary structure, and facilitate rapid nucleic acid annealing. Here we combine single molecule DNA stretching experiments with ensemble solution studies of protein-nucleic acid binding affinity, oligonucleotide annealing, and nucleic acid aggregation to measure the characteristics of wild-type (WT) and aromatic residue mutants of HIV-1 NC that are important for nucleic acid chaperone activity. These in vitro results are compared to in vivo HIV-1 replication for viruses containing the same mutations. This work allows us to directly relate HIV-1 NC structure with its function as a nucleic acid chaperone in vitro and in vivo. We show that replacement of either aromatic residue with another aromatic residue results in a protein that strongly resembles WT NC. In contrast, single amino acid substitutions of either Phe16Ala or Trp37Ala significantly slow down NC's DNA interaction kinetics, while retaining some helix-destabilization capability. A double Phe16Ala/Trp37Ala substitution further reduces the latter activity. Surprisingly, the ensemble nucleic acid binding, annealing, and aggregation properties are not significantly altered for any mutant except the double aromatic substitution with Ala. Thus, elimination of a single aromatic residue from either zinc finger strongly reduces NC's chaperone activity as determined by single molecule DNA stretching experiments without significantly altering its ensemble-averaged biochemical properties. Importantly, the substitution of aromatic residues with Ala progressively decreases NC's nucleic acid chaperone activity while also progressively inhibiting viral replication. Taken together, these data support the critical role of HIV-1 NC's aromatic residues, and establish a direct and statistically significant correlation between nucleic acid chaperone activity and viral replication.

Optical tweezers to study viruses

A virus is a complex molecular machine that propagates by channeling its genetic information from cell to cell. Unlike macroscopic engines, it operates in a nanoscopic world under continuous thermal agitation. Viruses have developed efficient passive and active strategies to pack and release nucleic acids. Some aspects of the dynamic behavior of viruses and their substrates can be studied using structural and biochemical techniques. Recently, physical techniques have been applied to dynamic studies of viruses in which their intrinsic mechanical activity can be measured directly. Optical tweezers are a technology that can be used to measure the force, torque and strain produced by molecular motors, as a function of time and at the single-molecule level. Thanks to this technique, some bacteriophages are now known to be powerful nanomachines; they exert force in the piconewton range and their motors work in a highly coordinated fashion for packaging the viral nucleic acid genome. Nucleic acids, whose elasticity and condensation behavior are inherently coupled to the viral packaging mechanisms, are also amenable to examination with optical tweezers. In this chapter, we provide a comprehensive analysis of this laser-based tool, its combination with imaging methods and its application to the study of viruses and viral molecules.

Tension induces a base-paired overstretched DNA conformation

Mixed-sequence DNA molecules undergo mechanical overstretching by approximately 70% at 60-70 pN. Since its initial discovery 15 y ago, a debate has arisen as to whether the molecule adopts a new form [Cluzel P, et al. (1996) Science 271:792-794; Smith SB, Cui Y, Bustamante C (1996) Science 271:795-799], or simply denatures under tension [van Mameren J, et al. (2009) Proc Natl Acad Sci USA 106:18231-18236]. Here, we resolve this controversy by using optical tweezers to extend small 60-64 bp single DNA duplex molecules whose base content can be designed at will. We show that when AT content is high (70%), a force-induced denaturation of the DNA helix ensues at 62 pN that is accompanied by an extension of the molecule of approximately 70%. By contrast, GC-rich sequences (60% GC) are found to undergo a reversible overstretching transition into a distinct form that is characterized by a 51% extension and that remains base-paired. For the first time, results proving the existence of a stretched basepaired form of DNA can be presented. The extension observed in the reversible transition coincides with that produced on DNA by binding of bacterial RecA and human Rad51, pointing to its possible relevance in homologous recombination.

A single zinc finger optimizes the DNA interactions of the nucleocapsid protein of the yeast retrotransposon Ty3

Reverse transcription in retroviruses and retrotransposons requires nucleic acid chaperones, which drive the rearrangement of nucleic acid conformation. The nucleic acid chaperone properties of the human immunodeficiency virus type-1 (HIV-1) nucleocapsid (NC) protein have been extensively studied, and nucleic acid aggregation, duplex destabilization and rapid binding kinetics have been identified as major components of its activity. However, the properties of other nucleic acid chaperone proteins, such as retrotransposon Ty3 NC, a likely ancestor of HIV-1 NC, are not well understood. In addition, it is unclear whether a single zinc finger is sufficient to optimize the properties characteristic of HIV-1 NC. We used single-molecule DNA stretching as a method for detailed characterization of Ty3 NC chaperone activity. We found that wild type Ty3 NC aggregates single- and double-stranded DNA, weakly stabilizes dsDNA, and exhibits rapid binding kinetics. Single-molecule studies in the presence of Ty3 NC mutants show that the N-terminal basic residues and the unique zinc finger at the C-terminus are required for optimum chaperone activity in this system. While the single zinc finger is capable of optimizing Ty3 NC's DNA interaction kinetics, two zinc fingers may be necessary in order to facilitate the DNA destabilization exhibited by HIV-1 NC.

B-DNA to zip-DNA: simulating a DNA transition to a novel structure with enhanced charge-transport characteristics

The forced extension of a DNA segment is studied in a series of steered molecular dynamics simulations, employing a broad range of pulling forces. Throughout the entire force range, the formation of a zipper-like (zip-) DNA structure is observed. In that structure, first predicted by Lohikoski et al., the bases of the DNA strands interdigitate with each other and form a single-base aromatic stack. Similar motifs, albeit only a few base pairs in extent, have been observed in experimental crystal structures. Analysis of the dynamics of structural changes in pulled DNA shows that S-form DNA, thought to be adopted by DNA under applied force, serves as an intermediate between B-DNA and zip-DNA. Therefore, the phase transition plateau observed in force-extension curves of DNA is suggested to reflect the B-DNA to zip-DNA structural transition. Electronic structure analysis of purine bases in zip-DNA indicates a several-fold to order of magnitude increase in the π-π electronic coupling among nearest-neighbor nucleobases, compared to B-DNA. We further observe that zip-DNA does not require base pair complementarity between DNA strands, and we predict that the increased electronic coupling in zip-DNA will result in a much higher rate of charge transfer through an all-purine zip-DNA compared to B-DNA of equal length.

Anchoring nascent RNA to the DNA template could interfere with transcription

During normal transcription, the nascent RNA product is released from the DNA template. However, in some cases, the RNA remains bound or can become reattached to the template DNA duplex (for example, through R-loop formation). We have analyzed the effect on transcription elongation of nascent RNA anchoring to the template DNA duplex. Because the RNA polymerase follows a helical path along DNA duplex during transcription, the anchoring would result in wrapping the nascent RNA around the DNA in the region between the anchoring point and the translocating polymerase. This wrapping would cause an unfavorable loss of conformation entropy of the nascent RNA. It consequently would create an apparent force to unwrap the RNA by disrupting either the transcription complex or the anchoring structure. We have estimated that this force would be comparable to those required to melt nucleic acid duplexes or to arrest transcription elongation in single-molecule experiments. We predict that this force would create negative supercoiling in the DNA duplex region between the anchoring point and the transcribing RNA polymerase: this can promote the formation of unusual DNA structures and facilitate RNA invasion into the DNA duplex. Potential biological consequences of these effects are discussed.

Competition for hydrogen-bond formation in the helix-coil transition and protein folding

The problem of the helix-coil transition of biopolymers in explicit solvents, such as water, with the ability for hydrogen bonding with a solvent is addressed analytically using a suitably modified version of the Generalized Model of Polypeptide Chains. Besides the regular helix-coil transition, an additional coil-helix or reentrant transition is also found at lower temperatures. The reentrant transition arises due to competition between polymer-polymer and polymer-water hydrogen bonds. The balance between the two types of hydrogen bonding can be shifted to either direction through changes not only in temperature, but also by pressure, mechanical force, osmotic stress, or other external influences. Both polypeptides and polynucleotides are considered within a unified formalism. Our approach provides an explanation of the experimental difficulty of observing the reentrant transition with pressure and underscores the advantage of pulling experiments for studies of DNA. Results are discussed and compared with those reported in a number of recent publications with which a significant level of agreement is obtained.

Paired mutations abolish and restore the balanced annealing and melting activities of ORF1p that are required for LINE-1 retrotransposition

Retrotransposition amplifies LINE-1 (L1) to high copy number in mammalian genomes. The L1 protein encoded by ORF1 (ORF1p) is required for retrotransposition. This dependence on ORF1p was investigated by mutating three highly conserved residues, R238, R284 and Y318 to alanine, thereby inactivating retrotransposition. R284A and Y318A were rescued by further substituting the alanine with the appropriate conservative amino acid, e.g. lysine or phenylalanine, respectively, whereas R238K remained inactive. Quantification of the steady-state levels of L1 RNA and ORF1p failed to discriminate active from inactive variants, indicating loss of L1 retrotransposition resulted from loss of function rather than reduced expression. The two biochemical properties known for ORF1p are high-affinity RNA binding and nucleic acid chaperone activity. Only R238A/K exhibited significantly reduced RNA affinities. The nucleic acid chaperone activities of the remaining paired mutants were assessed by single-molecule DNA stretching and found to mirror retrotransposition activity. To further examine ORF1p chaperone function, their energetic barriers to DNA annealing and melting were derived from kinetic work. When plotted against each other, the ratio of these two activities distinguished functional from non-functional ORF1p variants. These findings enhance our understanding of the requirements for ORF1p in LINE-1 retrotransposition and, more generally, nucleic acid chaperone function.

Overstretching DNA at 65 pN does not require peeling from free ends or nicks

DNA exhibits a remarkable mechanical transition where its extension increases by 70% at 65 pN. Notwithstanding more than a decade of experimental and theoretical studies, there remains a significant debate on the nature of overstretched DNA. We developed a topologically closed but rotationally unconstrained DNA assay, which contains no nicks or free ends. DNA in this assay exhibited the canonical overstretching transition at 65 pN but without hysteresis upon retraction (v(stage) = 5 μm/s). Introduction of a controlled nick led to hysteresis in the force-extension curve. Moreover, the degree of hysteresis increased with the number of nicks. Hence, the generation of single-stranded DNA from free ends or nicks is not an obligatory step in overstretching DNA, but rather a consequence.

Measuring DNA-protein binding affinity on a single molecule using optical tweezers

DNA-protein interactions may be observed on single molecules with a variety of techniques. However, quantifying the binding affinity is difficult and often requires many (∼100) individual events to characterize the interaction. We use a single λ DNA molecule that provides a lattice of binding sites for proteins. Extending and relaxing the tethered molecule reversibly melts DNA, allowing it to be converted between double-stranded (ds) and single-stranded (ss) forms. By monitoring changes in the properties of the DNA as a function of added protein concentration and fitting to binding models, the DNA-protein interaction may be characterized and quantified. As an example, the high mobility group protein HMGB1(box A  +  B) is observed to stabilize dsDNA. Measuring the strength of this effect allows us to determine the equilibrium association constant for HMGB1(box A  +  B) binding to dsDNA.

Single-molecule stretching studies of RNA chaperones

RNA chaperone proteins play significant roles in diverse biological contexts. The most widely studied RNA chaperones are the retroviral nucleocapsid proteins (NC), also referred to as nucleic acid (NA) chaperones. Surprisingly, the biophysical properties of the NC proteins vary significantly for different viruses, and it appears that HIV-1 NC has optimal NA chaperone activity. In this review we discuss the physical nature of the NA chaperone activity of NC. We conclude that the optimal NA chaperone must saturate NA binding, leading to strong NA aggregation and slight destabilization of all NA duplexes. Finally, rapid kinetics of the chaperone protein interaction with NA is another primary component of its NA chaperone activity. We discuss these characteristics of HIV-1 NC and compare them with those of other NA binding proteins and ligands that exhibit only some characteristics of NA chaperone activity, as studied by single molecule DNA stretching.

Study of force induced melting of dsDNA as a function of length and conformation

We measure the constant force required to melt double-stranded (ds) DNA as a function of length for lengths from 12 to 100,000 base pairs, where the force is applied to the 3'3' or 5'5' ends of the dsDNA. Molecules with 32 base pairs or fewer melt before overstretching. For these short molecules, the melting force is independent of the ends to which the force is applied and the shear force as a function of length is well described by de Gennes theory with a de Gennes length of less than 10 bp. Molecules with lengths of 500 base pairs or more overstretch before melting. For these long molecules, the melting force depends on the ends to which the force is applied. The melting force as a function of length increases even when the length exceeds 1000 bp, where the length dependence is inconsistent with de Gennes theory. Finally, we expand de Gennes melting theory to 3'5' pulling and compare the predictions with experimental results.

The nature of the force-induced conformation transition of dsDNA studied by using single molecule force spectroscopy

Single-stranded DNA binding proteins (SSB) interact with single-stranded DNA (ssDNA) specifically. Taking advantage of this character, we have employed Bacillus subtilis SSB protein to investigate the nature of force-induced conformation transition of double-stranded DNA (dsDNA) by using AFM-based single molecule force spectroscopy (SMFS) technique. Our results show that, when a dsDNA is stretched beyond its contour length, the dsDNA is partially melted, producing some ssDNA segments which can be captured by SSB proteins. We have also systematically investigated the effects of stretching length, waiting time, and salt concentration on the conformation transition of dsDNA and SSB-ssDNA interactions, respectively. Furthermore, the effect of proflavine, a DNA intercalator, on the SSB-DNA interactions has been investigated, and the results indicate that the proflavine-saturated dsDNA can be stabilized to the extent that the dsDNA will no longer melt into ssDNA under the mechanical force even up to 150 pN, and no SSB-DNA interactions are detectable.

Biophysical characterization of DNA binding from single molecule force measurements

Single molecule force spectroscopy is a powerful method that uses the mechanical properties of DNA to explore DNA interactions. Here we describe how DNA stretching experiments quantitatively characterize the DNA binding of small molecules and proteins. Small molecules exhibit diverse DNA binding modes, including binding into the major and minor grooves and intercalation between base pairs of double-stranded DNA (dsDNA). Histones bind and package dsDNA, while other nuclear proteins such as high mobility group proteins bind to the backbone and bend dsDNA. Single-stranded DNA (ssDNA) binding proteins slide along dsDNA to locate and stabilize ssDNA during replication. Other proteins exhibit binding to both dsDNA and ssDNA. Nucleic acid chaperone proteins can switch rapidly between dsDNA and ssDNA binding modes, while DNA polymerases bind both forms of DNA with high affinity at distinct binding sites at the replication fork. Single molecule force measurements quantitatively characterize these DNA binding mechanisms, elucidating small molecule interactions and protein function.

Different pulling modes in DNA overstretching: a theoretical analysis

We study the thermally driven denaturation of a double-stranded polymer in the presence of a stretching force via Monte-Carlo simulations. When one strand only is stretched, the denaturation transition is first order, while when both strands are stretched, melting is second order. By revisiting the Poland-Scheraga model for DNA melting, we show that at room temperature, the most likely scenario is that DNA melts as it overstretches. Our results are in general agreement with the most recent experiments and suggest how varying temperature and stretching mode may help settle the question whether S-DNA exists or not.

C-terminal domain modulates the nucleic acid chaperone activity of human T-cell leukemia virus type 1 nucleocapsid protein via an electrostatic mechanism

Retroviral nucleocapsid (NC) proteins are molecular chaperones that facilitate nucleic acid (NA) remodeling events critical in viral replication processes such as reverse transcription. Surprisingly, the NC protein from human T-cell leukemia virus type 1 (HTLV-1) is an extremely poor NA chaperone. Using bulk and single molecule methods, we find that removal of the anionic C-terminal domain (CTD) of HTLV-1 NC results in a protein with chaperone properties comparable with that of other retroviral NCs. Increasing the ionic strength of the solution also improves the chaperone activity of full-length HTLV-1 NC. To determine how the CTD negatively modulates the chaperone activity of HTLV-1 NC, we quantified the thermodynamics and kinetics of wild-type and mutant HTLV-1 NC/NA interactions. The wild-type protein exhibits very slow dissociation kinetics, and removal of the CTD or mutations that eliminate acidic residues dramatically increase the protein/DNA interaction kinetics. Taken together, these results suggest that the anionic CTD interacts with the cationic N-terminal domain intramolecularly when HTLV-1 NC is not bound to nucleic acids, and similar interactions occur between neighboring molecules when NC is NA-bound. The intramolecular N-terminal domain-CTD attraction slows down the association of the HTLV-1 NC with NA, whereas the intermolecular interaction leads to multimerization of HTLV-1 NC on the NA. The latter inhibits both NA/NC aggregation and rapid protein dissociation from single-stranded DNA. These features make HTLV-1 NC a poor NA chaperone, despite its robust duplex destabilizing capability.