Atomistic Investigation of the Effect of Incremental Modification of

May 21, 2014 - algorithm and its use in nucleic acid simulations. Struct. Fold Des. 1999,. 7, R55−R60. (60) Humphrey, W.; Dalke, A.; Schulten, K. VM...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCB

Atomistic Investigation of the Effect of Incremental Modification of Deoxyribose Sugars by Locked Nucleic Acid (β‑D‑LNA and α‑L‑LNA) Moieties on the Structures and Thermodynamics of DNA−RNA Hybrid Duplexes Gorle Suresh and U. Deva Priyakumar* Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India S Supporting Information *

ABSTRACT: Chemically modified oligonucleotides offer many possibilities in utilizing their special features for a vast number of applications in nucleic acid based therapies and synthetic molecular biology. Locked nucleic acid analogues (α-/β-LNA) are modifications having an extra ring of 2′-O,4′-C-methylene group in the furanose sugar. LNA strands have been shown to exhibit high binding affinity toward RNA and DNA strands, and the resultant duplexes show significantly high melting temperatures. In the present study, molecular dynamics (MD) simulations were performed on DNA−RNA hybrid duplexes by systematically modifying their deoxyribose sugars with locked nucleic acid analogues. Several geometrical and energetic analyses were performed using principal component (PCA) analysis and binding free energy methods to understand the consequence of incorporated isomeric LNA modifications on the structure, dynamics, and stability of DNA− RNA hybrid duplex. The β-modification systematically changes the conformation of the DNA−RNA hybrid duplex whereas drastic changes are observed for α-modification. The fully modified duplexes have distinct properties compared to partial and unmodified duplexes, and the partly modified duplexes have properties intermediate to full strand and unmodified duplexes. The distribution of BI versus BII populations suggests that backbone rearrangement is minimal for β-LNA modification in order to accommodate it in duplexes whereas extensive backbone rearrangement is necessary in order to incorporate α-LNA modification which subsequently alters the energetic and structural properties of the duplexes. The simulation results also suggest that the alteration of DNA−RNA hybrid properties depends on the position of modification and the gap between the modifications.



toward North-form (N-form is characterized by C3′-endo and C2′-exo conformation of the sugar moiety, which is often found in A-type nucleic acids) similar to RNA duplexes.5,6 Besides high binding affinity, the single-stranded antisense oligonucleotides acting on the mRNA should also exhibit high biostability and efficient cellular uptake.7−9 Previous studies have shown that the above-mentioned features are attainable up to some extent by introducing chemical modifications at sugar and phosphate moieties. Some of the chemically modified AOs include 2′-O-methyl,10 locked nucleic acids (LNA),11−13 peptide nucleic acids (PNA),14 tricyclo-DNA,15 morpholinoNAs,16 hexitol nucleic acids (HNA),17 arabino nucleic acids (ANA), and 2′-F-arabino nucleic acids (2′-F-ANA).18 Locked nucleic acid (LNA) duplexes are chemically modified nucleic acids in which an extra ring is introduced in sugar connecting to C2′ and C4′ atoms with −O−CH2− group

INTRODUCTION Being central to transmission, expression, and conservation of genetic information, nucleic acids are perceived as effective drug targets for the treatment of many diseases.1−4 The drug therapies which target the nucleic acids (genetic level) attribute high generality and rationale than traditional drug therapies. One alternative approach is antisense therapy, where the introduced synthetic oligonucleotide (AO) binds to the mRNA with high affinity and stops the transmission of the genetic information by either forming stable duplex or degrading the target mRNA with the help of RNase H enzyme (this process is known as RNase H activity).1−4 The basic principle behind the antisense therapy is observed in nature during the gene regulation process.1,2 Since the hybrid duplexes (formed by pure DNA and RNA strands) exhibit low stability and moderate affinity toward RNase H activity, there is an immense interest in studying the chemically modified nucleic acid strands (AOs) to possibly increase binding affinity with complementary nucleic acid strands. Recent studies have been focusing on modifications which can drive the furanose sugar conformation © 2014 American Chemical Society

Received: February 11, 2014 Revised: May 5, 2014 Published: May 21, 2014 5853

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

systematic incorporation of β-LNA modifications in DNA strand turns the duplex conformation toward N-type.37 Similarly, NMR studies performed on 9-mer and 10-mer DNA duplexes showed that the conformation is largely retained after introducing up to three to four α-LNA modifications.38 Structural studies had been performed on DNA−RNA hybrid duplex by introducing three α-LNA modifications in the DNA strand and observed an intermediate A/B-like conformation similar to native hybrid, suggesting a similar behavior of deoxyribose sugar to α-LNA.20 Hence, the α-LNA and β-LNA modifications were proposed as DNA and RNA mimics, respectively, in the context of helical structure. However, a recent study showed that the structure of fully modified β-LNA duplex differs significantly compared to pure RNA duplex.39 The full β-LNA modification results in decreasing twist and roll and narrowing minor groove compared to pure RNA duplex. Classical molecular dynamics (MD) simulations have been shown to be valuable tools to study nucleic acid properties.40−46 Pande and co-workers have done computational study on the fully modified β-LNA:β-LNA duplexes and explained several aspects like hydration, conformation, and internal dynamics in comparison to regular A-RNA duplexes.40 The structural differences between full β-LNA and RNA duplexes were also observed. The observed geometric differences indicate a nonregular A-type geometry for full β-LNA duplex. The recent MD simulations study on full α-LNA and β-LNA duplexes suggest that the full β-LNA duplex has a conformation similar to the regular RNA duplex but the full α-LNA duplex exhibits a conformation which is distinct from typical A- and Btype conformations despite having similar thermodynamic stability.41 Despite having several studies on the LNA-modified nucleic acid duplexes, the variations in structure, dynamics, and stability of nucleic acid duplexes upon introducing α-/β-LNA modifications are still not well understood. For example, the NMR study on a DNA−RNA hybrid duplex suggests that when three α-LNA modifications are introduced, a structure similar to a DNA−RNA hybrid with a higher melting temperature results.20 To gain a more detailed insight into the structures, dynamics, flexibility, and stability of DNA−RNA hybrids upon introducing LNA analogues (α-LNA and β-LNA) into its DNA strand, MD simulations have been performed on DNA−RNA hybrid duplex by systematically modifying its deoxyribose sugars (Figure 1). This kind of study is indispensable to understand how the chemical modifications affect the nature of the nucleic acid duplexes.

(Figure 1). Locked nucleic acid analogues such as 2′-O,4′-Cmethylene-β-D-ribofuranosyl (β-D-LNA) and 2′-O,4′-C-meth-

Figure 1. (a) Chemical representation of α-LNA and β-LNA modifications. (b) Numbering of the DNA and RNA strands of DNA−RNA hybrid used during the simulations and various analyses. (c) Sequences used to model the α/β-LNA-modified DNA−RNA hybrid duplexes studied in the present study. The modified site is indicated in orange color.

