Ab-inito Quantum Mechanical Calculations of NMR Chemical Shifts in Nucleic Acids Constituents II Conformational Dependence of the lH and 13C Chemical Shifts in the Ribose

Abstract

The magnetic shielding constant of the different ¹³C and ¹³H nuclei of a deoxyribose are calculated for the C2′ endo and C3′ endo puckerings of the furanose ring as a function of the conformation about the C4′C5′ bond. For the carbons the calculated variations are of several ppm, the C3′ endo puckering corresponding in most cases to a larger shielding than the C2′ endo one. For the protons the calculated variations of chemical shifts are all smaller than 1.3 ppm, that is of the order of magnitude of the variation of the geometrical shielding produced on these protons by the other units of a DNA double helix, with a change of the overall structure of the helix. The computations carried out on the deoxyribose ?3′ and 5′ phosphates for several conformations of the phosphate group tend to show that the changes of conformation of the charged group of atoms produce chemical shift variations smaller than the two conformational parameters of the deoxyribose itself. The calculations carried out for a ribose do give the general features of the differences between the carbon and proton spectra of deoxynucleosides and nucleosides.

The comparison of the measured and calculated phosphorylation shifts tend to show that the counterion contributes significantly, for some nuclei of the deoxyribose, to the shifts measured. The calculated magnitude of this polarization effect on carbon shifts suggests a tentative qualitative interpretation of carbon spectra of the ribose part of DNA double helices.     

Flexibility of the furanose ring in nucleic acids: a compilation of crystal data

A compilation of crystal structure data on deoxyribo- and ribonucleosides and their higher derivatives is presented. The aim of this paper is to highlight the flexibility of deoxyribose and ribose rings. So far, the conformational parameters of nucleic acids constituents of ribose and deoxyribose have not been analysed separately. This paper aims to correlate the conformational parameters with the nature and puckering of the sugar. Deoxyribose puckering occurs in the C2′ endo region while ribose puckering is observed both in the C3′ endo and C2′ endo regions. A few endocyclic and exocyclic bond angles depend on the puckering and the nature of the sugar. The majority of structures have an anti conformation about the glycosyl bond. There appears to be a puckering dependence on the torsion angle about the C4′C5′ bonds. Such stereochemical information is useful in model building studies of polynucleotides and nucleic acids.     

The occurrence of the syn-C3′ endo conformation and the distorted backbone conformations for C4′-C5′ and P-O5′ in oligo and polynucleotides

Abstract

The nucleoside constituents of nucleic acids prefer the anti conformation (1). When the sugar pucker is taken into account the nucleosides prefer the C2′endo-anti conformation. Of the nearly 300 nucleosides known, about 250 are in the anti conformation and 50 are in the syn-conformation, i.e., anti to syn conformation is 5:1. The nucleotide building blocks of nucleic acids show the same trend as nucleosides. Both the deoxy-guanosine and ribo- guanosine residues in nucleosides and nucleotides prefer the syn-C2′endo conformation with an intra-molecular hydrogen bond (for nucleosides) between the O5′- H and the N3 of the base and, a few syn-C3′endo conformations are also observed. Evidence is presented for the occurrence of the C3′endo-syn conformation for guanines in mis-paired double helical right-handed structures with the distorted sugar phosphate C4′-C5′ and P-O5′ bonds respectively, from g+ (gg) and g- to trans. Evidence is also provided for guanosine nucleotides in left-handed double-helical (Z-DNA) oligo and polynucleotides which has the same syn-C3′endo conformation and the distorted backbone sugar-phosphate bonds (C4′-C5′ and P- O5′) as in the earlier right-handed case.     

Why ribose was selected as the sugar component of nucleic acids

During evolution ribose was selected as the exclusive sugar component of nucleic acids. The selection is explained by using molecular models and by eliminating most of the other common sugars by looking at their chemical structure and envisioning how they would fit in a nucleic acid model. Comparisons of sugar pucker conformations and configurations of pentoses indicate that ribose was not randomly selected but the only choice, since beta-D-ribose fits best into the structure of physiological forms of nucleic acids. In other nucleotides containing arabinose, xylose, or lyxose, the C(2)'-OH and/or the C(3)'-OH are above the furanose ring, causing steric interference with the bulky base and the C(5)'-OH group.     

Studies of the B‐Z transition of DNA: The temperature dependence of the free‐energy difference,the composition of the counterion sheath in mixed salt,and the preparation of a sample of the 5′‐d[T‐(m5C‐G)12‐T] duplex in pure B‐DNA or Z‐DNA form

