Molecular modeling of nucleic acid structure

Defect-Induced Double-Stranded DNA Unwinding on Graphene

Several works have shown that graphene materials can effectively regulate the double-stranded DNA (dsDNA) structures and are used to remove antibiotic resistance genes in the environment, during which the morphology of the graphene surface plays a key role. However, the mechanism of how different graphene surfaces interact with dsDNA is poorly documented. Here, the interactions of dsDNA with defective graphene (D-Gra) and pristine graphene (P-Gra) have been explored by molecular dynamics simulations. Our data clearly showed that both D-Gra and P-Gra were able to attract dsDNA to form stable bindings. However, the structure evolutions of dsDNA are distinctly different. In detail, D-Gra can initiate quick unwinding of dsDNA and cause significant structural disruption. While for P-Gra, it demonstrated a much weaker capability to disrupt the dsDNA structure. This difference is due to the strong electrostatic interaction between defects and DNA nucleotides. Nucleotides can be highly restricted by the defect while the other parts of dsDNA could move along the transverse directions of D-Gra. This effectively introduces a "pulling force" from the defect that causes the breaking of the hydrogen bonds between dsDNA base pairs. Such force finally leads to the serious unwinding of dsDNA. Our present findings could help us to better understand the molecular mechanism of how the dsDNA canonical B-form was lost upon adsorption to graphene. The findings of the key roles of defects on graphene are beneficial for the design of functional graphenic materials for biological and medical applications through nanostructure engineering.

Properties and biological impact of RNA G-quadruplexes: from order to turmoil and back

Guanine-quadruplexes (G4s) are non-canonical four-stranded structures that can be formed in guanine (G) rich nucleic acid sequences. A great number of G-rich sequences capable of forming G4 structures have been described based on in vitro analysis, and evidence supporting their formation in live cells continues to accumulate. While formation of DNA G4s (dG4s) within chromatin in vivo has been supported by different chemical, imaging and genomic approaches, formation of RNA G4s (rG4s) in vivo remains a matter of discussion. Recent data support the dynamic nature of G4 formation in the transcriptome. Such dynamic fluctuation of rG4 folding-unfolding underpins the biological significance of these structures in the regulation of RNA metabolism. Moreover, rG4-mediated functions may ultimately be connected to mechanisms underlying disease pathologies and, potentially, provide novel options for therapeutics. In this framework, we will review the landscape of rG4s within the transcriptome, focus on their potential impact on biological processes, and consider an emerging connection of these functions in human health and disease.

Multiscale Conformational Sampling Reveals Excited-State Locality in DNA Self-Repair Mechanism

Ultraviolet (UV) irradiation is known to be responsible for DNA damage. However, experimental studies in DNA oligonucleotides have shown that UV light can also induce sequence-specific self-repair. Following charge transfer from a guanine adenine sequence adjacent to a cyclobutane pyrimidine dimer (CPD), the covalent bond between the two thymines could be cleaved, recovering the intact base sequence. Mechanistic details promoting the self-repair remained unclear, however. In our theoretical study, we investigated whether optical excitation could directly lead to a charge-transfer state, thereby initiating the repair, or whether the initial excited state remains localized on a single nucleobase. We performed conformational sampling of 200 geometries of the damaged DNA double strand solvated in water and used a hybrid quantum and molecular mechanics approach to compute excited states at the complete active space perturbation level of theory. Analysis of the conformational data set clearly revealed that the excited-state properties are uniformly distributed across the fluctuations of the nucleotide in its natural environment. From the electronic wavefunction, we learned that the electronic transitions remained predominantly local on either adenine or guanine, and no direct charge transfer occurred in the experimentally accessed energy range. The investigated base sequence is not only specific to the CPD repair mechanism but ubiquitously occurs in nucleic acids. Our results therefore give a very general insight into the charge locality of UV-excited DNA, a property that is regarded to have determining relevance in the structural consequences following absorption of UV photons.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

NMR Structure Determination for Oligonucleotides

NMR spectroscopy is a versatile tool for determining the structure and dynamics of nucleic acids under solution conditions. In this unit, we provide an overview and detail of the experiments and methods used in our laboratory to determine the structure of oligonucleotides at natural abundance, thus limiting our approach to ¹ H, ¹³ C, and ³¹ P NMR techniques. Isotopic labeling is heavily used in RNA NMR studies, however, labeling of DNA is still less common and, if modified nucleotides are investigated, is exceptionally expensive or not feasible. Each method described here is extensively documented and annotated with tips and observations to facilitate their application. Sections are devoted to sample preparation, NMR experiments and setup, resonance assignment, structure generation protocols, evaluation, tips that may be useful, and software sources. © 2018 by John Wiley & Sons, Inc.

Allosterism and signal transfer in DNA

We analysed the basic mechanisms of signal transmission in DNA and the origins of the allostery exhibited by systems such as the ternary complex BAMHI-DNA-GRDBD. We found that perturbation information generated by a primary protein binding event travels as a wave to distant regions of DNA following a hopping mechanism. However, such a structural perturbation is transient and does not lead to permanent changes in the DNA geometry and interaction properties at the secondary binding site. The BAMHI-DNA-GRDBD allosteric mechanism does not occur through any traditional models: direct (protein-protein), indirect (reorganization of the secondary site) readout or solvent-release. On the contrary, it is generated by a subtle and less common entropy-mediated mechanism, which might have an important role to explain other DNA-mediated cooperative effects.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Molecular modeling of nucleic acid structure

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. The unit provides an overview of the computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, the generation of an initial model, and characterization of the overall three-dimensional structure.

Molecular modeling of nucleic acid structure

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. This unit provides an overview of computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, generation of an initial model, and characterization of the overall three-dimensional structure.