ylene-α-L-ribofuranosyl (α-L-LNA) were shown to increase the thermodynamic stability of the DNA, RNA, and DNA−RNA hybrid duplexes when they are incorporated into strands and increase the melting temperatures up to ∼10 °C per modification.19,20 The high binding affinity of LNA analogues toward other nucleic acid strands make them useful in many nucleic acid research fields covering broad range of applications in therapeutics, biotechnology and molecular biology.19,21 LNAs are utilized in DNAzymes,22 siRNA approaches,23 and molecular beacons24 and in RNA in situ hybridization.25,26 LNA-modified oligonucleotides are shown to be capable of improving inhibitor efficiency for HIV-I27 and nucleotide cleavage for hammerhead ribozymes.28 β-LNA modification and its stereoisomer α-LNA had caught much attention compared to other LNA modifications.21−23,27,30 These two isomers (α-LNA and β-LNA) differ in the orientation of the −O−CH2− group (Figure 1) but share many features: both the modifications in nucleic acid strands result in strong binding to complementary strands, show high potency as antisense oligonucleotides as found in both in vitro and in vivo experiments, and have commercial availability and ease of synthesis by use of conventional phosphoroamidite chemistry. Given the wide range of applications, it is important to understand the structure, dynamics, and stability of LNAmodified duplexes with respect to pure nucleic acids. Several attempts have been made to synthesize and characterize LNAmodified nucleic acids. NMR studies performed on βLNA:RNA and β-LNA:DNA duplexes had suggested that these duplexes display helical features similar to native nucleic acids (like Watson−Crick base pairing, base stacking, helical twist, anti orientation etc.).31,32 The NMR studies performed on 9-mer DNA−RNA duplexes containing one, two, and three LNA modifications showed increasing A-type character with respect to the number of modifications.33−36 Previous 2D NMR study showed that the change in the conformation of the backbone is local especially if the number of β-LNA modifications is less, whereas the global structure resembles DNA−RNA hybrid. Furthermore, MD simulation studies by Ivanovo et al. on LNA−DNA hybrid have shown that the



COMPUTATIONAL METHODS Model Systems. To investigate the variations in the structure and dynamics of the DNA−RNA hybrid duplex, the sugars corresponding to one, two, three, and full strand deoxyribose nucleotides were modified with α-/β-LNA moieties (see Figure 1). Three different systems were considered for mono (LNA11, LNA12, LNA13) and di (LNA21, LNA22, LNA23) modifications according to the position and the number of nucleotides between them, respectively. The modifications were placed on the sugars of thymine present close to the 5′ terminal (LNA11), at the center (LNA12), and close to the 3′ terminal (LNA13) of the DNA strand in single-modified duplexes. The double-modified duplexes LNA21, LNA23, and LNA22 have one, two, and three nucleotide separation between two LNA modifications, respectively. The model systems were conceived in such a way

5854

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

tions were extended up to 20 ns in NPT ensemble by applying a weak harmonic restraint using a force constant 2 kcal/mol/Å2 on the hydrogen bonds of the terminal base pairs in order to prevent their chances of opening.40,41 The Leapfrog integrator was used to integrate Newton’s equations of motion that used a 2 fs integration time step. The coordinates were saved for every 5 ps for further analysis. The translational and rotational motions of the solute were removed by aligning all the frames of solute on its initial conformation. Final 10 ns simulation trajectories were used for all the structural and energetic analysis treating the first 10 ns simulation time as equilibration period based on the calculated root-mean-square deviations for duplexes with respect to the average structure obtained from final 10 ns (Figure S1 in the Supporting Information). VMD60 and Grace61 were used for image generation and making plots, respectively. The analysis corresponding to helical, base step parameters and groove widths were done with the help of Curves+ program.62

to understand the change in duplex properties with respect to the number of modifications, the position of modification and the gap between two modifications, which resulted in 16 model systems. The available DNA−RNA hybrid NMR experimental structure (1HG9)29 was used as template to generate the model systems by introducing appropriate modifications corresponding to α-/β-LNA nucleotides using CHARMM biomolecular simulation program.47 Available NMR structures were used as starting conformations for some of the LNA incorporated systems wherever available (Table 1).20,29,31,32 Control simulation on the DNA−RNA hybrid NMR experimental structure was also run using the same simulation protocol. Table 1. Protein Data Bank Ids of the Available α-/β-LNAModified DNA−RNA Hybrid Duplexes name

α-LNA

β-LNA

hybrid LNA11 LNA12 LNA13 LNA21 LNA22 LNA23 LNA3 FLNA

1HG929 − − − − − − 1OKF20 −

1HG929 − 1HHW31 − − − − 1HHX31 1H0Q32



RESULTS AND DISCUSSION A comprehensive study has been done by performing several structural and energetic calculations to understand how the incorporated isomeric LNA modifications affect the DNA− RNA hybrid duplex properties. This section is mainly divided into four subsections: first subsection discusses the structural variations in hybrid duplex upon the incorporation of one to many LNA modifications in terms of duplex conformation, thermodynamic stability, deformation, and global conformation. The second subsection discusses the similarities/differences in α-LNA:RNA and β-LNA:RNA duplexes in comparison to the pure DNA−RNA hybrid. The final two subsections cover how different the structures are when isomeric LNA modifications are introduced at different positions and with respect to the gap between two LNA modifications in DNA− RNA hybrid duplex. The present MD results have been compared with the available experimental results whenever possible. Structural Deviations and Dynamics. Root mean square (rms) deviation values calculated for all the duplexes studied here with respect to the average structure obtained from the final 10 ns simulations and averaged for the same simulation time are presented in Table 2. The rmsd values have been calculated by considering only the non-hydrogen atoms; the extent of deviation was also calculated for the atoms of the backbone and nucleobases individually. All the structures show deviation values within the range of 0.89−2.14 Å, indicating closeness to their initial conformation. The α-LNA-modified duplexes show high deviation values than β-LNA-modified, suggesting significant difference in the duplexes compared to regular DNA−RNA hybrid. The evolutions of rms deviations with simulation time are given in Figure S1a,b in the Supporting Information. Rms fluctuation (rmsf) calculations were done to understand the change in the DNA−RNA hybrid duplex flexibility with the incorporation of α-/β-LNA modifications. The rmsf values in Table S1 in the Supporting Information indicate that the α-LNA-modified duplexes show different flexibility pattern compared to the β-LNA-modified duplexes. The α-LNA-modified duplexes show significantly high rmsf values than β-LNA-modified duplexes. The observed changes are mainly due to the high fluctuations of the DNA− RNA hybrid duplex backbone. Noticeably, α-LNA modification increases the duplex flexibility whereas the β-LNA modification decreases the flexibility of the DNA−RNA hybrid duplex,

