Base Pair Fraying in Molecular Dynamics Simulations of DNA and RNA

The effect of the loop on the thermodynamic and kinetic of single base pair in pseudoknot

RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

Insights on the comparative affinity of ribonucleic acids with plant-based beta carboline alkaloid, harmine: Spectroscopic, calorimetric and computational evaluation

Small molecules as ligands target multifunctional ribonucleic acids (RNA) for therapeutic engagement. This study explores how the anticancer DNA intercalator harmine interacts various motifs of RNAs, including the single-stranded A-form poly (rA), the clover leaf tRNAphe, and the double-stranded A-form poly (rC)-poly (rG). Harmine showed the affinity to the polynucleotides in the order, poly (rA) > tRNAphe > poly (rC)·poly (rG). While no induced circular dichroism change was detected with poly (rC)poly (rG), significant structural alterations of poly (rA) followed by tRNAphe and occurrence of concurrent initiation of optical activity in the attached achiral molecule of alkaloid was reported. At 25 °C, the affinity further showed exothermic and entropy-driven binding. The interaction also highlighted heat capacity (ΔC o p ) and Gibbs energy contribution from the hydrophobic transfer (ΔG hyd) of binding with harmine. Molecular docking calculations indicated that harmine exhibits higher affinity for poly (rA) compared to tRNAphe and poly (rC)·poly (rG). Subsequent molecular dynamics simulations were conducted to investigate the binding mode and stability of harmine with poly(A), tRNAphe, and poly (rC)·poly (rG). The results revealed that harmine adopts a partial intercalative binding with poly (rA) and tRNAphe, characterized by pronounced stacking forces and stronger binding free energy observed with poly (rA), while a comparatively weaker binding free energy was observed with tRNAphe. In contrast, the stacking forces with poly (rC)·poly (rG) were comparatively less pronounced and adopts a groove binding mode. It was also supported by ferrocyanide quenching analysis. All these findings univocally provide detailed insight into the binding specificity of harmine, to single stranded poly (rA) over other RNA motifs, probably suggesting a self-structure formation in poly (rA) with harmine and its potential as a lead compound for RNA based drug targeting.

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

Dynamics of terminal fraying-peeling and hydrogen bonds dictates the sequential vs. cooperative melting pathways of nanoscale DNA and PNA triplexes

Peptide nucleic acids (PNAs) are charge-neutral synthetic DNA/RNA analogues. In many aspects of biology and biotechnology, the details of DNA and PNA melting reaction coordinates are crucial, and their associative/dissociative details remain inadequately understood. In the current study, we have attempted to gain insights into comparative melting pathways and binding affinity of iso-sequences of an 18-mer PNA-DNA-PNA triplex and the analogous DNA-DNA-DNA triplex, and DNA-DNA and PNA-DNA duplexes. It is intriguing that while the DNA-DNA-DNA triplex melts in two sequential steps, the PNA-DNA-PNA triplex melts in a single step and the mechanistic aspects for this difference are still not clear. We report an all-atom molecular dynamics simulation of both complexes in the temperature range of 300 to 500 K with 20 K intervals. Based on the trajectory analysis, we provide evidence that the association and dissociation are dictated by the differences in fraying-peeling effects from either terminus to the center in a zipper pattern among the PNA-DNA-PNA triplex and DNA-DNA-DNA triplexes. These are shown to be governed by the different characteristics of H-bonding, RMSD, and Free Energy Landscape (FEL) as analyzed by PCA, leading to the DNA-DNA-DNA triplex exhibiting sequential melting, while the PNA-DNA-PNA triplex shows cooperative melting of the whole fragment in a single-step. The PNA-DNA-PNA triplex base pairs are thermodynamically more stable than the DNA-DNA-DNA triplex, with the binding affinity of PNA-TFO to the PNA : DNA duplex being higher than that of DNA-TFO to the DNA : DNA duplex. The investigation of the association/dissociation of PNA-TFO to the PNA-DNA duplex has relevance and importance in the emerging effective applications of oligonucleotide therapy.

Simulations predict preferred Mg2+ coordination in a nonenzymatic primer-extension reaction center

The mechanism by which genetic information was copied prior to the evolution of ribozymes is of great interest because of its importance to the origin of life. The most effective known process for the nonenzymatic copying of an RNA template is primer extension by a two-step pathway in which 2-aminoimidazole-activated nucleotides first react with each other to form an imidazolium-bridged intermediate that subsequently reacts with the primer. Reaction kinetics, structure-activity relationships, and X-ray crystallography have provided insight into the overall reaction mechanism, but many puzzles remain. In particular, high concentrations of Mg2+ are required for efficient primer extension, but the mechanism by which Mg2+ accelerates primer extension remains unknown. By analogy with the mechanism of DNA and RNA polymerases, a role for Mg2+ in facilitating the deprotonation of the primer 3'-hydroxyl is often assumed, but no catalytic metal ion is seen in crystal structures of the primer-extension complex. To explore the potential effects of Mg2+ binding in the reaction center, we performed atomistic molecular dynamics simulations of a series of modeled complexes in which a Mg2+ ion was placed in the reaction center with inner-sphere coordination with different sets of functional groups. Our simulations suggest that coordination of a Mg2+ ion with both O3' of the terminal primer nucleotide and the pro-Sp nonbridging oxygen of the reactive phosphate of an imidazolium-bridged dinucleotide would help to pre-organize the structure of the primer/template substrate complex to favor the primer-extension reaction. Our results suggest that the catalytic metal ion may play an important role in overcoming electrostatic repulsion between a deprotonated O3' and the reactive phosphate of the bridged dinucleotide and lead to testable predictions of the mode of Mg2+ binding that is most relevant to catalysis of primer extension.

