Search

Simulating water accurately has been a major challenge in atomistic simulations for decades. Inclusion of electronic polarizability effects holds considerable promise, yet existing approaches suffer from significant computational overheads compared to the widely used nonpolarizable water models. We have developed a globally optimal polarizable water model, OPC3-pol, that explicitly accounts for electronic polarizability with minimal impact on the computational efficiency. OPC3-pol reproduces five key bulk water properties at room temperature with an average relative error of 0.6%. In atomistic simulations, OPC3-pol's computational efficiency is in between that of 3- and 4-point nonpolarizable models; the model supports increased (4 fs) integration time step. OPC3-pol is tested in simulations of globular protein ubiquitin and a B-DNA dodecamer with several AMBER force fields, ff99SB, ff14SB, ff19SB, and OL15, demonstrating structure stability close to reference on multi-microsecond time scale. Simulation of an intrinsically disordered amyloid β-peptide yields an ensemble with the radius of gyration of a random coil. The proposed water model can be trivially adopted by any package that supports standard nonpolarizable force fields and water models; its intended use is in long classical atomistic simulations where water polarization effects are expected to be important.

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

A fundamental property of DNA built from four informational nucleotide units (GCAT) is its ability to adopt different helical forms within the context of the Watson-Crick pair. Well-characterized examples include A-, B-, and Z-DNA. For this study, we created an isoinformational biomimetic polymer, built (like standard DNA) from four informational "letters", but with the building blocks being artificial. This ALternative Isoinformational ENgineered (ALIEN) DNA was hypothesized to support two nucleobase pairs, the

This report describes the chemical synthesis of aminotroponyl-conjugated deoxyuridine analog (at-dU) and its difluoroboron complex (dfbat-dU) and their phosphoramidites by using the versatile phosphorylating reagent 2-Cyanoethyl N,N-diisopropylchlorophosphoramidite. Tropolone is a non-benzenoid aromatic bioactive natural fluorescent molecule, possessing intramolecular charge transfer and metal chelating properties with transition metal ions such as Cu

An accurate molecular mechanics force field forms the basis of Molecular Dynamics simulations to obtain a realistic view of the structure and dynamics of biomolecules such as DNA. Although frequently updated to improve agreement with available experimental data, DNA force fields still rely in part on parameters introduced more than 20 years ago. We have developed an entirely new DNA force field, Tumuc1, derived from quantum mechanical calculations to obtain a consistent set of bonded parameters and partial atomic charges. The performance of the force field was extensively tested on a variety of DNA molecules. It excels in accuracy of B-DNA simulations but also performs very well on other types of DNA structures and structure formation processes such as hairpin folding, duplex formation, and dynamics of DNA-protein complexes. It can complement existing force fields in order to provide an increasingly accurate description of the structure and dynamics of DNA during simulation studies.

The thermal conductivity of double-stranded (ds) B-DNA was systematically investigated using classical molecular dynamics (MD) simulations. The effect of changing the base pair (bp) on the thermal conductivity of dsDNA needed investigation at a molecular level. Hence, four sequences, viz., poly(A), poly(G), poly(CG), and poly(AT), were initially analyzed in this work. First, the length of these sequences was varied from 4 to 40 bp at 300 K, and the respective thermal conductivity (κ) was computed. Second, the temperature-dependent thermal conductivities between 100 and 400 K were obtained in 50 K steps at 28 bp length. The Muller-Plathe reverse nonequilibrium molecular dynamics (RNEMD) was employed to set a thermal gradient and obtain all thermal conductivities in this work. Moreover, mixed sequences using AT and CG sequences, namely, A(CG)nT (n = 3-7), ACGC(AT)mGCGT (m = 0-5), and ACGC(AT)nAGCGT (n = 1-4), were investigated based on the hypothesis that these sequences could be better thermoelectrics. One-dimensional lattices are said to have diverging thermal conductivities at longer lengths, which violate the Fourier law. These follow the power law, where κ ∝ Lβ. At longer lengths, the exponent β needs to satisfy the condition β > 1/3 for divergent thermal conductivity. We find no such significant Fourier law violation through divergence of thermal conductivities at 80 bp lengths or 40 bp lengths. Also, in the case of the second study, the presence of short (m ≤ 2) encapsulated AT sequences within CG sequences shows an increasing trend. These results are important for engineering DNA-based thermal devices.