It is often envisioned that cations might coordinate at specific sites of nucleic acids and play an important structural role, for instance in the transition between B‐DNA and Z‐DNA. However, nucleic acid models explicitly devoid of specific sites may also exhibit features previously considered as evidence for specific binding. Such is the case of the “composite cylinder” (or CC) model which spreads out localized features of DNA structure and charge by cylindrical averaging, while sustaining the main difference between the B and Z structures, namely the better immersion of the B‐DNA phosphodiester charges in the solution. Here, we analyze the non‐electrostatic component of the free‐energy difference between B‐DNA and Z‐DNA. We also compute the composition of the counterion sheath in a wide range of mixed‐salt solutions and of temperatures: in contrast with the large difference of composition between the B‐DNA and Z‐DNA forms, the temperature dependence of sheath composition, previously unknown, is very weak. In order to validate the model, the mixed‐salt predictions should be compared to experiment. We design a procedure for future measurements of the sheath composition based on Anomalous Small‐Angle X‐ray Scattering and complemented by ³¹P NMR. With due consideration for the kinetics of the B‐Z transition and for the capacity of generating at will the B or Z form in a single sample, the 5′‐d[T‐(m⁵C‐G)₁₂‐T] 26‐mer emerges as a most suitable oligonucleotide for this study. Finally, the application of the finite element method to the resolution of the Poisson‐Boltzmann equation is described in detail. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 369–384, 2016.     

Analysis of the possible helical structures of nucleic acids and polynucleotides. Application of (n-h) plots.

