Coarse-grained Modeling of DNA Curvature

Coarse-grained Modeling of DNA Curvature

Gordon S. Freeman Department of Chemical and Biological Engineering
University of Wisconsin–Madison, Madison, WI, 53706 USA
   Daniel M. Hinckley Department of Chemical and Biological Engineering
University of Wisconsin–Madison, Madison, WI, 53706 USA
   Joshua P. Lequieu Institute for Molecular Engineering
University of Chicago, Chicago, IL 60637 USA
   Jonathan K. Whitmer Materials Science Division
Argonne National Laboratory, Argonne, IL 60439 USA
   Juan J. de Pablo Institute for Molecular Engineering
University of Chicago, Chicago, IL 60637 USA
Materials Science Division
Argonne National Laboratory, Argonne, IL 60439 USA
Corresponding Author: depablo@uchicago.edu
July 25, 2019
Abstract

Modeling of DNA–protein interactions is a complex process involving many important time and length scales. This can be facilitated through the use of coarse-grained models which reduce the number of degrees of freedom and allow efficient exploration of binding configurations. It is known that the local structure of DNA can significantly affect its protein-binding properties (i.e. intrinsic curvature in DNA-histone complexes). In a step towards comprehensive DNA–protein modeling, we expand the 3SPN.2 coarse-grained model to include intrinsic shape, and validate the refined model against experimental data including melting temperature, local flexibility, persistence length, and minor groove width profile.

I Introduction

Coarse-grained (CG) models provide access to time and length scales that are not generally accessible to all-atom (AA) molecular simulations. CG models have been applied to a wide variety of systems, including liquid crystals, block copolymers, proteins, and dexoxyribonucleic acid (DNA). Here we confine our discussion to CG models of DNA, which have been used to study various phenomena, including hybridizationHinckley, Lequieu, and de Pablo (2014), stretchingRomano et al. (2013), and bubble formationZeida et al. (2012). Several DNA CG models have been proposed in recent yearsDans et al. (2010); Ouldridge, Louis, and Doye (2011); Edens, Brozik, and Keller (2012); He et al. (2013); Hinckley et al. (2013); to the best of our knowledge, models designed specifically for studies of DNA structure and thermodynamic properties have not been applied to explore DNA–protein interactions.

Proteins do not have direct access to the functional groups that lead to the formation of the Watson–Crick (W–C) base pairs; however, the sequence of contiguous W–C base pairs results in local deviations from the mean double helix that facilitate the binding of a protein to a unique location along the double-stranded DNA (dsDNA). These changes can be quantified by the intrinsic flexibility and geometry of sequential nucleotide pairs (a base-stepOlson et al. (1998)) or quartetsLavery et al. (2010). Variations in the widths of the major and minor grooves, as well as the flexibility of individual base-steps, dictate the energetic benefit of binding. A CG model suitable for modeling DNA interacting with proteins should capture these sequence effects. In addition, the CG model should include explicit electrostatics, as charged amino acid side chains interact with the negatively–charged DNA backbone. Most CG models include sequence–dependence to some degree, usually in the CG topology and the energy parameters of the base stacking and base pairing interactions. Some models have demonstrated consistency with local base-step properties such as rise and twistDans et al. (2010); Edens, Brozik, and Keller (2012). However, no CG model to date has been demonstrated to capture the intrinsic curvature and global flexibility of DNA.

We present 3SPN.2C, a modification of the 3SPN.2 CG DNA modelHinckley et al. (2013). 3SPN.2 was previously shown to correctly model DNA biophysics by penalizing deviations from the ideal B-form DNA (B–DNA) crystallographic structureArnott, Campbell Smith, and Chandrasekaran (1976). The use of ideal B–DNA to define the minimum energy configuration suggests that sequence effects can be included by using sequence–specific base-step parameters to build the configuration. Additional sequence–dependence can be added by making the flexibility of bonded interactions dependent on the sequence context. Here we adopt both of these changes to extend the 3SPN.2 model to include sequence-dependent shape and flexibility.

