Melting of the solvent structure around a RNA duplex: a molecular dynamics simulation study

Roles of non-canonical structures of nucleic acids in cancer and neurodegenerative diseases

Cancer and neurodegenerative diseases are caused by genetic and environmental factors. Expression of tumour suppressor genes is suppressed by mutations or epigenetic silencing, whereas for neurodegenerative disease-related genes, nucleic acid-based effects may be presented through loss of protein function due to erroneous protein sequences or gain of toxic function from extended repeat transcripts or toxic peptide production. These diseases are triggered by damaged genes and proteins due to lifestyle and exposure to radiation. Recent studies have indicated that transient, non-canonical structural changes in nucleic acids in response to the environment can regulate the expression of disease-related genes. Non-canonical structures are involved in many cellular functions, such as regulation of gene expression through transcription and translation, epigenetic regulation of chromatin, and DNA recombination. Transcripts generated from repeat sequences of neurodegenerative disease-related genes form non-canonical structures that are involved in protein transport and toxic aggregate formation. Intracellular phase separation promotes transcription and protein assembly, which are controlled by the nucleic acid structure and can influence cancer and neurodegenerative disease progression. These findings may aid in elucidating the underlying disease mechanisms. Here, we review the influence of non-canonical nucleic acid structures in disease-related genes on disease onset and progression.

New Insights into the Functions of Nucleic Acids Controlled by Cellular Microenvironments

The right-handed double-helical B-form structure (B-form duplex) has been widely recognized as the canonical structure of nucleic acids since it was first proposed by James Watson and Francis Crick in 1953. This B-form duplex model has a monochronic and static structure and codes genetic information within a sequence. Interestingly, DNA and RNA can form various non-canonical structures, such as hairpin loops, left-handed helices, triplexes, tetraplexes of G-quadruplex and i-motif, and branched junctions, in addition to the canonical structure. The formation of non-canonical structures depends not only on sequence but also on the surrounding environment. Importantly, these non-canonical structures may exhibit a wide variety of biological roles by changing their structures and stabilities in response to the surrounding environments, which undergo vast changes at specific locations and at specific times in cells. Here, we review recent progress regarding the interesting behaviors and functions of nucleic acids controlled by molecularly crowded cellular conditions. New insights gained from recent studies suggest that nucleic acids not only code genetic information in sequences but also have unknown functions regarding their structures and stabilities through drastic structural changes in cellular environments.

Chemical biology of non-canonical structures of nucleic acids for therapeutic applications

DNA forms not only the canonical duplex structure but also non-canonical structures. Most potential sequences that induce the formation of non-canonical structures are present in disease-related genes. Interestingly, biological reactions are inhibited or dysregulated by non-canonical structure formation in disease-related genes. To control biological reactions, methods for inducing the formation of non-canonical structures have been developed using small molecules and oligonucleotides. In this feature article, we review biological reactions such as replication, transcription, and reverse transcription controlled by non-canonical DNA structures formed by disease-related genes. Furthermore, we discuss recent studies aimed at developing methods for regulating these biological reactions using drugs targeting the DNA structure.

Salt Dependence of A-Form RNA Duplexes: Structures and Implications

The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.

Molecular Dynamics Simulation of a RNA Aptasensor

Single-stranded RNA aptamers have emerged as novel biosensor tools. However, the immobilization procedure of the aptamer onto a surface generally induces a loss of affinity. To understand this molecular process, we conducted a complete simulation study for the Flavin mononucleotide aptamer for which experimental data are available. Several molecular dynamics simulations (MD) of the Flavin in complex with its RNA aptamer were conducted in solution, linked with six thymidines (T6) and, finally, immobilized on an hexanol-thiol-functionalized gold surface. First, we demonstrated that our MD computations were able to reproduce the experimental solution structure and to provide a meaningful estimation of the Flavin free energy of binding. We also demonstrated that the T6 linkage, by itself, does not generate a perturbation of the Flavin recognition process. From the simulation of the complete biosensor system, we observed that the aptamer stays oriented parallel to the surface at a distance around 36 Å avoiding, this way, interaction with the surface. We evidenced a structural reorganization of the Flavin aptamer binding mode related to the loss of affinity and induced by an anisotropic distribution of sodium cationic densities. This means that ionic diffusion is different between the surface and the aptamer than above this last one. We suggest that these findings might be extrapolated to other nucleic acids systems for the future design of biosensors with higher efficiency and selectivity.

Spectroscopic and Molecular Docking Study of the Interaction of DNA with a Morpholinium Ionic Liquid