Isotope-Edited Variable Temperature Infrared Spectroscopy for Measuring Transition Temperatures of Single A-T Watson-Crick Base Pairs in DNA Duplexes

Experimental methods to determine transition temperatures for individual base pair melting events in DNA duplexes are lacking despite intense interest in these thermodynamic parameters. Here, we determine the dimensions of the thymine (T) C2═O stretching vibration when it is within the DNA duplex via isotopic substitutions at other atomic positions in the structure. First, we determined that this stretching state was localized enough to specific atoms in the molecule to make submolecular scale measurements of local structure and stability in high molecular weight complexes. Next, we develop a new isotope-edited variable temperature infrared method to measure melting transitions at various locations in a DNA structure. As an initial test of this "sub-molecular scale thermometer", we applied our T13C2 difference infrared signal to measure location-dependent melting temperatures (TmL) in a DNA duplex via variable temperature attenuated total reflectance Fourier transform infrared (VT-ATR-FTIR) spectroscopy. We report that the TmL of a single Watson-Crick A-T base pair near the end of an A-T rich sequence (poly T) is ∼34.9 ± 0.7°C. This is slightly lower than the TmL of a single base pair near the middle position of the poly T sequence (TmL ∼35.6±0.2°C). In addition, we also report that the TmL of a single Watson-Crick A-T base pair near the end of a 50% G-C sequence (12-mer) is ∼52.5 ± 0.3°C, which is slightly lower than the global melting Tm of the 12-mer sequence (TmL ∼54.0±0.9°C). Our results provide direct physical evidence for end fraying in DNA sequences with our novel spectroscopic methods.

Site-Specific Investigation of DNA Holliday Junction Dynamics and Structure with 6-Methylisoxanthopterin, a Fluorescent Guanine Analog

DNA Holliday Junction (HJ) formation and resolution is requisite for maintaining genomic stability in processes such as replication fork reversal and double-strand break repair. If HJs are not resolved, chromosome disjunction and aneuploidy result, hallmarks of tumor cells. To understand the structural features that lead to processing of these four-stranded joint molecule structures, we seek to identify structural and dynamic features unique to the central junction core. We incorporate the fluorescent guanine analog 6-methylisoxanthopterin (6-MI) at ten different locations throughout a model HJ structure to obtain site-specific information regarding the structure and dynamics of bases relative to those in a comparable sequence context in duplex DNA. These comparisons were accomplished through measuring fluorescence lifetime, relative brightness, fluorescence anisotropy, and thermodynamic stability, along with fluorescence quenching assays. These time-resolved and steady-state fluorescence measurements demonstrate that the structural distortions imposed by strand crossing result in increased solvent exposure, less stacking of bases and greater extrahelical nature of bases within the junction core. The 6-MI base analogs in the junction reflect these structural changes through an increase in intensity relative to those in the duplex. Molecular dynamics simulations performed using a model HJ indicate the primary sources of deformation are in the shift and twist parameters of the bases at the central junction step. These results suggest that junction-binding proteins may use the unique structure and dynamics of the bases at the core for recognition.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD

Molecular dynamics simulations can be used in combination with experimental techniques to uncover the intricacies of biomolecular structure, dynamics, and the resulting interactions. However, many noncanonical nucleic acid structures have proven to be challenging to replicate in accurate agreement with experimental data, often attributed to known force field deficiencies. A common force field criticism is the handling of van der Waals (vdW) parameters, which have not been updated since the regular use of Ewald's methods became routine. This work dives into the effects of minute vdW radii shifts on RNA tetranucleotide, B-DNA, and Z-DNA model systems described by commonly used Amber force fields. Using multidimensional replica exchange molecular dynamics (M-REMD), the GACC RNA tetranucleotide demonstrated changes in the structural distribution between the NMR minor and anomalous structure populations based on the O2' vdW radii scanning. However, no significant change in the NMR Major conformation population was observed. There were minimal changes in the B-DNA structure but there were more substantial improvements in Z-DNA structural descriptions, specifically with the Tumuc1 force field. This occurred with both LJbb vdW radii adjustments and incorporation of the CUFIX nonbonded parameter modifications. Though the limited vdW modifications tested did not provide a universal fix to the challenge of simulating the various known nucleic acid structures, they do provide direction and a greater understanding for future force field development efforts.

Diversity of Self-Assembled RNA Complexes: From Nanoarchitecture to Nanomachines

New tool development for various nucleic acid applications is an essential task in RNA nanotechnology. Here, we determined the ability of self-limited complex formation by a pair of oligoribonucleotides carrying two pairwise complementary blocks connected by a linker of different lengths in each chain. The complexes were analyzed using UV melting, gel shift assay analysis, and molecular dynamics (MD) simulations. We have demonstrated the spontaneous formation of various self-limited and concatemer complexes. The linear concatemer complex is formed by a pair of oligomers without a linker in at least one of them. Longer linkers resulted in the formation of circular complexes. The self-limited complexes formation was confirmed using the toehold strand displacement. The MD simulations indicate the reliability of the complexes' structure and demonstrate their dynamics, which increase with the rise of complex size. The linearization of 2D circular complexes into 1D structures and a reverse cyclization process were demonstrated using a toehold-mediated approach. The approach proposed here for the construction and directed modification of the molecularity and shape of complexes will be a valuable tool in RNA nanotechnology, especially for the rational design of therapeutic nucleic acids with high target specificity and the programmable response of the immune system of organisms.

Twist-stretch relations in nucleic acids

