Folding of a DNA hairpin loop structure in explicit solvent using replica-exchange molecular dynamics simulations

Localization of Potential Energy in Hydrogen Bonds of the ATXN2 Gene

It is known that a number of neurodegenerative diseases, also called diseases of waiting, are associated with the expansion of the polyQ tract in the first exon of the ATXN2 gene. In the expanded polyQ tract, the probability of occurrence of non-canonical configurations (hairpins, G-quadruplexes, etc.) is significantly higher than in the normal one. Obviously, for their formation, the occurrence of open states (OSs) is necessary. Calculations were made for these processes using the angular mechanical model of DNA. It has been established that the probability of the large OS zones genesis in a DNA segment depends not only on the "strength" of the nucleotide sequence but also on the factors determining the dynamics of DNA; localization of the energy in the DNA molecule and the potential energy of interaction between pairs of nitrogenous bases also depend on environmental parameters. The potential energy of hydrogen bonds does not remain constant, and oscillatory movements lead to its redistribution and localization. In this case, OSs effectively dissipate the energy of oscillations. Thus, mathematical modeling makes it possible to calculate the localization of mechanical energy, which is necessary for the OSs formation, and to predict the places of their origin, taking into account the mechanical oscillations of the DNA molecule.

Free Energy-Based Refinement of DNA Force Field via Modification of Multiple Nonbonding Energy Terms

The amber-OL21 force field (ff) was developed to better describe noncanonical DNA, including Z-DNA. Despite its improvements for DNA simulations, this study found that OL21's scope of application was limited by embedded ff artifacts. In a benchmark set of seven DNA molecules, including two double-stranded DNAs transitioning between B- and Z-DNA and five single-stranded DNAs folding into mini-dumbbell or G-quadruplex structures, the free energy landscapes obtained using OL21 revealed several issues: Z-DNA was overly stabilized; misfolded states in mini-dumbbell DNAs were most stable; DNA GQ folding was consistently biased toward an antiparallel topology. To address these issues, a simple van der Waals (vdW) correction scheme, referred to as vdW5, was proposed for OL21. This involved revising multiple nonbonding energy terms to improve the overall quality of the free energy landscapes. The vdW5 correction effectively eliminated the artifacts in OL21, providing significantly improved free energy representations for DNAs tested. The vdW5-level revision substantially enhanced the predictive power of DNA simulations, thereby extending the scope of application for amber-OL21.

Fluctuations in Medium Viscosity May Affect the Stability of the CAG Tract in the ATXN2 Gene

Background: Trinucleotide repeats are the cause of many neurodegenerative diseases that are currently incurable. In this regard, the question of the causes of occurrence and methods of prevention or treatment of diseases caused by the expansion of repeats in the CAG tract of the ATXN2 gene remains relevant. Previously, it was shown that the frequency of occurrence of additional OS (open states) zones increases with increasing length of the CAG tract, and the value inverse to the frequency correlates with the age of disease onset. Methods: In this work, the influence of the viscosity of the medium and the external torque on the stability of the CAG tract in the ATXN2 gene was studied using mathematical modeling methods. Results: It has been established that the probability of the appearance of additional OS zones of significant size increases with an increase in the CAG of the tract (k > 40 CAG repeats) for all viscosity values, however, at k ≤ 40, the change in viscosity does not significantly affect the probability of additional OS zones in the tract. Conclusions: It was found that under normal conditions (absence of pathology), viscosity does not have a reliable effect on the stability of the DNA molecule, but when pathology appears, an increase in viscosity contributes to an increase in DNA stability, and, accordingly, a decrease has a negative effect on the stabilization of the DNA molecule. In the zone of close to incomplete penetrance of the disease, viscosity does not have a reliable effect on the stability of the CAG tract.

Aptamer-based targeted delivery systems for cancer treatment using DNA origami and DNA nanostructures

Due to the limitations of conventional cancer treatment methods, nanomedicine has appeared as a promising alternative, allowing improved drug targeting and decreased drug toxicity. In the development of cancer nanomedicines, among various nanoparticles (NPs), DNA nanostructures are more attractive because of their precisely controllable size, shape, excellent biocompatibility, programmability, biodegradability, and facile functionalization. Aptamers are introduced as single-stranded RNA or DNA molecules with recognize their corresponding targets. So, incorporating aptamers into DNA nanostructures led to influential vehicles for bioimaging and biosensing as well as targeted cancer therapy. In this review, the recent developments in the application of aptamer-based DNA origami and DNA nanostructures in advanced cancer treatment have been highlighted. Some of the main methods of cancer treatment are classified as chemo-, gene-, photodynamic- and combined therapy. Finally, the opportunities and problems for targeted DNA aptamer-based nanocarriers for medicinal applications have also been discussed.

Controlled destabilization of caged circularized DNA oligonucleotides predicted by replica exchange molecular dynamics simulations