MD Simulation Protocol. All the MD simulations were performed using the NAMD48 program employing the CHARMM36 all-atom nucleic acid force field.49−51 The topology and parameters for LNA modification were taken from the study of Pande and co-workers,40 developed by considering the point charges based on the RESP charges taken from the previous study of Nielsen et al.32 Previous studies that used these parameters on fully modified LNA duplexes obtained good agreement with the experimental studies.40,41 Missing hydrogen atoms were added to nucleic acids by using HBUILD command in CHARMM program, followed by 500step steepest decent minimization (SD) by applying massweighted harmonic constraints on heavy atoms. The minimized structures were then immersed in a pre-equilibrated water box with edge lengths 48 Å × 38 Å × 38 Å, and sodium ions were added randomly for neutralization. All water molecules making nonphysical contacts with the solute atoms were deleted based on a distance cutoff of 2 Å. A modified TIP3P52 model was used for the explicit solvent environment. These systems were subjected to a 500-step SD minimization followed by a 500-step adopted basis Newton−Raphson (ABNR) minimization. All the non-hydrogen atoms of the nucleic acid were restrained using harmonic potential with a force constant of 5 kcal/mol/ Å2. These systems were equilibrated for 100 ps under the same harmonic restraints on the nucleic acids in NVT ensemble. The covalent bonds involving hydrogens were constrained by applying SHAKE algorithm53 in all the simulations and calculations. CRYSTAL module54 was employed for setting up periodic boundary conditions. Temperature of the systems was maintained at 300 K using the Nose−Hoover thermostat,55 and constant pressure was maintained using the Langevin piston algorithm.56 The electrostatic real space and the Lennard-Jones (LJ) interactions were truncated at 12 Å with a LJ force switch smoothing function from 10 to12 Å.57 Particle mesh Ewald (PME) summation58,59 method was used to treat long-range electrostatic interactions. The production simula5855

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

The hydrogen bond distances corresponding to all the base pairs present in all the duplexes were calculated. The probability distributions corresponding to the distances of central hydrogen bond (N1−N3) (Figure S2a,b in the Supporting Information) show a sharp peak at 3 Å, indicating stability of the base pairs with the maintenance of Watson− Crick (WC) base pairing. But partial base pair opening is observed for certain base pairs in some duplexes indicated by the peaks corresponding to distances longer than 4.5 Å. Backbone and Sugar Conformations. The sugar puckering angles, characteristic backbone dihedral angles α, β, γ, δ, ε, ζ, and sugar−base dihedral angle (χ) are capable of differentiating the nucleic acid conformations. To understand the backbone conformation and dynamics in modified duplexes, an additional parameter ε − ζ is reported which is generally used for characterizing BI and BII conformations (Figure 2). The backbone conformation of the DNA is BI when ε − ζ < 0, and BII for ε − ζ > 0.63,64 The previous MD simulations on full LNA-modified duplexes suggest that the α-LNA-modified duplex has backbone conformation that is different from the typical DNA, RNA, and β-LNA-modified duplex.41 The probabilities corresponding to ε − ζ were computed by using the ε and ζ dihedral angles of all the nucleotides present in the duplex (Figure 2). The probability distributions in Figure 2 indicate a systematic increase in probability in the sampling region around 60−120° (BII) in α-LNA-modified duplexes whereas such kind of increase is not seen in β-LNA-modified duplexes. This suggests that the incorporation of β-LNA modification does not require much rearrangement in backbone in order to accommodate its sugar moiety in the context

Table 2. Average Values of Root Mean Square Deviations for Non-hydrogen Atoms of Overall, Backbone (Sugar and Phosphate) and Bases (Nucleobases) Corresponding to All the Duplexes Studied with Respect to the Averaged Structure Obtained from the Final 10 nsa duplex hybrid LNA11 LNA12 LNA13 LNA21 LNA22 LNA23 LNA3 FLNA a

overall 1.15 1.22 1.51 1.16 1.85 1.14 1.35 1.20 1.45 1.24 1.43 1.11 1.40 1.29 1.43 1.34 1.73

± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±

0.01 0.05 0.04 0.02 0.18 0.02 0.08 0.04 0.06 0.02 0.19 0.06 0.09 0.04 0.03 0.01 0.11

backbone 1.26 1.33 1.69 1.29 2.14 1.27 1.50 1.35 1.71 1.38 1.51 1.24 1.48 1.45 1.55 1.44 1.84

± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±

0.02 0.06 0.06 0.03 0.23 0.02 0.07 0.05 0.08 0.03 0.14 0.07 0.10 0.04 0.04 0.02 0.11

bases 0.95 1.02 1.07 0.92 1.28 0.89 1.06 0.93 1.00 0.97 1.21 1.38 1.19 1.01 1.20 1.11 1.45

± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±

0.01 0.03 0.02 0.01 0.08 0.01 0.08 0.03 0.02 0.02 0.26 0.03 0.07 0.03 0.02 0.01 0.10

Normal letters, β-LNA, italics, α-LNA).

agreeing with the previous MD simulation results on β-LNAmodified duplexes.40 Some additional discussions on the stiffness/flexibility of the DNA−RNA hybrid duplex will be discussed in the later sections.

Figure 2. Probability distributions of dihedral angles (ε − ζ) for all the α-LNA (left) and β-LNA (right) modified DNA−RNA hybrid duplexes. The difference (ε − ζ) corresponding to pure DNA−RNA hybrid duplex (black) and observed in experimental structures (green) is also included in all the sections for reference. 5856

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

Figure 3. (a) Pictorial representation of pseudorotation angles sampled by α-/β-LNA-modified, DNA, and RNA duplexes. (b) Percentage of 2E + 3E (green) and 3E (blue) conformations of DNA strand of the DNA−RNA hybrid duplex when α-/β-LNA modifications are incorporated. The 2E + 3E (black) and 3E (red) conformations of pure DNA−RNA hybrid duplex are also shown in both the panels for comparison.

slight transformation from high anti to anti regions as the number of LNA modifications (for both isomers) increases. There were no appreciable changes observed in backbone dihedral angles as the increase in number of β-LNA modifications whereas the backbone sample additional regions when incorporating α-LNA modifications and their probability increases with the increase in number of modifications. Generally, ribose sugars exhibit 3E conformation while deoxyribose sugar exhibits 2E conformation in addition to minor sampling in 3E and 0E conformations (Figure 3). Previous studies have shown that the ribose sugars present in DNA−RNA hybrid duplex exhibit conformations similar to the conformations exhibited by the ribose sugars of pure RNA duplex but a dynamic equilibrium exists between 2E and 3E

of WC base pairing and better stacking. However, the incorporated α-LNA moiety rearranges the nucleic acid backbone in order to adjust its sugar moiety into the WC pairing. This rearrangement sometimes causes destabilization in the DNA−RNA hybrid duplex through opening of the same and neighbor base pairs. Full strand modified duplexes are stable for both the modifications irrespective of the change in backbone conformation. To further understand the impact of incorporated isomeric LNA analogues on the DNA−RNA hybrid duplex conformation, the probability distributions corresponding to all the backbone dihedral angles are presented in Figure S3a,b in the Supporting Information. The distributions of glycosidic bond dihedral angle (χ angle) indicate that all the bases are in anti conformation despite a 5857

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