The manuscript begins with a description of the model and the data and methods used to assign sequence-dependent parameters. Results are then presented to demonstrate consistency with experimental melting temperatures and flexibilities. Lastly, we present a comparison of simulated minor groove widths to available experimental data.

Ii Methods

ii.1 3SPN.2C DNA Model

The 3SPN.2C model represents an extension of 3SPN.2, a third-generation CG modelKnotts et al. (2007); Sambriski, Schwartz, and de Pablo (2009); Hinckley et al. (2013) shown recently to capture the correct structural, mechanical, and thermodynamic properties of DNA. The bonded and non-bonded potentials of 3SPN.2 were constructed to penalize deviations from the B-DNA crystal structureArnott, Campbell Smith, and Chandrasekaran (1976). Force constants of the bond, angle, and dihedral potentials were independent of sequence, while base stacking energies and base pairing energies were sequence–dependent. Special emphasis was placed on capturing the correct flexibilities of both single–stranded DNA (ssDNA) and dsDNA. 3SPN.2C is intended for simulations of dsDNA interacting with proteins. Consequently, emphasis is placed on capturing the sequence-dependent shape and flexibility of dsDNA. The properties of ssDNA are not prioritized, making 3SPN.2 better suited for studies involving ssDNA.

In the 3SPN.2C model, the original version of 3SPN.2 is modified as follows: First, the reference configuration that defines the minimum energy structure of dsDNA is modified to include sequence-dependent shape. The methodology for so doing is described in Section II.2. This results in equilibrium distances and angles that are a function of the base-step111Code for generating initial topologies and simulating 3SPN.2C is available upon request. Second, each base step is assigned unique force constants for each bend angle, as described in Section II.3. Lastly, weak dihedral potentials are assigned to all dihedral angles formed by the three-site-per-nucleotide topology. These dihedrals, with the functional form

(1)

provide additional stability to the helix when deformed severely from the equilibrium structure, as is often the case when DNA binds to protein (i.e. DNA-histone binding). The magnitude of the torsion force constants, which are independent of sequence, is modified to provide qualitative agreement with experimental data from Ref. Geggier and Vologodskii, 2010.

The aforementioned modifications require the bonded and non-bonded energy parameters to be adjusted in order to preserve agreement with experimental data. Metadynamics simulations are used to assign both the intrastrand stacking energies and interstrand base pairing and cross stacking energies, as in the original parameterization of 3SPN.2Hinckley et al. (2013). The force constants and energy parameters for all interactions can be found in the Appendix. All simulations but the melting temperature calculations (cf. Fig. III.1) are performed at 300 K and 150 mM ionic strength in the NVT ensemble using a Langevin thermostat.

ii.2 Imparting Sequence-dependent Shape

Figure 1: Comparison of aligned minimum energy topologies used in 3SPN.2C (opaque) and 3SPN.2 (transparent). The dotted lines represents the helical axis of the topology without sequence-dependence. Sequence: 5’-CTGGAGAATCCCGGTGCCGAGGCCGCTCAATGGATCCTTGCAAGCTCTTGGTGCGCTTTTTCGGCTGTTGACGC-3’. This sequence was taken from a larger sequence shown to have high curvature (Sequence d1 in Ref. Segal et al., 2006).

Recent studies have suggested that DNA shape is an essential component of DNA–protein recognition. In particular, sequence attributes such as the minor groove width and the intrinsic curvature have been shown to play a roleWest et al. (2010); Rohs et al. (2009); Scipioni et al. (2004); Peters and Maher (2010). Building on such studies, we incorporate sequence-dependent shape into the 3SPN.2 model.

To implement sequence-dependent DNA shape, we followed the studies of Olson and coworkersOlson et al. (1998), who determined average base pair and base-step parameters for each of the ten unique DNA base-steps. For this work, we employed more recent average base-step parameters that were demonstrated to be well-suited for DNA–protein binding (Table 1)Morozov et al. (2009). Base pair parameters, shown in Table 2, are from Ref. Olson et al., 2006. It should be noted that employing base-step parameters (rather than longer range approaches such as trimeric parameters) has been shown to produce good agreement with experimental electrophoretic gel retardation studiesDe Santis et al. (1990).