Spatiotemporal control is a critical issue in the design of strategies for the photoregulation of oligonucleotide activity. Efficient uncaging, i.e., activation by removal of photolabile protecting groups (PPGs), often necessitates multiple PPGs. An alternative approach is based on circularization strategies, exemplified by intrasequential circularization, also denoted photo-tethering, as introduced in [Seyfried et al., Angew. Chem., Int. Ed., 2017, 56, 359]. Here, we develop a computational protocol, relying on replica exchange molecular dynamics (REMD), in order to characterize the destabilization of a series of circularized, caged DNA oligonucleotides addressed in the aforementioned study. For these medium-sized (32 nt) oligonucleotides, melting temperatures are computed, whose trend is in good agreement with experiment, exhibiting a large destabilization and, hence, reduction of the melting temperature of the order of ΔTm ∼ 30 K as compared with the native species. The analysis of free energy landscapes confirms the destabilization pattern experienced by the circularized oligonucleotides. The present study underscores that computational protocols that capture controlled destabilization and uncaging of oligonucleotides are promising as predictive tools in the tailored photocontrol of nucleic acids.

Force-induced melting of DNA hairpin: Unfolding pathways and phase diagrams

Using the exact enumeration technique, we have studied the force-induced melting of a DNA hairpin on the face centered cubic lattice for two different sequences which differ in terms of loop closing base pairs. The melting profiles obtained from the exact enumeration technique is consistent with the Gaussian network model and Langevin dynamics simulations. Probability distribution analysis based on the exact density of states revealed the microscopic details of the opening of the hairpin. We showed the existence of intermediate states near the melting temperature. We further showed that different ensembles used to model single-molecule force spectroscopy setups may give different force-temperature diagrams. We delineate the possible reasons for the observed discrepancies.

Toward Force Fields with Improved Base Stacking Descriptions

Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.

Large fragment Sanger sequencing identifies the newly encountered variant that caused null alleles in parentage testing

Short tandem repeat (STR) is regarded as a crucial tool for personal identification as well as parentage testing. Thus, genotyping errors of STRs could have negative effects on the reliability of forensic identification. A null allele at the combined DNA index system (CODIS) core loci D2S1338 was found in a father-daughter pair with the AGCU Expressmarker 22 kit which was a commonly used commercial kit during our daily laboratory work. This null allele caused the father and daughter to not conform to the laws of inheritance, thus potentially generating erroneous conclusions that excluded parentage. To figure out the reason for this phenomenon, re-amplification with new primers and then large fragment Sanger sequencing was conducted. We found a G to G/T variation at the position which is fifty-nine bases away from the 3' end of the core repeat in both samples. This probably could be considered a novel variant at the primer binding region which had not been reported that resulted in the emergence of the null allele. We also found that there was more than one single-nucleotide polymorphism (SNP) with minor allele frequency (MAF) greater than 0.1 in the upstream and downstream sequences of D2S1338. When designing primers for amplification of D2S1338, the possible adverse results of these SNPs should be taken into account and avoided.

Structural switching/polymorphism by sequential base substitution at quasi-palindromic SNP site (G → A) in LCR of human β-globin gene cluster

The human β-globin gene Locus Control Region (LCR), a dominant regulator of globin gene expression contains five tissue-specific DNase I-hypersensitive sites (HSs). A single nucleotide polymorphism (SNP) (A → G) present in HS4 region of locus control region (LCR), have shown a notable association between the G allele and the occurrence of β-thalassemia. This SNP site exhibiting a hairpin - duplex equilibrium manifested in A → B like DNA transition has previously been reported from this laboratory. Since, DNA is a dynamic and adaptable molecule, so any change of a single base within a primary DNA sequence can produce major biological consequences commonly manifested in genetic disorders such as sickle cell anemia and β-thalassemia. Herein, the differential behavior of sequential single base substitutions G → A on the quasi-palindromic sequence (d-TGGGGGCCCCA; HPG11) has been explored. A combination of native gel electrophoresis, circular dichroism (CD), and UV-thermal denaturation (Tm) techniques have been used to investigate the structural polymorphism associated with various variants of HPG11 i.e. HPG11A2 to HPG11A5. The CD spectra confirmed that all the HPG11 variants exhibit a hairpin - duplex equilibrium. Oligomer concentration dependence on CD spectra has been correlated with A → B DNA conformational transition. However, as revealed in gel electrophoresis, HPG11A2 → A5 exhibit the formation of a tetramolecular structure (four-way junction) at higher oligomer concentration. UV-melting studies also supported the melting of hairpin, duplex and four-way junction structure. This polymorphism pattern may possibly be significant for DNA-protein recognition, in the process of regulation of LCR in the β-globin gene.

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry

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

Influence of Environmental Parameters on the Stability of the DNA Molecule