The structural integrity of a nucleic acid under various conditions determines its utility in biocatalysis and biotechnology. Exploration of the ionic liquids (ILs) for extraction of DNA and other nucleic acid based applications requires an understanding of the nature of interaction between the IL and DNA. Considering these aspects, we have studied the interaction between calf-thymus DNA and a less toxic morpholinium IL, [Mor1,2][Br], employing fluorescence correlation spectroscopy (FCS), conventional steady state and time-resolved fluorescence, circular dichroism (CD) and molecular docking techniques. While the CD spectra indicate the stability of DNA and retention of its B-form in the presence of the morpholinium IL, the docking study reveals that [Mor1,2](+) binds to the minor groove of DNA with a binding energy of -4.57 kcal mol(-1). The groove binding of the cationic component of the IL is corroborated by the steady state fluorescence data, which indicated displacement of a known minor groove binder, DAPI, from its DNA-bound state on addition of [Mor1,2][Br]. The FCS measurements show that the hydrodynamic radius of DNA remains more or less constant in the presence of [Mor1,2][Br], thus suggesting that the structure of DNA is retained in the presence of the IL. DNA melting experiments show that the thermal stability of DNA is enhanced in the presence of morpholinium IL.

Structure, stability and behaviour of nucleic acids in ionic liquids

Nucleic acids have become a powerful tool in nanotechnology because of their conformational polymorphism. However, lack of a medium in which nucleic acid structures exhibit long-term stability has been a bottleneck. Ionic liquids (ILs) are potential solvents in the nanotechnology field. Hydrated ILs, such as choline dihydrogen phosphate (choline dhp) and deep eutectic solvent (DES) prepared from choline chloride and urea, are 'green' solvents that ensure long-term stability of biomolecules. An understanding of the behaviour of nucleic acids in hydrated ILs is necessary for developing DNA materials. We here review current knowledge about the structures and stabilities of nucleic acids in choline dhp and DES. Interestingly, in choline dhp, A-T base pairs are more stable than G-C base pairs, the reverse of the situation in buffered NaCl solution. Moreover, DNA triplex formation is markedly stabilized in hydrated ILs compared with aqueous solution. In choline dhp, the stability of Hoogsteen base pairs is comparable to that of Watson-Crick base pairs. Moreover, the parallel form of the G-quadruplex is stabilized in DES compared with aqueous solution. The behaviours of various DNA molecules in ILs detailed here should be useful for designing oligonucleotides for the development of nanomaterials and nanodevices.

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.

An insight into structure and stability of DNA in ionic liquids from molecular dynamics simulation and experimental studies

Molecular dynamics simulation and biophysical analysis were employed to reveal the characteristics and the influence of ionic liquids (ILs) on the structural properties of DNA. Both computational and experimental evidence indicate that DNA retains its native B-conformation in ILs. Simulation data show that the hydration shells around the DNA phosphate group were the main criteria for DNA stabilization in this ionic media. Stronger hydration shells reduce the binding ability of ILs' cations to the DNA phosphate group, thus destabilizing the DNA. The simulation results also indicated that the DNA structure maintains its duplex conformation when solvated by ILs at different temperatures up to 373.15 K. The result further suggests that the thermal stability of DNA at high temperatures is related to the solvent thermodynamics, especially entropy and enthalpy of water. All the molecular simulation results were consistent with the experimental findings. The understanding of the properties of IL-DNA could be used as a basis for future development of specific ILs for nucleic acid technology.

Control of stability and structure of nucleic acids using cosolutes

The stabilities, structures, and functions of nucleic acids are responsive to surrounding conditions. Living cells contain biomolecules, including nucleic acids, proteins, polysaccharides, and other soluble and insoluble low-molecular weight components, that occupy a significant fraction of the cellular volume (up to 40%), resulting in a highly crowded intracellular environment. We have proven that conditions that mimic features of this intra-cellular environment alter the physical properties affect the stability, structure, and function of nucleic acids. The ability to control structure of nucleic acids by mimicking intra-cellular conditions will be useful in nanotechnology applications of nucleic acids. This paper describes methods that can be used to analyze quantitatively the intra-cellular environment effects caused by cosolutes on nucleic acid structures and to regulate properties of nucleic acids using cosolutes.

Thermosensing to adjust bacterial virulence in a fluctuating environment

The lifecycle of most microbial pathogens can be divided into two states: existence outside and inside their hosts. The sudden temperature upshift experienced upon entry from environmental or vector reservoirs into a warm-blooded host is one of the most crucial signals informing the pathogens to adjust virulence gene expression and their host-stress survival program. This article reviews the plethora of sophisticated strategies that bacteria have evolved to sense temperature, and outlines the molecular signal transduction mechanisms used to modulate synthesis of crucial virulence determinants. The molecular details of thermal control through conformational changes of DNA, RNA and proteins are summarized, complex and diverse thermosensing principles are introduced and their potential as drug targets or synthetic tools are discussed.