In order to translate these base pair and base-step parameters into an actual structure for a given DNA sequence, we used X3DNALu and Olson (2003) to build an AA structure. We then coarse-grained the AA structure into 3SPN.2C by mapping each base, sugar, and phosphate to a bead placed at the respective center of mass. The result is a topology with sequence-dependent groove widths and intrinsic curvature, as shown in Fig. 1. The lengths and angles of each bond, bend, and dihedrals in this structure represent the minimum energy values of these interactions. Minimum energy distances and angles of base stacking, base pairing, and cross stacking interactions were also obtained from this structure.

Base-Step Twist Roll Tilt Shift Slide Rise
Å Å Å
AA 35.31 0.76 -1.84 -0.05 -0.21 3.27
AT 31.21 -1.39 0.00 0.00 -0.56 3.39
AC 31.52 0.91 -0.64 0.21 -0.54 3.39
AG 33.05 3.15 -1.48 0.12 -0.27 3.38
TA 36.20 5.25 0.00 0.00 0.03 3.34
TT 35.31 0.76 1.84 0.05 -0.21 3.27
TC 34.80 3.87 1.52 0.27 -0.03 3.35
TG 35.02 5.95 0.05 0.16 0.18 3.38
CA 35.02 5.95 -0.05 -0.16 0.18 3.38
CT 33.05 3.15 1.48 -0.12 -0.27 3.38
CC 33.17 3.86 0.40 0.02 -0.47 3.28
CG 35.30 4.29 0.00 0.00 0.57 3.49
GA 34.80 3.87 -1.52 -0.27 -0.03 3.35
GT 31.52 0.91 0.64 -0.21 -0.54 3.39
GC 34.38 0.67 0.00 0.00 -0.07 3.38
GG 33.17 3.86 -0.40 -0.02 -0.47 3.28
Table 1: Average base-step parameters for X3DNALu and Olson (2003) required to obtain the minimum energy configuration in 3SPN.2C. These values are from Ref. Morozov et al., 2009.
base pair Buckle Propeller Opening Shear Stretch Stagger
Å Å Å
A–T 1.8 -15.0 1.5 0.07 -0.19 0.07
T–A -1.8 -15.0 1.5 -0.07 -0.19 0.07
G–C 4.9 -8.7 -0.6 -0.16 -0.17 0.15
C–G -4.9 -8.7 -0.6 0.16 -0.17 0.15
Table 2: Average base pair parameters used to construct the 3SPN.2C model. These values are from Ref. Olson et al., 2006.

ii.3 Local Base-Step Flexibility

The local flexibility of DNA base-steps has been the subject of significant interest over the past several decadesOlson et al. (1998); Morozov et al. (2009); Olson et al. (2006). Olson and coworkers mined structural databases and calculated the relative flexibilities of DNA base-steps from these data Olson et al. (1998). In order to perform this analysis on the atomistic representation of DNA models, the many degrees of freedom (DOF) inherent in a DNA base-step (each atom has three spatial DOF) are reduced to six. Each base pair is represented by a plane; the six DOFs correspond to the three translational and three rotational transformations required to superimpose a base pair on its neighbor in the 3’ direction. By calculating the required transformations for a wide range of DNA structures in a variety of contexts (e.g. protein–bound, in solution, etc.), one can construct a 66 covariance matrix, , for each unique base-step. This covariance matrix can be used to estimate the conformational volume of the base-step, hereafter referred to as “,” as well as the force constants that penalize deformations of the six degrees of freedom for each base-stepOlson et al. (1998, 2006). In short, this analysis can give extensive insight into the local flexibility of each base-step in a global sense through or in great detail through information about the flexibility of each degree of freedom.

