DNA Structure Design Is Improved Using an Artificially Expanded Alphabet of Base Pairs Including Loop and Mismatch Thermodynamic Parameters.

We show that in silico design of DNA secondary structures is improved by extending the base pairing alphabet beyond A-T and G-C to include the pair between 2-amino-8-(1'-ß-d-2'-deoxyribofuranosyl)-imidazo-[1,2-a]-1,3,5-triazin-(8H)-4-one and 6-amino-3-(1'-ß-d-2'-deoxyribofuranosyl)-5-nitro-(1H)-pyridin-2-one, abbreviated as P and Z. To obtain the thermodynamic parameters needed to include P-Z pairs in the designs, we performed 47 optical melting experiments and combined the results with previous work to fit free energy and enthalpy nearest neighbor folding parameters for P-Z pairs and G-Z wobble pairs. We find G-Z pairs have stability comparable to that of A-T pairs and should therefore be included as base pairs in structure prediction and design algorithms. Additionally, we extrapolated the set of loop, terminal mismatch, and dangling end parameters to include the P and Z nucleotides. These parameters were incorporated into the RNAstructure software package for secondary structure prediction and analysis. Using the RNAstructure Design program, we solved 99 of the 100 design problems posed by Eterna using the ACGT alphabet or supplementing it with P-Z pairs. Extending the alphabet reduced the propensity of sequences to fold into off-target structures, as evaluated by the normalized ensemble defect (NED). The NED values were improved relative to those from the Eterna example solutions in 91 of 99 cases in which Eterna-player solutions were provided. P-Z-containing designs had average NED values of 0.040, significantly below the 0.074 of standard-DNA-only designs, and inclusion of the P-Z pairs decreased the time needed to converge on a design. This work provides a sample pipeline for inclusion of any expanded alphabet nucleotides into prediction and design workflows.

DNA Structure Design Is Improved Using an Artificially Expanded Alphabet of Base Pairs Including Loop and Mismatch Thermodynamic Parameters.

We show that in silico design of DNA secondary structures is improved by extending the base pairing alphabet beyond A-T and G-C to include the pair between 2-amino-8-(1'-ß-D-2'-deoxyribofuranosyl)-imidazo-[1,2- a ]-1,3,5-triazin-(8 H )-4-one and 6-amino-3-(1'-ß-D-2'-deoxyribofuranosyl)-5-nitro-(1 H )-pyridin-2-one, simply P and Z. To obtain the thermodynamic parameters needed to include P-Z pairs in the designs, we performed 47 optical melting experiments and combined the results with previous work to fit a new set of free energy and enthalpy nearest neighbor folding parameters for P-Z pairs and G-Z wobble pairs. We find that G-Z pairs have stability comparable to A-T pairs and therefore should be considered quantitatively by structure prediction and design algorithms. Additionally, we extrapolated the set of loop, terminal mismatch, and dangling end parameters to include P and Z nucleotides. These parameters were incorporated into the RNAstructure software package for secondary structure prediction and analysis. Using the RNAstructure Design program, we solved 99 of the 100 design problems posed by Eterna using the ACGT alphabet or supplementing with P-Z pairs. Extending the alphabet reduced the propensity of sequences to fold into off-target structures, as evaluated by the normalized ensemble defect (NED). The NED values were improved relative to those from the Eterna example solutions in 91 of 99 cases where Eterna-player solutions were provided. P-Z-containing designs had average NED values of 0.040, significantly below the 0.074 of standard-DNA-only designs, and inclusion of the P-Z pairs decreased the time needed to converge on a design. This work provides a sample pipeline for inclusion of any expanded alphabet nucleotides into prediction and design workflows.

Stability of the pH-Dependent Parallel-Stranded d(CGA) Motif.

Noncanonical DNA structures that retain programmability and structural predictability are increasingly being used in DNA nanotechnology applications, in which they offer versatility beyond traditional Watson-Crick interactions. The d(CGA) triplet repeat motif is structurally dynamic and can transition between parallel-stranded homo-base paired duplex and antiparallel unimolecular hairpin in a pH-dependent manner. Here, we evaluate the thermodynamic stability and nuclease sensitivity of oligonucleotides composed of the d(CGA) motif and several structurally related sequence variants. These results show that the structural transition resulting from decreasing the pH is accompanied by both a significant energetic stabilization and decreased nuclease sensitivity as unimolecular hairpin structures are converted to parallel-stranded homo-base paired duplexes. Furthermore, the stability of the parallel-stranded duplex form can be altered by changing the 5'-nucleobase of the d(CGA) triplet and the frequency and position of the altered triplets within long stretches of d(CGA) triplets. This work offers insight into the stability and versatility of the d(CGA) triplet repeat motif and provides constraints for using this pH-adaptive structural motif for creating DNA-based nanomaterials.