Fluctuations in viscosity within the cell nucleus have wide limits. When a DNA molecule passes from the region of high viscosity values to the region of low values, open states, denaturation bubbles, and unweaving of DNA strands can occur. Stabilization of the molecule is provided by energy dissipation-dissipation due to interaction with the environment. Separate sections of a DNA molecule in a twisted state can experience supercoiling stress, which, among other things, is due to complex entropic effects caused by interaction with a solvent. In this work, based on the numerical solution of a mechanical mathematical model for the interferon alpha 17 gene and a fragment of the Drosophila gene, an analysis of the external environment viscosity influence on the dynamics of the DNA molecule and its stability was carried out. It has been shown that an increase in viscosity leads to a rapid stabilization of the angular vibrations of nitrogenous bases, while a decrease in viscosity changes the dynamics of DNA: the rate of change in the angular deviations of nitrogenous bases increases and the angular deformations of the DNA strands increase at each moment of time. These processes lead to DNA instability, which increases with time. Thus, the paper considers the influence of the external environment viscosity on the dissipation of the DNA nitrogenous bases' vibrational motion energy. Additionally, the study on the basis of the described model of the molecular dynamics of physiological processes at different indicators of the rheological behavior of nucleoplasm will allow a deeper understanding of the processes of nonequilibrium physics of an active substance in a living cell to be obtained.

The structure and unzipping behavior of dumbbell and hairpin DNA revealed by real-time nanopore sensing

Hairpin structures play an essential role in DNA replication, transcription, and recombination. Single-molecule studies enable the real-time measurement and observation of the energetics and dynamics of hairpin structures, including folding and DNA-protein interactions. Nanopore sensing is emerging as a powerful tool for DNA sensing and sequencing, and previous research into hairpins using an α-hemolysin (α-HL) nanopore suggested that hairpin DNA enters from its stem side. In this work, the translocation and interaction of hairpin and dumbbell DNA samples with varying stems, loops, and toeholds were investigated systematically using a Mycobacterium smegmatis porin A (MspA) nanopore. It was found that these DNA constructs could translocate through the pore under a bias voltage above +80 mV, and blockage events with two conductance states could be observed. The events of the lower blockage were correlated with the loop size of the hairpin or dumbbell DNA (7 nt to 25 nt), which could be attributed to non-specific collisions with the pore, whereas the dwell time of events with the higher blockage were correlated with the stem length, thus indicating effective translocation. Furthermore, dumbbell DNA with and without a stem opening generated different dwell times when driven through the MspA nanopore. Finally, a new strategy based on the dwell time difference was developed to detect single nucleotide polymorphisms (SNPs). These results demonstrated that the unzipping behaviors and DNA-protein interactions of hairpin and dumbbell DNA could be revealed using nanopore technology, and this could be further developed to create sensors for the secondary structures of nucleic acids.

Multimodal imaging and genetic characteristics of Chinese patients with USH2A-associated nonsyndromic retinitis pigmentosa

BACKGROUND

To determine the clinical characteristics and molecular genetic background responsible for USH2A mutations associated with nonsyndromic retinitis pigmentosa (RP) in five Chinese families, a retrospective cross-sectional study was performed.

METHODS

Data on detailed history and comprehensive ophthalmological examinations were extracted from medical charts. Genomic DNA was sequenced by whole-exome sequencing. The pathogenicity predictions were evaluated by in silico analysis. The structural modeling of the wide-type and mutant USH2A proteins was displayed based on the I-Tasser software.

RESULTS

The ultra-wide-field fundus imaging showed a distinctive pattern of hyperautofluorescence in the parafoveal ring with macular sparing. Ten USH2A variants were detected, including seven missense mutations, two splicing mutations, and one insertion mutation. Six of these variants have already been reported, and the remaining four were novel. Of the de novo mutations, the p.C931Y and p.G4489S mutations were predicted to be deleterious or probably damaging; the p.M4853V mutation was predicted to be neutral or benign; and the IVS22+3A>G mutation was a splicing mutation that could influence mRNA splicing and affect the formation of the hairpin structure of the USH2A protein.

CONCLUSIONS

Our data further confirm that USH2A protein plays a pivotal role in the maintenance of photoreceptors and expand the spectrum of USH2A mutations that are associated with nonsyndromic RP in Chinese patients.

Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot

The assembly mechanism of RNA, vital to describing its functions, depends on both the sequence and the metal ion concentration. How the latter influences the folding trajectories remains an important unsolved problem. Here, we examine the folding pathways of an RNA pseudoknot (PK) with key functional roles in transcription and translation, using a combination of experiments and simulations. We demonstrate that the PK, consisting of two hairpins with differing stabilities, folds by parallel pathways. Surprisingly, the flux between them is modulated by monovalent salt concentration. Our work shows that the order of assembly of PKs is determined by the relative stability of the hairpins, implying that the folding landscape can be controlled by sequence and ion concentration.

The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.

