Assessing the Current State of Amber Force Field Modifications for DNA

Investigation of the Interaction between Aloe vera Anthraquinone Metabolites and c-Myc and C-Kit G-Quadruplex DNA Structures

G-quadruplexes are nucleotide sequences present in the promoter region of numerous oncogenes, having a key role in the suppression of gene transcription. Recently, the binding of anthraquinones from Aloe vera to G-quadruplex structures has been studied through various physico-chemical techniques. Intrigued by the reported results, we investigated the affinity of aloe emodin, aloe emodin-8-glucoside, and aloin to selected G-quadruplex nucleotide sequences by NMR spectroscopy. The structural determinants for the formation of the ligand/nucleotide complexes were elucidated and a model of the interactions between the tested compounds and C-Kit and c-Myc G-quadruplex DNA structures was built by integrated NMR and molecular modeling studies. Overall, the obtained results confirmed and implemented the previously reported findings, pointing out the complementarity of the different approaches and their contribution to a more detailed overview of the ligand/nucleotide complex formation. Furthermore, the proposed models of interaction could pave the way to the design of new nature-derived compounds endowed with increased G-quadruplex stabilizing activity.

Replicative bypass studies of l-deoxyribonucleosides in Vitro and in E. coli cell

L-nucleosides were the most important antiviral lead compounds because they can inhibit viral DNA polymerase and DNA synthesis of many viruses, whereas they may lead to mutations in DNA replication and cause genomic instability. In this study, we reported the replicative bypass of L-deoxynucleosides in recombinant DNA by restriction enzyme-mediated assays to examine their impact on DNA replication in vitro and in E. coli cells. The results showed that a template L-dC inhibited Taq DNA polymerase reaction, whereas it can be bypassed by Vent (exo^-) DNA polymerase as well as in cell replication, inserting correct nucleotides opposite L-dC. L-dG can be bypassed by Taq DNA polymerase and in E. coli cells, maintaining insertion of correct incoming nucleotides, and L-dG induced mutagenic replication by Vent (exo^-) DNA polymerase. In contrast, L-dA can induced mutagenic replication in vitro and in E. coli cells. MD simulations were performed to investigate how DNA polymerase affected replicative bypass and mutations when D-nucleosides replaced with L-nucleosides. This study will provide a basis for the ability to assess the cytotoxic and mutagenic properties of the L-nucleoside drugs.

The Essential Role of Prolines and Their Conformation in Allosteric Regulation of Kaiso Zinc Finger DNA-Binding Activity by the Adjacent C-Terminal Loop

Kaiso is a methyl-DNA-binding protein containing three C2H2 zinc fingers with a C-terminal extension that participates in DNA binding. The linker between the last zinc finger and the DNA-binding portion of the extension contains two prolines that are highly conserved in vertebrates and in cognate ZBTB4 and ZBTB38 proteins. Prolines provide chain rigidity and can exist in cis and trans conformations that can be switched by proline isomerases, affecting protein function. We found that substitution of the conserved proline P588, but not of P577, to alanine, negatively affected KaisoDNA-binding according to molecular dynamics simulation and in vitro DNA-binding assays. Molecular dynamics simulations of the Kaiso DNA-binding domain with P588 either substituted to alanine or switched to the cis-conformation revealed similar alterations in the H-bonding network and uncovered allosteric effects leading to structural rearrangements in the entire domain that resulted in the weakening of DNA-binding affinity. The substitution of proline with a large hydrophobic residue led to the same negative effects despite its ability to partially rescue the intrinsic DNA-binding activity of the C-terminal loop. Thus, the presence of the C-terminal extension and cis-conformation of proline residues are essential for efficient Kaiso-DNA binding, which likely involves intramolecular tension squeezing the DNA chain.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

Duplex DNA Retains the Conformational Features of Single Strands: Perspectives from MD Simulations and Quantum Chemical Computations

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.

Translesion Synthesis across the N2-Ethyl-deoxyguanosine Adduct by Human PrimPol

Primase-DNA polymerase (PrimPol) is involved in reinitiating DNA synthesis at stalled replication forks. PrimPol also possesses DNA translesion (TLS) activity and bypasses several endogenous nonbulky DNA lesions in vitro. Little is known about the TLS activity of PrimPol across bulky carcinogenic adducts. We analyzed the DNA polymerase activity of human PrimPol on DNA templates with seven N²-dG lesions of different steric bulkiness. In the presence of Mg²⁺ ions, bulky N²-isobutyl-dG, N²-benzyl-dG, N²-methyl(1-naphthyl)-dG, N²-methyl(9-anthracenyl)-dG, N²-methyl(1-pyrenyl)-dG, and N²-methyl(1,3-dimethoxyanthraquinone)-dG adducts fully blocked PrimPol activity. At the same time, PrimPol incorporated complementary deoxycytidine monophosphate (dCMP) opposite N²-ethyl-dG with moderate efficiency but did not extend DNA beyond the lesion. We also demonstrated that mutation of the Arg288 residue abrogated dCMP incorporation opposite the lesion in the presence of Mn²⁺ ions. When Mn²⁺ replaced Mg²⁺, PrimPol carried out DNA synthesis on all DNA templates with N²-dG adducts in standing start reactions with low efficiency and accuracy, possibly utilizing a lesion "skipping" mechanism. The TLS activity of PrimPol opposite N²-ethyl-dG but not bulkier adducts was stimulated by accessory proteins, polymerase delta-interacting protein 2 (PolDIP2), and replication protein A (RPA). Molecular dynamics studies demonstrated the absence of stable interactions with deoxycytidine triphosphate (dCTP), large reactions, and C1'-C1' distances for the N²-isobutyl-dG and N²-benzyl-dG PrimPol complexes, suggesting that the size of the adduct is a limiting factor for efficient TLS across minor groove adducts by PrimPol.