The two helical parameters n and h where n is the number of nucleotide residues per turn and h is the height per nucleotide residue have been evaluated for single stranded helical polynucleotide chains comprising C(3') -endo and C(2') endo class of nucleotides. The helical parameters are found to be especially sensitive to the C(4')-C(3') (sugar pucker) and the C(4')-C(5') torsions. The (n-h) plots display only one important helix forming domain for each class of nucleotides characterized by the sugar pucker and the C(4')-C(5') torsion. A correlation between the (n-h) plots and the known RNA (A,A') and DNA (A,B,C) helical forms has been established. It is found that all forms of helices except the C-DNA possess a favorable combination of P-O torsions. The analysis of the (n-h) plots suggests that C-DNA can have a conformation very similar to B-DNA. Although the (n-h) plots predict the stereochemical possibility of both right-handed and left-handed helices, nucleic acids apparently prefer right-handed conformation because of the energetics associated with the sugar-phosphate backbone and the base.     

Oligonucleotide N3′→P5′ Phosphoramidates and Thio‐Phoshoramidates as Potential Therapeutic Agents

Nucleic acids analogues, i.e., oligonucleotide N3′→P5′ phosphoramidates and N3′→P5′ thio‐phosphoramidates, containing 3′‐amino‐3′‐deoxy nucleosides with various 2′‐substituents were synthesized and extensively studied. These compounds resist nuclease hydrolysis and form stable duplexes with complementary native phosphodiester DNA and, particularly, RNA strands. An increase in duplexes' melting temperature, ΔT_m, relative to their phosphodiester counterparts, reaches 2.2–4.0° per modified nucleoside. 2′‐OH‐ (RNA‐like), 2′‐O‐Me‐, and 2′‐ribo‐F‐nucleoside substitutions result in the highest degree of duplex stabilization. Moreover, under close to physiological salt and pH conditions, the 2′‐deoxy‐ and 2′‐fluoro‐phosphoramidate compounds form extremely stable triple‐stranded complexes with either single‐ or double‐stranded phosphodiester DNA oligonucleotides. Melting temperature, T_m, of these triplexes exceeds T_m values for the isosequential phosphodiester counterparts by up to 35°. 2′‐Deoxy‐N3′→P5′ phosphoramidates adopt RNA‐like C3′‐endo or N‐type nucleoside sugar‐ring conformations and hence can be used as stable RNA mimetics. Duplexes formed by 2′‐deoxy phosphoramidates with complementary RNA strands are not substrates for RNase H‐mediated cleavage in vitro. Oligonucleotide phosphoramidates and especially thio‐phosphoramidates conjugated with lipid groups are cell‐permeable and demonstrate high biological target specific activity in vitro. In vivo, these compounds show good bioavailability and efficient biodistribution to all major organs, while exerting acceptable toxicity at therapeutically relevant doses. Short oligonucleotide N3′→P5′ thio‐phosphoramidate conjugated to 5′‐palmitoyl group, designated as GRN163L (Imetelstat), was recently introduced as a potent human telomerase inhibitor. GRN163L is not an antisense agent; it is a direct competitive inhibitor of human telomerase, which directly binds to the active site of the enzyme and thus inhibits its activity. This compound is currently in multiple Phase‐I and Phase‐I/II clinical trials as potential broad‐spectrum anticancer agent.     

Conformational Characteristics and Correlations in Crystal Structures of Nucleic Acid Oligonucleotides: Evidence for Sub-states

Abstract

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible ?, ζ, α, β and γ backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain α/γ changes are accompanied by C3′ endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual ?/ζ combinations occur with C2′ (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all con- formational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

47 Incorporation of alternative nucleic acid chemistries into G-quadruplex structure

Nucleic acids that form G-quadruplex (G4) structure have found applications in a host of research and technology regimes. Numerous G4 based aptamer drugs have been identified with pharmacological activity against cancer, HIV, prions, and blood coagulation (1). In the field of nanotechnology, G4 based sensors and nano-machines have also received much attention. The ability to synthesize nucleic acid ex-vivo allows for the site-specific incorporation of non-natural chemistries into nucleic acids that can be used to tune their physical and pharmacological properties. We summarize the results of a series of studies investigating the effective incorporation of alternative nucleic acid chemistries into G4 DNA. These modified chemistries include C8-modified guanine bases, as well as 2′-F, 2′-F-ANA, and Locked nucleic acid (LNA) modifications to the ribose sugar. We report primarily on the effect of these modifications on G-quadruplex folding topology, thermal stability, and structure. The substitution of LNA-guanosine into the core guanine tetrads disrupts structure in specific structural environments. On the other hand, 2′-F- and 2′-F-ANA guanosine can generally be incorporated without disrupting the structure when substituted into guanine bases in certain structural conformations. We find that 2′-F-ANA-guanosine and 2′-F-guanosine are powerful tools for controling the conformation of G4 structures (2). Functionalization at the C8 of the guanine base stabilizes in a manner dependent on the glycosidic conformation of the base, with different modification chemistries stabilizing to varying extents (3). The results of these studies provide useful insight on how to effectively incorporate some useful chemical tools from the growing toolbox of modified nucleic acid chemistries into G-quadruplex nucleic acid.     

Ab-inito quantum mechanical calculations of NMR chemical shifts in nucleic acids constituents. II. Conformational dependence of the 1H and 13C chemical shifts in the ribose

The magnetic shielding constant of the different 13C and 1H nuclei of a deoxyribose are calculated for the C2' endo and C3' endo puckerings of the furanose ring as a function of the conformation about the C4'C5' bond. For the carbons the calculated variations are of several ppm, the C3' endo puckering corresponding in most cases to a larger shielding than the C2' endo one. For the protons the calculated variations of chemical shifts are all smaller than 1.3 ppm, that is of the order of magnitude of the variation of the geometrical shielding produced on these protons by the other units of a DNA double helix, with a change of the overall structure of the helix. The computations carried out on the deoxyribose-3' and 5' phosphates for several conformations of the phosphate group tend to show that the changes of conformation of the charged group of atoms produce chemical shift variations smaller than the two conformational parameters of the deoxyribose itself. The calculations carried out for a ribose do give the general features of the differences between the carbon and proton spectra of deoxynucleosides and nucleosides. The comparison of the measured and calculated phosphorylation shifts tend to show that the counterion contributes significantly, for some nuclei of the deoxyribose, to the shifts measured. The calculated magnitude of this polarization effect on carbon shifts suggests a tentative qualitative interpretation of carbon spectra of the ribose part of DNA double helices.     

On the conformational dependence of the proton chemical shifts in nucleosides and nucleotides III. Proton chemical shifts of 5′-nucleotides as a function of different conformational parameters

The variation of the chemical shift of the protons of 5′-UMP and 5′-AMP is calculated as a function of χ_CN, ψ and ? torsion angles. The shift of H8 of 5′-AMP and H6 of 5′-UMP is found to be very sensitive to the value of χ_CN. For the anti conformations the shift of these protons is more sensitive to the value of the rotation about CS′-05′ than about C4′-CS′. For the protons of the ribose the calculations show that for the C2′-endo pucker H3′ and H2′ undergo the largest chemical shift variations when ? and ψ vary. The calculated variations are considered in relation with the role of the conformation of the nucleotides in the chemical shift variation between mono and polynucleotides and between the different helical structures of polynucleotides.     

Stereochemistry of 2′-5′ Linked Nucleic Acids: Crystal and Molecular Structure of Ammonium Adenylyl-2′, 5′-Adenosine Tetrahydrate: A Core Fragment of 2′-5′ Oligo A′s Produced by Interferon Induced Adenylate Synthetased

Abstract

The preponderance of 3′-5′ phosphodiester links in nucleic acids is well known. Albeit less prevalent, the 2′-5′ links are specifically utilised in the formation of ‘lariat’ in group II introns and in the msDNA-RNA junction in myxobacterium. As a sequel to our earlier study on cytidylyl-2′,5′-adenosine we have now obtained the crystal structure of adenylyl-2′,5′-adenosine (A2′p5′A) at atomic resolution. This dinucleoside monophosphate crystallises in the orthorhombic space group P2₁2₁2₁ with a= 7.956(3)Å, b = 12.212(3)Å and c = 36.654 (3) Å. CuKα intensity data were collected on a diffractometer. The structure was sloved by direct methods and refined by full matrix least squares methods to R = 10.8 %. The 2′ terminal adenine is in the commonly observed anti (χ₂ =?161°) conformation and the 5′ terminal base has a syn (χ₁ = 55°) conformation more often seen in purine nucleotides. A noteworthy feature of A2′p5′ A is the intranucleotide hydrogen bond between N3 and 05′ atoms of the 5′ adenine base. The two furanose rings in A2′ p5′ A show different conformations-C2′ endo, C3′ endo puckering for the 5′ and 2′ ends respectively. In this structure too there is a stacking of the purine base on the ribose 04′ just as in other 2′-5′ dinucleoside structures, a feature characteristically seen in the left handed ZDNA. In having syn, anti conformation about the glycosyl bonds, C2′ endo, C3′ endo mixed sugar puckering and N3–05′ intramolecular hydrogen bond A2′p5′ A resembles its 3′-5′ analogue and several other 2′-5′ dinucleoside monophosphate structures solved so far. Striking similarities between the 2′-5′ dinucleoside monophosphate structures suggest that the conformation of the 5′-end nucleoside dictates the conformation of the 2′ end nucleoside. Also, the 2′-5′ dimers do not favour formation of miniature classical double helical structures like the 3′-5′ dimers. It is conceivable, 2–5(A) could be using the stereochemical features of A2′p5′ A which accounts for its higher activity.     

Z-RNA: the solution NMR structure of r(CGCGCG).

Left-handed double-helical Z-RNA has been studied using the ribohexanucleotide pentaphosphate r(CpGpCpGpCpG). One-dimensional and two-dimensional proton nmr experiments were used to probe the structural details of the left-handed helix in concentrated sodium perchlorate solution. In 1M NaClO4 the RNA adopts the normal A-form double helix, and in 6M NaClO4 it is nearly all in the Z form. In 4M NaClO4 it exists as nearly equal parts of A form and Z form. Resonances corresponding to both A and Z form appear in the nmr spectrum, indicating that the duplex exchanges slowly between forms. Spin-spin coupling constants between protons in the ribose rings were used to determine the sugar-pucker conformations of the individual nucleotides. Quantitative nuclear Overhauser experiments were used to determine proton-proton distances within the nucleoside, and from these distances values for the glycosidic torsion angle were determined. The results show that the cytidines adopt C2'-endo sugar puckers (S type) with pseudo-rotation phase values (P) of approximately 165 degrees. The bases are in the anti conformation, with chi values of approximately -140 degrees. The internal guanosines adopt C3'-endo sugar puckers (N type) with P approximately 18 degrees, while the 3'-terminal guanosine ribose exists in an equilibrium between S- and N-type conformations. All three guanosine bases adopt the syn conformation, with chi approximately 70 degrees. The results indicate that the solution structure of Z-RNA is very similar to that of Z-DNA.     

BII nucleotides in the B and C forms of natural-sequence polymeric DNA: A new model for the C form of DNA

A combination of solid-state (31)P and (13)C NMR, X-ray diffraction, and model building is used to show that the B and C forms of fibrous macromolecular DNA consist of two distinct nucleotide conformations, which correspond closely to the BI and BII nucleotide conformations known from oligonucleotide crystals. The proportion of the BII conformation is higher in the C form than in the B form. We show structural models for a 10(1) double helix involving BI nucleotides and a 9(1) double helix involving BII nucleotides. The 10(1) BI model is similar to a previous model of B-form DNA, while the 9(1) BII model is novel. The BII model has a very deep and narrow minor groove, a shallow and wide major groove, and highly inclined bases. This work shows that the B to C transition in fibers corresponds to BI to BII conformational changes of the individual nucleotides.     

Molecular dynamics simulations of xDNA

xDNA is a modified DNA, which contains natural as well as expanded bases. Expanded bases are generated by the addition of a benzene spacer to the natural bases. A set of AMBER force‐field parameters were derived for the expanded bases and the structural dynamics of the xDNA decamer ( xT5 ′ G xT A xC xG C xA xG T3′ ) · ( xA5′ C T xG C G xT A xC A3′) was explored using a 22 ns molecular dynamics simulation in explicit solvent. During the simulation, the duplex retained its Watson‐Crick base‐pairing and double helical structure, with deviations from the starting B‐form geometry towards A‐form; the deviations are mainly in the backbone torsion angles and in the helical parameters. The sugar pucker of the residues were distributed among a variety of modes; C2′ endo, C1′ exo, O4′ endo, C4′ exo, C2′ exo, and C3′ endo. The enhanced stacking interactions on account of the modification in the bases could help to retain the duplex nature of the helix with minor deviations from the ideal geometry. In our simulation, the xDNA showed a reduced minor groove width and an enlarged major groove width in comparison with the NMR structure. Both the grooves are larger than that of standard B‐DNA, but major groove width is larger than that of A‐DNA with almost equal minor groove width. The enlarged groove widths and the possibility of additional hydration in the grooves makes xDNA a potential molecule for various applications. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 351–360, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Rapid communication capturing the destabilizing effect of dihydrouridine through molecular simulations

The structural effects of the commonly occurring modified nucleoside dihydrouridine (D) observed experimentally in model oligonucleotides include a strong destabilization of the C3′‐endo sugar conformation of D, the disruption of stacking interactions of neighboring residues with D and a possible destabilization of the C3′‐endo sugar pucker of the 5′‐neighboring nucleoside. Our simulations with a combination of a set of parameters for modified RNA residues with the recently developed AMBER FF99χ force field having reoptimized glycosidic torsion angle parameters for standard nucleosides was found to reproduce the destabilizing effect of dihydrouridine better than with the AMBER FF99 force field for nucleic acids for which the parameters for the modified residues were originally developed. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 985–991, 2014.     

Crystal structure of the protein from Arabidopsis thaliana gene At5g06450, a putative DnaQ‐like exonuclease domain‐containing protein with homohexameric assembly

Arabidopsis thaliana gene At5g06450 encodes a putative DnaQ‐like 3′‐5′ exonuclease domain‐containing protein (AtDECP). The DnaQ‐like 3′‐5′ exonuclease domain is often found as a proofreading domain of DNA polymerases. The overall structure of AtDECP adopts an RNase H fold that consists of a mixed β‐sheet flanked by α‐helices. Interestingly, AtDECP forms a homohexameric assembly with a central six fold symmetry, generating a central cavity. The ring‐shaped structure and comparison with WRN‐exo, the best structural homologue of AtDECP, suggest a possible mechanism for implementing its exonuclease activity using positively charged patch on the N‐terminal side of the homohexameric assembly. The homohexameric structure of AtDECP provides unique information about the interaction between the DnaQ‐like 3′‐5′ exonuclease and its substrate nucleic acids.Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Properties of Antigenomic Hepatitis Delta Virus Ribozyme Cis‐ and Trans‐ Analogs

A series of permuted variants of antigenomic HDV ribozyme and trans‐acting variants were constructed. The catalytic activity study of the ribozymes has shown that all the variants were capable of self‐cleaving with equally biphasic kinetics. Ribonuclease and Fe(II)‐EDTA cleavage have provided evidence that all designed ribozymes fold according to the pseudoknot model and the conformations of the initial and cleaved ribozyme are different. A scheme of HDV ribozyme self‐cleavage reaction was suggested. The role of hydrogen bonds in the reaction was evaluated by substitution of ribose in the ribozyme for deoxyribose. It was found that the 2′‐OH group of U23 and C27 is critical for the reaction to occur; the 2′‐OH group of U32 and U39 is important, while 2′‐OH groups of other nucleotides of loop 3, stem 4 and stem 1 are unimportant for the cleavage activity.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016