In 3SPN.2C, a local orthogonal coordinate system or “triad” is defined for each base pair, as shown in Fig. 2. The six DOFs mentioned previously are then the three translational and rotational transformations required to move from one base pair to the next base pair that constitutes a base-step. The triad for each base pair is defined as follows: first, a vector is drawn from the base site of the sense strand to the base site of the antisense strand. This constitutes the “” axis of the triad. Second, the “” axis is constructed by drawing a vector from the center-of-mass of the two base sites that constitute the base pair that points in the 3’ direction along the helical axis. In practice this is done by finding the projection of the “” axis on the helical axis; the “” axis is then a unit vector parallel to the vector rejection of the “” axis on the helical axis. Lastly, the “” axis is the cross product of the previous two vectors () and points in the direction of the major groove. This is done for each base pair in a sequence and the six parameters for each base-step are determined using the method described by Calladine, making use of the so-called “mid-step triad”El Hassan and Calladine (1995). Using an ensemble of configurations from direct simulations, the covariance matrix is determined for each base-step. The conformational volume, , of each base-step is calculated by taking the square root of the determinant of . provides a measure of the flexibility of a base-step (more accurately, it provides an estimate of the variability in the six degrees of freedom for the base-step), with a larger value indicating greater flexibility. (Interested readers are referred to Ref. Olson et al., 2006 for additional details).

Figure 2: Base-step triads for two base pairs comprising a base-step. T and T correspond to base pair triads at the 5’ and 3’ ends of a base-step, respectively. The conformational volume, , of a base-step is determined from the covariance matrix of the six degrees of freedom of a base-step, as described by Olson et al.Olson et al. (2006). The degrees of freedom are determined from the base-step triads according to the methodology of Calladine et al.El Hassan and Calladine (1995). The T triad is colored blue and the T is colored red.

In order to incorporate the experimental values of in 3SPN.2C, we set the force constants of the Base–Sugar–Phosphate (B–S–P, where B={A, T, G, or C}) and Sugar–Phosphate–Base (P–S–B) bending angles in the following manner:

(2)

where the index indicates a specific base-step, is the configurational volume of that step as characterized by Olson and coworkersOlson et al. (2006), and 0.4 and 4.1 are the minimum and maximum values of , normalized by , given by Olson et al. This limits values of between and . The choice of this function and the range of force constants are arbitrary but were found to give semi-quantitative agreement between simulation and experiment as shown in Fig. 3b. Table 3 in the Appendix provides the values of bending penalties used in the model.

Iii Results

Figure 3: (a) Consistency between melting temperatures calculated using the 3SPN.2C CG DNA model. The black error bar represents the predicted temperature of the sequence for which the model parameters were adjusted (5’-TACTAACATTAACTA-3’; = 69 mM). The red error bars represent the melting temperatures predicted for other sequences and ionic strengths. (b) Semi-quantitative agreement between experimentalOlson et al. (2006) and simulated local flexibility . To facilitate comparison of our data to that of Olson et al., all values of in the figure are normalized by , as done in Ref. Olson et al., 2006. The solid black line represents a best-fit to . (c) Comparison of the ability of each CG model to capture the effect of sequence on collective flexibility via the persistence length . Experimental persistence lengths are from Ref. Geggier and Vologodskii, 2010. The solid lines represent linear fits to the data and provide a guide to the eye. In all figures, the dotted line represents exact correspondence between simulation and experiment.

iii.1 Melting Temperature

While the 3SPN.2C model is intended to primarily simulate dsDNA, it is important to preserve the ability of dsDNA to melt. Melting temperature calculations were performed using multiple walker metadynamicsRaiteri et al. (2006), as done previouslyHinckley et al. (2013). Interstrand base pair interactions were scaled uniformly until good agreement was achieved between experimental and simulated melting temperatures for a reference sequence. The melting temperatures of several other sequences and ionic strengths where then predicted. Figure 3a demonstrates good agreement, as expected given the relatively minor modifications to the 3SPN.2 model. This result also highlights that, for a single validation metric, there exist many sets of CG parameters that provide satisfactory performance.