Helicene Structure between DNA and Cyanuric Acid: The Role of Noncovalent Interactions

The self-assembly between DNA and small organic molecules can expand the structural space of and introduce novel functionalities to DNA nanomaterials. In particular, it was demonstrated that poly(adenosine) DNA self-assembles with cyanuric acid (CA) to form a triplex helical structure. Previous molecular dynamics simulations showed that the DNA-CA assemblies adopt a novel noncovalent helicene structure that has a continuous helical hydrogen bond network. This article explores why the assemblies adopt the helicene geometry instead of an alternative planar hexameric rosette geometry. Analysis of the hydrogen bonding and stacking interaction energies indicates that constraining the system to the hexameric rosette geometry strains the hydrogen bonds without significantly improving the interaction energy. Molecular dynamics simulations for the assemblies between adenosine nucleosides and CA confirm that the formation of helicene structure is primarily driven by base-pair interactions and not because of the DNA backbone.

Quantum-based machine learning and AI models to generate force field parameters for drug-like small molecules

Force fields for drug-like small molecules play an essential role in molecular dynamics simulations and binding free energy calculations. In particular, the accurate generation of partial charges on small molecules is critical to understanding the interactions between proteins and drug-like molecules. However, it is a time-consuming process. Thus, we generated a force field for small molecules and employed a machine learning (ML) model to rapidly predict partial charges on molecules in less than a minute of time. We performed density functional theory (DFT) calculation for 31770 small molecules that covered the chemical space of drug-like molecules. The partial charges for the atoms in a molecule were predicted using an ML model trained on DFT-based atomic charges. The predicted values were comparable to the charges obtained from DFT calculations. The ML model showed high accuracy in the prediction of atomic charges for external test data sets. We also developed neural network (NN) models to assign atom types, phase angles and periodicities. All the models performed with high accuracy on test data sets. Our code calculated all the descriptors that were needed for the prediction of force field parameters and produced topologies for small molecules by combining results from ML and NN models. To assess the accuracy of the predicted force field parameters, we calculated solvation free energies for small molecules, and the results were in close agreement with experimental free energies. The AI-generated force field was effective in the fast and accurate generation of partial charges and other force field parameters for small drug-like molecules.

Insights into the binding mode of AS1411 aptamer to nucleolin

AS1411 aptamer can function as a recognition probe to detect the cell surface nucleolin overexpressed in cancer cells, however, little is known about their binding process. This study proposed a feasible binding mode for the first time and provided atomic-level descriptions for the high affinity and specific binding of AS1411. The binding pose predicted by docking was screened using knowledge-based criteria, and a microsecond molecular dynamics (MD) simulation showed the stable existence of the predicted structure in the solution. Structural analysis shows that the unique capping of the 5′ end of AS1411 provides the specific binding with RBD1, and the interactions of hydrogen bond, salt bridge, and water-mediated network between AS1411 and RBD1,2 stabilize the binding. The calculation of per-residue decomposition emphasizes the dominant contribution of van der Waals energy and critical residues are screened. Our study provides the molecular basis of this specific binding and can guide rational AS1411-based aptamers design. Further insights require tight collaborations between the experiments and in silico studies.

Molecular dynamics simulations of cyanine dimers attached to DNA Holliday junctions

Dye aggregates and their excitonic properties are of interest for their applications to organic photovoltaics, non-linear optics, and quantum information systems. DNA scaffolding has been shown to be effective at promoting the aggregation of dyes in a controllable manner. Specifically, isolated DNA Holliday junctions have been used to achieve strongly coupled cyanine dye dimers. However, the structural properties of the dimers and the DNA, as well as the role of Holliday junction isomerization are not fully understood. To study the dynamics of cyanine dimers in DNA, molecular dynamics simulations were carried out for adjacent and transverse dimers attached to Holliday junctions in two different isomers. It was found that dyes attached to adjacent strands in the junction exhibit stronger dye-DNA interactions and larger inter-dye separations compared to transversely attached dimers, as well as end-to-end arrangements. Transverse dimers exhibit lower inter-dye separations and more stacked configurations. Furthermore, differences in Holliday junction isomer are analyzed and compared to dye orientations. For transverse dyes exhibiting the smaller inter-dye separations, excitonic couplings were calculated and shown to be in agreement with experiment. Our results suggested that dye attachment locations on DNA Holliday junctions affect dye-DNA interactions, dye dynamics, and resultant dye orientations which can guide the design of DNA-templated cyanine dimers with desired properties.

Molecular dynamics simulations reveal dye attachment and DNA Holliday junction isomer effects on dye dimer orientations and excitonic couplings. These simulations can guide synthesis and experiments of dye-DNA structures for excitonic applications.

Efficient fluorescent recognition of ATP/GTP by a water-soluble bisquinolinium pyridine-2,6-dicarboxamide compound. Crystal structures, spectroscopic studies and interaction mode with DNA

