Urea Mimics Nucleobases by Preserving the Helical Integrity of B-DNA Duplexes via Hydrogen Bonding and Stacking Interactions.

Urea lesions are formed in DNA because of free radical damage of the thymine base, and their occurrence in DNA blocks DNA polymerases, which has deleterious consequences. Recently, it has been shown that urea is capable of forming hydrogen bonding and stacking interactions with nucleobases, which are responsible for the unfolding of RNA in aqueous urea. Base pairing and stacking are inherent properties of nucleobases; because urea is able to form both, this study attempts to investigate if urea can mimic nucleobases in the context of nucleic acid structures by examining the effect of introducing urea lesions complementary to the four different nucleobases on the overall helical integrity of B-DNA duplexes and their thermodynamic stabilities using molecular dynamics (MD) simulations. The MD simulations resulted in stable duplexes without significant changes in the global B-DNA conformation. The urea lesions occupy intrahelical positions by forming hydrogen bonds with nitrogenous nucleobases, in agreement with experimental results. Furthermore, these urea lesions form hydrogen bonding and stacking interactions with other nucleobases of the same and partner strands, analogous to nucleobases in typical B-DNA duplexes. Direct hydrogen bond interactions are observed for the urea-purine pairs within DNA duplexes, whereas two different modes of pairing, namely, direct hydrogen bonds and water-mediated hydrogen bonds, are observed for the urea-pyrimidine pairs. The latter explains the complexities involved in interpreting the experimental nuclear magnetic resonance data reported previously. Binding free energy calculations were further performed to confirm the thermodynamic stability of the urea-incorporated DNA duplexes with respect to pure duplexes. This study suggests that urea mimics nucleobases by pairing opposite all four nucleobases and maintains the overall structure of the B-DNA duplexes.

Ligand-Induced Stabilization of a Duplex-like Architecture Is Crucial for the Switching Mechanism of the SAM-III Riboswitch.

Riboswitches are structured RNA motifs that control gene expression by sensing the concentrations of specific metabolites and make up a promising new class of antibiotic targets. S-Adenosylmethionine (SAM)-III riboswitch, mainly found in lactic acid bacteria, is involved in regulating methionine and SAM biosynthetic pathways. SAM-III riboswitch regulates the gene expression by switching the translation process on and off with respect to the absence and presence of the SAM ligand, respectively. In this study, an attempt is made to understand the key conformational transitions involved in ligand binding using atomistic molecular dynamics (MD) simulations performed in an explicit solvent environment. G26 is found to recognize the SAM ligand by forming hydrogen bonds, whereas the absence of the ligand leads to opening of the binding pocket. Consistent with experimental results, the absence of the SAM ligand weakens the base pairing interactions between the nucleobases that are part of the Shine-Dalgarno (SD) and anti-Shine-Dalgarno (aSD) sequences, which in turn facilitates recognition of the SD sequence by ribosomes. Detailed analysis reveals that a duplex-like structure formed by nucleotides from different parts of the RNA and the adenine base of the ligand is crucial for the stability of the completely folded state in the presence of the ligand. Previous experimental studies have shown that the SAM-III riboswitch exists in equilibrium between the unfolded and partially folded states in the absence of the ligand, which completely folds upon binding of the ligand. Comparison of the results presented here to the available experimental data indicates the structures obtained using the MD simulations resemble the partially folded state. Thus, this study provides a detailed understanding of the fully and partially folded structures of the SAM-III riboswitch in the presence and absence of the ligand, respectively. This study hypothesizes a dual role for the SAM ligand, which facilitates conformational switching between partially and fully folded states by forming a stable duplex-like structure and strengthening the interactions between SD and aSD nucleotides.

Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers.

A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (Rg), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrimer as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study.

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.

Double zipper helical assembly of deoxyoligonucleotides: mutual templating and chiral imprinting to form hybrid DNA ensembles.

Herein, the conventional and unconventional hydrogen bonding potential of adenine in APA for double zipper helical assembly of deoxyoligonucleotides is demonstrated under ambient conditions. The quantum mechanical calculations supported the formation of hybrid DNA ensembles.

DNA-RNA hybrid duplexes with decreasing pyrimidine content in the DNA strand provide structural snapshots for the A- to B-form conformational transition of nucleic acids.

DNA-RNA hybrids are heterogeneous nucleic acid duplexes consisting of a DNA strand and a RNA strand, and are formed as key intermediates in many important biological processes. They serve as substrates for the RNase H enzymatic activity, which has been exploited for several biomedical technologies such as antiviral and antisense therapies. To understand the relation of structural properties with the base composition in DNA-RNA hybrids, molecular dynamics (MD) simulations were performed on selected model systems by systematically varying the deoxypyrimidine (dPy) content from 0 to 100% in the DNA strand. The results suggest that the hybrid duplex properties are highly dependent on their deoxypyrimidine content of the DNA strand. However, such variations are not seen in their corresponding pure DNA and RNA duplex counterparts. It is also noticed that the systematic variation in deoxypyrimidine content of hybrids leads to gradual transformation between B- and A-form nucleic acid structures. Binding free energy calculations explain the previous experimental findings that the hybrids with high deoxypyrimidine content (>50%) are more stable than their respective pure counterparts. Pseudorotation angles, minor groove widths, phosphodiester angles, and glycosidic dihedral angles exhibit gradual A- to A/B-like conformation with decreasing deoxypyrimidine content. Based on extensive analysis, possible factors that affect RNase H enzymatic activity on hybrid duplexes with high dPy composition are proposed.

