DNA/TNA mesoscopic modeling of melting temperatures suggests weaker hydrogen bonding of CG than in DNA/RNA

Bubble Relaxation Dynamics in Homopolymer DNA Sequences

Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work, we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond time scale. By studying energy autocorrelation functions of relatively large bubbles inserted into thermalised DNA molecules, we extract characteristic relaxation times from the equilibration process for both adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of different amplitudes and widths are investigated through extensive statistics and appropriate fittings of their relaxation. Characteristic relaxation times increase with bubble amplitude and width. We show that, within the model, relaxation times are two orders of magnitude longer in GC sequences than in AT sequences. Overall, our results confirm that large bubbles leave a lasting impact on the molecule's dynamics, for times between 0.5-500 ns depending on the homopolymer type and bubble shape, thus clearly affecting long-time evolutions of the molecule.

Nucleic Acids and Their Analogues for Biomedical Applications

Nucleic acids are emerging as powerful and functional biomaterials due to their molecular recognition ability, programmability, and ease of synthesis and chemical modification. Various types of nucleic acids have been used as gene regulation tools or therapeutic agents for the treatment of human diseases with genetic disorders. Nucleic acids can also be used to develop sensing platforms for detecting ions, small molecules, proteins, and cells. Their performance can be improved through integration with other organic or inorganic nanomaterials. To further enhance their biological properties, various chemically modified nucleic acid analogues can be generated by modifying their phosphodiester backbone, sugar moiety, nucleobase, or combined sites. Alternatively, using nucleic acids as building blocks for self-assembly of highly ordered nanostructures would enhance their biological stability and cellular uptake efficiency. In this review, we will focus on the development and biomedical applications of structural and functional natural nucleic acids, as well as the chemically modified nucleic acid analogues over the past ten years. The recent progress in the development of functional nanomaterials based on self-assembled DNA-based platforms for gene regulation, biosensing, drug delivery, and therapy will also be presented. We will then summarize with a discussion on the advanced development of nucleic acid research, highlight some of the challenges faced and propose suggestions for further improvement.

Complete Mesoscopic Parameterization of Single LNA Modifications in DNA Applied to Oncogene Probe Design

The use of mesoscopic models to describe the thermodynamic properties of locked nucleic acid (LNA)-modified nucleotides can provide useful insights into their properties, such as hydrogen-bonding and stacking interactions. In addition, the mesoscopic parameters can be used to optimize LNA insertion in probes, to achieve accurate melting temperature predictions, and to obtain duplex opening profiles at the base-pair level. Here, we applied this type of model to parameterize a large set of melting temperatures for LNA-modified sequences, from published sources, covering all possible nearest-neighbor configurations. We have found a very large increase in Morse potentials, which indicates very strong hydrogen bonding as the main cause of improved LNA thermodynamic stability. LNA-modified adenine-thymine (AT) was found to have similar hydrogen bonding to unmodified cytosine-guanine (CG) base pairs, while for LNA CG, we found exceptionally large hydrogen bonding. In contrast, stacking interactions, which were thought to be behind the stability of LNA, were similar to unmodified DNA in most cases. We applied the new LNA parameters to the design of BRAF, KRAS, and EGFR oncogene variants by testing all possible LNA modifications. Selected sequences were then synthesized and had their hybridization temperatures measured, achieving a prediction accuracy within 1 °C. We performed a detailed base-pair opening analysis to discuss specific aspects of these probe hybridizations that may be relevant for probe design.

Distributions of bubble lifetimes and bubble lengths in DNA

We investigate the distribution of bubble lifetimes and bubble lengths in DNA at physiological temperature, by performing extensive molecular dynamics simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended version (ePBD) having a sequence-dependent stacking interaction, emphasizing the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content on these distributions. For both models we find that base pair-dependent (GC vs AT) thresholds for considering complementary nucleotides to be separated are able to reproduce the observed dependence of the melting temperature on the GC content of the DNA sequence. Using these thresholds for base pair openings, we obtain bubble lifetime distributions for bubbles of lengths up to ten base pairs as the GC content of the sequences is varied, which are accurately fitted with stretched exponential functions. We find that for both models the average bubble lifetime decreases with increasing either the bubble length or the GC content. In addition, the obtained bubble length distributions are also fitted by appropriate stretched exponential functions and our results show that short bubbles have similar likelihoods for any GC content, but longer ones are substantially more likely to occur in AT-rich sequences. We also show that the ePBD model permits more, longer-lived, bubbles than the PBD system.