The new dicationic pyridine-2,6-dicarboxamide-based compound 1 bearing two N-alkylquinolinium units was synthesized, structurally determined by single-crystal X-ray diffraction, and studied in-depth as a fluorescent receptor for nucleotides and inorganic phosphorylated anions in pure water. The addition of nucleotides to 1 at pH = 7.0 quenches its blue emission with a selective affinity towards adenosine 5'-triphosphate (ATP) and guanosine 5'-tripohosphate (GTP) over other nucleotides such CTP, UTP, ADP, AMP, dicarboxylates and inorganic anions. On the basis of the spectroscopic tools (1H, 31P NMR, UV-vis, fluorescence), MS measurements and DFT calculations, receptor 1 binds ATP with high affinity (log K = 5.04) through the simultaneous formation of strong hydrogen bonds and π-π interactions between the adenosine fragment and quinolinium ring with binding energy calculated in 8.7 kcal mol-1. High affinity for ATP/GTP is attributed to the high acidity of amides and preorganized rigid structure of 1. Receptor 1 is an order of magnitude more selective for ATP than GTP. An efficient photoinduced electron transfer quenching mechanism with simultaneous receptor-ATP complexation in both the excited and ground states is proposed. Additionally, multiple spectroscopic studies and molecular dynamics simulations showed that 1 can intercalate into DNA base pairs.

Microsecond MD Simulations to Explore the Structural and Energetic Differences between the Human RXRα-PPARγ vs. RXRα-PPARγ-DNA

The heterodimeric complex between retinoic X receptor alpha (RXRα) and peroxisome proliferator-activated receptor gamma (PPARγ) is one of the most important and predominant regulatory systems, controlling lipid metabolism by binding to specific DNA promoter regions. X-ray and molecular dynamics (MD) simulations have revealed the average conformation adopted by the RXRα-PPARγ heterodimer bound to DNA, providing information about how multiple domains communicate to regulate receptor properties. However, knowledge of the energetic basis of the protein-ligand and protein-protein interactions is still lacking. Here we explore the structural and energetic mechanism of RXRα-PPARγ heterodimer bound or unbound to DNA and forming complex with co-crystallized ligands (rosiglitazone and 9-cis-retinoic acid) through microsecond MD simulations, molecular mechanics generalized Born surface area binding free energy calculations, principal component analysis, the free energy landscape, and correlated motion analysis. Our results suggest that DNA binding alters correlated motions and conformational mobility within RXRα–PPARγ system that impact the dimerization and the binding affinity on both receptors. Intradomain correlated motions denotes a stronger correlation map for RXRα-PPARγ-DNA than RXRα-PPARγ, involving residues at the ligand binding site. In addition, our results also corroborated the greater role of PPARγ in regulation of the free and bound DNA state.

Is an Inductive Effect Explicit Account Required for Atomic Charges Aimed at Use within the Force Fields?

Polarization and inductive effects are the concepts that have been widely used in qualitative and even quantitative descriptions of experimentally observed properties in chemistry. The polarization effect has proven to be important in cases of biomolecular modeling though still the vast majority of molecular simulations use the classical non-polarizable force fields. In the last few decades, a lot of effort has been put into promoting the polarization effect and incorporating it into modern force fields and charge calculation methods. In contrast, the inductive effect has not attracted such attention and is effectively absent in both classic and modern force fields. Thus, a question is whether this difference corresponds to the difference in the physical significance of the effects and their explicit account, or is an artifact that should be corrected in the next generation of force fields. The significance of the electronic effects is studied in this paper through the prism of performance of specific models for atomic charge calculation that take into explicit account a nested set of effects: the formal charge, the nearest neighbors, the inductive effect, and finally the model, which takes into account all effects, which are possible to account for using atomic charges. The specific choice for the methods is the following: formal charges, MMFF94 bond charge increments, Dynamic Electronegativity Relaxation (DENR), and RESP. We propose a special scheme for the separate estimation of each particular effect contribution. By pairwise comparing the residual molecular electrostatic potential (MEP) errors of those charge models (aimed at best reproducing the quantum chemical reference MEP), we sequentially revealed how the account of each effect contributes to the better-quality MEP reproduction. The following relative importance of effects was estimated; thus, the natural hierarchy of the effects was established. First, the account of formal charges is of primordial importance. Second, the nearest neighbors account is the next in significance. Third, the explicit account of inductive effect in empirical charge calculation schemes was shown to significantly─both qualitatively and quantitatively─improve the quality of MEP reproduction. Fourth, the contribution of polarization is indirectly assessed. Surprisingly, it is of the order of magnitude of the inductive effect even for the molecular systems, for which it is anticipated to be more significant. Finally, the relative importance of anisotropic effects in neutral molecules was additionally reviewed.

H2B Lys34 Ubiquitination Induces Nucleosome Distortion to Stimulate Dot1L Activity

Ubiquitination-dependent histone crosstalk plays critical roles in chromatin-associated processes and is highly associated with human diseases. Mechanism studies of the crosstalk have been of the central focus. Here our study on the crosstalk between H2BK34ub and Dot1L-catalyzed H3K79me suggests a novel mechanism of ubiquitination-induced nucleosome distortion to stimulate the activity of an enzyme. We determined the cryo-electron microscopy structures of Dot1L-H2BK34ub nucleosome complex and the H2BK34ub nucleosome alone. The structures reveal that H2BK34ub induces an almost identical orientation and binding pattern of Dot1L on nucleosome as H2BK120ub, which positions Dot1L for the productive conformation through direct ubiquitin-enzyme contacts. However, H2BK34-anchored ubiquitin does not directly interact with Dot1L as occurs in the case of H2BK120ub, but rather induces DNA and histone distortion around the modified site. Our findings establish the structural framework for understanding the H2BK34ub-H3K79me trans-crosstalk and highlight the diversity of mechanisms for histone ubiquitination to activate chromatin-modifying enzymes.

Applications of Molecular Dynamics Simulation in Protein Study

Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.

Structural Analysis of the Poly(thymidine)-Melamine Assembly

Hydrogen bonding between the DNA nucleobases and small organic molecules, such as melamine, is a new strategy for the design of novel DNA materials. Poly(thymidine) DNA and melamine self-assemble into a duplex structure containing two antiparallel DNA strands hydrogen bonded to central melamine units. In this Article, molecular dynamics simulations rationalize the observed antiparallel duplex structure. Alternative duplex and triplex structures with parallel and antiparallel strand orientations are shown to be unstable because of the increase in unfavorable interactions between the DNA backbones.