iii.2 Sequence-dependent Flexibility

iii.2.1 Local Flexibility

To assess the behavior of the model with respect to local flexibility, , we compare 3SPN.2C to the data of Olson and coworkersOlson et al. (2006). Simulations were performed for each possible DNA tetramer and the local flexibility of each base step was calculated by averaging over the sixteen tetramers centered on the base step. The agreement between 3SPN.2C and experimental values of , shown in Fig. 3b, is semi-quantitative, with a Pearson Correlation Coefficient of 0.88. Quantitative agreement is difficult to achieve because of the approximate methods used to define the base-step triads and the number of parameters that act in concert to affect local flexibility. It is possible to alter the parameters governing the anisotropy of the stacking and cross-stacking interactions in 3SPN.2C and so tune the local flexibility; however, this requires extensive parameterization via multiple iterations in order to recover simultaneously stacking, melting, and bulk persistence length behaviors. We view the semi-quantitative level of agreement between model and experiment achieved here as satisfactory, given the difficulty in determining the six relevant degrees of freedom in the coarse-grained model.

iii.2.2 Collective Flexibility

The persistence length of the 3SPN.2C representation of dsDNA is predominantly influenced by the strength of the bonded angle and dihedral potentials along the backbone. The strengths of the angle potentials were assigned previously, with the other energy parameters as described in Ref. Hinckley et al., 2013. The force constants of the Sugar–Phosphate–Sugar (S–P–S) and Phosphate–Sugar–Phosphate (P–S–P) angles along the DNA backbone, as well as the backbone dihedral force constants, were assigned using the experimental persistence length data of Geggier and VologodskiiGeggier and Vologodskii (2010). They determined the persistence length of 200 base pair DNA segments using cyclization assays and used the results to assign a bending penalty to each base-step. The S–P–S bending penalties in our model are given a force constant proportional to the base step bending penalties reported by Geggier and Vologodskii (Table 3). As with B–S–P and P–S–B bending penalties, the range of the force constant was chosen such that there was reasonable agreement with experimental persistence length data. The force constants of the P–S–P bend and the dihedral potentials were modified until the persistence length of dsDNA came into good agreement with Ref. Geggier and Vologodskii, 2010.

Persistence length calculations were performed using 75 base pair sequences taken from the middle of the sequences listed in Table S3 of Ref. Geggier and Vologodskii, 2010. Simulations were performed using both the 3SPN.2 and 3SPN.2C CG models. The persistence length was calculated from the resulting trajectories using the helical axis autocorrelation function, as in previous workHinckley et al. (2013). As demonstrated in Fig. 3c, 3SPN.2C is able to correctly capture trends in persistence length as a function of sequence. In contrast, 3SPN.2 is not able to capture the trends in persistence length, despite the dependence of intrastrand stacking energies on sequence. Such inability to capture trends in persistence length has also been observed in another coarse-grained DNA model with sequence-dependent base stacking energiesŠulc et al. (2012).

iii.3 Sequence-Dependent Minor Groove Width

3SPN.2C simulations were performed in order to determine the minor groove width profiles of unbound DNA sequences. The minor groove width was determined from the resulting trajectories employing the method of El Hassan and CalladineEl Hassan and Calladine (1995). In this method the minor groove width is defined by the distance between phosphate groups on opposing strands and refined to capture only the minimum separation between the two backbone curves.

To characterize the performance of 3SPN.2C with respect to sequence-dependent DNA shape, we employed three sources of experimental data. The first two rely on algorithms that predict the shape of the DNA minor groove based on experimental data employing photo-chemical cleavage of double-stranded DNA by the uranyl(IV) ionLindemose et al. (2011) and hydroxyl radical cleavage (ORChID2)Bishop et al. (2011), respectively. The third source of data is a library of DNA tetramers mined from the PDB databank by Rohs et al.Rohs et al. (2009). In all cases, a comparison is made to the minor groove width determined through simulations of all 256 possible tetramers. Figure 4 shows the correlation between 3SPN.2C and these three sources of experimental data. In panels (a) and (b), the middle 13 base-steps from tetrameric sequences are compared to model predictions (). In panel (c), only the minor groove width at the central base-step of the tetramer is compared to experimental data (). In this third comparison, each experimental tetrameric minor groove width represents the average from all PDB structures containing the tetramer. These are the same data given in Ref. Rohs et al., 2009. In all cases, the correlation is statistically significant ( in the worst case, the comparison to PDB data). However, the best correlation is between 3SPN.2C and the uranyl photo-cleavage modelLindemose et al. (2011).