Molecular dynamics simulations and geometry optimizations carried out at the quantum level as well as by quantum mechanical/molecular mechanics methods predict that short, single-stranded DNA oligonucleotides adopt conformations very similar to those observed in crystallographic double-stranded B-DNA, with rise coordinates close to ≈3.3 Å. In agreement with the experimental evidence, the computational results show that DNA single strands rich in adjacent purine nucleobases assume more regular arrangements than poly-thymine. The preliminary results suggest that single-stranded poly-cytosine DNA should also retain a substantial helical order in solution. A comparison of the structures of single and double helices confirms that the B-DNA motif is a favorable arrangement also for single strands. Indeed, the optimal geometry of the complementary single helices is changed to a very small extent in the formation of the duplex.

DNA-protein interactions are ubiquitous in cellular processes. Impeding unwanted nucleic acid interactions and protein recognition have therapeutic implications. Therefore, new chemical scaffolds and studies related to their structural basis of nucleic acid recognition are essential for developing high-affinity DNA binders. In this study, we have employed a fragment-based strategy to design and synthesize benzimidazole-guanidinium hybrid compounds and study the individual fragment's role in imparting selectivity and specificity in DNA recognition. The fragments were extensively studied using thermal denaturation, circular dichroism, UV-vis absorption spectroscopy, and molecular docking techniques. The results indicate an interdependent role of the benzimidazole core, polar ends, and the DNA composition in imparting sequence-selective binding of the benzimidazole-guanidinium hybrid compounds in the DNA minor groove. Circular dichroism and molecular docking studies indicated minor groove binding analogous to classical minor groove binders such as DAPI and Hoechst 33258. Thermal denaturation studies indicated that the best binder (compound

Hole transfer along the axis of duplex DNA has been the focus of physical chemistry research for decades, with implications in diverse fields, from nanotechnology to cell oxidative damage. Computational approaches are particularly amenable for this problem, to complement experimental data for interpretation of transfer mechanisms. To be predictive, computational results need to account for the inherent mobility of biological molecules during the time frame of experimental measurements. Here, we address the structural variability of B-DNA and its effects on hole transfer in a combined molecular dynamics (MD) and real-time time-dependent density functional theory (RT-TDDFT) study. Our results show that quantities that characterize the charge transfer process, such as the time-dependent dipole moment and hole population at a specific site, are sensitive to structural changes that occur on the nanosecond time scale. We extend the range of physical properties for which such a correlation has been observed, further establishing the fact that quantitative computational data on charge transfer properties should include statistical averages. Furthermore, we use the RT-TDDFT results to assess an efficient tight-binding method suitable for high-throughput predictions. We demonstrate that charge transfer, although affected by structural variability, on average, remains strong in AA and GG dimers.

Branched DNA (bDNA) nanostructures have emerged as self-assembled biomaterials and are being considered for biomedical applications. Herein, we report the biophysical interaction between self-assembled bDNA nanostructure with circulating protein bovine serum albumin (BSA) and cellular enzyme bovine liver catalase (BLC). The binding between bDNA and BSA or BLC was confirmed through the decrease in fluorescence spectra. The Stern-Volmer data supports for non-covalent bonding with ~1 binding site in case of BSA and BLC thus advocating a static binding. Furthermore, FTIR and ITC study confirmed the binding of bDNAs with proteins through hydrogen bonding and van der Waals interaction. The negative free energy observed in ITC represent spontaneous reaction for BLC-bDNA interaction. The biophysical interaction between bDNA nanostructures and proteins was also supported by DLS and zeta potential measurement. With an increase in bDNA concentrations up to 100 nM, no significant change in absorbance and CD spectra was observed for both BLC and BSA which suggests structural stability and unaffected secondary conformation of proteins in presence of bDNA. Furthermore, the catalytic activity of BLC was unaltered in presence of bDNAscr even with increasing the incubation period from 1 h to 24 h. Interestingly, the time-dependent decrease in activity of BLC was protected by bDNAmix. The thermal melting study suggests a higher Tm value for proteins in presence of bDNAmix which demonstrates that interaction with bDNAmix increases the thermal stability of proteins. Collectively these data suggest that self-assembled DNA nanostructure may bind to BSA for facilitating circulation in plasma or binding to intracellular proteins like BLC for stabilization, however the secondary conformation of protein or catalytic activity of enzyme is unaltered in presence of bDNA nanostructure. Thus, the newly established genomic sequence-driven self-assembled DNA nanostructure can be explored for in vitro or in vivo experimental work in recent future.