Early Stages of RNA-Mediated Conversion of Human Prions

Prion diseases are characterized by the conversion of prion proteins from a PrPC fold into a disease-causing PrPSC form that is self-replicating. A possible agent to trigger this conversion is polyadenosine RNA, but both mechanism and pathways of the conversion are poorly understood. Using coarse-grained molecular dynamic simulations we study the time evolution of PrPC over 600 μs. We find that both the D178N mutation and interacting with polyadenosine RNA reduce the helicity of the protein and encourage formation of segments with strand-like motifs. We conjecture that these transient β-strands nucleate the conversion of the protein to the scrapie conformation PrPSC.

RelEB3 toxin-antitoxin system of Salmonella Typhimurium with a ribosome-independent toxin and a mutated non-neutralising antitoxin

The RelEB3 toxin-antitoxin (TA) system of Salmonella enterica subsp. enterica serovar Typhimurium consists of a RelE3 toxin which suppresses bacterial growth, but its RelB3 antitoxin does not neutralise the toxin. The relEB3 operon is widespread in Proteobacteria and is related to higBA2 from Vibrio cholerae. In contrast to the ribosome-dependent HigB2 toxin, however, the RelE3 toxin degraded free RNA independently of the ribosome. A basic loop possibly involved in HigB2's binding to the ribosome is shortened in RelE3, which instead contains a uniquely conserved R51 important for RelE3's toxicity. The RelB3 antitoxin, meanwhile, specifically recognised the CACCTGGTG palindromic motif in the promoter site. RelB3 contains P14 which is conserved as Ala in most homologues, and mutating P14 to Ala enabled the antitoxin to bind to RelE3 and restored bacterial growth. The P14 RelB3 variant, which most likely arose by a point mutation in a recent ancestor of S. Typhimurium and closely related serovars, could have possibly provided the bacteria with a faster response to stress, and might have spread to other serovars through homologous recombination.

Preserving the Integrity of Empirical Force Fields

Generalized force fields (FFs) act as extensions to biomolecular FFs to provide a wide coverage of organic molecules. However, their precise application to an arbitrary molecule presents a separate challenge. We show that MATCH assigns different atom types and bonded and nonbonded parameters than CGenFF, and the AM1-BCC charge model, commonly used with GAFF/GAFF2, does not exactly reproduce the performance of the RESP charge model. The results indicate the need for caution when employing FFs to ensure their integrity with respect to their implementation and validation.

Stability and Existence of Noncanonical I-motif DNA Structures in Computer Simulations Based on Atomistic and Coarse-Grained Force Fields

Cytosine-rich DNA sequences are able to fold into noncanonical structures, in which semi-protonated cytosine pairs develop extra hydrogen bonds, and these bonds are responsible for the overall stability of a structure called the i-motif. The i-motif can be formed in many regions of the genome, but the most representative is the telomeric region in which the CCCTAA sequences are repeated thousands of times. The ability to reverse folding/unfolding in response to pH change makes the above sequence and i-motif very promising components of nanomachines, extended DNA structures, and drug carriers. Molecular dynamics analysis of such structures is highly beneficial due to direct insights into the microscopic structure of the considered systems. We show that Amber force fields for DNA predict the stability of the i-motif over a long timescale; however, these force fields are not able to predict folding of the cytosine-rich sequences into the i-motif. The reason is the kinetic partitioning of the folding process, which makes the transitions between various intermediates too time-consuming in atomistic force field representation. Application of coarse-grained force fields usually highly accelerates complex structural transitions. We, however, found that three of the most popular coarse-grained force fields for DNA (oxDNA, 3SPN, and Martini) were not able to predict the stability of the i-motif structure. Obviously, they were not able to accelerate the folding of unfolded states into an i-motif. This observation must be strongly highlighted, and the need to develop suitable extensions of coarse-grained force fields for DNA is pointed out. However, it will take a great deal of effort to successfully solve these problems.

Exploring the Interaction of G-quadruplex Binders with a (3 + 1) Hybrid G-quadruplex Forming Sequence within the PARP1 Gene Promoter Region

The enzyme PARP1 is an attractive target for cancer therapy, as it is involved in DNA repair processes. Several PARP1 inhibitors have been approved for clinical treatments. However, the rapid outbreak of resistance is seriously threatening the efficacy of these compounds, and alternative strategies are required to selectively regulate PARP1 activity. A noncanonical G-quadruplex-forming sequence within the PARP1 promoter was recently identified. In this study, we explore the interaction of known G-quadruplex binders with the G-quadruplex structure found in the PARP gene promoter region. The results obtained by NMR, CD, and fluorescence titration, also confirmed by molecular modeling studies, demonstrate a variety of different binding modes with small stabilization of the G-quadruplex sequence located at the PARP1 promoter. Surprisingly, only pyridostatin produces a strong stabilization of the G-quadruplex-forming sequence. This evidence makes the identification of a proper (3+1) stabilizing ligand a challenging goal for further investigation.

Transient Hoogsteen Base Pairs Observed in Unbiased Molecular Dynamics Simulations of DNA