Figure 4: Minor groove width from Langevin Dynamics simulations of all 256 possible tetramers compared to (a) Electronegativity score from uranyl(IV) ionLindemose et al. (2011), (b) ORChID2 predictionsBishop et al. (2011) and (c) average minor groove widths from PDB structuresRohs et al. (2009);

Representative examples of agreement between 3SPN.2C and the algorithms based on experimental cleavage data are given in Fig. 5. It is immediately clear that the qualitative agreement between 3SPN.2C and uranyl cleavage is better than that between 3SPN.2C and ORChID2. This is in large part due to the noise in the ORChID2 signal. This is especially apparent in Fig. 5e, in which both 3SPN.2C and uranyl cleavage predictions show an unchanging minor groove width for the sequence (CG) while ORChID2 predicts an oscillating cleavage signal.

Figure 5: Comparison of 3SPN.2C minor groove width profile (black lines) for six representative sequences to minor groove width prediction algorithms based on photo-chemical cleavage by uranyl (red lines) Lindemose et al. (2011) and hydroxyl radical cleavage (ORChID2, blue lines)Bishop et al. (2011). Panels (a-c) show AA, AT, and TA motifs while panels (d-f) show CC, GC, and CG motifs.

Iv Conclusion

We have presented an extension of an existing coarse-grained DNA model for simulating DNA interactions with proteins. The model uses experimental measurements of base-step parameters and mobility to inform the parameters imparting sequence-dependent shape and flexibility to dsDNA. We have demonstrated that the model correctly predicts the effect of sequence on the persistence length of dsDNA. We have also shown that the model is consistent with experimental measurements of minor groove width.

This model should find applications in scenarios where the structure and flexibility of DNA are important. A forthcoming publication will demonstrate the successful use of the model to examine the origin of DNA affinity in the nucleosome. The 3SPN.2C model has been implemented in LAMMPSPlimpton (1995); LAM () and is available online3SP ().

Acknowledgements.
The University of Wisconsin–Madison Center for High Throughput Computing is gratefully acknowledged for providing computational resources and computer expertise. We further acknowledge computational resources provide by the Midway computing cluster at the University of Chicago. G.S.F. and J.P.L. gratefully acknowledge partial support of this research by the NSF-funded University of Wisconsin Nanoscale Science and Engineering Center (NSEC). D.M.H. was funded by a Graduate Research Fellowship from the National Science Foundation (Grant No. DGE-1256259). J.K.W. and J.J.D.P. acknowledge support from UChicago Argonne, LLC, Operator of Argonne National Laboratory (“Argonne”). Argonne, a U.S. Department of Energy Office of Science laboratory, is operated under Contract No. DE-AC02-06CH11357.

*

Appendix A Shape-Dependent DNA Model Parameters

Changing the topology of the 3SPN.2C requires modification to the equilibrium bond lengths, bend angles, and dihedral angles. As stated in Section II.2, equilibrium DNA configurations are generated from X3DNALu and Olson (2003) using parameters from Refs. Morozov et al., 2009 and Olson et al., 2006. 3SPN.2C is mapped onto the resulting atomistic structure according to the method outlined in Ref. Hinckley et al., 2013. From this equilibrium structure, every bond, bend, and torsion in the structure is characterized and the relevant lengths and angles are written to files used as inputs for actual 3SPN.2C simulations. The code is available3SP () for generating 3SPN.2C topologies and performing simulations in LAMMPS using the modified energy parameters discussed below.