Structural Switch from Hairpin to Duplex/Antiparallel G-Quadruplex at Single-Nucleotide Polymorphism (SNP) Site of Human Apolipoprotein E (APOE) Gene Coding Region

A gradual dementia, which leads to the loss of memory and intellectual abilities, is the main characteristics of Alzheimer's disease. Amyloid-β (Aβ) plaques are the main components that accumulate and form clumps in the brains of people suffering from Alzheimer's disease. Apolipoprotein E (APOE), an amyloid-binding protein is considered as one of the main genetic risk factor of the late-onset Alzheimer's disease. Different isoforms of APOE gene named APOE2, APOE3, and APOE4 are known to exist, which differ from each other at certain positions involving single-nucleotide polymorphisms (SNPs). Out of these isoforms, APOE4 increases the risk of developing late-onset Alzheimer's disease, whereas APOE3 is the most common among the general population. APOE4 differs from the common APOE3 by only one nucleotide at position +2985 (T to C), which results in immense alteration in the structure and function of the APOE gene. A combination of gel electrophoresis (polyacrylamide gel electrophoresis, PAGE), circular dichroism (CD), CD melting, thermal difference spectra and UV-thermal denaturation (T_M) techniques was used to investigate the structural polymorphism associated with T → C single-nucleotide polymorphism (SNP) at the GC-rich sequence (d-TGGAGGACGTGTGCGGCCGCCT; APOE22T). Herein, we report that APOE22T DNA sequence switches between hairpin to antiparallel quadruplex from low to high oligomer concentration. On the contrary, its C-counterpart (APOE22C) forms hairpin as well as intermolecular antiparallel duplex structure. This structural change may possibly contribute to the protein recognition pattern, which, in turn, might control the APOE gene expression.

Origins of biological function in DNA and RNA hairpin loop motifs from replica exchange molecular dynamics simulation

Temperature REMD reveals how local chemical changes can result in markedly differing conformational landscapes for DNA and RNA hairpin loops.

Insilico direct folding of thrombin-binding aptamer G-quadruplex at all-atom level

The reversible folding of the thrombin-binding DNA aptamer G-quadruplexes (GQs) (TBA-15) starting from fully unfolded states was demonstrated using a prolonged time scale (10-12 μs) parallel tempering metadynamics (PTMetaD) simulation method in conjunction with a modified version of the AMBER bsc1 force field. For unbiased descriptions of the folding free energy landscape of TBA-15, this force field was minimally modified. From this direct folding simulation using the modified bsc1 force field, reasonably converged free energy landscapes were obtained in K+-rich aqueous solution (150 mM), providing detailed atomistic pictures of GQ folding mechanisms for TBA-15. This study found that the TBA folding occurred via multiple folding pathways with two major free energy barriers of 13 and 15 kcal/mol in the presence of several intermediate states of G-triplex variants. The early formation of these intermediates was associated with a single K+ ion capturing. Interestingly, these intermediate states appear to undergo facile transitions among themselves through relatively small energy barriers.

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.

Ab initio RNA folding

RNA molecules are essential cellular machines performing a wide variety of functions for which a specific three-dimensional structure is required. Over the last several years, the experimental determination of RNA structures through x-ray crystallography and NMR seems to have reached a plateau in the number of structures resolved each year, but as more and more RNA sequences are being discovered, the need for structure prediction tools to complement experimental data is strong. Theoretical approaches to RNA folding have been developed since the late nineties, when the first algorithms for secondary structure prediction appeared. Over the last 10 years a number of prediction methods for 3D structures have been developed, first based on bioinformatics and data-mining, and more recently based on a coarse-grained physical representation of the systems. In this review we are going to present the challenges of RNA structure prediction and the main ideas behind bioinformatic approaches and physics-based approaches. We will focus on the description of the more recent physics-based phenomenological models and on how they are built to include the specificity of the interactions of RNA bases, whose role is critical in folding. Through examples from different models, we will point out the strengths of physics-based approaches, which are able not only to predict equilibrium structures, but also to investigate dynamical and thermodynamical behavior, and the open challenges to include more key interactions ruling RNA folding.

Molecular Dynamics Simulations of Nucleic Acids. From Tetranucleotides to the Ribosome

We present a brief overview of explicit solvent molecular dynamics (MD) simulations of nucleic acids. We explain physical chemistry limitations of the simulations, namely, the molecular mechanics (MM) force field (FF) approximation and limited time scale. Further, we discuss relations and differences between simulations and experiments, compare standard and enhanced sampling simulations, discuss the role of starting structures, comment on different versions of nucleic acid FFs, and relate MM computations with contemporary quantum chemistry. Despite its limitations, we show that MD is a powerful technique for studying the structural dynamics of nucleic acids with a fast growing potential that substantially complements experimental results and aids their interpretation.