Duplex DNA is modeled as canonical B-DNA displaying the characteristic Watson-Crick base pairs. A less common and short-lived pairing of the nucleobases is the Hoogsteen (HG) conformation. The low population of the HG base pairs (<1%) necessitates extended sampling times in order to analyze through unbiased molecular dynamics (MD) simulations. We have discovered that with extended sampling times using multiple independent copies of an 18-mer sequence, the MD trajectories reproduce the expected and transient HG base pairing. Consistent with experiment, the percentage of the HG events are within the range of ∼0.1-1.0% over the combined aggregate sampling time of more than 3.6 ms. We present the reliability of the current AMBER set of nucleic acid force fields and tools to accurately simulate naturally occurring base-pairing and opening events without any bias or restraints. The mechanism consists of base pair fraying, flipping of the purine, and reformation with HG base pairs.

Strategies for Controlling the Spatial Orientation of Single Molecules Tethered on DNA Origami Templates Physisorbed on Glass Substrates: Intercalation and Stretching

Nanoarchitectural control of matter is crucial for next-generation technologies. DNA origami templates are harnessed to accurately position single molecules; however, direct single molecule evidence is lacking regarding how well DNA origami can control the orientation of such molecules in three-dimensional space, as well as the factors affecting control. Here, we present two strategies for controlling the polar (θ) and in-plane azimuthal (ϕ) angular orientations of cyanine Cy5 single molecules tethered on rationally-designed DNA origami templates that are physically adsorbed (physisorbed) on glass substrates. By using dipolar imaging to evaluate Cy5's orientation and super-resolution microscopy, the absolute spatial orientation of Cy5 is calculated relative to the DNA template. The sequence-dependent partial intercalation of Cy5 is discovered and supported theoretically using density functional theory and molecular dynamics simulations, and it is harnessed as our first strategy to achieve θ control for a full revolution with dispersion as small as ±4.5°. In our second strategy, ϕ control is achieved by mechanically stretching the Cy5 from its two tethers, being the dispersion ±10.3° for full stretching. These results can in principle be applied to any single molecule, expanding in this way the capabilities of DNA as a functional templating material for single-molecule orientation control. The experimental and modeling insights provided herein will help engineer similar self-assembling molecular systems based on polymers, such as RNA and proteins.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

Deciphering the functional mechanism of zinc ions of PARP1 binding with single strand breaks and double strand breaks

Poly(ADP-ribose)polymerase 1 (PARP1) is a key target for the treatment of cancer-related diseases, and plays an important role in biological processes such as DNA repair, regulating a variety of metabolic and signal transduction processes. Understanding the dynamic binding mechanisms between each domain of PARP1 and DNA is of great significance to deepen the understanding on the function of PARP1 and to facilitate the design of inhibitors. Herein, strategies such as classical molecular dynamics simulation, conformational analysis, binding free energy calculation and energy decomposition were used to shed light on the binding mechanisms of different DNA binding domains (DBDs, including ZnF1, ZnF2 and ZnF3) in PARP1 with DNA and on the influences of zinc ions on the binding process. On one hand, during binding with DNA, ZnF2 tends to expand its space to identify the DNA damage sites and ZnF1/ZnF2 recognizes the interfaces on both sides of DNA damage rather than one side during the process of DNA repair. More importantly, the stable secondary structure of L₂ of ZnF2 (PRO146 to MET153) is the key conformational change for ZnF1 and ZnF2 to recognize DNA damage. Meanwhile, ZnF3 has little effect on the binding mechanisms of PARP1. On the other hand, for the structural differences of DBD domains, zinc ions in ZnF1 and ZnF2 (Zn1 and Zn2) have an impact not only on the conformational changes of PARP1, but also on the conformational changes brought by the interaction of double strand breaks (DSB) and single strand breaks (SSB). And meanwhile, Zn3 also has little effect on ZnF3 for the system of ZnF3/DSB. The findings presented in this work deepen the understanding on the functional mechanism of PARP1 and provide a theoretical basis for further study on the interaction between different inhibitors and DBD domains to design more potential inhibitors.

Poly(ADP-ribose)polymerase 1 (PARP1) is a key target for treatment of cancer-related diseases. Detailed structural changes DBD in PARP1 during the binding process with DNA were investigated and the dynamic conformational differences of DBD caused by zinc ions were revealed.

Guiding the folding of G-quadruplexes through loop residue interactions

A G-rich sequence was designed to allow folding into either a stable parallel or hybrid-type topology. With the parent sequence featuring coexisting species, various related sequences with single and double mutations and with a shortened central propeller loop affected the topological equilibrium. Two simple modifications, likewise introduced separately to all sequences, were employed to lock folds into one of the topologies without noticeable structural alterations. The unique combination of sequence mutations, high-resolution NMR structural information, and the thermodynamic stability for both topological competitors identified critical loop residue interactions. In contrast to first loop residues, which are mostly disordered and exposed to solvent in both propeller and lateral loops bridging a narrow groove, the last loop residue in a lateral three-nucleotide loop is engaged in stabilizing stacking interactions. The propensity of single-nucleotide loops to favor all-parallel topologies by enforcing a propeller-like conformation of an additional longer loop is shown to result from their preference in linking two outer tetrads of the same tetrad polarity. Taken together, the present studies contribute to a better structural and thermodynamic understanding of delicate loop interactions in genomic and artificially designed quadruplexes, e.g. when employed as therapeutics or in other biotechnological applications.

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes

