Physical rationale behind the nonlinear enthalpy-entropy compensation in DNA duplex stability

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

A Rational Study of the Origin and Generality of Anti-Enthalpy–Entropy Compensation (AEEC) Phenomenon

In addition to enthalpy–entropy compensation (EEC), anti-enthalpy–entropy compensation (AEEC) phenomenon is also found in literature. The reports on the latter are limited, and analyses and justifications are so far unclear. Herein we present demonstration on the nature and possibility of the AEEC phenomenon. Although literature reports so far have mostly shown linear AEEC, we have found both linear and non-linear dependences. The non-linearity we consider arises from large “free energy window” (FEW) like EEC, recently presented and discussed. Attempts have been made to rationalize the observations in terms of solvation–desolvation phenomenon of the involved processes anticipated in previous studies. We have found that enthalpy and entropy of formation of gaseous and solid materials may also exhibit AEEC, it can be thus classified as a global phenomenon. In addition to AEEC, the phenomenon of no-enthalpy–entropy compensation (NEEC) is also reported. Thus, the phenomena EEC, AEEC and NEEC are thermodynamic puzzles that require close attention and analysis. With the help of wide range of physical chemical processes, an elaborate general understanding of the linear and non-linear AEEC phenomena has been attempted. The experimental data herein used are collected from literature reports, new measurements were not done. A variety of examples have supported possible generality of AEEC.

Relating DNA base-pairing in aqueous media to DNA polymerase fidelity

Controversy surrounds the perceived absence of a relationship between DNA polymerase fidelity (kinetic discrimination) and free energy changes determined from DNA melting studies (thermodynamic discrimination). Thermodynamic discrimination together with aqueous solvent effects can account for kinetic fidelities on the order of those observed experimentally.

Evaluation of the Gibbs Free Energy Changes and Melting Temperatures of DNA/DNA Duplexes Using Hybridization Enthalpy Calculated by Molecular Dynamics Simulation

A molecular dynamics simulation approach was applied for the prediction of the thermal stability of oligonucleotide duplexes. It was shown that the enthalpy of the DNA/DNA complex formation could be calculated using this approach. We have studied the influence of various simulation parameters on the secondary structure and the hybridization enthalpy value of Dickerson-Drew dodecamer. The optimal simulation parameters for the most reliable prediction of the enthalpy values were determined. The thermodynamic parameters (enthalpy and entropy changes) of a duplex formation were obtained experimentally for 305 oligonucleotides of various lengths and GC-content. The resulting database was studied with molecular dynamics (MD) simulation using the optimized simulation parameters. Gibbs free energy changes and the melting temperatures were evaluated using the experimental correlation between enthalpy and entropy changes of the duplex formation and the enthalpy values calculated by the MD simulation. The average errors in the predictions of enthalpy, the Gibbs free energy change, and the melting temperature of oligonucleotide complexes were 11%, 10%, and 4.4 °C, respectively. We have shown that the molecular dynamics simulation gives a possibility to calculate the thermal stability of native DNA/DNA complexes a priori with an unexpectedly high accuracy.

Kinetics of DNA duplex formation: A-tracts versus AT-tracts

A-tracts (AAAA…:TTTT…) form much faster than AT-tracks (ATAT…:TATA…).

Hydrogen bonding or stacking interactions in differentiating duplex stability in oligonucleotides containing synthetic nucleoside probes for alkylated DNA

Oligonucleotides that hybridize to modified DNA are useful chemical tools to probe the noncovalent interactions that stabilize DNA duplexes. In an effort to better understand the interactions that influence the specificity of hybridization probes for O(6)-alkylguanine lesions, we examined a series of synthetic nucleoside analogues (BIM, Benzi, and Peri) with respect to their ability to stabilize duplex DNA comprised of native or damaged DNA oligonucleotides. The base-modified nucleoside analogues contained systematically varied hydrogen-bonding and π-stacking properties. The nucleoside probes were incorporated into DNA and paired opposite canonical bases (A, T, C, or G), O(6) -methylguanine (O(6)-MeG), O(6)-benzylguanine (O(6)-BnG), or a stable abasic site analogue (tetrahydrofuran, THF). On the basis of the free energy of duplex formation, the highest degree of stabilization was observed when Peri was paired opposite O(6)-MeG. The thermodynamic data suggest that the smaller probes stabilize DNA duplexes more through hydrogen bonding, whereas the larger probes, with a greater capacity to π stack, contribute to duplex stabilization more on the basis of base stacking. These results demonstrate that increased helix stability could be achieved when BIM, Benzi, or Peri were paired opposite damage-containing DNA rather than unmodified DNA (that is, O(6)-MeG rather than G). This knowledge is expected to be useful in the design and development of nucleoside analogues for uses in DNA-based technologies.