Nucleic acids are highly deformable helical molecules constantly stretched, twisted and bent in their biological functioning. Single molecule experiments have shown that double stranded (ds)-RNA and standard ds-DNA have opposite twist-stretch patterns and stretching properties when overwound under a constant applied load. The key structural features of the A-form RNA and B-form DNA helices are here incorporated in a three-dimensional mesoscopic Hamiltonian model which accounts for the radial, bending and twisting fluctuations of the base pairs. Using path integral techniques which sum over the ensemble of the base pair fluctuations, I compute the average helical repeat of the molecules as a function of the load. The obtained twist-stretch relations and stretching properties, for short A- and B-helical fragments, are consistent with the opposite behaviors observed in kilo-base long molecules.

19F electron nuclear double resonance (ENDOR) spectroscopy for distance measurements using trityl spin labels in DNA duplexes

The combination of fluorine labeling and pulsed electron-nuclear double resonance (ENDOR) is emerging as a powerful technique for obtaining structural information about proteins and nucleic acids. In this work, we explored the capability of Mims 19F ENDOR experiments on reporting intermolecular distances in trityl- and 19F-labeled DNA duplexes at three electron paramagnetic resonance (EPR) frequencies (34, 94, and 263 GHz). For spin labeling, we used the hydrophobic Finland trityl radical and hydrophilic OX063 trityl radical. Fluorine labels were introduced into two positions of a DNA oligonucleotide. The results indicated that hyperfine splittings visible in the ENDOR spectra are consistent with the most populated interspin distances between 19F and the trityl radical predicted from molecular dynamic (MD) simulations. Moreover, for some cases, ENDOR spectral simulations based on MD results were able to reproduce the conformational distribution reflected in the experimental ENDOR line broadening. Additionally, MD simulations provided more detailed information about the melting of terminal base pairs of the oligonucleotides and about the configuration of the trityls relative to a DNA end.

Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation

RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.

Can molecular dynamics be used to simulate biomolecular recognition?

There are many problems in biochemistry that are difficult to study experimentally. Simulation methods are appealing due to direct availability of atomic coordinates as a function of time. However, direct molecular simulations are challenged by the size of systems and the time scales needed to describe relevant motions. In theory, enhanced sampling algorithms can help to overcome some of the limitations of molecular simulations. Here, we discuss a problem in biochemistry that offers a significant challenge for enhanced sampling methods and that could, therefore, serve as a benchmark for comparing approaches that use machine learning to find suitable collective variables. In particular, we study the transitions LacI undergoes upon moving between being non-specifically and specifically bound to DNA. Many degrees of freedom change during this transition and that the transition does not occur reversibly in simulations if only a subset of these degrees of freedom are biased. We also explain why this problem is so important to biologists and the transformative impact that a simulation of it would have on the understanding of DNA regulation.

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Cation valence dependence of hydrogen bond and stacking potentials in DNA mesoscopic models

Monovalent and divalent cations play a crucial role in living cells and for molecular techniques such as PCR. Here we evaluate DNA melting temperatures in magnesium (Mg²⁺) and magnesium‑potassium (Mg²⁺+ K⁺) buffers with a mesoscopic model that allows us to estimate hydrogen bonds and stacking interaction potentials. The Mg²⁺ and Mg²⁺+ K⁺ results are compared to previous calculations for sodium ions (Na⁺), in terms of equivalent sodium concentration and ionic strength. Morse potentials, related to hydrogen bonding, were found to be essentially constant and unaffected by cation conditions. However, for stacking interactions we find a clear dependence with ionic strength and cation valence. The highest ionic strength variations, for both hydrogen bonds and stacking interactions, was found at the sequence terminals. This suggests that end-to-end interactions in DNA will be strongly dependent on cation valence and ionic strength.

Thermodynamics and kinetics of an A-U RNA base pair under force studied by molecular dynamics simulations

Mechanical force has been widely used to study RNA folding and unfolding. Understanding how the force affects the opening and closing of a single base pair, which is a basic step for RNA folding and unfolding and a fundamental behavior in some important biological activities, is crucial to understanding the mechanism of RNA folding and unfolding under mechanical force. In this work, we investigated the opening and closing process of an RNA base pair under mechanical force with constant-force stretching molecular dynamics simulations. It was found that high mechanical force results in overstretching, and the open state is a high-energy state. The enthalpy and entropy change of the base-pair opening-closing transition were obtained and the results at low forces were in good agreement with the nearest-neighbor model. The temperature and force dependence of the opening and closing rates were also obtained. The position of the transition state for the base-pair opening-closing transition under mechanical force was determined. The free energy barrier of opening a base pair without force is the enthalpy increase, and the work done by the force from the closed state to the transition state decreases the barrier and increases the opening rate. The free energy barrier of closing the base pair without force results from the entropy loss, and the work done by the force from the open state to the transition state increases the barrier and decreases the closing rate. The transition rates are strongly dependent on the temperature and force, while the transition path times are weakly dependent on force and temperature.

Maturation of siRNA by strand separation: Steered molecular dynamics study

RNA interference, particularly siRNA induced gene silencing is becoming an important avenue of modern therapeutics. The siRNA is delivered to the cells as short double helical RNA which becomes single stranded for forming the RISC complex. Significant experimental evidence is available for most of the steps except the process of the separation of the two strands. We have attempted to understand the pathway for double stranded siRNA (dsRNA) to single stranded (ssRNA) molecules using steered molecular dynamics simulations. As the process is completely unexplored we have applied force from all possible directions restraining all possible residues to convert dsRNA to ssRNA. We found pulling one strand along the helical axis direction restraining the far end of the other strand demands excessive force for ssRNA formation. Pulling a central residue of one strand, in a direction perpendicular to the helix axis, while keeping the base paired residue fixed requires intermediate force for strand separation. Moreover, we found that in this process the force requirement is quite high for the first bubble formation (nucleation energy) and the bubble propagation energies are quite small. We believe the success rate of the design of siRNA sequences for gene silencing may increase if this mechanistic knowledge is utilized for such a design process.Communicated by Ramaswamy H. Sarma.