Computational exploration of mobile ion distributions around RNA duplex

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 ns of simulation using 64 replicas). The simulations show distinct behavior of single and double charged cations. Binding of Mg(2+) ion includes tight binding to specific sites while Na(+) binds only diffusively. The tight binding of Mg(2+) is with a solvation shell while Na(+) can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16 A from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) has a significant impact on the pair correlation functions. Predicted binding sites of Mg(2+) at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.

Thermodynamics of RNA melting, one base pair at a time

The melting of base pairs is a ubiquitous feature of RNA structural transitions, which are widely used to sense and respond to cellular stimuli. A recent study employing solution nuclear magnetic resonance (NMR) imino proton exchange spectroscopy provides a rare base-pair-specific view of duplex melting in the Salmonella FourU RNA thermosensor, which regulates gene expression in response to changes in temperature at the translational level by undergoing a melting transition. The authors observe "microscopic" enthalpy-entropy compensation--often seen "macroscopically" across a series of related molecular species--across base pairs within the same RNA. This yields variations in base-pair stabilities that are an order of magnitude smaller than corresponding variations in enthalpy and entropy. A surprising yet convincing link is established between the slopes of enthalpy-entropy correlations and RNA melting points determined by circular dichroism (CD), which argues that unfolding occurs when base-pair stabilities are equalized. A single AG-to-CG mutation, which enhances the macroscopic hairpin thermostability and folding cooperativity and renders the RNA thermometer inactive in vivo, spreads its effect microscopically throughout all base pairs in the RNA, including ones far removed from the site of mutation. The authors suggest that an extended network of hydration underlies this long-range communication. This study suggests that the deconstruction of macroscopic RNA unfolding in terms of microscopic unfolding events will require careful consideration of water interactions.

Direct observation of the temperature-induced melting process of the Salmonella fourU RNA thermometer at base-pair resolution

In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5'-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine-Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14-C25 base pair. The mismatch base pair in the wt RNA thermometer (A8-G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.

The pathway of oligomeric DNA melting investigated by molecular dynamics simulations

Details of the reaction coordinate for DNA melting are fundamental to much of biology and biotechnology. Recently, it has been shown experimentally that there are at least three states involved. To clarify the reaction mechanism of the melting transition of DNA, we perform 100-ns molecular dynamics simulations of a homo-oligomeric, 12-basepair DNA duplex, d(A(12)).d(T(12)), with explicit salt water at 400 K. Analysis of the trajectory reveals the various biochemically important processes that occur on different timescales. Peeling (including fraying from the ends), searching for Watson-Crick complements, and dissociation are recognizable processes. However, we find that basepair searching for Watson-Crick complements along a strand is not mechanistically tied to or directly accessible from the dissociation steps of strand melting. A three-step melting mechanism is proposed where the untwisting of the duplex is determined to be the major component of the reaction coordinate at the barrier. Though the observations are limited to the characteristics of the system being studied, they provide important insight into the mechanism of melting of other more biologically relevant forms of DNA, which will certainly differ in details from those here.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

High performance computing in biology: multimillion atom simulations of nanoscale systems

Computational methods have been used in biology for sequence analysis (bioinformatics), all-atom simulation (molecular dynamics and quantum calculations), and more recently for modeling biological networks (systems biology). Of these three techniques, all-atom simulation is currently the most computationally demanding, in terms of compute load, communication speed, and memory load. Breakthroughs in electrostatic force calculation and dynamic load balancing have enabled molecular dynamics simulations of large biomolecular complexes. Here, we report simulation results for the ribosome, using approximately 2.64 million atoms, the largest all-atom biomolecular simulation published to date. Several other nano-scale systems with different numbers of atoms were studied to measure the performance of the NAMD molecular dynamics simulation program on the Los Alamos National Laboratory Q Machine. We demonstrate that multimillion atom systems represent a 'sweet spot' for the NAMD code on large supercomputers. NAMD displays an unprecedented 85% parallel scaling efficiency for the ribosome system on 1024 CPUs. We also review recent targeted molecular dynamics simulations of the ribosome that prove useful for studying conformational changes of this large biomolecular complex in atomic detail.

Simulation and modeling of nucleic acid structure, dynamics and interactions

In moving towards the simulation of larger nucleic acid assemblies over longer timescales that include more accurate representations of the environment, we are nearing the end of an era characterized by single nanosecond molecular dynamics simulation of nucleic acids. We are excited by the promise and predictability of the modeling methods, yet remain prudently cautious of sampling and force field limitations. Highlights include the accurate representation of subtle drug-DNA interactions, the detailed study of modified and unusual nucleic acid structures, insight into the influence of dynamics on the structure of DNA, and exploration of the interaction of solvent and ions with nucleic acids.