Thermodynamic studies of ionic hydration and interactions for amino acid ionic liquids in aqueous solutions at 298.15 K

Amino acid ionic liquids are a special class of ionic liquids due to their unique acid-base behavior, biological significance, and applications in different fields such as templates in synthetic chemistry, stabilizers for biological macromolecules, etc. The physicochemical properties of these ionic liquids can easily be altered by making the different combinations of amino acids as anion along with possible cation modification which makes amino acid ionic liquids more suitable to understand the different kinds of molecular and ionic interactions with sufficient depth so that they can provide fruitful information for a molecular level understanding of more complicated biological processes. In this context, volumetric and osmotic coefficient measurements for aqueous solutions containing 1-ethyl-3-methylimidazolium ([Emim]) based amino acid ionic liquids of glycine, alanine, valine, leucine, and isoleucine are reported at 298.15 K. From experimental osmotic coefficient data, mean molal activity coefficients of ionic liquids were estimated and analyzed using the Debye-Hückel and Pitzer models. The hydration numbers of ionic liquids in aqueous solutions were obtained using activity data. Pitzer ion interaction parameters are estimated and compared with other electrolytes reported in the literature. The nonelectrolyte contribution to the aqueous solutions containing ionic liquids was studied by calculating the osmotic second virial coefficient through an application of the McMillan-Mayer theory of solution. It has been found that the second osmotic virial coefficient which includes volume effects correlates linearly with the Pitzer ion interaction parameter estimated independently from osmotic data as well as the hydrophobicity of ionic liquids. The enthalpy-entropy compensation effect, explained using the Starikov-Nordén model of enthalpy-entropy compensation, and partial molar entropy analysis for aqueous [Emim][Gly] solutions are made by using experimental Gibb's free energy data and literature enthalpy data. This study highlights that the hydrophobic interaction persists even in the limit of infinite dilution where the hydration effects are usually dominant, implying importance of hydrophobic hydration. Analysis of the results further shows that the hydration of amino acid ionic liquids occurs through the cooperative H-bond formation with the kosmotropic effect in contrast to the usual inorganic salts or hydrophobic salts like tetraalkylammonium halides.

Stability and mismatch discrimination of locked nucleic acid-DNA duplexes

Locked nucleic acids (LNA; symbols of bases, +A, +C, +G, and +T) are introduced into chemically synthesized oligonucleotides to increase duplex stability and specificity. To understand these effects, we have determined thermodynamic parameters of consecutive LNA nucleotides. We present guidelines for the design of LNA oligonucleotides and introduce free online software that predicts the stability of any LNA duplex oligomer. Thermodynamic analysis shows that the single strand–duplex transition is characterized by a favorable enthalpic change and by an unfavorable loss of entropy. A single LNA modification confines the local conformation of nucleotides, causing a smaller, less unfavorable entropic loss when the single strand is restricted to the rigid duplex structure. Additional LNAs adjacent to the initial modification appear to enhance stacking and H-bonding interactions because they increase the enthalpic contributions to duplex stabilization. New nearest-neighbor parameters correctly forecast the positive and negative effects of LNAs on mismatch discrimination. Specificity is enhanced in a majority of sequences and is dependent on mismatch type and adjacent base pairs; the largest discriminatory boost occurs for the central +C·C mismatch within the +T+C+C sequence and the +A·G mismatch within the +T+A+G sequence. LNAs do not affect specificity in some sequences and even impair it for many +G·T and +C·A mismatches. The level of mismatch discrimination decreases the most for the central +G·T mismatch within the +G+G+C sequence and the +C·A mismatch within the +G+C+G sequence. We hypothesize that these discrimination changes are not unique features of LNAs but originate from the shift of the duplex conformation from B-form to A-form.