A Simulation Study of the Effect of Naturally Occurring Point Mutations on the SRY-DNA Complex

Molecular dynamics (MD) simulations were conducted in order to investigate the effect of the naturally occurring point mutations of the transcription factor (TF) sex-determining region Y (SRY) on the structure and dynamics of the SRY-DNA complex. The normal SRY, along with the two mutants I13T and G40R, comprising point mutations on the SRY chain, which have been clinically identified in patients with sex developmental disorders, were modeled as DNA complexes. Our modeling work aims at elucidating atomic-level structural determinants of the aberrant SRY-DNA complexation by means of μs-long MD. The results suggest that the observed disorders brought about by the G40R-DNA and I13T-DNA may arise predominantly from the destabilization of the complex being in accord with in vitro assays found elsewhere and from modifications of the DNA bending as revealed in this study. Comparative potential of mean force computations, over a sequence of short separation distances for the three complexes, verified a higher stability of the normal SRY-DNA. Examining the way the SRY mutations modulate the SRY-DNA complex dynamics at the microscopic level is important also toward elucidating molecular determinants of function for proteins capable of binding to DNA.

Probabilistic Interpretation of NMR J-Couplings Determines BI-BII State Equilibria in DNA

Interpretation of ³J_P,H3' NMR scalar spin-spin coupling constants in DNA becomes more reliable by including distinct structural states such as BI and BII, using the weighted-static or, better still, the recently implemented adiabatic-MD (Ad-MD) method. The calculation method employs an adiabatic ("Ad") dependence of ³J_P,H3' coupling on NMR-assigned torsion angle, ε, weighted by P(ε) probability distribution calculated by molecular dynamics (MD). Ad-MD calculations enable cross-validation of the bsc1, OL15, and OL21 force fields and various parametrizations of the Karplus equation describing the dependence of ³J_P,H3' coupling on ε torsion (KE). The mean absolute deviation of Ad-MD ³J_P,H3' couplings from the experimental values in Dickerson-Drew DNA is comparable to the scatter of ³J_P,H3' couplings among four separate NMR experiments. A commonly accepted assumption of homogeneity of one kind of structure-dynamic state within DNA (BI or BII) is questionable because the principal characteristics of relevant P(ε) probabilities (shapes and positioning) vary with DNA sequence. The theory outlined in the present work sets limits to future reparameterization of MD force fields, as relevant to NMR data.

A single-molecule approach to unravel the molecular mechanism of the action of Deinococcus radiodurans RecD2 and its interaction with SSB and RecA in DNA repair

Helicases are ATP-driven molecular machines that directionally remodel nucleic acid polymers in all three domains of life. They are responsible for resolving double-stranded DNA (dsDNA) into single-strands, which is essential for DNA replication, nucleotide excision repair, and homologous recombination. RecD2 from Deinococcus radiodurans (DrRecD2) has important contributions to the organism's unusually high tolerance to gamma radiation and hydrogen peroxide. Although the results from X-ray Crystallography studies have revealed the structural characteristics of the protein, direct experimental evidence regarding the dynamics of the DNA unwinding process by DrRecD2 in the context of other accessory proteins is yet to be found. In this study, we have probed the exact binding event and processivity of DrRecD2 at single-molecule resolution using Protein-induced fluorescence enhancement (smPIFE) and Forster resonance energy transfer (smFRET). We have found that the protein prefers to bind at the 5' terminal end of the single-stranded DNA (ssDNA) by Drift and has helicase activity even in absence of ATP. However, a faster and iterative mode of DNA unwinding was evident in presence of ATP. The rate of translocation of the protein was found to be slower on dsDNA compared to ssDNA. We also showed that DrRecD2 is recruited at the binding site by the single-strand binding protein (SSB) and during the unwinding, it can displace RecA from ssDNA.

Accelerating the Ensemble Convergence of RNA Hairpin Simulations with a Replica Exchange Structure Reservoir

RNA is a key participant in many biological processes, but studies of RNA using computer simulations lag behind those of proteins, largely due to less-developed force fields and the slow dynamics of RNA. Generating converged RNA ensembles for force field development and other studies remains a challenge. In this study, we explore the ability of replica exchange molecular dynamics to obtain well-converged conformational ensembles for two RNA hairpin systems in an implicit solvent. Even for these small model systems, standard REMD remains computationally costly, but coupling to a pre-generated structure library using the reservoir REMD approach provides a dramatic acceleration of ensemble convergence for both model systems. Such precise ensembles could facilitate RNA force field development and validation and applications of simulation to more complex RNA systems. The advantages and remaining challenges of applying R-REMD to RNA are investigated in detail.

Self-Assembly of All-DNA Rods with Controlled Patchiness

Double-stranded DNA (dsDNA) fragments exhibit noncovalent attractive interactions between their tips. It is still unclear how DNA liquid crystal self-assembly is affected by such blunt-end attractions. It is demonstrated that stiff dsDNA fragments with moderate aspect ratio can specifically self-assemble in concentrated aqueous solutions into different types of smectic mesophases on the basis of selectively screening of blunt-end DNA stacking interactions. To this end, this type of attractions are engineered at the molecular level by constructing DNA duplexes where the attractions between one or both ends are screened by short hairpin caps. All-DNA bilayer and monolayer smectic-A type of phases, as well as a columnar phase, can be stabilized by controlling attractions strength. The results imply that the so far elusive smectic-A in DNA rod-like liquid crystals is a thermodynamically stable phase. The existence of the bilayer smectic phase is confirmed by Monte-Carlo simulations of hard cylinders decorated with one attractive terminal site. This work demonstrates that DNA blunt-ends behave as well-defined monovalent attractive patches whose strength and position can be potentially precisely tuned, highlighting unique opportunities concerning the stabilization of nonconventional DNA-based lyotropic liquid crystal phases assembled by all-DNA patchy particles with arbitrary geometry and composition.