A variety of sensitive and specific molecular diagnostic assays has been described for detecting nucleic acids in biological samples that may harbor pathogens of interest. These methods include very rapid, isothermal nucleic acid amplification methods that can be deployed outside of the laboratory environment, such as loop-mediated isothermal DNA amplification (LAMP) and recombinase-polymerase amplification (RPA). However, all molecular diagnostic assays must be preceded by nucleic acid extraction from the biological samples of interest, which provides suitable template molecules for the assays. To exploit the features of the amplification assays and be utilized outside of the lab, these methods must be rapid and avoid the need for typical laboratory chemicals and equipment. We describe a protocol for the extraction of DNA from field-collected insects that can be implemented at the point of collection and used to detect the presence of DNA sequences from potential plant pathogens that may be vectored by the insects. This protocol provides template DNA that is suitable for PCR, LAMP, and RPA. The FTA PlantSaver card-based DNA extraction product was also confirmed to amplify the mitochondrial cytochrome oxidase 1 (CO1) universal barcode that could later be sequenced to identify any insect. Lastly, we provide an example using field-collected insects, Neokolla (Graphocephala) heiroglyphica, and demonstrate the detection of the plant pathogen Xylella fastidiosa in carrier insects using PCR, RPA, and LAMP.

Guanine is the most susceptible base to oxidation damage induced by reactive oxygen species including singlet oxygen (1O2, 1Δg). We clarify whether the first step of guanine oxidation in B−DNA proceeds via either a zwitterionic or a diradical intermediate. The free energy profiles are calculated by means of a combined quantum mechanical and molecular mechanical (QM/MM) method coupled with the adaptive biasing force (ABF) method. To describe the open‐shell electronic structure of 1O2 correctly, the broken‐symmetry spin‐unrestricted density functional theory (BS−UDFT) with an approximate spin projection (AP) correction is applied to the QM region. We find that the effect of spin contamination on the activation and reaction free energies is up to ∼8 kcal mol−1, which is too large to be neglected. The QM(AP−ULC−BLYP)/MM‐based free energy calculations also reveal that the reaction proceeds through a diradical transition state, followed by a conversion to a zwitterionic intermediate. Our computed activation energy of 5.2 kcal mol−1 matches experimentally observed range (0∼6 kcal mol−1)., Realistic (free energy) profile: We demonstrate how singlet oxygen reacts with guanine in B−DNA by means of the spin‐projected QM(AP−ULC−BLYP)/MM free energy simulations. The present approach gives a realistic free energy profile, which agrees well with experimental findings. The rate‐limiting syn‐addition of 1O2 is found to proceed via a diradical transition state, followed by the smooth conversion from the open‐shell surface to the closed‐shell one.

Eleven copper chloride coordination compounds (1-11) with 4'-(4'-substituted-phenyl)-2,2':6',2''-terpyridine ligands bearing hydrogen (L1), cyano (L2), p-hydroxyl (L3), m-hydroxyl (L4), o-hydroxyl (L5), methoxyl (L6), iodo (L7), bromo (L8), chloro (L9), fluoro (L10) or methylsulfonyl (L11) were prepared and characterized by IR spectroscopy, elemental analysis and single crystal X-ray diffraction. Antiproliferative activities against tumor cells were investigated and DNA interactions were studied by circular dichroism spectroscopy and molecular modeling methods. In vitro data demonstrate that all the compounds exhibit higher antiproliferative activities as compared to cisplatin against five human carcinoma cell lines: A549, Bel-7402, Eca-109, HeLa and MCF-7. Compound 6 with methoxyl shows the best anti-proliferation activity. Spectrophotometric results reveal the strong affinity of the compounds for binding with DNA as intercalators and induce DNA conformational transitions. The results of molecular docking studies show that the compounds interact with DNA through π-π stacking, van der Waals forces, hydrophobic interactions and hydrogen bonds. The binding energies between compound 11 and three macromolecules, including DNA duplex, oligonucleotide and DNA-Topo I complex, are the lowest. The binding stability of compounds containing hydroxyl, methoxy and methylsulfonyl groups with biological macromolecules mainly relies on the hydrogen bonds. The ability of a compound to form hydrogen bonds can promote its binding to biological targets, thereby exhibiting high antiproliferative activity.