of relation is not observed for α-LNA-modified duplexes. Some of the α-LNA-modified duplexes (LNA11 and LNA12) exhibit wider minor groove width values than pure DNA−RNA hybrid duplex. The internal dynamics of the nucleic acid duplexes depends on the stacking of nitrogenous bases present in the same strand and complementary strand and these are crucial for the duplex stability. This base stacking can be explained by using helicoidal and step parameters. Helicoidal parameters represent the orientation of base pair relative to its global helical axis, and step parameters give the orientation and position of base pair with respect to next base pair. The average values of the helical parameters shown in Figures S4 and S5 in the Supporting Information indicate that the β-LNA modification changes the duplex structure in a systematic manner whereas α-LNA modification shows drastic changes in the duplex structure. All the helical parameters, especially shift, roll, and twist, are more sensitive to the structural modifications similar to the previous MD studies on modified nucleic acid duplexes.37,69 The β-LNA modification imparts a systematic variation in most of the parameters whereas drastic changes were observed for α-LNA modifications. For example, the twist of the helix decreases as the increase in number of α-LNA modifications compared to βLNA modifications. Helical unwinding in the DNA−RNA hybrid duplex is observed for partial α-LNA modifications. The duplexes with full α-/β-LNA-modified strands exhibit average values that are distinct from partially modified duplexes (Figure S5). The average values corresponding to individual base pair step parameters indicate that partial modifications corresponding to these α/β-modifications perturb the duplex locally (Figures S6 and S7 in the Supporting Information). But the effect is significantly higher for α-LNA modification resulting in change in the global structure of the DNA−RNA hybrid duplex. The variations observed in base pair step parameters indicate that the LNA modifications show different impact on the duplex internal dynamics. These dynamics are independent of the number, position, and gap between two modifications. But the duplexes with full LNA strand exhibit properties that differ from partly modified duplexes and the observed changes adequately reflect this. Some of the internal rotations (twist) and translational (shift) motions are considerably dependent on the number of modifications and support the fact that the partly modified duplexes have intermediate properties than pure and fully modified duplexes. Though significant variations are observed in some of the helical parameters of partly β-LNAmodified duplexes, the overall effect is still the same, implying that the changes due to partial β-LNA modification are local, supporting the experimental results.34−36 Change in Thermodynamic Stability. The stability of helix depends primarily on the interactions among WC base pairs and stacking interactions with respect to other bases present in the same strand and its complementary strand. Binding free energy calculations were performed using the MM-GBSA method to understand the change in thermodynamic stability of DNA−RNA hybrid when incorporating various numbers of isomeric LNA modifications. In general, the total free energy can be determined using the following equation: G = EMM + Gsolv − TSMM, where G represents total free energy, EMM represents the total molecular mechanical energy, Gsolv gives contribution from solvent, T is absolute temperature, and SMM represents the molecular mechanical entropy. The binding free energy between two complementary strands can be calculated as ΔGbinding = Ecoul(DNA−RNA) +

conformations for the deoxyribose sugars, resulting in an intermediate A/B-like conformation to the DNA−RNA hybrid duplex.45,65 Experimental studies have suggested that the α- and β-LNA analogues act as DNA and RNA mimics, respectively.20 The percentages of puckering angles corresponding to DNA strand were calculated for all the duplexes in order to understand the variation in conformation of the pure and modified DNA−RNA hybrid duplex (Figure 3). The changes observed in percentages in Figure 3 indicate that the systematic incorporation of α- and β-LNA modifications progressively turn the DNA strand toward 3E and 3E conformations, respectively. But the impact of α-LNA modification is high compared to its analogue β-LNA. Percentages corresponding to the individual sugar conformations in all the duplexes are provided in the Supporting Information (Table S2). The sugar puckering amplitudes, the maximum torsion angle with respect to other atoms in the plane, were calculated to check the differences in sugar puckering due to LNA modifications (data not shown). It is noticed that the incorporated LNA increases the sugar moiety amplitude with the increase in sugar amplitude of the complementary nucleotide. However, these two modifications act differently on the neighbor sugar moieties; i.e., the α-LNA modification increases their amplitudes and the β-LNA modification decreases it. This could be due to the orientation of the added −O−CH2− group. The observed changes in conformation and amplitudes of the sugars are due to incorporation of LNA modifications and they indicate instability of the DNA−RNA hybrid duplex in terms of the conformation. Minor Groove Widths and Helical Parameters. The previous studies have suggested that the intermediate minor groove width of DNA−RNA hybrids act as one of the key factors in the RNase H activity.66 Average values of minor groove width regions present in all the pure and modified DNA−RNA hybrid duplexes were calculated and are shown in Figure 4. These values indicate that the β-LNA modification

Figure 4. Average values of minor groove width regions for β-LNA (red) and α-LNA (black) modified duplexes calculated for the final 10 ns simulation time.

does not show appreciable change in the minor groove width whereas the α-LNA modification narrows the minor groove width. The extent of impact on the minor groove width is high for α-LNA modification than β-LNA modification. The minor groove width regions of the partial β-LNA-modified DNA− RNA hybrid duplexes are very close to the minor groove width values observed for the full β-LNA duplex41 whereas such kind 5858

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

Evdw(DNA−RNA) + Gsolv(DNA−RNA) − Gsolv(DNA) − Gsolv(RNA); here the double strand represents duplex and DNA, RNA represents individual strands. The solvation free energy was estimated as the sum of the electrostatic solvation energy (GGB) and the nonpolar solvation energy (Gnonpolar). Gsolv = GGB + Gnonpolar; the electrostatic contribution to the solvation free energy was estimated by using the generalized Born molecular volume (GBMV) method. The nonpolar contribution toward solvation free energy was estimated as Gnonpolar = γ × SASA, where γ = 0.0072 kcal/mol/Å2 and SASA is the solvent-accessible surface area,67,68 calculated by using a probe particle with a radius of 1.4 Å on DNA−RNA surface. The contribution from entropy is not considered here and such approach is widely used in the literature and has been found to be adequate for obtaining qualitative trends.43,69−71 Within the approximations used here, the entropies of the individual strands are similar when they are part of the duplex and when they are free/considered individually. Hence, including the entropic factors does not change the qualitative trends (see Figure S8 in the Supporting Information). Relative binding free energies (ΔΔG) with respect to the DNA−RNA hybrid were calculated (Figure 5).

to understand the helical stiffness of the modified DNA−RNA hybrid duplexes. Diagnolization of positional covariance matrices provide eigenvectors and eigenvalues where eigenvectors describe the nature of the essential movements and eigenvalues give the contribution of each movement to the positional variance of the trajectory. The force constants calculated for essential modes (in this case, first 15 modes) represent the rigidity of the duplex. It has been shown before that the 10 lowest modes account for most of the fluctuations of nucleic acid duplexes.45,65,72,73 The force constants, within the harmonic limit, can be calculated from the eigenvalues with the following relation ki =

kBT λi