The Ad-MD method to calculate NMR shift including effects due to conformational dynamics: The 31 P NMR shift in DNA

A method for averaging of NMR parameters by molecular dynamics (MD) has been derived from the method of statistical averaging in MD snapshots, benchmarked and applied to structurally dynamic interpretation of the ³¹ P NMR shift (δ_31P ) in DNA phosphates. The method employs adiabatic dependence of an NMR parameter on selected geometric parameter(s) that is weighted by MD-calculated probability distribution(s) for the geometric parameter(s) (Ad-MD method). The usage of Ad-MD for polymers is computationally convenient when one pre-calculated structural dependence of an NMR parameter is employed for all chemically equivalent units differing only in dynamic behavior. The Ad-MD method is benchmarked against the statistical averaging method for δ_31P in the model phosphates featuring distinctively different structures and dynamic behavior. The applicability of Ad-MD is illustrated by calculating ³¹ P NMR spectra in the Dickerson-Drew DNA dodecamer. δ_31P was calculated with the B3LYP/IGLO-III/PCM(water) and the probability distributions for the torsion angles adjacent to the phosphorus atoms in the DNA phosphates were calculated using the OL15 force field.

RNA-Specific Force Field Optimization with CMAP and Reweighting

RNA plays a key role in a variety of cell activities. However, it is difficult to capture its structure dynamics by the traditional experimental methods because of the inherent limitations. Molecular dynamics simulation has become a valuable complement to the experimental methods. Previous studies have indicated that the current force fields cannot accurately reproduce the conformations and structural dynamics of RNA. Therefore, an RNA-specific force field was developed to improve the conformation sampling of RNA. The distribution of ζ/α dihedrals of tetranucleotides was optimized by a reweighting method, and the grid-based energy correction map (CMAP) term was first introduced into the Amber RNA force field of ff99bsc0χOL3, named ff99OL3_CMAP1. Extensive validations of tetranucleotides and tetraloops show that ff99OL3_CMAP1 can significantly decrease the population of an incorrect structure, increase the consistency between the simulation results and experimental values for tetranucleotides, and improve the stability of tetraloops. ff99OL3_CMAP1 can also precisely reproduce the conformation of a duplex and riboswitches. These findings confirm that the newly developed force field ff99OL3_CMAP1 can improve the conformer sampling of RNA.

Chemical reactivity and binding interactions in ribonucleic acid-peptide complexes

The covalent and noncovalent backbone binding interactions in RNA-peptide complexes were studied by DFT methods. Four RNA structures R1(GGCUAGCC), R2(AAUCGAUU), R3(GGGAUCCC), and R4(AAAGCUUU) has been selected for eight protonated peptides (DR, ER, GR, KR, NGR, RR, tmeGnd (tme), and VR) interactions based on an experimental study (Anal Chem. 2019; 91:1659-1664). Chemical reactivity theory is used to study the reactivity of eight peptides with global descriptors. Lower hardness values reflected low stability and high reactivity for the protonated peptides. DR, ER, GR, KR, NGR, RR, and VR show lower value of ω, μ while tme has high value of ω, μ. Larger highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap for ER, GR, and KR showed greater structural stability for peptides. AutoDock and PatchDock results indicated that R1, R2, and R4 retain hairpin structures while interacting with peptide complexes. The calculated binding energies of (R1-R4)-peptide complexes from AutoDock tools are (1.49-11.12) kcal/mol. Results showed that the noncovalent interactions are stronger than the covalent interactions for R1-peptide complexes. The reason might be the transfer of proton from protonated ligand to deprotonated RNA, which has initiated the loss of the ligand. Also it has been observed that proton transfer has become energetically unfavorable in presence of additional hydrogen bonds which is predicted in the experimental results.

Determining Sequence-Dependent DNA Oligonucleotide Hybridization and Dehybridization Mechanisms Using Coarse-Grained Molecular Simulation, Markov State Models, and Infrared Spectroscopy

A robust understanding of the sequence-dependent thermodynamics of DNA hybridization has enabled rapid advances in DNA nanotechnology. A fundamental understanding of the sequence-dependent kinetics and mechanisms of hybridization and dehybridization remains comparatively underdeveloped. In this work, we establish new understanding of the sequence-dependent hybridization/dehybridization kinetics and mechanism within a family of self-complementary pairs of 10-mer DNA oligomers by integrating coarse-grained molecular simulation, machine learning of the slow dynamical modes, data-driven inference of long-time kinetic models, and experimental temperature-jump infrared spectroscopy. For a repetitive ATATATATAT sequence, we resolve a rugged dynamical landscape comprising multiple metastable states, numerous competing hybridization/dehybridization pathways, and a spectrum of dynamical relaxations. Introduction of a G:C pair at the terminus (GATATATATC) or center (ATATGCATAT) of the sequence reduces the ruggedness of the dynamics landscape by eliminating a number of metastable states and reducing the number of competing dynamical pathways. Only by introducing a G:C pair midway between the terminus and the center to maximally disrupt the repetitive nature of the sequence (ATGATATCAT) do we recover a canonical "all-or-nothing" two-state model of hybridization/dehybridization with no intermediate metastable states. Our results establish new understanding of the dynamical richness of sequence-dependent kinetics and mechanisms of DNA hybridization/dehybridization by furnishing quantitative and predictive kinetic models of the dynamical transition network between metastable states, present a molecular basis with which to understand experimental temperature jump data, and furnish foundational design rules by which to rationally engineer the kinetics and pathways of DNA association and dissociation for DNA nanotechnology applications.

Z-DNA as a Touchstone for Additive Empirical Force Fields and a Refinement of the Alpha/Gamma DNA Torsions for AMBER