Left-handed Z-DNA is radically different from the most common right-handed B-DNA and can be stabilized by interactions with the Zα domain, which is found in a group of proteins, such as human ADAR1 and viral E3L proteins. It is well-known that most Zα domains bind to Z-DNA in a conformation-specific manner and induce rapid B–Z transition in physiological conditions. Although many structural and biochemical studies have identified the detailed interactions between the Zα domain and Z-DNA, little is known about the molecular basis of the B–Z transition process. In this study, we successfully converted the B–Z transition-defective Zα domain, vvZαE3L, into a B–Z converter by improving B-DNA binding ability, suggesting that B-DNA binding is involved in the B–Z transition. In addition, we engineered the canonical B-DNA binding protein GH5 into a Zα-like protein having both Z-DNA binding and B–Z transition activities by introducing Z-DNA interacting residues. Crystal structures of these mutants of vvZαE3L and GH5 complexed with Z-DNA confirmed the significance of conserved Z-DNA binding interactions. Altogether, our results provide molecular insight into how Zα domains obtain unusual conformational specificity and induce the B–Z transition.

Extrachromosomal DNA (ecDNA) amplification promotes intratumoral genetic heterogeneity and accelerated tumor evolution1–3, but its frequency and clinical impact are unclear. Here we show, using computational analysis of whole-genome sequencing data from 3,212 cancer patients, that ecDNA amplification frequently occurs in most cancer types, but not in blood or normal tissue. Oncogenes were highly enriched on amplified ecDNA and the most common recurrent oncogene amplifications arise on ecDNA. EcDNA amplifications resulted in higher levels of oncogene transcription compared to copy number matched linear DNA, coupled with enhanced chromatin accessibility and more frequently resulted in transcript fusions. Patients whose cancers carry ecDNAs have significantly shorter survival, even when controlled for tissue type, than do patients whose cancers are not driven by ecDNA-based oncogene amplification. The results presented here demonstrate that ecDNA-based oncogene amplification is common in cancer, is different from chromosomal amplification and drives poor outcome for patients across many cancer types.

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.

Left-handed G quadruplexes (LHG4) have been recently discovered as a new class of G quadruplexes. The biological functions of LHG4s are still unknown, but they share a striking resemblance to Z-DNA in their helicity and jagged phosphate backbone. To further understand structural features of the LHG4s that define their left handedness, we have employed human-interpretable machine-learning methods to classify right- and left-handed G4s purely based on torsional angle analysis. Our results reveal the importance of the α, β, δ, and χ angles in left-handed structuring across both Z-DNAs and LHG4s. Our analysis may serve as the first step to understanding the conditions of formation for LHG4s and their potential biological relevance.

Existence of numerous biomolecules results in biological fluids to be extremely crowded. Thus, Macromolecular crowding is an essential phenomenon to sustain active conformation of proteins in biological systems. Herein, double helical deoxyribonucleic acid (B-DNA) is presented for the first time as a biomacromolecular crowding system for sustainable packaging of cytochrome c (Cyt C). The peroxidase activity of Cyt C was investigated in the presence of various concentrations of B-DNA (from salmon milt). At an optimized concentration of 0.125 mg/mL B-DNA, an 11-fold higher catalytic activity was found than in native Cyt C with improved stability. Molecular docking and spectroscopic analyses revealed that electrostatic and H-bonding are the main interactions between DNA and Cyt C that affect the structural stability and activity of the protein. Moreover, the catalytic activity and stability of the protein were further investigated in the presence of severe process conditions by UV-visible, circular dichroism, and Fourier-transform infrared spectroscopies. Molecularly crowded Cyt C showed significantly higher activity and stability under severe environments such as high temperature (110 °C), oxidative stress, high pH (pH 10) and biological (trypsin) and chemical denaturants (urea) compared to bare Cyt C. The observed results support the suitability of DNA-based macromolecular crowding media as a viable and effective stabilizer of proteins against multiple stresses.

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double‐stranded B‐DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson–Crick base pairing, using dispersion‐corrected DFT at the BLYP‐D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better‐known effects of partial desolvation. Furthermore, we show that the sequence‐dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so‐called “diagonal interactions” or “cross terms”. Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases., All about that sequence! Quantum chemical analyses reveal that the stability of B‐DNA triplets depends not only on the nucleotide composition, but also on the order in which the nucleotides occur. This sequence dependence in the duplex stability results primarily from attractive or repulsive cross terms between bases in adjacent Watson–Crick pairs.