where λi represents the eigenvalue (in Å2) corresponding to mode i, kB is the Boltzmann constant, and T is the absolute temperature. The calculated force constants indicate similar plasticity of the duplexes at small deformation modes (Figure 6). But the variations in force constants become different at higher modes for α-/β-LNA-modified duplexes, indicating their different duplex nature in terms of the helix stiffness. It is also noticed that the duplexes with fully modified α-/β-LNA strand are softer than partially modified and pure DNA−RNA hybrid, agreeing with the experimental results. As the number of α-/βLNA modifications increases in the DNA−RNA hybrid, the duplex becomes softer and less stiff than the regular DNA− RNA hybrid. The duplex with three LNA modifications is less stiff than the full strand modified duplex, further supporting the experimental results and validating the MD simulations. The intramolecular entropies were calculated employing the Schlitter method by using the eigenvalues obtained from the mass-weighted covariance matrix.45 These plots indicate that the DNA−RNA hybrid duplex is the most disordered structure at smaller essential modes. But this situation changes at higher modes for both α-/β-LNA modifications (Figure 6). Even though the stiffness analysis using essential motions provides useful information regarding deformability, the modes of such motions are not easily understandable.72,74 The elastic force constants corresponding to the helix deformations along their translational and rotational motions in helical space (shift, rise, slide, roll, tilt, and twist) were calculated to examine the nature of fluctuations. The inversion of the helical covariance matrix (C = ⟨XiXj⟩) of helical parameters (Xi,Xj) gives the stiffness matrix (F = F[kij]), and the diagonal elements (kij, i = j) give force constants corresponding to pure steps and offdiagonal elements (ki,j, i ≠ j) representing the force constant for cross terms.45,72

Figure 5. Relative binding free energies (in kcal/mol) for DNA−RNA hybrids with increasing incorporated LNA modifications.

These ΔΔG values indicate that the full strand modified DNA−RNA hybrid duplexes are more stable than pure hybrid duplex, in agreement with the experimental results.32 Some of the partial β-LNA-modified duplexes are also more stable than pure DNA−RNA hybrid, agreeing with the previous studies. However, the partial incorporation of α-LNA modifications results in duplexes that are unstable compared to pure DNA− RNA hybrids. Further, it is observed that the enhanced duplex stability via incorporating α-/β-LNA modifications is dependent on the position and gap (number of nucleotides) between two successive LNA modifications. It is also noticed that the increased LNA modifications in DNA strand show strong affinity toward complementary RNA strand and favorable LNA−RNA duplex formation for β-LNA modification whereas the reverse is observed for partial α-LNA modification. Helical Stiffness/Deformability and Effect of LNA Modifications. Previous experimental studies showed that the LNA−RNA hybrid is more rigid than pure duplex and the fully modified LNA strand makes the duplex more rigid than partially modified LNA strand.18 Force constants and intramolecular entropies in harmonic limit (along essential modes) were computed by using principal component analysis (PCA)

F = kBTC −1

where F represents the stiffness matrix, kB is the Boltzmann constant, T is absolute temperature, and C represents the covariance matrix. To further explain the stiffness at the global translational (ktrans), rotational (krotat), and global (kglobal) levels, the respective force constants are calculated by using the following equations. k trans = kshiftk risekslide

k rotat = k rollk twistk tilt kglobal = k transk rotat 5859

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

Figure 6. Harmonic force constants (in kcal/mol/Å2) and intramolecular entropies (in kcal/mol/K) corresponding to the first 15 modes for all the duplexes studied.

unwinding than the β-LNA:RNA duplex and different thermodynamic stability. The extent of structural and energetic perturbation for α-LNA modification is larger compared to βLNA modification observed in base pair step parameters (Figures S12 and S13 in the Supporting Information). The binding of α-/β-LNA strands complementary to RNA strand results in a stable LNA:RNA duplex but the extent of stabilization is higher for β-LNA modification, suggesting a better selection of β-LNA modification over α-LNA modification. A dissimilar backbone rearrangement is observed for α-/β-LNA-modified DNA strands which further support their distinct nature. Do the Hybrid Properties Depend on the Position of Modification and the Gap between Two Modifications? Previous studies on LNA-modified duplexes suggest that the changes in duplexes are dependent on the position (5′ or 3′ or center) of modification.33−35 In this section, the variations observed in hybrid properties when introducing the LNA analogues at different positions, i.e., close to 5′, center, and close to the 3′ end, are discussed. The single-modified systems LNA11, LNA12, and LNA13 are discussed in this section. The LNA11 and LNA13 have modifications near the terminal ends, and LNA12 has modification at middle position (see Figure 1). The observed variations in structural and energetic features indicate that the dissimilar nature of α-/β-LNA modifications toward affecting the DNA−RNA hybrid duplex properties. The modification at central position results in a duplex different from the one with terminal modifications for both α-/β-LNA modifications. Qualitative trends of the relative binding free energies suggest that α-/β-LNA modification at the central position stabilizes the duplex compared to the terminal modifications. In other words, the influence of LNA modification on the DNA−RNA hybrid duplex is position dependent. Similar differences can also be found in other duplexes like pure DNA and pure RNA duplexes. This subsection covers the impact of isomeric LNA modifications when modifying the deoxyribose sugars at two positions while maintaining different gaps between them. The

It is noticed that the calculated force constants are higher for βLNA-modified duplexes and smaller for α-LNA-modified duplexes than DNA−RNA hybrid duplex (Table S3 in the Supporting Information). This indicates that the β-LNAmodified duplexes are more rigid than the DNA−RNA hybrid and α-LNA-modified duplexes in all the local and global helical deformations. In other words, the deformation of the DNA− RNA hybrid duplex along its translational and rotational motions will become moderately easy after introducing α-LNA modifications, and it will be slightly stiff when introducing βLNA modifications. This indicates the opposite nature of the isomeric α-/β-LNA modifications on the stiffness of the DNA− RNA hybrid duplex, in support of the harmonic force constants and intramolecular entropies. Comparison of 3α-LNA−DNA:RNA versus 3β-LNA− DNA:RNA and α-LNA:RNA versus β-LNA:RNA duplexes. The previous studies suggest that more than three modifications are sufficient to change the duplex conformation into Atype. To understand this, the three α/β-modified duplexes are compared here. Appreciable changes are not found from the calculated intra-phosphate−phosphate (P−P) distances (Figure S9 in the Supporting Information). The partial α-LNAmodified DNA−RNA hybrid duplexes show higher helical unwinding than partial β-LNA-modified nucleic acid duplexes as observed in full LNA-modified duplexes41 (Figures S10 and S11 in the Supporting Information). The present MD simulation resulted in a 3 α-LNA-modified duplex (LNA3) structure that exhibits a larger helical unwinding compared to DNA−RNA hybrid and 3 β-LNA-modified duplex. The variations in global and local base pair step parameters suggest an altered geometrical nature for both the LNA3 duplexes. The fully modified DNA−RNA hybrid duplexes indicate that the full α-LNA modification slightly increases the intrastrand P−P distance of the modified strand and decreases the intrastrand P−P distances of the complementary strand (Figure S9). But the β-modified strand has small P−P distances compared to pure DNA−RNA hybrid and full α-LNA-modified duplex. The α-LNA:RNA duplex exhibits higher helical 5860

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