Nucleic acid-cleaving catalytic DNA for sensing and therapeutics.

DNAzymes with nucleic acid-cleaving catalytic activity are increasing in versatility through concerted efforts to discover new sequences with unique functions, and they are generating excitement in the sensing community as cheap, stable, amplifiable detection elements. This review provides a comprehensive list and detailed descriptions of the DNAzymes identified to date, classified by their associated small molecule or ion needed for catalysis; of note, this classification clarifies conserved regions of various DNAzymes that are not obvious in the literature. Furthermore, we detail the breadth of functionality of these DNA sequences as well as the range of reaction conditions under which they are useful. In addition, the utility of the DNAzymes in a variety of sensing and therapeutic applications is presented, detailing both their advantages and disadvantages.

Biophysics of Artificially Expanded Genetic Information Systems. Thermodynamics of DNA Duplexes Containing Matches and Mismatches Involving 2-Amino-3-nitropyridin-6-one (Z) and Imidazo[1,2-a]-1,3,5-triazin-4(8H)one (P).

Synthetic nucleobases presenting non-Watson-Crick arrangements of hydrogen bond donor and acceptor groups can form additional nucleotide pairs that stabilize duplex DNA independent of the standard A:T and G:C pairs. The pair between 2-amino-3-nitropyridin-6-one 2'-deoxyriboside (presenting a {donor-donor-acceptor} hydrogen bonding pattern on the Watson-Crick face of the small component, trivially designated Z) and imidazo[1,2-a]-1,3,5-triazin-4(8H)one 2'-deoxyriboside (presenting an {acceptor-acceptor-donor} hydrogen bonding pattern on the large component, trivially designated P) is one of these extra pairs for which a substantial amount of molecular biology has been developed. Here, we report the results of UV absorbance melting measurements and determine the energetics of binding of DNA strands containing Z and P to give short duplexes containing Z:P pairs as well as various mismatches comprising Z and P. All measurements were done at 1 M NaCl in buffer (10 mM Na cacodylate, 0.5 mM EDTA, pH 7.0). Thermodynamic parameters (ΔH°, ΔS°, and ΔG°37) for oligonucleotide hybridization were extracted. Consistent with the Watson-Crick model that considers both geometric and hydrogen bonding complementarity, the Z:P pair was found to contribute more to duplex stability than any mismatches involving either nonstandard nucleotide. Further, the Z:P pair is more stable than a C:G pair. The Z:G pair was found to be the most stable mismatch, forming either a deprotonated mismatched pair or a wobble base pair analogous to the stable T:G mismatch. The C:P pair is less stable, perhaps analogous to the wobble pair observed for C:O6-methyl-G, in which the pyrimidine is displaced into the minor groove. The Z:A and T:P mismatches are much less stable. Parameters for predicting the thermodynamics of oligonucleotides containing Z and P bases are provided. This represents the first case where this has been done for a synthetic genetic system.

Applying Thymine Isostere 2,4-Difluoro-5-Methylbenzene as a NMR Assignment Tool and Probe of Homopyrimidine/Homopurine Tract Structural Dynamics.

Proton assignment of nuclear magnetic resonance (NMR) spectra of homopyrimidine/homopurine tract oligonucleotides becomes extremely challenging with increasing helical length due to severe cross-peak overlap. As an alternative to the more standard practice of (15)N and (13)C labeling of oligonucleotides, here, we describe a method for assignment of highly redundant DNA sequences that uses single-site substitution of the thymine isostere 2,4-difluoro-5-methylbenzene (dF). The impact of this approach in facilitating the assignment of intractable spectra and analyzing oligonucleotide structure and dynamics is demonstrated using A-tract and TATA box DNA and two polypurine tract-containing RNA:DNA hybrids derived from HIV-1 and the Saccharomyces cerevisiae long-terminal repeat-containing retrotransposon Ty3. Only resonances proximal to the site of dF substitution exhibit sizable chemical shift changes, providing spectral dispersion while still allowing chemical shift mapping of resonances from unaffected residues distal to the site of modification directly back to the unmodified sequence. It is further illustrated that dF incorporation can subtly alter the conformation and dynamics of homopyrimidine/homopurine tract oligonucleotides, and how these NMR observations can be correlated, in the cases of the TATA box DNA, with modulation in the TATA box-binding protein interaction using an orthogonal gel assay.

Direct observation of dynamic mechanical regulation of DNA condensation by environmental stimuli.