Reliable oligonucleotide conformational ensemble generation in explicit solvent for force field assessment using reservoir replica exchange molecular dynamics simulations

Molecular dynamics force field development and assessment requires a reliable means for obtaining a well-converged conformational ensemble of a molecule in both a time-efficient and cost-effective manner. This remains a challenge for RNA because its rugged energy landscape results in slow conformational sampling and accurate results typically require explicit solvent which increases computational cost. To address this, we performed both traditional and modified replica exchange molecular dynamics simulations on a test system (alanine dipeptide) and an RNA tetramer known to populate A-form-like conformations in solution (single-stranded rGACC). A key focus is on providing the means to demonstrate that convergence is obtained, for example, by investigating replica RMSD profiles and/or detailed ensemble analysis through clustering. We found that traditional replica exchange simulations still require prohibitive time and resource expenditures, even when using GPU accelerated hardware, and our results are not well converged even at 2 μs of simulation time per replica. In contrast, a modified version of replica exchange, reservoir replica exchange in explicit solvent, showed much better convergence and proved to be both a cost-effective and reliable alternative to the traditional approach. We expect this method will be attractive for future research that requires quantitative conformational analysis from explicitly solvated simulations.

Frontiers in molecular dynamics simulations of DNA

It has been known for decades that DNA is extremely flexible and polymorphic, but our knowledge of its accessible conformational space remains limited. Structural data, primarily from X-ray diffraction studies, is sparse in comparison to the manifold configurations possible, and direct experimental examinations of DNA's flexibility still suffer from many limitations. In the face of these shortcomings, molecular dynamics (MD) is now an essential tool in the study of DNA. It affords detailed structural and dynamical insights, which explains its recent transition from a small number of highly specialized laboratories to a large variety of groups dealing with challenging biological problems. MD is now making an irreversible journey to the mainstream of research in biology, with the attendant opportunities and challenges. But given the speed with which MD studies of DNA have spread, the roots remain somewhat shallow: in many cases, there is a lack of deep knowledge about the foundations, strengths, and limits of the technique. In this Account, we discuss how MD has become the most important source of structural and flexibility data on DNA, focusing on advances since 2007 of atomistic MD in the description of DNA under near-physiological conditions and highlighting the possibilities and shortcomings of the technique. The evolution in the field over the past four years is a prelude to the ongoing revolution. The technique has gained in robustness and predictive power, which when coupled with the spectacular improvements in software and hardware has enabled the tackling of systems of increasing complexity. Simulation times of microseconds have now been achieved, with even longer times when specialized hardware is used. As a result, we have seen the first real-time simulation of large conformational transitions, including folding and unfolding of short DNA duplexes. Noteworthy advances have also been made in the study of DNA-ligand interactions, and we predict that a global thermodynamic and kinetic picture of the binding landscape of DNA will become available in a few years. MD will become a crucial tool in areas such as biomolecular engineering and synthetic biology. MD has also been shown to be an excellent source of parameters for mesoscopic models of DNA flexibility. Such models can be refined through atomistic MD simulations on small duplexes and then applied to the study of entire chromosomes. Recent evidence suggests that MD-derived elastic models can successfully predict the position of regulatory regions in DNA and can help advance our understanding of nucleosome positioning and chromatin plasticity. If these results are confirmed, MD simulations can become the ultimate tool to decipher a physical code that can contribute to gene regulation. We are entering the golden age of MD simulations of DNA. Undoubtedly, the expectations are high, but the challenges are also enormous. These include the need for more accurate potential energy functionals and for longer and more complex simulations in more realistic systems. The joint research effort of several groups will be crucial for adapting the technique to the requirements of the coming decade.

A kinetic zipper model with intrachain interactions applied to nucleic acid hairpin folding kinetics

Single-stranded DNA and RNA hairpin structures with 4-10 nucleotides (nt) in the loop and 5-8 basepairs (bp) in the stem fold on 10-100 μs timescale. In contrast, theoretical estimate of first contact time of two ends of an ideal semiflexible polymer of similar lengths (with persistence length ~2-nt) is 10-100 ns. We propose that this three-orders-of-magnitude difference between these two timescales is a result of roughness in the folding free energy surface arising from intrachain interactions. We present a statistical mechanical model that explicitly includes all misfolded microstates with nonnative Watson-Crick (WC) and non-WC contacts. Rates of interconversion between different microstates are described in terms of two adjustable parameters: the strength of the non-WC interactions (ΔG(nWC)) and the rate at which a basepair is formed adjacent to an existing basepair (k(bp)(+)). The model accurately reproduces the temperature and loop-length dependence of the measured relaxation rates in temperature-jump studies of a 7-bp stem, single-stranded DNA hairpin with 4-20-nt-long poly(dT) loops, with ΔG(nWC) ≈ -2.4 kcal/mol and k(bp)(+) ≥ (1 ns)(-1), in 100 mM NaCl. Thus, our model provides a microscopic interpretation of the slow hairpin folding times as well as an estimate of the strength of intrachain interactions.