of hydrogen bond lengths in the base pairs, probability distributions of backbone dihedral angles, average values of base pair step and helical parameters, time series of helical parameters, and intrastrand phosphate−phosphate distances. This material is available free of charge via the Internet at http://pubs.acs.org.

systems LNA23, LNA21, and LNA22 have 1, 2, and 3 intermediate nucleotides (gap) between two successive LNA modifications. The variations observed in the properties of these systems indicate that the effect of each modification is similar to the respective single modification, suggesting an additive nature of the modifications. The various energy and structural calculations suggest that these two isomeric α-/βLNA modifications affect the duplex peculiarly when the gap is increased between them. The β-LNA modifications stabilize the duplex with large gap between them whereas the α-LNA modification destabilizes it. It may be true that the stability of the duplex increases as the increase in gap between two modifications increases for β-LNA modification but not for αLNA modification. Future studies with various number of sequence lengths are required to confirm this. The various trends observed in sugar conformations and helicoidal parameters suggest that the gap between two modifications play a significant role in enhancing the duplex properties.



Corresponding Author

*E-mail: [email protected]. Phone: +91-40-6653 1161. Fax: +9140-6653 1413. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS G.S. thanks Council of Scientific and Industrial Research (CSIR), India, for a senior research fellowship. U.D.P. thanks the Department of Biotechnology (Government of India) for the IYBA award. Department of Atomic Energy (Govt. of India) is acknowledged for financial assistance.



CONCLUSIONS MD simulations have been performed to understand the change in structures, dynamics, and stability of DNA−RNA hybrid duplexes when systematically modifying its deoxyribose sugars with isomeric LNA (α-/β-LNA) modifications. The changes in structure and energetics due to these modifications were examined using various analyses, including principal component (PCA) and free energy calculations. The observed trends suggest that these two isomeric modifications (α-/βLNA) significantly affect the stability/stiffness of the DNA− RNA hybrid duplex. The binding free energies indicate that the full strand modified α-LNA:RNA and β-LNA:RNA duplexes are quite stable compared to pure DNA−RNA hybrid duplex. Some of the partial modified duplexes are also stable. The structural perturbation in the DNA−RNA hybrid duplex is high for α-LNA modification compared to β-LNA. For example, there is an unusual disparity in helical twist for α-LNA-modified duplexes compared to β-LNA-modified duplexes. Combination of the present results and previous MD results on full α-/βLNA-modified duplexes41 suggests that the β-LNA modification can mimic ribose sugar, but the α-LNA does not mimic the deoxyribose sugar. The present MD results also suggest that 3 β-LNA modifications were not sufficient to change the DNA− RNA duplex conformation into A-like. But the increasing βLNA modifications slowly turn the DNA−RNA duplex into Alike but the α-LNA modification changes the conformation instead of preserving its A/B-like conformation. The full strand modified α-LNA:RNA and β-LNA:RNA duplexes show properties distinct from pure DNA−RNA hybrid duplex. The current MD simulation results also suggest that the perturbations in DNA−RNA hybrid duplex properties are dependent on the position of modification and the gap between two modifications. The population distribution of backbone dihedral angles BI versus BII suggests that the β-LNA modification can be incorporated into the DNA strand of the DNA−RNA hybrid duplex without backbone rearrangement whereas backbone rearrangement is required in order to accommodate the α-LNA modification, resulting in unstable duplexes for some of the partial modifications.



AUTHOR INFORMATION



REFERENCES

(1) Stephenson, M. L.; Zamecnik, P. C. Inhibition of rous-sarcoma viral-RNA translation by a specific oligodeoxyribonucleotide. Proc. Natl. Acad. Sci. U.S.A. 1978, 75, 285−288. (2) Zamecnik, P. C.; Stephenson, M. L. Inhibition of rous-sarcoma virus-replication and cell transformation by a specific oligodeoxynucleotide. Proc. Natl. Acad. Sci. U.S.A. 1978, 75, 280−284. (3) Crooke, S. T. Molecular mechanisms of action of antisense drugs. Biochim. Biophys. Acta 1999, 1489, 31−44. (4) Kurreck, J. Antisense technologies: improvement through novel chemical modifications. Eur. J. Biochem. 2003, 270, 1628−1644. (5) Helene, C.; Toulme, J. J. Specific regulation of gene expression by antisense, sense and antigene nucleic acids. Biochim. Biophys. Acta 1990, 1049, 99−125. (6) Dias, N.; Stein, C. A. Antisense Oligonucleotides: Basic Concepts and Mechanisms. Mol. Cancer Ther. 2002, 1, 347−355. (7) Bennett, C. F.; Swayze, E. E. RNA targeting therapeutics: molecular mechanisms of antisense oligonucleotides as a therapeutic platform. Annu. Rev. Pharmacol. Toxicol. 2010, 50, 259−293. (8) Uhlmann, E.; Peyman, A. Antisense oligonucleotides: a new therapeutic principle. Chem. Rev. 1990, 90, 543−584. (9) Freier, S. M.; Altmann, K. H. The ups and downs of nucleic acid duplex stability: structure stability studies on chemically modified DNA:RNA duplexes. Nucleic Acids Res. 1997, 25, 4429−4443. (10) Manoharan, M. 2′-Carbohydrate modifications in antisense oligonucleotide therapy: importance of conformation, configuration and conjugation. Biochim. Biophys. Acta 1999, 1489, 117−130. (11) Braasch, D. A.; Corey, D. R. Locked nucleic acid (LNA):finetuning the recognition of DNA and RNA. Chem. Biol. 2001, 8, 1−7. (12) Petersen, M.; Wengel, J. LNA: a versatile tool for therapeutics and genomics. Trends Biotechnol. 2003, 21, 74−81. (13) Imanishi, T.; Obika, S. BNAs: novel nucleic acid analogs with a bridged sugar moiety. Chem. Commun. 2002, 1653−1659. (14) Nielsen, P. E. Peptide nucleic acids (PNA) in chemical biology and drug discovery. Chem. Biodivers. 2010, 7, 786−804. (15) Renneberg, D.; Schumperli, D.; Leumann, C. J. Biological and antisense properties of tricyclo-DNA. Nucleic Acids Res. 2002, 30, 2751−2757. (16) Moulton, J. D.; Jiang, S. Gene knockdowns in adult animals: PPMOs and vivo-morpholinos. Molecules 2009, 14, 1304−1323. (17) Kang, H.; Fisher, M. H.; Xu, D.; Miyamoto, Y. J.; Marchand, A.; Van Aerschot, A.; Herdewijn, P.; Juliano, R. L. Inhibition of MDR1 gene expression by chimeric HNA antisense oligonucleotides. Nucleic Acids Res. 2004, 32, 4411−4419.

ASSOCIATED CONTENT