Gene delivery is a promising way to treat hereditary diseases and cancer; however, there is little understanding of DNA:carrier complex mechanical properties, which may be critical for the protection and release of nucleic acids. We applied optical tweezers to directly measure single-molecule mechanical properties of DNA condensed using 19-mer poly-L-lysine (PLL) or branched histidine-lysine (HK) peptides. Force-extension profiles indicate that both carriers condense DNA actively, showing force plateaus during stretching and relaxation cycles. As the environment such as carrier concentration, pH, and the presence of zinc ions changes, DNA:HK complexes showed dynamically regulated mechanical properties at multiple force levels. The fundamental knowledge from this study can be applied to design a mechanically tailored complex which may enhance transfection efficiency by controlling the stability of the complex temporally and spatially.

DNA, flexibly flexible.

Investigators have constructed dsDNA molecules with several different base modifications and have characterized their bending and twisting flexibilities using atomic force microscopy, DNA ring closure, and single-molecule force spectroscopy with optical tweezers. The three methods provide persistence length measurements that agree semiquantitatively, and they show that the persistence length is surprisingly similar for all of the modified DNAs. The circular dichroism spectra of modified DNAs differ substantially. Simple explanations based on base stacking strength, polymer charge, or groove occupancy by functional groups cannot explain the results, which will guide further high-resolution theory and experiments.

Tethered particle motion reveals that LacI·DNA loops coexist with a competitor-resistant but apparently unlooped conformation.

The lac repressor protein (LacI) efficiently represses transcription of the lac operon in Escherichia coli by binding to two distant operator sites on the bacterial DNA and causing the intervening DNA to form a loop. We employed single-molecule tethered particle motion to observe LacI-mediated loop formation and breakdown in DNA constructs that incorporate optimized operator binding sites and intrinsic curvature favorable to loop formation. Previous bulk competition assays indirectly measured the loop lifetimes in these optimized DNA constructs as being on the order of days; however, we measured these same lifetimes to be on the order of minutes for both looped and unlooped states. In a range of single-molecule DNA competition experiments, we found that the resistance of the LacI-DNA complex to competitive binding is a function of both the operator strength and the interoperator sequence. To explain these findings, we present what we believe to be a new kinetic model of loop formation and DNA competition. In this proposed new model, we hypothesize a new unlooped state in which the unbound DNA-binding domain of the LacI protein interacts nonspecifically with nonoperator DNA adjacent to the operator site at which the second LacI DNA-binding domain is bound.

Enhanced silencing and stabilization of siRNA polyplexes by histidine-mediated hydrogen bonds.

Branched peptides containing histidines and lysines (HK) have been shown to be effective carriers for DNA and siRNA. We anticipate that elucidation of the binding mechanism of HK with siRNA will provide greater insight into the self-assembly and delivery of the HK:siRNA polyplex. Non-covalent bonds between histidine residues and nucleic acids may enhance the stability of siRNA polyplexes. We first compared the polyplex biophysical properties of a branched HK with those of branched asparagine-lysine peptide (NK). Consistent with siRNA silencing experiments, gel electrophoresis demonstrated that the HK siRNA polyplex maintained its integrity with prolonged incubation in serum, whereas siRNA in complex with NK was degraded in a time-dependent manner. Isothermal titration calorimetry of various peptides binding to siRNA at pH 7.3 showed that branched polylysine, interacted with siRNA was initially endothermic, whereas branched HK exhibited an exothermic reaction at initial binding. The exothermic interaction indicates formation of non-ionic bonds between histidines and siRNA; purely electrostatic interaction is entropy-driven and endothermic. To investigate the type of non-ionic bond, we studied the protonation state of imidazole rings of a selectively (15)N labeled branched HK by heteronuclear single quantum coherence NMR. The peak of Nδ1-H tautomers of imidazole shifted downfield (in the direction of deprotonation) by 0.5-1.0 ppm with addition of siRNA, providing direct evidence that histidines formed hydrogen bonds with siRNA at physiological pH. These results establish that histidine-rich peptides form hydrogen bonds with siRNA, thereby enhancing the stability and biological activity of the polyplex in vitro and in vivo.

Rationally designed coiled-coil DNA looping peptides control DNA topology.