Oncogene amplification on extrachromosomal DNA (ecDNA) is prevalent in human cancer and is associated with poor outcomes. Clonal, megabase-sized circular ecDNAs in cancer are distinct from nonclonal, small sub-kilobase-sized DNAs that may arise during normal tissue homeostasis. ecDNAs enable profound changes in gene regulation beyond copy-number gains. An emerging principle of ecDNA regulation is the formation of ecDNA hubs: micrometer-sized nuclear structures of numerous copies of ecDNAs tethered by proteins in spatial proximity. ecDNA hubs enable cooperative and intermolecular sharing of DNA regulatory elements for potent and combinatorial gene activation. The 3D context of ecDNA shapes its gene expression potential, selection for clonal heterogeneity among ecDNAs, distribution through cell division, and reintegration into chromosomes. Technologies for studying gene regulation and structure of ecDNA are starting to answer long-held questions on the distinct rules that govern cancer genes beyond chromosomes.

The mechanisms of interaction of DNA with pharmacological molecules are critical to understanding their therapeutic actions on physiological systems. Piperlongumine is widely studied for its anticancer potential. Multi-spectrometry, calorimetry and in silico studies were employed to study the interaction of piperlongumine and calf thymus DNA. UV-Vis spectroscopy illustrated a hyperchromic pattern in spectra of the calf thymus DNA-piperlongumine complex, while fluorescent quenching was observed in emission spectral studies. Competitive displacement assay demonstrated higher displacement and binding constant for DNA-rhodamine B complex by piperlongumine than DNA-methylene blue complex. Differential scanning calorimetry presented non-significant changes in melting temperature and molecular docking presented the precise interaction site of piperlongumine with calf thymus DNA at minor groove. Further, piperlongumine treatment did not result in pBluescript KS plasmid DNA cleavage as revealed from the DNA topology assay. All these experiments confirmed the binding of piperlongumine with DNA through minor groove binding mode.

Understanding the role of base pairing and stacking displayed by polynucleotide chains interwind together resulting in a double-helical B-DNA type structure is crucial to gaining access to the sophisticated structural arrangement of DNA. Several computational and experimental studies hinted towards the dominance of base pairing over stacking for duplex stability. To find out how significant the individual Watson-Crick hydrogen bonds are in maintaining the double-helical integrity of the DNA, in the present article, we selectively switched off the hydrogen bonds (one specific bond or their combinations in all the base pairs at a time) via manipulating the force fields for A-T and G-C base pairs. We studied 12 systems in total via all-atom explicit-solvent molecular dynamics simulations (200 ns each). The MD output structures were compared with the control system by means of structural, dynamic, and energetic properties to monitor the overall consequences of removing H-bond(s) on the B-DNA characteristics of the model systems. Our findings suggest that all the individual hydrogen bonds involved in base pairing are vital for maintaining the DNA structural integrity as any possible alteration in Watson-Crick hydrogen bond(s) leads to the disintegration/collapse of DNA strands resulting in unfolded states.

Methylcytosines in mammalian genomes are the main epigenetic molecular codes that switch off the repertoire of genes in cell-type and cell-stage dependent manners. DNA methyltransferases (DMT) are dedicated to managing the status of cytosine methylation. DNA methylation is not only critical in normal development, but it is also implicated in cancers, degeneration, and senescence. Thus, the chemicals to control DMT have been suggested as anticancer drugs by reprogramming the gene expression profile in malignant cells. Here, we report a new optical technique to characterize the activity of DMT and the effect of inhibitors, utilizing the methylation-sensitive B-Z transition of DNA without bisulfite conversion, methylation-sensing proteins, and polymerase chain reaction amplification. With the high sensitivity of single-molecule FRET, this method detects the event of DNA methylation in a single DNA molecule and circumvents the need for amplification steps, permitting direct interpretation. This method also responds to hemi-methylated DNA. Dispensing with methylation-sensitive nucleases, this method preserves the molecular integrity and methylation state of target molecules. Sparing methylation-sensing nucleases and antibodies helps to avoid errors introduced by the antibody’s incomplete specificity or variable activity of nucleases. With this new method, we demonstrated the inhibitory effect of several natural bio-active compounds on DMT. All taken together, our method offers quantitative assays for DMT and DMT-related anticancer drugs.