Although current AMBER force fields are relatively accurate for canonical B-DNA, many noncanonical structures are still described incorrectly. As noncanonical motifs are attracting increasing attention due to the role they play in living organisms, further improvement is desirable. Here, we have chosen the Z-DNA molecule, which can be considered a touchstone of the universality of empirical force fields, since the noncanonical α and γ backbone conformations native to Z-DNA are also found in protein-DNA complexes, i-motif DNA, and other noncanonical DNAs. We show that spurious α/γ conformations occurring in simulations with current AMBER force fields, OL15 and bsc1, are largely due to inaccurate α/γ parametrization. Moreover, stabilization of native Z-DNA substates involving γ = trans conformations appears to be in conflict with the correct description of the canonical B-DNA structure. Because the balance of the native and spurious conformations is influenced by nonadditive effects, this is a difficult case for an additive dihedral energy scheme such as AMBER. We propose new α/γ parameters, denoted OL21, and show that they improve the stability of native α/γ Z-DNA substates while keeping the canonical DNA description virtually unchanged, thus representing a reasonable compromise within the additive force field framework. Although further extensive testing is needed, the new modification appears to be a promising step toward a more reliable description of noncanonical DNA motifs and provides the best performance for Z-DNA molecules among current AMBER force fields.

The Time Scale of Electronic Resonance in Oxidized DNA as Modulated by Solvent Response: An MD/QM-MM Study

The time needed to establish electronic resonant conditions for charge transfer in oxidized DNA has been evaluated by molecular dynamics simulations followed by QM/MM computations which include counterions and a realistic solvation shell. The solvent response is predicted to take ca. 800-1000 ps to bring two guanine sites into resonance, a range of values in reasonable agreement with the estimate previously obtained by a kinetic model able to correctly reproduce the observed yield ratios of oxidative damage for several sequences of oxidized DNA.

Base pair fluctuations in helical models for nucleic acids

A statistical method is developed to estimate the maximum amplitude of the base pair fluctuations in a three dimensional mesoscopic model for nucleic acids. The base pair thermal vibrations around the helix diameter are viewed as a Brownian motion for a particle embedded in a stable helical structure. The probability to return to the initial position is computed, as a function of time, by integrating over the particle paths consistent with the physical properties of the model potential. The zero time condition for the first-passage probability defines the constraint to select the integral cutoff for various macroscopic helical conformations, obtained by tuning the twist, bending, and slide motion between adjacent base pairs along the molecule stack. Applying the method to a short homogeneous chain at room temperature, we obtain meaningful estimates for the maximum fluctuations in the twist conformation with ∼10.5 base pairs per helix turn, typical of double stranded DNA helices. Untwisting the double helix, the base pair fluctuations broaden and the integral cutoff increases. The cutoff is found to increase also in the presence of a sliding motion, which shortens the helix contour length, a situation peculiar of dsRNA molecules.

The N-terminal Helix-Turn-Helix Motif of Transcription Factors MarA and Rob Drives DNA Recognition

DNA-binding proteins play an important role in gene regulation and cellular function. The transcription factors MarA and Rob are two homologous members of the AraC/XylS family that regulate multidrug resistance. They share a common DNA-binding domain, and Rob possesses an additional C-terminal domain that permits binding of low-molecular weight effectors. Both proteins possess two helix-turn-helix (HTH) motifs capable of binding DNA; however, while MarA interacts with its promoter through both HTH-motifs, prior studies indicate that Rob binding to DNA via a single HTH-motif is sufficient for tight binding. In the present work, we perform microsecond time scale all-atom simulations of the binding of both transcription factors to different DNA sequences to understand the determinants of DNA recognition and binding. Our simulations characterize sequence-dependent changes in dynamical behavior upon DNA binding, showcasing the role of Arg40 of the N-terminal HTH-motif in allowing for specific tight binding. Finally, our simulations demonstrate that an acidic C-terminal loop of Rob can control the DNA binding mode, facilitating interconversion between the distinct DNA binding modes observed in MarA and Rob. In doing so, we provide detailed molecular insight into DNA binding and recognition by these proteins, which in turn is an important step toward the efficient design of antivirulence agents that target these proteins.

Stable DNA Sequence Over Close-Ending and Pairing Sequences Constraint

DNA computing is a new method based on molecular biotechnology to solve complex problems. The design of DNA sequences is a multi-objective optimization problem in DNA computing, whose objective is to obtain optimized sequences that satisfy multiple constraints to improve the quality of the sequences. However, the previous optimized DNA sequences reacted with each other, which reduced the number of DNA sequences that could be used for molecular hybridization in the solution and thus reduced the accuracy of DNA computing. In addition, a DNA sequence and its complement follow the principle of complementary pairing, and the sequence of base GC at both ends is more stable. To optimize the above problems, the constraints of Pairing Sequences Constraint (PSC) and Close-ending along with the Improved Chaos Whale (ICW) optimization algorithm were proposed to construct a DNA sequence set that satisfies the combination of constraints. The ICW optimization algorithm is added to a new predator–prey strategy and sine and cosine functions under the action of chaos. Compared with other algorithms, among the 23 benchmark functions, the new algorithm obtained the minimum value for one-third of the functions and two-thirds of the current minimum value. The DNA sequences satisfying the constraint combination obtained the minimum of fitness values and had stable and usable structures.

W-RESP: Well-Restrained Electrostatic Potential-Derived Charges. Revisiting the Charge Derivation Model