Artificial DNA looping peptides were engineered to study the roles of protein and DNA flexibility in controlling the geometry and stability of protein-mediated DNA loops. These LZD (leucine zipper dual-binding) peptides were derived by fusing a second, C-terminal, DNA-binding region onto the GCN4 bZip peptide. Two variants with different coiled-coil lengths were designed to control the relative orientations of DNA bound at each end. Electrophoretic mobility shift assays verified formation of a sandwich complex containing two DNAs and one peptide. Ring closure experiments demonstrated that looping requires a DNA-binding site separation of 310 bp, much longer than the length needed for natural loops. Systematic variation of binding site separation over a series of 10 constructs that cyclize to form 862-bp minicircles yielded positive and negative topoisomers because of two possible writhed geometries. Periodic variation in topoisomer abundance could be modeled using canonical DNA persistence length and torsional modulus values. The results confirm that the LZD peptides are stiffer than natural DNA looping proteins, and they suggest that formation of short DNA loops requires protein flexibility, not unusual DNA bendability. Small, stable, tunable looping peptides may be useful as synthetic transcriptional regulators or components of protein-DNA nanostructures.

LacI-DNA-IPTG loops: equilibria among conformations by single-molecule FRET.

The E. coli Lac repressor (LacI) tetramer binds simultaneously to a promoter-proximal DNA binding site (operator) and an auxiliary operator, resulting in a DNA loop, which increases repression efficiency. Induction of the lac operon by allolactose reduces the affinity of LacI for DNA, but induction does not completely prevent looping in vivo. Our previous work on the conformations of LacI loops used a hyperstable model DNA construct, 9C14, that contains a sequence directed bend flanked by operators. Single-molecule fluorescence resonance energy transfer (SM-FRET) on a dual fluorophore-labeled LacI-9C14 loop showed that it adopts a single, stable, high-FRET V-shaped LacI conformation. Ligand-induced changes in loop geometry can affect loop stability, and the current work assesses loop population distributions for LacI-9C14 complexes containing the synthetic inducer IPTG. SM-FRET confirms that the high-FRET LacI-9C14 loop is only partially destabilized by saturating IPTG. LacI titration experiments and FRET fluctuation analysis suggest that the addition of IPTG induces loop conformational dynamics and re-equilibration between loop population distributions that include a mixture of looped states that do not exhibit high-efficiency FRET. The results show that repression by looping even at saturating IPTG should be considered in models for regulation of the operon. We propose that persistent DNA loops near the operator function biologically to accelerate rerepression upon exhaustion of inducer.

Surface-modified HK:siRNA nanoplexes with enhanced pharmacokinetics and tumor growth inhibition.

We characterized in this study the pharmacokinetics and antitumor efficacy of histidine-lysine (HK):siRNA nanoplexes modified with PEG and a cyclic RGD (cRGD) ligand targeting αvß3 and αvß5 integrins. With noninvasive imaging, systemically administered surface-modified HK:siRNA nanoplexes showed nearly 4-fold greater blood levels, 40% higher accumulation in tumor tissue, and 60% lower luciferase activity than unmodified HK:siRNA nanoplexes. We then determined whether the surface-modified HK:siRNA nanoplex carrier was more effective in reducing MDA-MB-435 tumor growth with an siRNA targeting Raf-1. Repeated systemic administration of the selected surface modified HK:siRNA nanoplexes targeting Raf-1 showed 35% greater inhibition of tumor growth than unmodified HK:siRNA nanoplexes and 60% greater inhibition of tumor growth than untreated mice. The improved blood pharmacokinetic results and tumor localization observed with the integrin-targeting surface modification of HK:siRNA nanoplexes correlated with greater tumor growth inhibition. This investigation reveals that through control of targeting ligand surface display in association with a steric PEG layer, modified HK: siRNA nanoplexes show promise to advance RNAi therapeutics in oncology and potentially other critical diseases.

FRET studies of a landscape of Lac repressor-mediated DNA loops.

DNA looping mediated by the Lac repressor is an archetypal test case for modeling protein and DNA flexibility. Understanding looping is fundamental to quantitative descriptions of gene expression. Systematic analysis of LacI•DNA looping was carried out using a landscape of DNA constructs with lac operators bracketing an A-tract bend, produced by varying helical phasings between operators and the bend. Fluorophores positioned on either side of both operators allowed direct Förster resonance energy transfer (FRET) detection of parallel (P1) and antiparallel (A1, A2) DNA looping topologies anchored by V-shaped LacI. Combining fluorophore position variant landscapes allows calculation of the P1, A1 and A2 populations from FRET efficiencies and also reveals extended low-FRET loops proposed to form via LacI opening. The addition of isopropyl-ß-D-thio-galactoside (IPTG) destabilizes but does not eliminate the loops, and IPTG does not redistribute loops among high-FRET topologies. In some cases, subsequent addition of excess LacI does not reduce FRET further, suggesting that IPTG stabilizes extended or other low-FRET loops. The data align well with rod mechanics models for the energetics of DNA looping topologies. At the peaks of the predicted energy landscape for V-shaped loops, the proposed extended loops are more stable and are observed instead, showing that future models must consider protein flexibility.