S Supporting Information *

Tables of rmsf and elastic force constants, and figures of time series of root-mean-square distances, probability distributions 5861

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

(18) Kawasaki, A. M.; Casper, M. D.; Freier, S. M.; Lesnik, E. A.; Zounes, M. C.; Cummins, L. L.; Gonzalez, C.; Cook, P. D. Uniformly modified 2′-deoxy-2′-fluoro phosphorothioate oligonucleotides as nuclease-resistant antisense compounds with high affinity and specificity for RNA targets. J. Med. Chem. 1993, 36, 831−841. (19) Kauppinen, S.; Vester, B.; Wengel, J. Locked nucleic acid (LNA): High affinity targeting of RNA for diagnostics and therapeutics. Drug Discovery Today: Technol. 2005, 2, 287−290. (20) Nielsen, J. T.; Stein, P. C.; Petersen, M. NMR structure of an αL-LNA:RNA hybrid: structural implications for RNase H recognition. Nucleic Acids Res. 2003, 31, 5858−5867. (21) Jepsen, J. S.; Sorensen, M. D.; Wengel, J. Locked Nucleic Acid: A Potent Nucleic Acid Analog in Therapeutics and Biotechnology. Oligonucleotides 2004, 14, 130−146. (22) Vester, B.; Lundberg, L. B.; Sorensen, M. D.; Babu, B. R.; Douthwaite, S.; Wengel, J. LNAzymes: incorporation of LNA-type monomers into DNAzymes markedly increases RNA cleavage. J. Am. Chem. Soc. 2002, 124, 13682−13683. (23) Elmen, J.; Thonberg, H.; Ljungberg, K.; Frieden, M.; Westergaard, M.; Xu, Y. H.; Wahren, B.; Liang, Z. C.; Urum, H.; Koch, T.; Wahlestedt, C. Locked nucleic acid (LNA) mediated improvements in siRNA stability and functionality. Nucleic Acids Res. 2005, 33, 439−447. (24) Wang, L.; Yang, C. Y. J.; Medley, C. D.; Benner, S. A.; Tan, W. H. Locked nucleic acid molecular beacons. J. Am. Chem. Soc. 2005, 127, 15664−15665. (25) Thomsen, R.; Nielsen, P. S.; Jensen, T. H. Dramatically improved RNA in situ hybridization signals using LNA-modified probes. RNA 2005, 11, 1745−1748. (26) Kloosterman, W. P.; Wienholds, E.; de Bruijn, E.; Kauppinen, S.; Plasterk, R. H. A. In situ detection of miRNAs in animal embryos using LNA-modified oligonucteotide probes. Nat. Methods 2006, 3, 27−29. (27) Jakobsen, M. R.; Haasnoot, J.; Wengel, J.; Berkhout, B.; Kjems, J. Efficient inhibition of HIV-1 expression by LNA-modified antisense oligonucleotides and DNAzymes targeted to functionally selected binding sites. Retrovirology 2007, 4, 1−13. (28) Christiansen, J. K.; Lobedanz, S.; Arar, K.; Wengel, J.; Vester, B. LNA nucleotides improve cleavage efficiency of singular and binary hammerhead ribozymes. Bioorg. Med. Chem. 2007, 15, 6135−6143. (29) Bondensgaard, K.; Petersen, M.; Singh, S. K.; Rajwanshi, V. K.; Kumar, R.; Wengel, J.; Jacobsen, J. P. Structural studies of LNA: RNA duplexes by NMR: conformations and implications for RNase H activity. Chemistry 2000, 6, 2687−2695. (30) Sorensen, M. D.; Kvaerno, L.; Bryld, T.; Hakansson, A. E.; Verbeure, B.; Gaubert, G.; Herdewijn, P.; Wengel, J. α-L-riboConfigured locked nucleic acid (α-L-LNA): synthesis and properties. J. Am. Chem. Soc. 2002, 124, 2164−2176. (31) Petersen, M.; Bondensgaard, K.; Wengel, J.; Jacobsen, J. P. Locked nucleic acid (LNA) recognition of RNA: NMR solution structures of LNA:RNA hybrids. J. Am. Chem. Soc. 2000, 124, 5974− 5982. (32) Nielsen, K. E.; Rasmussen, J.; Kumar, R.; Wengel, J.; Jacobsen, J. P.; Petersen, M. NMR studies of fully modified locked nucleic acid (LNA) hybrids: solution structure of an LNA: RNA hybrid and characterization of an LNA:DNA hybrid. Bioconjugate Chem. 2004, 15, 449−457. (33) Nielsen, K. E.; Singh, S. K.; Wengel, J.; Jacobsen, J. P. Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC):d(GCAGAAGCAG) duplex, containing four locked nucleotides. Bioconjugate Chem. 2000, 11, 228− 238. (34) McTingue, P. M.; Peterson, R. J.; Kahn, J. D. SequenceDependent Thermodynamic Parameters for Locked Nucleic Acid (LNA)-DNA Duplex Formation. Biochemistry 2004, 43, 5388−5405. (35) Kaur, H.; Arora, A.; Wengel, J.; Maiti, S. Thermodynamic, Counterion, and Hydration Effects for the Incorporation of Locked Nucleic Acid Nucleotides into DNA Duplexes. Biochemistry 2006, 45, 7347−7355.