Role of the closing base pair for d(GCA) hairpin stability: free energy analysis and folding simulations

Hairpin loops belong to the most important structural motifs in folded nucleic acids. The d(GNA) sequence in DNA can form very stable trinucleotide hairpin loops depending, however, strongly on the closing base pair. Replica-exchange molecular dynamics (REMD) were employed to study hairpin folding of two DNA sequences, d(gcGCAgc) and d(cgGCAcg), with the same central loop motif but different closing base pairs starting from single-stranded structures. In both cases, conformations of the most populated conformational cluster at the lowest temperature showed close agreement with available experimental structures. For the loop sequence with the less stable G:C closing base pair, an alternative loop topology accumulated as second most populated conformational state indicating a possible loop structural heterogeneity. Comparative-free energy simulations on induced loop unfolding indicated higher stability of the loop with a C:G closing base pair by ~3 kcal mol(-1) (compared to a G:C closing base pair) in very good agreement with experiment. The comparative energetic analysis of sampled unfolded, intermediate and folded conformational states identified electrostatic and packing interactions as the main contributions to the closing base pair dependence of the d(GCA) loop stability.

Conformational changes of non-B DNA

In contrast to B-DNA that has a right-handed double helical structure with Watson-Crick base pairing under the ordinary physiological conditions, repetitive DNA sequences under certain conditions have the potential to fold into non-B DNA structures such as hairpin, triplex, cruciform, left-handed Z-form, tetraplex, A-motif, etc. Since the non-B DNA-forming sequences induce the genetic instability and consequently can cause human diseases, the molecular mechanism for their genetic instability has been extensively investigated. On the contrary, non-B DNA can be widely used for application in biotechnology because many DNA breakage hotspots are mapped in or near the sequences that have the potential to adopt non-B DNA structures. In addition, they are regarded as a fascinating material for the nanotechnology using non-B DNAs because they do not produce any toxic byproducts and are robust enough for the repetitive working cycle. This being the case, an understanding on the mechanism and dynamics of their structural changes is important. In this critical review, we describe the latest studies on the conformational dynamics of non-B DNAs, with a focus on G-quadruplex, i-motif, Z-DNA, A-motif, hairpin and triplex (189 references).

Lattice model of oligonucleotide hybridization in solution. I. Model and thermodynamics

A coarse-grained lattice model of DNA oligonucleotides is proposed to investigate the general mechanisms by which single-stranded oligonucleotides hybridize to their complementary strands in solution. The model, based on a high-coordination cubic lattice, is simple enough to allow the direct simulation of DNA solutions, yet capturing how the fundamental thermodynamic processes are microscopically encoded in the nucleobase sequences. Physically relevant interactions are considered explicitly, such as interchain excluded volume, anisotropic base-pairing and base-stacking, and single-stranded bending rigidity. The model is studied in detail by a specially adapted Monte Carlo simulation method, based on parallel tempering and biased trials, which is designed to overcome the entropic and enthalpic barriers associated with the sampling of hybridization events of multiple single-stranded chains in solution. This methodology addresses both the configurational complexity of bringing together two complementary strands in a favorable orientation (entropic barrier) and the energetic penalty of breaking apart multiple associated bases in a double-stranded state (enthalpic barrier). For strands with sequences restricted to nonstaggering association and homogeneous pairing and stacking energies, base-pairing is found to dominate the hybridization over the translational and conformational entropy. For strands with sequence-dependent pairing corresponding to that of DNA, the complex dependence of the model's thermal stability on concentration, sequence, and degree of complementarity is shown to be qualitatively and quantitatively consistent both with experiment and with the predictions of statistical mechanical models.

Structural, mechanical, and thermodynamic properties of a coarse-grained DNA model

We explore in detail the structural, mechanical, and thermodynamic properties of a coarse-grained model of DNA similar to that recently introduced in a study of DNA nanotweezers [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, Phys. Rev. Lett. 134, 178101 (2010)]. Effective interactions are used to represent chain connectivity, excluded volume, base stacking, and hydrogen bonding, naturally reproducing a range of DNA behavior. The model incorporates the specificity of Watson-Crick base pairing, but otherwise neglects sequence dependence of interaction strengths, resulting in an "average base" description of DNA. We quantify the relation to experiment of the thermodynamics of single-stranded stacking, duplex hybridization, and hairpin formation, as well as structural properties such as the persistence length of single strands and duplexes, and the elastic torsional and stretching moduli of double helices. We also explore the model's representation of more complex motifs involving dangling ends, bulged bases and internal loops, and the effect of stacking and fraying on the thermodynamics of the duplex formation transition.

Binding of the bacteriophage P22 N-peptide to the boxB RNA motif studied by molecular dynamics simulations