Representation of electrostatic interactions by a Coulombic pairwise potential between atom-centered partial charges is a fundamental and crucial part of empirical force fields used in classical molecular dynamics simulations. The broad success of the AMBER force-field family originates mainly from the restrained electrostatic potential (RESP) charge model, which derives partial charges to reproduce the electrostatic field around the molecules. However, the description of the electrostatic potential around molecules by standard RESP may be biased for some types of molecules. In this study, we modified the RESP charge derivation model to improve its description of the electrostatic potential around molecules and thus electrostatic interactions in the force field. In particular, we reoptimized the atomic radii for definition of the grid points around the molecule, redesigned the restraining scheme, and included extra point (EP) charges. The RESP fitting was significantly improved for aromatic heterocyclic molecules. Thus, the suggested W-RESP(-EP) charge derivation model shows some potential for improving the performance of the nucleic acid force fields, for which the poor description of nonbonded interactions, such as the underestimated stability of base pairing, is well-established. We also report some preliminary simulation tests (around 1 ms of simulation data) on A-RNA duplexes, tetranucleotides, and tetraloops. The simulations reveal no adverse effects, while the description of base-pairing interactions might be improved. The new charges can thus be used in future attempts to improve the nucleic acid simulation force fields, in combination with reparametrization of the other terms.

Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations

Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)₂ also known as the Drew–Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 μs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.

Impact of DNA Adduct Size, Number, and Relative Position on the Toxicity of Aromatic Amines: A Molecular Dynamics Case Study of ANdG- and APdG-Containing DNA Duplexes