Quantitative methods for measuring DNA flexibility in vitro and in vivo.

The double-helical DNA biopolymer is particularly resistant to bending and twisting deformations. This property has important implications for DNA folding in vitro and for the packaging and function of DNA in living cells. Among the outstanding questions in the field of DNA biophysics are the underlying origin of DNA stiffness and the mechanisms by which DNA stiffness is overcome within cells. Exploring these questions requires experimental methods to quantitatively measure DNA bending and twisting stiffness both in vitro and in vivo. Here, we discuss two classical approaches: T4 DNA ligase-mediated DNA cyclization kinetics and lac repressor-mediated DNA looping in Escherichia coli. We review the theoretical basis for these techniques and how each can be applied to quantitate biophysical parameters that describe the DNA polymer. We then show how we have modified these methods and applied them to quantitate how apparent DNA physical properties are altered in vitro and in vivo by sequence-nonspecific architectural DNA-binding proteins such as the E. coli HU protein and eukaryotic HMGB proteins.

Gene repression by minimal lac loops in vivo.

The inflexibility of double-stranded DNA with respect to bending and twisting is well established in vitro. Understanding apparent DNA physical properties in vivo is a greater challenge. Here, we exploit repression looping with components of the Escherichia coli lac operon to monitor DNA flexibility in living cells. We create a minimal system for testing the shortest possible DNA repression loops that contain an E. coli promoter, and compare the results to prior experiments. Our data reveal that loop-independent repression occurs for certain tight operator/promoter spacings. When only loop-dependent repression is considered, fits to a thermodynamic model show that DNA twisting limits looping in vivo, although the apparent DNA twist flexibility is 2- to 4-fold higher than in vitro. In contrast, length-dependent resistance to DNA bending is not observed in these experiments, even for the shortest loops constraining <0.4 persistence lengths of DNA. As observed previously for other looping configurations, loss of the nucleoid protein heat unstable (HU) markedly disables DNA looping in vivo. Length-independent DNA bending energy may reflect the activities of architectural proteins and the structure of the DNA topological domain. We suggest that the shortest loops are formed in apical loops rather than along the DNA plectonemic superhelix.

Computational analysis of looping of a large family of highly bent DNA by LacI.

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.

Eukaryotic HMGB proteins as replacements for HU in E. coli repression loop formation.

DNA looping is important for gene repression and activation in Escherichia coli and is necessary for some kinds of gene regulation and recombination in eukaryotes. We are interested in sequence-nonspecific architectural DNA-binding proteins that alter the apparent flexibility of DNA by producing transient bends or kinks in DNA. The bacterial heat unstable (HU) and eukaryotic high-mobility group B (HMGB) proteins fall into this category. We have exploited a sensitive genetic assay of DNA looping in living E. coli cells to explore the extent to which HMGB proteins and derivatives can complement a DNA looping defect in E. coli lacking HU protein. Here, we show that derivatives of the yeast HMGB protein Nhp6A rescue DNA looping in E. coli lacking HU, in some cases facilitating looping to a greater extent than is observed in E. coli expressing normal levels of HU protein. Nhp6A-induced changes in the DNA length-dependence of repression efficiency suggest that Nhp6A alters DNA twist in vivo. In contrast, human HMGB2-box A derivatives did not rescue looping.

Effects of nucleoid proteins on DNA repression loop formation in Escherichia coli.

The intrinsic stiffness of DNA limits its ability to be bent and twisted over short lengths, but such deformations are required for gene regulation. One classic paradigm is DNA looping in the regulation of the Escherichia coli lac operon. Lac repressor protein binds simultaneously to two operator sequences flanking the lac promoter. Analysis of the length dependence of looping-dependent repression of the lac operon provides insight into DNA deformation energetics within cells. The apparent flexibility of DNA is greater in vivo than in vitro, possibly because of host proteins that bind DNA and induce sites of flexure. Here we test DNA looping in bacterial strains lacking the nucleoid proteins HU, IHF or H-NS. We confirm that deletion of HU inhibits looping and that quantitative modeling suggests residual looping in the induced operon. Deletion of IHF has little effect. Remarkably, DNA looping is strongly enhanced in the absence of H-NS, and an explanatory model is proposed. Chloroquine titration, psoralen crosslinking and supercoiling-sensitive reporter assays show that the effects of nucleoid proteins on looping are not correlated with their effects on either total or unrestrained supercoiling. These results suggest that host nucleoid proteins can directly facilitate or inhibit DNA looping in bacteria.