(36) Kaur, H.; Wengel, J.; Maiti, S. Thermodynamics of DNA-RNA Heteroduplex Formation: Effects of Locked Nucleic Acid Nucleotides Incorporated into the DNA Strand. Biochemistry 2008, 47, 1218− 1227. (37) Ivanova, A.; Rosch, N. The Structure of LNA:DNA Hybrids from Molecular Dynamics Simulations: The Effect of Locked Nucleotides. J. Phys. Chem. A 2007, 111, 9307−9319. (38) Nielson, K. M. E.; Petersen, M.; Hakansson, A. E.; Wengel, J.; Jacobsen, J. P. α-L-LNA (α-L-ribo configured locked nucleic acid) recognition of DNA: an NMR spectroscopic study. Chem.Eur. J. 2002, 8, 3001−3009. (39) Eichert, A.; Behling, K.; Betzel, C.; Erdmann, V. A.; Furste, J. P.; Forster, C. The crystal structure of an ‘All Locked’ nucleic acid duplex. Nucleic Acids Res. 2010, 38, 6729−6736. (40) Pande, V.; Nilsson, L. Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations. Nucleic Acids Res. 2008, 36, 1508− 1516. (41) Suresh, G.; Priyakumar, U. D. Structures, Dynamics and Stabilities of Fully Modified Locked Nucleic Acid (β-D-LNA and α-LLNA) Duplexes in Comparison to Pure DNA and RNA Duplexes. J. Phys. Chem. B 2013, 117, 5556−5564. (42) Nielsen, C. B.; Singh, S. K.; Wengel, J.; Jacobsen, J. P. The Solution Structure of a Locked Nucleic Acid (LNA) Hybridized to DNA. J. Biomol. Struct. Dyn. 1999, 17, 175−191. (43) Petersen, M.; Hakansson, A. E.; Wengel, J.; Jacobsen, J. P. α-LLNA (α-l-ribo Configured Locked Nucleic Acid) Recognition of RNA. A Study by NMR Spectroscopy and Molecular Dynamics Simulations. J. Am. Chem. Soc. 2001, 123, 7431−7432. (44) Condon, D. E.; Yildirim, I.; Kennedy, S. D.; Mort, B. C.; Kierzek, R.; Turner, D. H. Optimization of an AMBER Force Field for the Artificial Nucleic Acid, LNA, and Benchmarking with NMR of L(CAAU). J. Phys. Chem. B 2014, 118, 1216−1228. (45) Noy, A.; Luque, F. J.; Orozco, M. Theoretical Analysis of Antisense Duplexes: Determinants of the RNase H Susceptibility. J. Am. Chem. Soc. 2008, 130, 3486−3496. (46) Nina, M.; Fonne-Pfister, R.; Beaudegnies, R.; Chekatt, H.; Jung, P. M. J.; Murphy-Kessabi, F.; Mesmaeker, A. D.; Wendeborn, S. Recognition of RNA by Amide Modified Backbone Nucleic Acids: Molecular Dynamics Simulations of DNA-RNA Hybrids in Aqueous Solution. J. Am. Chem. Soc. 2005, 127, 6027−6038. (47) Brooks, B. R.; Brooks, C. L., III; Mackerell, A. D., Jr.; Nilsson, L.; Karplus, M. CHARMM: The biomolecular simulation program. J. Comput. Chem. 2009, 30, 1545−1614. (48) Phillips, J. C.; Braun, R.; Wang, Wei; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781− 1802. (49) Foloppe, N.; MacKerell, A. D., Jr. All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data. J. Comput. Chem. 2000, 21, 86−104. (50) MacKerell, A. D., Jr.; Banavali, N. K. All-atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution. J. Comput. Chem. 2000, 21, 105−120. (51) Denning, E. J.; Priyakumar, U. D.; Nilsson, L.; MacKerell, A. D., Jr. Impact of 2′-hydroxyl sampling on the conformational properties of RNA: Update of the CHARMM all-atom additive force field for RNA. J. Comput. Chem. 2011, 32, 1929−1943. (52) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926−936. (53) Ryckaert, J. P.; Ciccotti, G.; Berendsen, H. J. C. Numerical integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes. J. Comput. Phys. 1997, 23, 327−341. (54) Field, M. J.; Karplus, M. CRYSTAL: Program for Crystal Calculations in CHARMM; Harvard University: Cambridge, MA, 1992. 5862

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863

The Journal of Physical Chemistry B

Article

(55) Hoover, W. G. Canonical dynamics: equilibrium phase-space distributions. Phys. Rev., A 1985, 31, 1695−1697. (56) Feller, S. E.; Zhang, Y.; Pastor, R. W.; Brooks, R. W. Constant pressure molecular dynamics simulation: the Langevin piston method. J. Chem. Phys. 1995, 103, 4613−4621. (57) Steinbach, P. J.; Brooks, B. R. New spherical-cutoff methods for long-range forces in macromolecular simulation. J. Comput. Chem. 1994, 15, 667−683. (58) Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T. A.; Lee, H.; Pedersen, L. G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577−8593. (59) Darden, T.; Perera, L.; Li, L. P.; Pedersen, L. New tricks for modelers from the crystallography toolkit: the particle mesh ewald algorithm and its use in nucleic acid simulations. Struct. Fold Des. 1999, 7, R55−R60. (60) Humphrey, W.; Dalke, A.; Schulten, K. VMDVisual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33−38. (61) http://plasma-gate.weizmann.ac.il/Grace/. (62) Lavery, R.; Moakher, M.; Maddocks, J. H.; Petkeviciute, D.; Zakrzewska, K. Conformational analysis of nucleic acids revisited: Curves+. Nucleic Acids Res. 2009, 37, 5917−5929. (63) Hartmann, B.; Piazzola, D.; Lavery, R. BI−BII transitions in BDNA. Nucleic Acids Res. 1993, 21, 561−568. (64) Heddi, B.; Foloppe, N.; Bouchemal, N.; Hantz, E.; Hartmann, B. Quantification of DNA BI/BII backbone states in solution. Implications for DNA overall structure and recognition. J. Am. Chem. Soc. 2006, 128, 9170−9177. (65) Priyakumar, U. D.; MacKerell, A. D., Jr. Atomic detail investigation of structure and dynamics of DNA-RNA hybrids: A Molecular dynamics simulations study. J. Phys. Chem. B 2008, 112, 1515−1524. (66) Nowotny, M.; Gaidamakov, S. A.; Crouch, R. J.; Yang, W. Crystal Structures of RNase H Bound to an RNA/DNA Hybrid: Substrate specificity and Metal Dependent Catalysis. Cell 2005, 121, 1005−1016. (67) Lee, B.; Richards, F. M. The interpretation of protein structures: estimation of static accessibility. J. Mol. Biol. 1971, 55, 3799−400. (68) Still, W. C.; Tempczyk, A.; Hawley, R. C.; Hendrickson, T. Semianalytical treatment of solvation for molecular mechanics and dynamics. J. Am. Chem. Soc. 1990, 112, 10606−10614. (69) Doshi, U.; McGowan, L. C.; Ladani, T. S.; Hamelberg, D. Resolving the complex role of enzyme conformational dynamics in catalytic function. Proc. Natl. Acad. Sci. U.S.A. 2012, 15, 5699−5704. (70) Fulle, S.; Withers-Martinez, C.; Blackman, M. J.; Morris, G. M.; Finn, P. W. Molecular Determinants of Binding to the Plasmodium Subtilisin-like Protease 1. J. Chem. Info. Model. 2013, 53, 573−583. (71) Robaldo, L.; Pointggia, R.; Di Lella, S.; Estrin, E. A.; Engels, J. W.; Iribarren, A. M. Montserrat, J. M. Conformational states of 2′-CMethylpyrimidine Nucleosides in Single and Double Nucleic Acid Stranded Structures. J. Phys. Chem. B 2012, 117, 57−69. (72) Noy, A.; Perez, A.; Lankas, F.; Luque, F. J.; Orozco, M. Relative flexibility of DNA and RNA: a Molecular dynamics simulations study. J. Mol. Biol. 2004, 343, 627−638. (73) Noy, A.; Perez, A.; Marquez, M.; Luque, F. J.; Orozco, M. Structure, recognition properties, and flexibility of the DNA-RNA hybrid. J. Am. Chem. Soc. 2005, 127, 4910−4920. (74) Olson, W. K.; Gorin, A. A.; Lu, X. J.; Hock, L. M.; Zhurkin, V. B. DNA sequence-dependent deformability deduced from protein-DNA crystal complexes. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 11163−11168.

5863

dx.doi.org/10.1021/jp5014779 | J. Phys. Chem. B 2014, 118, 5853−5863