Although molecular dynamics (MD) simulations have been used extensively to study the structural dynamics of proteins, the role of MD simulation in studies of nucleic acid based systems has been more limited. One contributing factor to this disparity is the historically lower level of accuracy of the physical models used in such simulations to describe interactions involving nucleic acids. By modifying nonbonded and torsion parameters of a force field from the Amber family of models, we recently developed force field parameters for RNA that achieve a level of accuracy comparable to that of state-of-the-art protein force fields. Here we report force field parameters for DNA, which we developed by transferring nonbonded parameters from our recently reported RNA force field and making subsequent adjustments to torsion parameters. We have also modified the backbone charges in both the RNA and DNA parameter sets to make the treatment of electrostatics compatible with our recently developed variant of the Amber protein and ion force field. We name the force field resulting from the union of these three parameter sets (the new DNA parameters, the revised RNA parameters, and the existing protein and ion parameters) DES-Amber. Extensive testing of DES-Amber indicates that it can describe the thermal stability and conformational flexibility of single- and double-stranded DNA systems with a level of accuracy comparable to or, especially for disordered systems, exceeding that of state-of-the-art nucleic acid force fields. Finally, we show that, in certain favorable cases, DES-Amber can be used for long-timescale simulations of protein–nucleic acid complexes.

Twisting DNA by salt

The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.

Studying the Dynamics of a Complex G-Quadruplex System: Insights into the Comparison of MD and NMR Data

Molecular dynamics (MD) simulations are coming of age in the study of nucleic acids, including specific tertiary structures such as G-quadruplexes. While being precious for providing structural and dynamic information inaccessible to experiments at the atomistic level of resolution, MD simulations in this field may still be limited by several factors. These include the force fields used, different models for ion parameters, ionic strengths, and water models. We address various aspects of this problem by analyzing and comparing microsecond-long atomistic simulations of the G-quadruplex structure formed by the human immunodeficiency virus long terminal repeat (HIV LTR)-III sequence for which nuclear magnetic resonance (NMR) structures are available. The system is studied in different conditions, systematically varying the ionic strengths, ion numbers, and water models. We comparatively analyze the dynamic behavior of the G-quadruplex motif in various conditions and assess the ability of each simulation to satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints and structural parameters. The conditions taking into account K⁺-ions to neutralize the system charge, mimicking the intracellular ionic strength, and using the four-atom water model are found to be the best in reproducing the experimental NMR constraints and data. Our analysis also reveals that in all of the simulated environments residues belonging to the duplex moiety of HIV LTR-III exhibit the highest flexibility.

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M−1cm−1) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

DNA Sequence Specificity Reveals a Role of the HLTF HIRAN Domain in the Recognition of Trinucleotide Repeats

DNA damage tolerance (DDT) pathways enable cells to cope with a variety of replication blocks that threaten their ability to complete DNA replication. Helicase-like transcription factor (HLTF) plays a central role in the error-free DDT pathway, template switching (TS), by serving as a ubiquitin ligase to polyubiquitinate the DNA sliding clamp PCNA, which promotes TS initiation. HLTF also serves as an ATP-dependent DNA translocase facilitating replication fork remodeling. The HIP116, Rad5p N-terminal (HIRAN) domain of HLTF specifically recognizes the unmodified 3'-end of single-stranded DNA (ssDNA) at stalled replication forks to promote fork regression. Several crystal structures of the HIRAN domain in complex with ssDNA have been reported; however, optimal ssDNA sequences for high-affinity binding with the domain have not been described. Here we elucidated DNA sequence preferences of HLTF HIRAN through systematic studies of its binding to ssDNA substrates using fluorescence polarization assays and a computational analysis of the ssDNA:HIRAN interaction. These studies reveal that the HLTF HIRAN domain preferentially recognizes a (T/C)TG sequence at the 3'-hydroxyl ssDNA end, which occurs in the CTG trinucleotide repeat (TNR) regions that are susceptible to expansion and deletion mutations identified in neuromuscular and neurodegenerative disorders. These findings support a role for HLTF in maintaining the stability of difficult to replicate TNR microsatellite regions.

Novel eGZ-motif formed by regularly extruded guanine bases in a left-handed Z-DNA helix as a major motif behind CGG trinucleotide repeats

The expansion of d(CGG) trinucleotide repeats (TRs) lies behind several important neurodegenerative diseases. Atypical DNA secondary structures have been shown to trigger TR expansion: their characterization is important for a molecular understanding of TR disease. CD spectroscopy experiments in the last decade have unequivocally demonstrated that CGG runs adopt a left-handed Z-DNA conformation, whose features remain uncertain because it entails accommodating GG mismatches. In order to find this missing motif, we have carried out molecular dynamics (MD) simulations to explore all the possible Z-DNA helices that potentially form after the transition from B- to Z-DNA. Such helices combine either CpG or GpC Watson-Crick steps in Z-DNA form with GG-mismatch conformations set as either intrahelical or extrahelical; and participating in BZ or ZZ junctions or in alternately extruded conformations. Characterization of the stability and structural features (especially overall left-handedness, higher-temperature and steered MD simulations) identified two novel Z-DNA helices: the most stable one displays alternately extruded Gs, and is followed by a helix with symmetrically extruded ZZ junctions. The G-extrusion favors a seamless stacking of the Watson-Crick base pairs; extruded Gs favor syn conformations and display hydrogen-bonding and stacking interactions. Such conformations could have the potential to hijack the MMR complex, thus triggering further expansion.

Never Cared for What They Do: High Structural Stability of Guanine-Quadruplexes in the Presence of Strand-Break Damage

DNA integrity is an important factor that assures genome stability and, more generally, the viability of cells and organisms. In the presence of DNA damage, the normal cell cycle is perturbed when cells activate their repair processes. Although efficient, the repair system is not always able to ensure complete restoration of gene integrity. In these cases, mutations not only may occur, but the accumulation of lesions can either lead to carcinogenesis or reach a threshold that induces apoptosis and programmed cell death. Among the different types of DNA lesions, strand breaks produced by ionizing radiation are the most toxic due to the inherent difficultly of repair, which may lead to genomic instability. In this article we show, by using classical molecular simulation techniques, that compared to canonical double-helical B-DNA, guanine-quadruplex (G4) arrangements show remarkable structural stability, even in the presence of two strand breaks. Since G4-DNA is recognized for its regulatory roles in cell senescence and gene expression, including oncogenes, this stability may be related to an evolutionary cellular response aimed at minimizing the effects of ionizing radiation.