The energies of bonded and non-bonded interactions in 3SPN.2C differ from those in 3SPN.2. While the bond energies in 3SPN.2C are the same as 3SPN.2, the bending angles are different. The sequence-dependent bending energies are shown in Table 3. The magnitude of the dihedral energies are modified as follows: , the force constant for the Gaussian well-potential is modified from 6.0 kJ/mol/rad to 7.0 kJ/mol/rad and , the force constant of the additional dihedral potential (see Eq. 1) is 2.0 kJ/mol/rad.

The 3SPN.2C model also includes changes to the base pairing, cross-stacking, and base stacking energies. Base pairing and cross stacking interactions were scaled by a factor of 0.861. The base pairing energies were 14.41 kJ/mol for A–T and 18.24 kJ/mol for G–C. The base stacking energies were calculated using metadynamics simulations (see Ref. Hinckley et al., 2013 for additional details). The energies of base stacking and cross stacking interactions are provided in Table 4. The other parameters that appear in the angle-dependent potentials (, ; see Ref. Hinckley et al., 2013 for additional details) are unchanged from 3SPN.2.

Angle Base-Step Angle Base-Step Angle Base-Step
kJ/mol/rad kJ/mol/rad kJ/mol/rad
A-S-P AA 460 C-S-P CA 206 S-P-S AA 355
A-S-P AT 370 C-S-P CT 358 S-P-S AT 147
A-S-P AC 442 C-S-P CC 278 S-P-S AC 464
A-S-P AG 358 C-S-P CG 278 S-P-S AG 368
P-S-A AA 460 P-S-C AC 442 S-P-S TA 230
P-S-A TA 120 P-S-C TC 383 S-P-S TT 355
P-S-A CA 206 P-S-C CC 278 S-P-S TC 442
P-S-A GA 383 P-S-C GC 336 S-P-S TG 273
T-S-P TA 120 G-S-P GA 383 S-P-S CA 273
T-S-P TT 460 G-S-P GT 442 S-P-S CT 368
T-S-P TC 383 G-S-P GC 336 S-P-S CC 165
T-S-P TG 206 G-S-P GG 278 S-P-S CG 478
P-S-T AT 370 P-S-G AG 358 S-P-S GA 442
P-S-T TT 460 P-S-G TG 206 S-P-S GT 464
P-S-T CT 358 P-S-G CG 278 S-P-S GC 228
P-S-T GT 442 P-S-G GG 278 S-P-S GG 165
P-S-P all 300
Table 3: Bending angle energy constants employed in the 3SPN.2C DNA model. Single letter names for sites are as follows: A, T, C, and G denote the four DNA bases, S denotes a sugar moiety, and P denotes a phosphate moiety. B–S–P and P–S–B bend energies were assigned using Eq. 2; S–P–S bend energies were assigned using base-step bend energies from Ref. Geggier and Vologodskii, 2010.
(a)
Base
(kJ/mol)
Base
A T G C
A 13.82 15.05 13.32 15.82
T 9.15 12.44 9.58 13.11
G 13.76 14.59 14.77 15.17
C 9.25 12.42 8.83 14.01
(b)
Base
(kJ/mol)
Base
A T G C
A 1.882 2.388 2.439 1.680
T 2.388 1.882 2.187 2.566
G 2.439 2.187 3.250 0.972
C 1.680 2.566 0.972 4.135
(c)
Base
(kJ/mol)
Base
A T G C
A 1.882 2.388 2.566 2.187
T 2.388 1.882 1.680 2.439
G 2.566 1.680 4.135 0.972
C 2.187 2.439 0.972 3.250
Table 4: Base-stacking and cross-stacking energies for 3SPN.2C. Section (a) describes base-stacking energy scales. Sections (b) and (c) describe cross-stacking energy scales. Variables are as defined in Ref. Hinckley et al., 2013. Upward-pointing arrows denote the sense strand while downward-pointing arrows denote the anti-sense strand (for cross-stacking interactions).

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
208391
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description