Hidden modes of DNA binding by human nuclear receptors.

Human nuclear receptors (NRs) are a superfamily of ligand-responsive transcription factors that have central roles in cellular function. Their malfunction is linked to numerous diseases, and the ability to modulate their activity with synthetic ligands has yielded 16% of all FDA-approved drugs. NRs regulate distinct gene networks, however they often function from genomic sites that lack known binding motifs. Here, to annotate genomic binding sites of known and unexamined NRs more accurately, we use high-throughput SELEX to comprehensively map DNA binding site preferences of all full-length human NRs, in complex with their ligands. Furthermore, to identify non-obvious binding sites buried in DNA-protein interactomes, we develop MinSeq Find, a search algorithm based on the MinTerm concept from electrical engineering and digital systems design. The resulting MinTerm sequence set (MinSeqs) reveal a constellation of binding sites that more effectively annotate NR-binding profiles in cells. MinSeqs also unmask binding sites created or disrupted by 52,106 single-nucleotide polymorphisms associated with human diseases. By implicating druggable NRs as hidden drivers of multiple human diseases, our results not only reveal new biological roles of NRs, but they also provide a resource for drug-repurposing and precision medicine.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif.

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

Flexibility and structure of flanking DNA impact transcription factor affinity for its core motif.

Spatial and temporal expression of genes is essential for maintaining phenotype integrity. Transcription factors (TFs) modulate expression patterns by binding to specific DNA sequences in the genome. Along with the core binding motif, the flanking sequence context can play a role in DNA-TF recognition. Here, we employ high-throughput in vitro and in silico analyses to understand the influence of sequences flanking the cognate sites in binding of three most prevalent eukaryotic TF families (zinc finger, homeodomain and bZIP). In vitro binding preferences of each TF toward the entire DNA sequence space were correlated with a wide range of DNA structural parameters, including DNA flexibility. Results demonstrate that conformational plasticity of flanking regions modulates binding affinity of certain TF families. DNA duplex stability and minor groove width also play an important role in DNA-TF recognition but differ in how exactly they influence the binding in each specific case. Our analyses further reveal that the structural features of preferred flanking sequences are not universal, as similar DNA-binding folds can employ distinct DNA recognition modes.

Dynamics of physiologically relevant noncanonical DNA structures: an overview from experimental and theoretical studies.

DNA is a complex molecule with phenomenal inherent plasticity and the ability to form different hydrogen bonding patterns of varying stabilities. These properties enable DNA to attain a variety of structural and conformational polymorphic forms. Structurally, DNA can exist in single-stranded form or as higher-order structures, which include the canonical double helix as well as the noncanonical duplex, triplex and quadruplex species. Each of these structural forms in turn encompasses an ensemble of dynamically heterogeneous conformers depending on the sequence composition and environmental context. In vivo, the widely populated canonical B-DNA attains these noncanonical polymorphs during important cellular processes. While several investigations have focused on the structure of these noncanonical DNA, studying their dynamics has remained nontrivial. Here, we outline findings from some recent advanced experimental and molecular simulation techniques that have significantly contributed toward understanding the complex dynamics of physiologically relevant noncanonical forms of DNA.

Ionic liquids make DNA rigid.

Persistence length of double-stranded DNA (dsDNA) is known to decrease with an increase in ionic concentration of the solution. In contrast to this, here we show that the persistence length of dsDNA increases dramatically as a function of ionic liquid (IL) concentration. Using all atom explicit solvent molecular dynamics simulations and theoretical models, we present, for the first time, a systematic study to determine the mechanical properties of dsDNA in various hydrated ILs at different concentrations. We find that dsDNA in 50 wt % ILs have lower persistence length and stretch modulus in comparison to 80 wt % ILs. We further observe that both the persistence length and stretch modulus of dsDNA increase as we increase the concentration of ILs. The present trend of the stretch modulus and persistence length of dsDNA with IL concentration supports the predictions of the macroscopic elastic theory, in contrast to the behavior exhibited by dsDNA in monovalent salt. Our study further suggests the preferable ILs that can be used for maintaining DNA stability during long-term storage.

Dehydrated DNA in B-form: ionic liquids in rescue.

The functional B-conformation of DNA succumbs to the A-form at low water activity. Methods for room temperature DNA storage that rely upon 'anhydrobiosis', thus, often encounter the loss of DNA activity due to the B→A-DNA transition. Here, we show that ionic liquids, an emerging class of green solvents, can induce conformational transitions in DNA and even rescue the dehydrated DNA in the functional B-form. CD spectroscopic analyses not only reveal rapid transition of A-DNA in 78% ethanol medium to B-conformation in presence of ILs, but also the high resistance of IL-bound B-form to transit to A-DNA under dehydration. Molecular dynamics simulations show the unique ability of ILs to disrupt Na+ ion condensation and form 'IL spine' in DNA minor groove to drive the A→B transition. Implications of these findings range from the plausible use of ILs as novel anhydrobiotic DNA storage medium to a switch for modulating DNA conformational transitions.

Ion-water wires in imidazolium-based ionic liquid/water solutions induce unique trends in density.