Accurate Sequence-Dependent Coarse-Grained Model for Conformational and Elastic Properties of Double-Stranded DNA

We introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar, at the base moiety, and at the phosphate group, respectively. The sequence dependence is included by considering a step-dependent parametrization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from the literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist–stretch and twist–bend couplings, the persistence length, and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of single-molecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines. Particularly, an implementation of the model in LAMMPS is made available on an online repository to ease its usage within the DNA research community.

Diameter Dependent Melting and Softening of dsDNA Under Cylindrical Confinement

Carbon nanotubes (CNTs) are considered promising candidates for biomolecular confinement, including DNA encapsulation for gene delivery. Threshold values of diameters have been reported for double-stranded DNA (dsDNA) encapsulation inside CNTs. We have performed all-atom molecular dynamics (MD) simulations of dsDNAs confined inside single-walled CNTs (SWCNTs) at the physiologically relevant temperature of 300 K. We found that the dsDNA can be confined without being denatured only when the diameter of the SWCNT exceeds a threshold value. Below this threshold diameter, the dsDNA gets denatured and melts even at the temperature of 300 K. Our simulations using SWCNTs with chirality indices (20,20) to (30,30) at 300 K found the critical diameter to be 3.25 nm (corresponding to (24,24) chirality). Analyses of the hydrogen bonds (H-bonds), Van der Walls (VdW) energy, and other inter-base interactions show drastic reduction in the number of H-bonds, VdW energy, and electrostatic energies between the bases of dsDNA when it is confined in narrower SWCNTs (up to diameter of 3.12 nm). On the other hand, the higher interaction energy between the dsDNA and the SWCNT surface in narrower SWCNTs assists in the melting of the dsDNA. Electrostatic mapping and hydration status analyses show that the dsDNA is not adequately hydrated and the counter ion distribution is not uniform below the critical diameter of the SWCNT. As properly hydrated counter ions provide stability to the dsDNA, we infer that the inappropriate hydration of counter ions and their non-uniform distribution around the dsDNA cause the melting of the dsDNA inside SWCNTs of diameter below the critical value of 3.25 nm. For confined dsDNAs that do not get denatured, we computed their elastic properties. The persistence length of dsDNA was found to increase by a factor of about two and the torsional stiffness by a factor of 1.5 for confinement inside SWCNTs of diameters up to 3.79 nm, the stretch modulus also following nearly the same trend. Interestingly, for higher diameters of SWCNT, 3.79 nm and above, the dsDNA becomes more flexible, demonstrating that the mechanical properties of the dsDNA under cylindrical confinement depend non-monotonically on the confinement diameter.

Discovery, Topo I inhibitory activity and mechanism evaluation of two novel cytisine-type alkaloid dimers from the seeds of Sophora alopecuroides L

Alopecurosines A and B (CMs 1 and 2, respectively) are two novel cytisine-type alkaloid dimers first isolated from the aerial parts of Sophora alopecuroides L. CMs 1 and 2 are new dimeric alkaloids whose piperidine matrine ring is cleaved and connected via the N'-1 bond. Their chemical structures have been confirmed by IR, UV, HR-ESI-MS, and NMR. Preliminary screening shows that they have topoisomerase I (Topo I)-based anti-tumor activity. Their Topo I inhibitory activities and mechanism have been evaluated by agarose gel electrophoresis assay and a molecular docking study. The results show that the inhibition rate of CM 1 is 82.26% at 1 mM concentration and that it exhibits significantly Topo I inhibitory activity. Further research has illustrated that CMs 1 and 2 exert inhibitory activity by stabilising the Topo I-DNA cleavage complex, implying that they have the potential to be developed as novel Topo I inhibitors.

Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields

Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

Rad5 HIRAN domain: Structural insights into its interaction with ssDNA through molecular modeling approaches

The Rad5 protein is an SWI/SNF family ubiquitin ligase that contains an N-terminal HIRAN domain and a RING C3HC4 motif. The HIRAN domain is critical for recognition of the stalled replication fork during the replication process and acts as a sensor to initiate the damaged DNA checkpoint. It is a conserved domain widely distributed in eukaryotic organisms and is present in several DNA-binding proteins from all kingdoms. Here we showed that distant species have important differences in key residues that affect affinity for ssDNA. Based on these findings, we hypothesized that different HIRAN domains might affect fork reversal and translesion synthesis through different metabolic processes. To address this question, we predicted the tertiary structure of both yeast and human HIRAN domains using molecular modeling. Structural dynamics experiments showed that the yeast HIRAN domain exhibited higher structural denaturation than its human homolog, although both domains became stable in the presence of ssDNA. Analysis of atomic contacts revealed that a greater number of interactions between the ssDNA nucleotides and the Rad5 domain are electrostatic. Taken together, these results provide new insights into the molecular mechanism of the HIRAN domain of Rad5 and may guide us to further elucidate differences in the ancient eukaryotes HIRAN sequences and their DNA affinity. Communicated by Ramaswamy H. Sarma.

Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E

Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.

The molecular origin of the electrostatic gating of single-molecule field-effect biosensors investigated by molecular dynamics simulations