Extrachromosomal DNA (ecDNA) is a type of circular and tumor specific genetic element. EcDNA has been reported to display open chromatin structure, facilitate oncogene amplification and genetic material unequal segregation, and is associated with poor cancer patients’ prognosis. The ability of immune evasion is a typical feature for cancer progression, however the tumor intrinsic factors that determine immune evasion remain poorly understood. Here we show that the presence of ecDNA is associated with markers of tumor immune evasion, and obtaining ecDNA could be one of the mechanisms employed by tumor cells to escape immune surveillance. Tumors with ecDNA usually have comparable TMB and neoantigen load, however they have lower immune cell infiltration and lower cytotoxic T cell activity. The microenvironment of tumors with ecDNA shows increased immune-depleted, decreased immune-enriched fibrotic types. Both MHC class I and class II antigen presentation genes’ expression are decreased in tumors with ecDNA, and this could be the underlying mechanism for ecDNA associated immune evasion. This study provides evidence that ecDNA formation is an immune escape mechanism for cancer cells.

G-quadruplex existence was proved in cells by using both antibodies and small molecule fluorescent probes. However, the G-quadruplex probes designed thus far are structure- but not conformation-specific. Recently, a core-extended naphthalene diimide (cex-NDI) was designed and found to provide fluorescent signals of markedly different intensities when bound to G-quadruplexes of different conformations or duplexes. Aiming at evaluating how the fluorescence behaviour of this compound is associated with specific binding modes to the different DNA targets, cex-NDI was here studied in its interaction with hybrid G-quadruplex, parallel G-quadruplex, and B-DNA duplex models by biophysical techniques, molecular docking, and biological assays. cex-NDI showed different binding modes associated with different amounts of stacking interactions with the three DNA targets. The preferential binding sites were the groove, outer quartet, or intercalative site of the hybrid G-quadruplex, parallel G-quadruplex, and B-DNA duplex, respectively. Interestingly, our data show that the fluorescence intensity of DNA-bound cex-NDI correlates with the amount of stacking interactions formed by the ligand with each DNA target, thus providing the rationale behind the conformation-sensitive properties of cex-NDI and supporting its use as a fluorescent probe of G-quadruplex structures. Notably, biological assays proved that cex-NDI mainly localizes in the G-quadruplex-rich nuclei of cancer cells.

Accurate parametrization of force fields (FFs) is of ultimate importance for computer simulations to be reliable and to possess a predictive power. In this work, we analyzed, in multi-microsecond simulations of a 40-base-pair DNA fragment, the performance of four force fields, namely, the two recent major updates of CHARMM and two from the AMBER family. We focused on a description of double-helix DNA flexibility and dynamics both at atomistic and at mesoscale level in coarse-grained (CG) simulations. In addition to the traditional analysis of different base-pair and base-step parameters, we extended our analysis to investigate the ability of the force field to parametrize a CG DNA model by structure-based bottom-up coarse-graining, computing DNA persistence length as a function of ionic strength. Our simulations unambiguously showed that the CHARMM36 force field is unable to preserve DNA's structural stability at over-microsecond time scale. Both versions of the AMBER FF, parmbsc0 and parmbsc1, showed good agreement with experiment, with some bias of parmbsc0 parameters for intermediate A/B form DNA structures. The CHARMM27 force field provides stable atomistic trajectories and overall (among the considered force fields) the best fit to experimentally determined DNA flexibility parameters both at atomistic and at mesoscale level. Ministry of Education (MOE) Accepted version This work was supported by the Singapore Ministry of Education Academic Research Fund (AcRF) Tier 2 (MOE2014-T2-1-123 (ARC51/14)) and Tier 3 (MOE2012- T3-1-001) grants (to L.N.) and by the Swedish Research Council (grant 2017-03950 to A.P.L.).

GapR is a nucleoid-associated protein required for the cell cycle of Caulobacter cresentus. We have determined new crystal structures of GapR to high resolution. As in a recently published structure, a GapR monomer folds into one long N-terminal α helix and two shorter α helices, and assembles into a tetrameric ring with a closed, positively charged, central channel. In contrast to the conclusions drawn from the published structures, we observe that the central channel of the tetramer presented here could freely accommodate B-DNA. Mutation of six conserved lysine residues lining the cavity and electrophoretic mobility gel shift experiments confirmed their role in DNA binding and the channel as the site of DNA binding. Although present in our crystals, DNA could not be observed in the electron density maps, suggesting that DNA binding is non-specific, which could be important for tetramer-ring translocation along the chromosome. In conjunction with previous GapR structures we propose a model for DNA binding and translocation that explains key published observations on GapR and its biological functions.