Ionic liquid/water binary mixtures are rapidly gaining popularity as solvents for dissolution of cellulose, nucleobases, and other poorly water-soluble biomolecules. Hence, several studies have focused on measuring the thermophysical properties of these versatile mixtures. Among these, 1-ethyl-3-methylimidazolium ([emim]) cation-based ILs containing different anions exhibit unique density behaviours upon addition of water. While [emim][acetate]/water binary mixtures display an unusual rise in density with the addition of low-to-moderate amounts of water, those containing the [trifluoroacetate] ([Tfa]) anion display a sluggish decrease in density. The density of [emim][tetrafluoroborate] ([emim][BF4])/water mixtures, on the other hand, declines rapidly in close accordance with the experimental reports. Here, we unravel the structural basis underlying this unique density behavior of [emim]-based IL/water mixtures using all-atom molecular dynamics (MD) simulations. The results revealed that the distinct nature of anion-water hydrogen bonded networks in the three systems was a key in modulating the observed unique density behaviour. Vast expanses of uninterrupted anion-water-anion H-bonded stretches, denoted here as anion-water wires, induced significant structuring in [emim][Ac]/water mixtures that resulted in the density rise. Conversely, the presence of intermittent large water clusters disintegrated the anion-water wires in [emim][Tfa]/water and [emim][BF4]/water mixtures to cause a monotonic density decrease. The differential nanostructuring affected the dynamics of the solutions proportionately, with the H-bond making and breaking dynamics found to be greatly retarded in [emim][Ac]/water mixtures, while it exhibited a faster relaxation in the other two binary solutions.

High Nucleobase-Solubilizing Ability of Low-Viscous Ionic Liquid/Water Mixtures: Measurements and Mechanism.

Research on nucleobases has always been on the forefront owing to their ever-increasing demand in the pharmaceutical, agricultural, and other industries. The applications, however, became limited due to their poor solubility in water. Recently, ionic liquids (ILs) have emerged as promising solvents for nucleobase dissolution, as they exhibit >100-fold increased solubility compared to water. But the high viscosity of ILs remains as a bottleneck in the field. Here, by solubility and viscosity measurements, we show that addition of low-to-moderate quantity of water preserves the high solubilizing capacity of ILs, while reducing the viscosity of the solution by several folds. To understand the mechanism of nucleobase dissolution, molecular dynamics simulations were carried out, which showed that at low concentrations water incorporates into the IL-nucleobase network without much perturbing of the nucleobase-IL interactions. At higher concentrations, increasing numbers of IL anion-water hydrogen bonds replace IL-nucleobase interactions, which have been confirmed by (1)H- and (13)C NMR chemical shifts of the IL ions.

Nanostructural Reorganization Manifests in Sui-Generis Density Trend of Imidazolium Acetate/Water Binary Mixtures.

Ionic liquids (ILs) are emerging as a novel class of solvents in chemical and biochemical research. Their range of applications further expands when a small quantity of water is added. Thus, the past decade has seen extensive research on IL/water binary mixtures. While the thermophysical properties of most of these mixtures exhibited the expected trend, few others have shown deviations from the general course. One such example is the increase in density of the 1-alkyl-3-methyl imidazolium acetate ([Rn mim][Ac])-based ILs with the addition of low to moderate concentrations of water. Although such a unique trend was observed for imidazolium cations of different tail lengths and also from independent experiments, the molecular basis of this unique behavior remains unknown. In this study, we examine the nanostructural reordering in [Rn mim][Ac] (n = 2-6) ILs due to added water by means of molecular dynamics simulations, and correlate the observed changes to the sui-generis density trend. Results suggest that the initial rise in density in these ILs mainly pertains to the water-induced increased spatial correlation among the polar components, where high basicity of the acetate anion plays a key role. At moderate water concentration, the density can rise further for ILs with longer cation tails due to hydrophobic clustering. Thus, while [emim][Ac]/water mixtures exhibit the density turnover at Xw = 0.5, [bmim][Ac] and [hmim][Ac] show the turnover at Xw = 0.7. The detailed understanding provided here could help the preparation of optimal IL/water binary mixtures for various biochemical applications.

Groove binding mechanism of ionic liquids: a key factor in long-term stability of DNA in hydrated ionic liquids?

Nucleic acid sample storage is of paramount importance in biotechnology and forensic sciences. Very recently, hydrated ionic liquids (ILs) have been identified as ideal media for long-term DNA storage. Hence, understanding the binding characteristics and molecular mechanism of interactions of ILs with DNA is of both practical and fundamental interest. Here, we employ molecular dynamics simulations and spectroscopic experiments to unravel the key factors that stabilize DNA in hydrated ILs. Both simulation and experimental results show that DNA maintains the native B-conformation in ILs. Simulation results further suggest that, apart from the electrostatic association of IL cations with the DNA backbone, groove binding of IL cations through hydrophobic and polar interactions contributes significantly to DNA stability. Circular dichroism spectral measurements and fluorescent dye displacement assay confirm the intrusion of IL molecules into the DNA minor groove. Very interestingly, the IL ions were seen to disrupt the water cage around DNA, including the spine of hydration in the minor groove. This partial dehydration by ILs likely prevents the hydrolytic reactions that denature DNA and helps stabilize DNA for the long term. The detailed understanding of IL-DNA interactions provided here could guide the future development of novel ILs, specific for nucleic acid solutes.