Protein-RNA interactions are important for many cellular processes. The Nut-utilization site (N)-protein of bacteriophages contains an N-terminal arginine-rich motif that undergoes a folding transition upon binding to the boxB RNA hairpin loop target structure. Molecular dynamics simulations were used to investigate the dynamics of the P22 N-peptide-boxB complex and to elucidate the energetic contributions to binding. In addition, the free-energy changes of RNA and peptide conformational adaptation to the bound forms, as well as the role of strongly bound water molecules at the peptide-RNA interface, were studied. The influence of peptide amino acid substitutions and the salt dependence of interaction were investigated and showed good agreement with available experimental results. Several tightly bound water molecules were found at the RNA-binding interface in both the presence and absence of N-peptide. Explicit consideration of the waters resulted in shifts of calculated contributions during the energetic analysis, but overall similar binding energy contributions were found. Of interest, it was found that the electrostatic field of the RNA has a favorable influence on the coil-to-alpha-helix transition of the N-peptide already outside of the peptide-binding site. This result may have important implications for understanding peptide-RNA complex formation, which often involves coupled folding and association processes. It indicates that electrostatic interactions near RNA molecules can lead to a shift in the equilibrium toward the bound form of an interacting partner before it enters the binding pocket.

Elbow flexibility of the kt38 RNA kink-turn motif investigated by free-energy molecular dynamics simulations

Kink-turns (K-turns) are common structural motifs that can introduce sharp kinks into double-stranded RNA, and have been proposed to mediate large-scale motions in the ribosome. K-turns consist of a bulge loop region flanked by trans sugar-Hoogsteen G:A pairs, and the sharp kink conformation is stabilized by A-minor interactions (adenine contacting a G:C basepair in the minor groove). Umbrella-sampling molecular dynamics simulations were used to disrupt an A-minor interaction in the ribosomal kt38 turn and to calculate the associated free-energy change. Coupling of umbrella sampling with replica exchanges between neighboring umbrella-sampling intervals could further improve the convergence of the free-energy calculations. The simulations revealed a coupled A-minor disruption and global opening of the K-turn motif, and allowed us to characterize several intermediate A-minor conformations. The calculated free-energy profile indicated a meta-stable, semi-open structure of slightly higher free energy ( approximately 1 kcal mol(-1)), separated by a small free-energy barrier ( approximately 1.5 kcal mol(-1)) from the closed (highly kinked) form. Both K-turn states are stabilized by distinct variants of the A-minor interaction. Further opening of the K-turn structure required significantly larger free-energy changes. The semi-open form had a reduced kink angle compatible with experimental data on K-turn solution structures, and opening was coupled to a continuous global unwinding of the K-turn motif. The range of free-energy changes associated with kt38 opening and unwinding are compatible with the idea that K-turns may facilitate biologically relevant motions during large-scale ribosome dynamics.

Arc-repressor dimerization on DNA: folding rate enhancement by colocalization

Multimeric proteins are ubiquitous in many cellular processes that require high levels of regulation. Eukaryotic gene expression is often regulated by a mechanism of combinatorial control that involves the binding of dimeric transcription factors to DNA together with the coordinated activity of additional proteins. In this study, we investigated the dimerization of the Arc-repressor on DNA with the aim of achieving microscopic insight into the possible advantages of interacting with DNA as a complex rather than as a monomeric single-domain protein. We used a computational coarse-grained model in which the protein dynamics was governed by native interactions and protein-DNA interactions were dictated by electrostatic forces. Inspired by previous experimental work that showed an enhanced refolding rate for the Arc-repressor in the presence of DNA and other polyanions, we focused on the mechanism and kinetics of the assembly of Arc monomers in the presence of single- (ssDNA) and double-stranded DNA (dsDNA) molecules in a low-salt concentration environment. The electrostatic interactions that attract the protein to the dsDNA were shown to be fundamental in colocalizing the unfolded Arc chains and in accelerating refolding. Arc monomers bind the dsDNA efficiently and nonspecifically, and search for each other via one-dimensional diffusion. The fastest folding of Arc is observed for DNA of 30 bp. Longer DNA is significantly less efficient in accelerating the Arc refolding rate, since the two subunits search distinct regions of the one-dimensional DNA and are therefore much less colocalized. The probability that the two unfolded chains will meet on 200 bp DNA is similar to that in the bulk. The colocalization of Arc subunits on ssDNA results in much faster folding compared to that obtained on dsDNA of the same length. Differences in the rate of Arc refolding, cooperativity, and the structure of its transition state ensemble introduced by ssDNA and dsDNA molecules demonstrate the important role of colocalization in biological self-assembly processes.

Enhanced conformational sampling of nucleic acids by a new Hamiltonian replica exchange molecular dynamics approach