P element transposase catalyzes the mobility of P element DNA transposons within the Drosophila genome. P element transposase exhibits several unique properties, including the requirement for a guanosine triphosphate cofactor and the generation of long staggered DNA breaks during transposition. To gain insights into these features, we determined the atomic structure of the Drosophila P element transposase strand transfer complex using cryo-EM. The structure of this post-transposition nucleoprotein complex reveals that the terminal single-stranded transposon DNA adopts unusual A-form and distorted B-form helical geometries that are stabilized by extensive protein-DNA interactions. Additionally, we infer that the bound guanosine triphosphate cofactor interacts with the terminal base of the transposon DNA, apparently to position the P element DNA for catalysis. Our structure provides the first view of the P element transposase superfamily, offers new insights into P element transposition and implies a transposition pathway fundamentally distinct from other cut-and-paste DNA transposases.

Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here we report solution structures and dynamics of duplex DNA containing two all-natural C-HgII-T metallo base pairs separated by six canonical base pairs. NMR experiments reveal a 3:1 ratio of well-resolved structures in dynamic equilibrium. The major species contains two (N3)T-HgII-(N3)C base pairs in a predominantly B-form helix. The minor species contains (N3)T-HgII-(N4)C base pairs and greater A-form characteristics. Ten-fold different 1J coupling constants (15N,199Hg) are observed for (N3)C-HgII (114 Hz) versus (N4)C-HgII (1052 Hz) connectivities, reflecting differences in cytosine ionization and metal-bonding strengths. Dynamic interconversion between the two types of C-HgII-T base pairs are coupled to a global conformational exchange between the helices. These observations inspired the design of a repetitive DNA sequence capable of undergoing a global B-to-A-form helical transition upon adding HgII, demonstrating that C-HgII-T has unique switching potential in DNA-based materials and devices., Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here, the authors prepared a metallo-DNA duplex including two C-Hg(II)-T base pairs separated by six normal Watson-Crick base pairs and investigated its solution structure and dynamics using NMR spectroscopy.

Significance The main stabilizer of the DNA double helix is not the base-pair hydrogen bonds but coin-pile stacking of base pairs, whose hydrophobic cohesion, requiring abundant water, indirectly makes the DNA interior dry so that hydrogen bonds can exert full recognition power. We report that certain semihydrophobic agents depress the stacking energy (measurable in single-molecule experiments), leading to transiently occurring holes in the base-pair stack (monitorable via binding of threading intercalators). Similar structures observed in DNA complexes with RecA and Rad51, and previous observations of spontaneous strand exchange catalyzed in semihydrophobic model systems, make us propose that some hydrophobic protein residues may have roles in catalyzing homologous recombination. We speculate that hydrophobic catalysis is a general phenomenon in DNA enzymes., Hydrophobic base stacking is a major contributor to DNA double-helix stability. We report the discovery of specific unstacking effects in certain semihydrophobic environments. Water-miscible ethylene glycol ethers are found to modify structure, dynamics, and reactivity of DNA by mechanisms possibly related to a biologically relevant hydrophobic catalysis. Spectroscopic data and optical tweezers experiments show that base-stacking energies are reduced while base-pair hydrogen bonds are strengthened. We propose that a modulated chemical potential of water can promote “longitudinal breathing” and the formation of unstacked holes while base unpairing is suppressed. Flow linear dichroism in 20% diglyme indicates a 20 to 30% decrease in persistence length of DNA, supported by an increased flexibility in single-molecule nanochannel experiments in poly(ethylene glycol). A limited (3 to 6%) hyperchromicity but unaffected circular dichroism is consistent with transient unstacking events while maintaining an overall average B-DNA conformation. Further information about unstacking dynamics is obtained from the binding kinetics of large thread-intercalating ruthenium complexes, indicating that the hydrophobic effect provides a 10 to 100 times increased DNA unstacking frequency and an “open hole” population on the order of 10−2 compared to 10−4 in normal aqueous solution. Spontaneous DNA strand exchange catalyzed by poly(ethylene glycol) makes us propose that hydrophobic residues in the L2 loop of recombination enzymes RecA and Rad51 may assist gene recombination via modulation of water activity near the DNA helix by hydrophobic interactions, in the manner described here. We speculate that such hydrophobic interactions may have catalytic roles also in other biological contexts, such as in polymerases.