Human exposure to aromatic amines (AAs) can result in carcinogenic DNA adducts. To complement previous work geared toward understanding the mutagenicity of AA-derived adducts, which has almost exclusively studied (monoadducted) DNA containing a single lesion, the present work provides the first in-depth comparison of the structure of monoadducted and diadducted DNA duplexes. Specifically, molecular dynamics (MD) simulations were initially performed on DNA containing the nonmutagenic single-ringed N-(deoxyguanosin-8-yl)-aniline (^ANdG) or the mutagenic four-ringed N-(deoxyguanosin-8-yl)-1-aminopyrene (^APdG) lesion at G₁, G₂, or G₃ in the AA deletion hotspot (5'-G₁G₂CG₃CC) in the anti or syn glycosidic orientation (B/S duplex conformation). Subsequently, diadducted strands were assessed that span each combination of damaged sites (G₁G₂ (nearest neighbors), G₂G₃ (next-nearest neighbors), and G₁G₃ (two intervening nucleotides)) and anti/syn lesion glycosidic orientations. Despite other N-linked C8-dG adducts exhibiting sequence dependence conformational heterogeneity, a single ^ANdG or ^APdG lesion induces helical conformational homogeneity that is exclusively controlled by aryl moiety size. However, the preferred damaged DNA conformation can change upon the addition of a second adduct depending on lesion separation, with neighboring lesions stabilizing a nonmutagenic conformation and next-nearest damaged sites stabilizing a promutagenic conformation regardless of adduct size. As a result, diadducted DNA is found to adopt conformations that are unfavored for the corresponding monoadducted system, pointing to differential replication and repair outcomes for diadducted DNA compared to those for monoadducted DNA. Thus, although the toxicity of monoadducted DNA is most significantly dictated by lesion size, the toxicity can increase or decrease upon a second damaging event depending on lesion size and relative position. Overall, our work adds the number of lesions and their spatial separation to the growing list of factors that determine the structure and biological outcomes of adducted DNA.

Delocalization-Assisted Transport through Nucleic Acids in Molecular Junctions

The flow of charge through molecules is central to the function of supramolecular machines, and charge transport in nucleic acids is implicated in molecular signaling and DNA repair. We examine the transport of electrons through nucleic acids to understand the interplay of resonant and nonresonant charge carrier transport mechanisms. This study reports STM break junction measurements of peptide nucleic acids (PNAs) with a G-block structure and contrasts the findings with previous results for DNA duplexes. The conductance of G-block PNA duplexes is much higher than that of the corresponding DNA duplexes of the same sequence; however, they do not display the strong even-odd dependence conductance oscillations found in G-block DNA. Theoretical analysis finds that the conductance oscillation magnitude in PNA is suppressed because of the increased level of electronic coupling interaction between G-blocks in PNA and the stronger PNA-electrode interaction compared to that in DNA duplexes. The strong interactions in the G-block PNA duplexes produce molecular conductances as high as 3% G0, where G0 is the quantum of conductance, for 5 nm duplexes.

First passage time study of DNA strand displacement

DNA strand displacement, in which a single-stranded nucleic acid invades a DNA duplex, is pervasive in genomic processes and DNA engineering applications. The kinetics of strand displacement have been studied in bulk; however, the kinetics of the underlying strand exchange were obfuscated by a slow bimolecular association step. Here, we use a novel single-molecule fluorescence resonance energy transfer approach termed the "fission" assay to obtain the full distribution of first passage times of unimolecular strand displacement. At a frame time of 4.4 ms, the first passage time distribution for a 14-nucleotide displacement domain exhibited a nearly monotonic decay with little delay. Among the eight different sequences we tested, the mean displacement time was on average 35 ms and varied by up to a factor of 13. The measured displacement kinetics also varied between complementary invaders and between RNA and DNA invaders of the same base sequence, except for T → U substitution. However, displacement times were largely insensitive to the monovalent salt concentration in the range of 0.25-1 M. Using a one-dimensional random walk model, we infer that the single-step displacement time is in the range of ∼30-300 μs, depending on the base identity. The framework presented here is broadly applicable to the kinetic analysis of multistep processes investigated at the single-molecule level.

Secondary structural choice of DNA and RNA associated with CGG/CCG trinucleotide repeat expansion rationalizes the RNA misprocessing in FXTAS

CGG tandem repeat expansion in the 5'-untranslated region of the fragile X mental retardation-1 (FMR1) gene leads to unusual nucleic acid conformations, hence causing genetic instabilities. We show that the number of G…G (in CGG repeat) or C…C (in CCG repeat) mismatches (other than A…T, T…A, C…G and G…C canonical base pairs) dictates the secondary structural choice of the sense and antisense strands of the FMR1 gene and their corresponding transcripts in fragile X-associated tremor/ataxia syndrome (FXTAS). The circular dichroism (CD) spectra and electrophoretic mobility shift assay (EMSA) reveal that CGG DNA (sense strand of the FMR1 gene) and its transcript favor a quadruplex structure. CD, EMSA and molecular dynamics (MD) simulations also show that more than four C…C mismatches cannot be accommodated in the RNA duplex consisting of the CCG repeat (antisense transcript); instead, it favors an i-motif conformational intermediate. Such a preference for unusual secondary structures provides a convincing justification for the RNA foci formation due to the sequestration of RNA-binding proteins to the bidirectional transcripts and the repeat-associated non-AUG translation that are observed in FXTAS. The results presented here also suggest that small molecule modulators that can destabilize FMR1 CGG DNA and RNA quadruplex structures could be promising candidates for treating FXTAS.

Targeting structural features of viral genomes with a nano-sized supramolecular drug

RNA targeting is an exciting frontier for drug design. Intriguing targets include functional RNA structures in structurally-conserved untranslated regions (UTRs) of many lethal viruses. However, computational docking screens, valuable in protein structure targeting, fail for inherently flexible RNA. Herein we harness MD simulations with Markov state modeling to enable nanosize metallo-supramolecular cylinders to explore the dynamic RNA conformational landscape of HIV-1 TAR untranslated region RNA (representative for many viruses) replicating experimental observations. These cylinders are exciting as they have unprecedented nucleic acid binding and are the first supramolecular helicates shown to have anti-viral activity in cellulo: the approach developed in this study provides additional new insight about how such viral UTR structures might be targeted with the cylinder binding into the heart of an RNA-bulge cavity, how that reduces the conformational flexibility of the RNA and molecular details of the insertion mechanism. The approach and understanding developed represents a new roadmap for design of supramolecular drugs to target RNA structural motifs across biology and nucleic acid nanoscience.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments

Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Investigation into the mechanism and dynamics of DNA association and dissociation utilizing kinetic Monte Carlo simulations

In this work, we present a kinetic Markov state Monte Carlo model designed to complement temperature-jump (T-jump) infrared spectroscopy experiments probing the kinetics and dynamics of short DNA oligonucleotides. The model is designed to be accessible to experimental researchers in terms of both computational simplicity and expense while providing detailed insights beyond those provided by experimental methods. The model is an extension of a thermodynamic lattice model for DNA hybridization utilizing the formalism of the nucleation-zipper mechanism. Association and dissociation trajectories were generated utilizing the Gillespie algorithm and parameters determined via fitting the association and dissociation timescales to previously published experimental data. Terminal end fraying, experimentally observed following a rapid T-jump, in the sequence 5'-ATATGCATAT-3' was replicated by the model that also demonstrated that experimentally observed fast dynamics in the sequences 5'-C(AT)_nG-3', where n = 2-6, were also due to terminal end fraying. The dominant association pathways, isolated by transition pathway theory, showed two primary motifs: initiating at or next to a G:C base pair, which is enthalpically favorable and related to the increased strength of G:C base pairs, and initiating in the center of the sequence, which is entropically favorable and related to minimizing the penalty associated with the decrease in configurational entropy due to hybridization.

Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer

Some DNA sequences in crystals and in complexes with proteins can exist in the forms intermediate between the B- and A-DNA. Based on this, it was implied that the B-to-A transition for any DNA molecule should go through these intermediate forms also in kinetics. More precisely, the helix parameter Slide has to change first, and the molecule should take the E-form. After that, the Roll parameter changes. In the present work, we simulated the kinetics of the B-A transition in the Drew-Dickerson dodecamer, a known B-philic DNA oligomer. We used the "sugar" coarse-grained model that reproduces ribose flexibility, preserves sequence specificity, employs implicit water and explicit ions, and offers the possibility to vary friction. As the control parameter of the transition, we chose the volume available for a counterion and considered the change from a large to a small volume. In the described system, the B-to-A conformational transformation proved to correspond to a first-order phase transition. The molecule behaves like a small cluster in the region of such a transition, jumping between the A- and B-forms in a wide range of available volumes. The viscosity of the solvent does not affect the midpoint of the transition but only the overall mobility of the system. All helix parameters change synchronously on average, we have not observed the sequence "Slide first, Roll later" in kinetics, and the E-DNA is not a necessary step for the transition between the B- and A-forms in the studied system. So, the existence of the intermediate DNA forms requires specific conditions, shifting the common balance of interactions: certain nucleotide sequence in specific solution or/and the interaction with some protein.

Molecular dynamics simulations of acyclic analogs of nucleic acids for antisense inhibition

For antisense applications, oligonucleotides must be chemically modified to be resistant to endogenous nucleases. Until now, antisense oligonucleotide (ASO) analogs have been synthesized and then tested for their ability to duplex with a target nucleic acid, usually RNA. In this work, using molecular dynamics calculations simulations, we systematically tested a series of chemically modified analogs in which the 2-deoxyribose was substituted for by one or two methylene groups on each side of the phosphate backbone, producing four compounds, of which three were previously unknown. We used a 9-mer sequence of which the solution structure has been determined by NMR spectroscopy and tested the ability to form stable duplexes of these acyclic analogs to both DNA and RNA. In only one case out of eight, we unexpectedly found the formation of a stable duplex with complementary RNA. We also applied limitations on end fraying because of the terminal AT base pairs, in order to eliminate this as a factor in the comparative results. We consider this a predictive method to potentially identify target ASO analogs for synthesis and testing for antisense drug development.