Although molecular dynamics (MD) simulations have been applied frequently to study flexible molecules, the sampling of conformational states separated by barriers is limited due to currently possible simulation time scales. Replica-exchange (Rex)MD simulations that allow for exchanges between simulations performed at different temperatures (T-RexMD) can achieve improved conformational sampling. However, in the case of T-RexMD the computational demand grows rapidly with system size. A Hamiltonian RexMD method that specifically enhances coupled dihedral angle transitions has been developed. The method employs added biasing potentials as replica parameters that destabilize available dihedral substates and was applied to study coupled dihedral transitions in nucleic acid molecules. The biasing potentials can be either fixed at the beginning of the simulation or optimized during an equilibration phase. The method was extensively tested and compared to conventional MD simulations and T-RexMD simulations on an adenine dinucleotide system and on a DNA abasic site. The biasing potential RexMD method showed improved sampling of conformational substates compared to conventional MD simulations similar to T-RexMD simulations but at a fraction of the computational demand. It is well suited to study systematically the fine structure and dynamics of large nucleic acids under realistic conditions including explicit solvent and ions and can be easily extended to other types of molecules.

Brownian dynamics simulations of single-stranded DNA hairpins

We present a Brownian dynamics model which we use to study the kinetics and thermodynamics of single-stranded DNA hairpins, gaining insights into the role of stem mismatches and the kinetics rates underlying the melting transition. The model is a base-backbone type in which the DNA bases and sugar-phosphate backbone are represented as single units (beads) in the context of the Brownian dynamics simulations. We employ a minimal number of bead-bead interactions, leading to a simple computational scheme. To demonstrate the veracity of our model for DNA hairpins, we show that the model correctly captures the effects of base stacking, hydrogen bonding, and temperature on both the thermodynamics and the kinetics of hairpin formation and melting. When cast in dimensionless form, the thermodynamic results obtained from the present model compare favorably with default predictions of the m-fold server, although the present model is not sufficiently robust to provide dimensional results. The kinetic data at low temperatures indicate frequent but short-lived opening events, consistent with the measured chain end-to-end probability distribution. The model is also used to study the effect of base mismatches in the stem of the hairpin. With the parameters used here, the model overpredicts the relative shift in the melting temperature due to mismatches. The melting transition can be primarily attributed to a rapid increase in the hairpin opening rate rather than an equivalent decrease in the closing rate, in agreement with single-molecule experimental data.

Free energy landscape of a DNA-carbon nanotube hybrid using replica exchange molecular dynamics

The DNA-carbon nanotube hybrid (DNA-CN) consists of a single-wall carbon nanotube (SWCN) coated with a self-assembled monolayer of single-stranded DNA (ssDNA). Recent experiments have demonstrated that this nanomaterial is ideal for numerous nanotechnological applications. Despite this importance, the structure of this material remains poorly understood. Molecular dynamics (MD) simulations have provided information about the self-assembly mechanisms and ssDNA conformations that characterize DNA-CN. However, MD simulations of biopolymers at low temperatures (T approximately 300 K) result in kinetic trapping that limits conformational sampling. Here, we present results of a large-scale replica exchange molecular dynamics (REMD) simulation that provides extensive sampling of the entire ensemble of oligonucleotide conformations in a (GT)(7)-SWCN hybrid. We calculate the free energy landscape and find minima corresponding to six distinct conformations, with a nonhelical loop structure as the global minimum. The hybrid contains significant structural disorder, with desorbed bases as an important structural feature. These results expand our understanding of DNA-CN and indicate the relevance of REMD for explorations of the physical properties of organic-inorganic multifunctional nanomaterials.

Unfolding and melting of DNA (RNA) hairpins: the concept of structure-specific 2D dynamic landscapes

A 2D free-energy landscape model is presented to describe the (un)folding transition of DNA/RNA hairpins, together with molecular dynamics simulations and experimental findings. The dependence of the (un)folding transition on the stem sequence and the loop length is shown in the enthalpic and entropic contributions to the free energy. Intermediate structures are well defined by the two coordinates of the landscape during (un)zipping. Both the free-energy landscape model and the extensive molecular dynamics simulations totaling over 10 mus predict the existence of temperature-dependent kinetic intermediate states during hairpin (un)zipping and provide the theoretical description of recent ultrafast temperature-jump studies which indicate that hairpin (un)zipping is, in general, not a two-state process. The model allows for lucid prediction of the collapsed state(s) in simple 2D space and we term it the kinetic intermediate structure (KIS) model.

Structural insight into RNA hairpin folding intermediates

Hairpins are a ubiquitous secondary structure motif in RNA molecules. Despite their simple structure, there is some debate over whether they fold in a two-state or multi-state manner. We have studied the folding of a small tetraloop hairpin using a serial version of replica exchange molecular dynamics on a distributed computing environment. On the basis of these simulations, we have identified a number of intermediates that are consistent with experimental results. We also find that folding is not simply the reverse of high-temperature unfolding and suggest that this may be a general feature of biomolecular folding.