Protein-DNA interactions are of great biological importance. The specificity and strength of these intimate contacts are crucial in the proper functioning of a cell, wherein the role of DNA dynamic bendability has been a matter of discussion. We relate DNA bendability to protein binding by introducing some simple modifications in the DNA structure. We removed C5′ carbon in first modified structure and the second has an additional carbon between C3′ and 3′-OH, hereby pronounced as C(−) and C(+) nucleic acids respectively. We observed that C(+) nucleic acid retains B-DNA duplex as seen by means of 500 ns long molecular dynamics (MD) simulations, structural and energetic calculations, while C(−) nucleic acid attains a highly bend structure. We transferred these observations to a protein-DNA system in order to monitor as to what extent the bendability enhances the protein binding. The energetics of binding is explored by performing 100 ns long MD simulations on control and modified DNA-protein complexes followed by running MM-PBSA/GBSA calculations on the resultant structures. It is observed that C(+) nucleic acid has protein binding in close correspondence to the control system (∼−14 kcal/mol) due to their relatable structure, while the C(−) nucleic acid displayed high binding to the protein (∼−18 kcal/mol). DelPhi based calculations reveal that the high binding could be the result of enhanced electrostatic interactions caused by exposed bases in the bend structure for protein recognition. Such modified oligonucleotides, due to their improved binding to protein and resistance to nuclease degradation, have a great therapeutic value.

The integration of the HIV-1 genome into the host genome is an essential step in the life cycle of the virus and it plays a critical role in the expression, long-term persistence, and reactivation of HIV expression. To better understand the local genomic environment surrounding HIV-1 proviruses, we assessed the influence of non-canonical B-form DNA (non-B DNA) on the HIV-1 integration site selection. We showed that productively and latently infected cells exhibit different integration site biases towards non-B DNA motifs. We identified a correlation between the integration sites of the latent proviruses and non-B DNA features known to potently influence gene expression (e.g., cruciform, guanine-quadruplex (G4), triplex, and Z-DNA). The reactivation potential of latent proviruses with latency reversal agents also correlated with their proximity to specific non-B DNA motifs. The perturbation of G4 structures in vitro using G4 structure-destabilizing or -stabilizing ligands resulted in a significant reduction in integration within 100 base pairs of G4 motifs. The stabilization of G4 structures increased the integration within 300–500 base pairs from G4 motifs, increased integration near transcription start sites, and increased the proportion of latently infected cells. Moreover, we showed that host lens epithelium-derived growth factor (LEDGF)/p75 and cleavage and polyadenylation specificity factor 6 (CPSF6) influenced the distribution of integration sites near several non-B DNA motifs, especially G4 DNA. Our findings identify non-B DNA motifs as important factors that influence productive and latent HIV-1 integration and the reactivation potential of latent proviruses.

Although current AMBER force fields are relatively accurate for canonical B-DNA, many non-canonical structures are still described incorrectly. As non-canonical motifs are attracting increasing attention due to the role they play in living organisms, further improvement is desirable. Here, we have chosen Z-DNA molecule, can be considered a touchstone of the universality of empirical force fields, since the non-canonical α and γ backbone conformations native to Z-DNA are also found in protein-DNA complexes, i-motif DNA and other non-canonical DNAs. We show that spurious α/γ conformations occurring in simulations with current AMBER force fields, OL15 and bsc1, are largely due to inaccurate α/γ parameterization. Moreover, stabilization of native Z-DNA substates involving γ = trans conformations appears to be in conflict with the correct description of the canonical B-DNA structure. Because the balance of the native and spurious conformations is influenced by non-additive effects, this is a difficult case for an additive dihedral energy scheme such as AMBER. We propose new α/γ parameters, denoted OL21, and show that they improve the stability of native α/γ Z-DNA substates while keeping the canonical DNA description virtually unchanged, and thus represent a reasonable compromise within the additive force field framework. Although further extensive testing is needed, the new modification appears to be a promising step towards a more reliable description of non-canonical DNA motifs and provides the best performance for Z-DNA molecules among current AMBER force fields.