Field-effect biosensors (bioFETs) offer a novel way to measure the kinetics of biomolecular events such as protein function and DNA hybridization at the single-molecule level on a wide range of time scales. These devices generate an electrical current whose fluctuations are correlated to the kinetics of the biomolecule under study. BioFETs are indeed highly sensitive to changes in the electrostatic potential (ESP) generated by the biomolecule. Here, using all-atom solvent explicit molecular dynamics simulations, we further investigate the molecular origin of the variation of this ESP for two prototypical cases of proteins or nucleic acids attached to a carbon nanotube bioFET: the function of the lysozyme protein and the hybridization of a 10-nt DNA sequence, as previously done experimentally. Our results show that the ESP changes significantly on the surface of the carbon nanotube as the state of these two biomolecules changes. More precisely, the ESP distributions calculated for these molecular states explain well the magnitude of the conductance fluctuations measured experimentally. The dependence of the ESP with salt concentration is found to agree with the reduced conductance fluctuations observed experimentally for the lysozyme, but to differ for the case of DNA, suggesting that other mechanisms might be at play in this case. Furthermore, we show that the carbon nanotube does not impact significantly the structural stability of the lysozyme, corroborating that the kinetic rates measured using bioFETs are similar to those measured by other techniques. For DNA, we find that the structural ensemble of the single-stranded DNA is significantly impacted by the presence of the nanotube, which, combined with the ESP analysis, suggests a stronger DNA-device interplay. Overall, our simulations strengthen the comprehension of the inner working of field-effect biosensors used for single-molecule kinetics measurements on proteins and nucleic acids.

DNA methylation cues in nucleosome geometry, stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

Searching for Low Probability Opening Events in a DNA Sliding Clamp

The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.

Nucleoside Analogues with a Seven-Membered Sugar Ring: Synthesis and Structural Compatibility in DNA-RNA Hybrids

Herein we describe results on the pairing properties of synthetic DNA and RNA oligonucleotides that contain nucleotide analogues with a 7-membered sugar ring (oxepane nucleotides). Specifically, we describe the stereoselective synthesis of a set of three oxepane thymine nucleosides (OxT), their conversion to phosphoramidite derivatives, and their use in solid-phase synthesis to yield chimeric OxT-DNA and OxT-RNA strands. The different regioisomeric OxT phosphoramidites allowed for positional variations of the phosphate bridge and assessment of duplex stability when the oxepane nucleotides were incorporated in dsDNA, dsRNA, and DNA-RNA hybrids. Little to no destabilization was observed when two of the three regioisomeric OxT units were incorporated in the DNA strand of DNA-RNA hybrids, a remarkable result considering the dramatically different structure of oxepanes in comparison to 2'-deoxynucleosides. Extensive molecular modeling and dynamics studies further revealed the various structural features responsible for the tolerance of both OxT modifications in DNA-RNA duplexes, such as base-base stacking and sugar-phosphate H-bond interactions. These studies suggest that oxepane nucleotide analogues may find applications in synthetic biology, where synthetic oligonucleotides can be used to create new tools for biotechnology and medicine.

Riboflavin Stabilizes Abasic, Oxidized G-Quadruplex Structures

The G-quadruplex is a noncanonical fold of DNA commonly found at telomeres and within gene promoter regions of the genome. These guanine-rich sequences are highly susceptible to damages such as base oxidation and depurination, leading to abasic sites. In the present work, we address whether a vacancy, such as an abasic site, in a G-quadruplex serves as a specific ligand recognition site. When the G-tetrad is all guanines, the vacant (abasic) site is recognized and bound by free guanine nucleobase. However, we aim to understand whether the preference for a specific ligand recognition changes with the presence of a guanine oxidation product 8-oxo-7,8-dihydroguanine (OG) adjacent to the vacancy in the tetrad. Using molecular dynamics simulation, circular dichroism, and nuclear magnetic resonance, we examined the ability for riboflavin to stabilize abasic site-containing G-quadruplex structures. Through structural and free energy binding analysis, we observe riboflavin’s ability to stabilize an abasic site-containing G-quadruplex only in the presence of an adjacent OG-modified base. Further, when compared to simulation with the vacancy filled by free guanine, we observe that the free guanine nucleobase is pushed outside of the tetrad by OG to interact with other parts of the structure, including loop residues. These results support the preference of riboflavin over free guanine to fill an OG-adjacent G-quadruplex abasic vacancy.

Backbone Hydrocarbon-Constrained Nucleic Acids Modulate Hybridization Kinetics for RNA

The binding affinity of therapeutic oligonucleotides (ONs) for their cognate RNA is determined by the rates of association (k_a) and dissociation (k_d). Single-stranded ONs are highly flexible and can adopt multiple conformations in solution, some of which may not be conducive for hybridization. We investigated if restricting rotation around the sugar-phosphate backbone, by tethering two adjacent backbone phosphonate esters using hydrocarbon bridges, can modulate hybridization kinetics of the modified ONs for complementary RNA. Given the large number of possible analogues with different tether lengths and configurations at the phosphorus atoms, we employed molecular dynamic simulations to optimize the size of the hydrocarbon bridge to guide the synthetic efforts. The backbone-constrained nucleotide trimers with stereodefined configurations at the contiguous backbone phosphorus atoms were assembled using a ring-closing metathesis reaction, then incorporated into oligonucleotides by an in situ synthesis of the phosphoramidites followed by coupling to solid supports. Evaluation of the modified oligonucleotides revealed that 15-membered macrocyclic-constrained analogues displayed similar or slightly improved on-rates but significantly increased off-rates compared to unmodified DNA ONs, resulting in reduced duplex stability. In contrast, LNA ONs with conformationally preorganized furanose rings showed similar on-rates to DNA ONs but very slow off-rates, resulting in net improvement in duplex stability. Furthermore, the experimental data generally supported the molecular dynamics simulation results, suggesting that this strategy can be used as a predictive tool for designing the next generation of constrained backbone ON analogues with improved hybridization properties.