Atomistic investigation of the effect of incremental modification of deoxyribose sugars by locked nucleic acid (ß-D-LNA and α-L-LNA) moieties on the structures and thermodynamics of DNA-RNA hybrid duplexes.

Chemically modified oligonucleotides offer many possibilities in utilizing their special features for a vast number of applications in nucleic acid based therapies and synthetic molecular biology. Locked nucleic acid analogues (α-/ß-LNA) are modifications having an extra ring of 2'-O,4'-C-methylene group in the furanose sugar. LNA strands have been shown to exhibit high binding affinity toward RNA and DNA strands, and the resultant duplexes show significantly high melting temperatures. In the present study, molecular dynamics (MD) simulations were performed on DNA-RNA hybrid duplexes by systematically modifying their deoxyribose sugars with locked nucleic acid analogues. Several geometrical and energetic analyses were performed using principal component (PCA) analysis and binding free energy methods to understand the consequence of incorporated isomeric LNA modifications on the structure, dynamics, and stability of DNA-RNA hybrid duplex. The ß-modification systematically changes the conformation of the DNA-RNA hybrid duplex whereas drastic changes are observed for α-modification. The fully modified duplexes have distinct properties compared to partial and unmodified duplexes, and the partly modified duplexes have properties intermediate to full strand and unmodified duplexes. The distribution of BI versus BII populations suggests that backbone rearrangement is minimal for ß-LNA modification in order to accommodate it in duplexes whereas extensive backbone rearrangement is necessary in order to incorporate α-LNA modification which subsequently alters the energetic and structural properties of the duplexes. The simulation results also suggest that the alteration of DNA-RNA hybrid properties depends on the position of modification and the gap between the modifications.

Structures, dynamics, and stabilities of fully modified locked nucleic acid (ß-D-LNA and α-L-LNA) duplexes in comparison to pure DNA and RNA duplexes.

Locked nucleic acid (LNA) is a chemical modification which introduces a -O-CH2- linkage in the furanose sugar of nucleic acids and blocks its conformation in a particular state. Two types of modifications, namely, 2'-O,4'-C-methylene-ß-D-ribofuranose (ß-D-LNA) and 2'-O,4'-C-methylene-α-L-ribofuranose (α-L-LNA), have been shown to yield RNA and DNA duplex-like structures, respectively. LNA modifications lead to increased melting temperatures of DNA and RNA duplexes, and have been suggested as potential therapeutic agents in antisense therapy. In this study, molecular dynamics (MD) simulations were performed on fully modified LNA duplexes and pure DNA and RNA duplexes sharing a similar sequence to investigate their structure, stabilities, and solvation properties. Both LNA duplexes undergo unwinding of the helical structure compared to the pure DNA and RNA duplexes. Though the α-LNA substituent has been proposed to mimic deoxyribose sugar in its conformational properties, the fully modified duplex was found to exhibit unique structural and dynamic properties with respect to the other three nucleic acid structures. Free energy calculations accurately capture the enhanced stabilization of the LNA duplex structures compared to DNA and RNA molecules as observed in experiments. π-stacking interaction between bases from complementary strands is shown to be one of the contributors to enhanced stabilization upon LNA substitution. A combination of two factors, namely, nature of the -O-CH2- linkage in the LNAs vs their absence in the pure duplexes and similar conformations of the sugar rings in DNA and α-LNA vs the other two, is suggested to contribute to the stark differences among the four duplexes studied here in terms of their structural, dynamic, and energetic properties.

Crenarchaeal chromatin proteins Cren7 and Sul7 compact DNA by inducing rigid bends.

Archaeal chromatin proteins share molecular and functional similarities with both bacterial and eukaryotic chromatin proteins. These proteins play an important role in functionally organizing the genomic DNA into a compact nucleoid. Cren7 and Sul7 are two crenarchaeal nucleoid-associated proteins, which are structurally homologous, but not conserved at the sequence level. Co-crystal structures have shown that these two proteins induce a sharp bend on binding to DNA. In this study, we have investigated the architectural properties of these proteins using atomic force microscopy, molecular dynamics simulations and magnetic tweezers. We demonstrate that Cren7 and Sul7 both compact DNA molecules to a similar extent. Using a theoretical model, we quantify the number of individual proteins bound to the DNA as a function of protein concentration and show that forces up to 3.5 pN do not affect this binding. Moreover, we investigate the flexibility of the bending angle induced by Cren7 and Sul7 and show that the protein-DNA complexes differ in flexibility from analogous bacterial and eukaryotic DNA-bending proteins.