Simple Elastic Network Models for Exhaustive Analysis of Long Double-Stranded DNA Dynamics with Sequence Geometry Dependence
Simple elastic network models of DNA were developed to reveal the structure-dynamics relationships for several nucleotide sequences. First, we propose a simple all-atom elastic network model of DNA that can explain the profiles of temperature factors for several crystal structures of DNA. Second, we propose a coarse-grained elastic network model of DNA, where each nucleotide is described only by one node. This model could effectively reproduce the detailed dynamics obtained with the all-atom elastic network model according to the sequence-dependent geometry. Through normal-mode analysis for the coarse-grained elastic network model, we exhaustively analyzed the dynamic features of a large number of long DNA sequences, approximately bp in length. These analyses revealed positive correlations between the nucleosome-forming abilities and the inter-strand fluctuation strength of double-stranded DNA for several DNA sequences.
Elastic network models of proteins, including all-atom models(1); (2) and coarse-grained models(3); (4); (5); (6); (7); (8); (9) represent some of the simplest and most powerful types of theoretical models that can accurately reveal structure-dynamics relationships and the mechanisms underlying a protein’s functional activities (10); (11); (12); (13); (14). Such models have also been widely employed to accurately reproduce the temperature factors on the crystal structure of a protein via normal-mode analysis. These models can thus demonstrate the large and slow deformations of proteins that are essential for protein functions, but remain difficult to demonstrate via all-atom molecular dynamics simulations.
Along with proteins, DNA is the most important biomolecule for the activities of all living organisms. Recent developments in molecular biology have revealed that DNA contains several functional regions. Therefore DNA is no longer considered to have the sole function in the storage of genetic information, but is now known to be actively involved in gene regulation(15); (16); (17); (18), insulator activity(19); (20); (21); (22); (23); (24), and in the construction of chromosomal architectures(26); (27); (25); (28) such as heterochromatins and topologically associated domains through nucleosome formation and protein bindings(29); (30); (31); (32); (33). The functional behavior of each strand of DNA is determined not only by chemical aspects of the nucleotides and base pairs but also by its physical characteristics such as the structure and dynamics of the nucleotide sequences in each functional region. However, comprehensive understanding of the physical properties of nucleotide sequences lags far behind the knowledge accumulated of their biochemical aspects (15); (16); (26).
Since the last century, much progress has been made in revealing the physical aspects of DNA using all-atom normal-mode analysis(34); (35) and molecular dynamics simulations(36); (37); (38); (39); (40); (41); (42). Several coarse-grained models of DNA (and RNA) have also been proposed. Some of these describe the detailed shape of each nucleotide using three or more particles(43); (44); (45); (46); (47); (48); (49), whereas others describe each base pair by simply one or two particles(50); (52); (51); (53); (54). Molecular dynamics simulations and normal-mode analysis of these models have identified the flexibilities, nucleosome-forming abilities, and zip-unzip transitions of the double helices of some specific DNA sequences from tens to a few hundred base pairs in length. Although these methods have proven to be very powerful for the analysis of the physical aspects of DNA, molecular dynamics simulations are not suitable for conducting exhaustive analyses of several sequences simultaneously owing to the high computational costs of such extensive simulations; for example, analyses of whole genomic and whole possible sequences. Moreover, the normal-mode analysis of all-atomic models of long DNA sequences is also computationally heavy.
Alternatively, data-driven methods have been proposed for determining the mechanical properties of DNA with respect to the helical parameters and local flexibilities of base pairs from X-ray crystal structure analysis and all-atom molecular dynamics (55); (56); (57); (58); (59); (60). These methods also seem to be powerful and can be applied to the analysis of the mechanical properties of several DNA sequences simultaneously. However, the methods thus far proposed have focused on static and local mechanical properties. Therefore, they are not particularly useful for the study of the functional contribution of the dynamic and correlated motions of DNA.
The objective of this study was to construct a simple coarse-grained elastic network model of DNA to allow for exhaustive analysis of the dynamic correlated motions of several long DNA sequences. For this purpose, we first constructed a simple all-atom elastic network model of short DNA sequences based on the method introduced by Tirion(1) for the modeling of protein dynamics. We confirmed that this model could accurately reproduce the temperature factors of some DNA fragments from data obtained with X-ray crystal structure analysis.
Second, we developed a simple coarse-grained elastic network model of long DNA sequences, where each nucleotide is described by only one node. We confirmed that this simplified model has lower computational costs but can nonetheless reproduce the nucleotide sequence-dependent dynamics revealed by the all-atom elastic network model.
Finally, through the normal-mode analysis of this coarse-grained model, we conducted an exhaustive analysis of the general features of the sequence-dependent structure-dynamics relationships among several DNA sequences. We specifically focused on the dynamic properties of a large number of long DNA sequences, approximately bp in length (with respect to the length of the nucleosome-forming regions), for the genomes of some model organisms as well as for random sequences with several A, T, C, G ratios. Through these analyses, we found that the dynamic aspects of DNA are highly influenced by their sequences, and found positive correlations between the nucleosome-forming abilities and inter-strand fluctuations of double-stranded DNA.
In particular, we focus on the geometry dependencies of the dynamics of several DNA sequences, since a recent study demonstrated that the geometry of DNA sequences has a dominant contribution to their mechanical features(47). In the following arguments, for simplicity, we use “sequence-dependent” to mean “dynamics that depend on sequence geometry”.
Ii Models and Methods
ii.1 Basic structures of Elastic Network Models of DNA
In order to construct elastic network models of DNA, the basic DNA structure first needs to be determined. In the following arguments, we construct two types of models: i) a simple all-atom elastic network model (AAENM) to reproduce the characteristics of the crystal structures of DNA, and ii) a simple coarse-grain elastic network model (CGENM) to reproduce the characteristics of the AAENM. We obtained the basic DNA structures in the following two ways for the respective purposes of constructing each model.
For construction of the AAENM, we employed the atom coordinate sets of several naked DNA crystal structures included in the Protein Data Bank (PDB) as the given basic structures of the model. We used DNA crystal structures containing only DNA and water molecules under different conditions of crystallization, where all temperature factors are given as positive values (Table 1). The suitability of the model was evaluated by comparisons of the temperature factors obtained between the model and those obtained from X-ray crystal structure analysis.
The objective of our constructed arguments was to unveil the sequence-dependent dynamic features of several long DNA sequences simultaneously. In recent crystal structure analysis, only shorter DNA sequences (i.e., less than bp in length) were studied. On the other hand, several types of helical parameter sets, base pair parameters, and base step parameter sets have been proposed through experiments or molecular dynamics simulations(36); (61); (62); (63); (64); (65); (66)(Table 2, S1 Table). Here, we used the application X3DNA(67) to obtain the coordinates of each atom for any sequence from these helical parameters, and used these parameters to construct the AAENMs and CGENMs of longer DNA sequences (for detailed methods of the generation of the atom coordinates, see Lu and Olson (67)). In the following arguments, we compare the characteristics of the AAENM and CGENM constructed by X3DNA from three different helical parameter sets: i) a parameter set obtained by an in vitro experiment and X-ray crystal structure analysis(48); (61); (62); (63) (Table 2), ii) a parameter set obtained only by the X-ray crystal structure analysis, and iii) a parameter set obtained by all-atom molecular dynamics simulations(58); (65); (66) (S1 Table).
ii.2 All-atom Elastic Network Model of DNA
A simple all-atom elastic network model of double-stranded DNA was constructed based on the model proposed by Tirion(1). In this model, we regard all of the atoms of given DNA sequences as the nodes comprising the elastic network. For the analysis of the crystal structures of DNA involving water molecules, we also regard the oxygen atoms of water around DNA as the nodes of the elastic network. We define the mass and position of atom as and (), respectively. The potential of all atoms is given as
Here, is the position of atom of the basic DNA structure, as described above.
The first term indicates the interaction potential among atoms that are spatially closed in the basic DNA structure (Fig. 1(a)). Here, refers to the Van der Waals radius of atom , is an arbitrary cut-off parameter that models the decay of interactions with distance, and indicates the Heaviside function, where () for (). We assume , which is empirically considered to be an appropriate range for biomolecules, at least for proteins(1). The results of the arguments presented in this paper were qualitatively unchanged for the appropriate range of . In the crystal structures of DNA, the motions of water molecules are also restricted by the crystal packing. Thus, for all atoms, we assume that spatially closed pairs of atoms are connected by linear springs with their respective natural lengths. The elastic coefficient () is a phenomenological constant, which is assumed to be the same for all interacting pairs.
The second term indicates the boundary effects of each DNA and water molecule in each DNA crystal structure. This term plays a crucial role for the analysis of the fluctuations of the crystal structure of DNA, such as temperature factors, since the fluctuations of the nucleotides at the edges of DNA and water molecules are restricted due to the following facts.
In the crystal of DNA, the motion of edges at upper and lower streams of each DNA segment (left and right edges in Fig. 1(a)) is influenced by atoms of other adjacent segments in the long-axis direction of DNA segments. Moreover, the motions of water molecules around each DNA segment are influenced by atoms of other adjacent water molecules or DNA segments in the direction perpendicular to the longitudinal axis of the DNA segment. Thus, we need to consider the second term of Eq. (1) to implement such effects, where indicates the strength of such effects for atoms in the edge nucleotides of each DNA and water molecule, respectively.
Remarkably, as shown in the Results section, the distributions of the temperature factor of atoms exhibited various patterns from the same sequence of crystallized DNA (S1 Fig and S2 Fig) owing to the dependency on the conditions of crystallization. Therefore, in order to compare the characteristics of the present DNA model to those of the crystal structure of DNA, appropriate values of need to be assigned to the atom that belongs to the edge base pairs and water molecules. For simplicity, we assign and to the atoms at the edge base pairs and the water molecules, respectively, and otherwise.
It is noted that the internal structure of the crystal of DNA is spatially anisotropic. Thus, it is reasonable to assume that the interactions among local parts of the crystal of DNA exhibit different strengths in different directions. Accordingly, in general, the strengths of the restrictions of atoms belonging to the edge of DNA differ from those of water molecules. Thus, we assume that and are different values (Table 1).
For simplicity, the mass of each atom is assumed as a constant value (, is Avogadro’s number). However, we confirmed that the results were almost identical when using the precise masses of the atoms.
ii.3 Coarse-grained Elastic Network Model of DNA
A simple coarse-grained elastic network model of double-stranded DNA, where each nucleotide is described as one node, was constructed as follows (Fig. 1(a)). We define the coordinate of the C1’ carbon of nucleotide , (), as the position of nucleotide , and regard the motion of the C1’ carbon as that of the nucleotide. Here, we assume that the mass of the C1’ carbon obeys . The potential of all nucleotides is given as
Here, is the position of nucleotide of the basic structure of DNA, as defined above. We assume that nucleotides and belong to the same base pair. For nucleotide , the sum is restricted to the pair of nucleotides in the same base pair (), in the neighboring base pairs (), in the next neighboring bases pairs (), and in the next-next neighboring bases in another strand () (Fig. 1(b)). The elastic coefficient () is a phenomenological constant that is assumed to be the same for all interacting pairs. This model is considered as a simplified version of previously proposed one-site-per-nucleotide models(50); (53).
ii.4 Normal-mode Analysis
An overview of the theory of normal-mode analysis is provided in several recent studies(1); (2); (3); (4); (10); (11); (12); (13); (5); (6); (7). Thus, we here briefly show the results of this analysis. For this analysis, we define (, ) as a -dimensional position vector, and as the position vector of the basic structure. Here, for the AAENM and for the CGENM. The motions of small deviations of from , obey
where and () are the -th largest eigenvalue and its eigenvector of the Hessian matrix as
We assume that the system is in thermodynamic equilibrium with temperature . Thus, the amplitude is given as
with Boltzmann constant . Using this solution, the mean square fluctuation of the -th atom in the AAENM () is obtained as
with Boltzmann constant, and the temperature factor is displayed as . Here, represents the temporal average.
For the CGENM, we define the mean square fluctuation of the -th nucleotide () as . To consider the motion of the -th nucleotide in the AAENM, we define the average nucleotide motions as . Using this vector, we define the mean square fluctuation of the -th nucleotide as . For the motif of the -th nucleotide ( sugar, base, and phosphoric acid), the average temperature factor of the motif () is defined as the average of belonging to each motif of each nucleotide, .
ii.5 Treatment of the Temperature Factor in X-ray crystal Structure Analysis
In order to evaluate the validity of the AAENM, we measured the correlation coefficient between the profile of the temperature factor obtained from PDB data (via X-ray crystal structure analysis) and that obtained from the AAENM based on this crystal structure. It is noted that the temperature factor profiles for some of the PDB data often include unnaturally large or small values. Thus, the correlation coefficients were estimated using data excluding such outliers. In the present evaluations, the value was considered as an outlier if , where and are the average and standard deviation of , respectively, and is used based on the standard arguments of statistics.
ii.6 Evaluations of Anisotropic Fluctuations of DNA
We also focused on the relationships between the fluctuations of each nucleotide in the AAENMs and CGENMs in the directions parallel to the base pair axis, parallel to the helix axis, and vertical to both the base pair and helix axes.
Here, we name the nucleotides in one and the other strand constructing the -th base pair as the -th and -th nucleotide. We define the position vectors of the C1’ atoms belonging to the -th and -th nucleotides as and , and consider
(Fig. 1(c)). It is noted that and are not orthogonal in general; however, we confirmed that the angles of these vectors were always sufficiently close to rad for each .
Using these vectors, the mean square fluctuations of the -th and -th nucleotides of the AAENM in the directions parallel to the base pair axis are defined as and , those parallel to the helix axis are defined as and , and those in the torsional direction are defined as and , respectively. The mean square fluctuations of the -th and -th nucleotides of the CGENM in these same directions are obtained by and , and , and and , respectively. Moreover, we consider the mean square fluctuation of the relative base position of each base pair of the CGENM in the three directions listed above given by , , and , respectively.
ii.7 Evaluations of the Overall Geometry of DNA
The overall geometry of each modeled DNA molecule is characterized by the ratios among the square root of the three principal components of the populations of atoms, , , and (). Here, , , and are obtained as eigenvalues of the covariant matrix
where , is the position of -th atom, is the position of the center of mass of a given DNA molecule, and indicates the average for all s. We evaluate the overall geometry of a given DNA molecule using the linearity and the line symmetry with respect to the -axis . Here, is large when the DNA is straightened, and is large (small) when the DNA forms a three (two)-dimensional curve with wide (flat) envelope.
ii.8 Evaluations of Correlations Among the Results of AAENM, CGENM, and Experiments
We employ Pearson’s correlation coefficients, , to evaluate the correlations among the profiles of temperature factors and several anisotropic fluctuations of atoms obtained by AAENM, CGENM, and experiments.
Iii Results and Discussion
iii.1 Comparisons of Fluctuations between the AAENM and the Crystal Structure of DNA
The fluctuations of atoms of the AAENMs of several short DNA sequences were measured with normal-mode analysis. Here, the basic structures of the DNAs are given by the crystal structures of the naked DNAs with the following PDB IDs: 1BNA, 7BNA, 9BNA, 1D91, 1DC0, 122D, 123D, 181D, and 330D (Table 1(68); (69); (70); (71); (72); (73); (74); (75)). To confirm the validities of the AAENMs, the correlations between the distribution profiles of the temperature factor of atoms () and the average temperature factor of the motifs () in the crystal structures and those in the corresponding AAENMs were measured. In the following arguments, we employ the optimal values of , , and for each crystal structure (Table 1), which were manually found to maximize of between the results of the crystal structure analysis and those of the AAENM.
By choosing the appropriate conditions for the atoms of the two edge base pairs and water molecules ( and ) for each PDB ID, of each AAENM exhibited a similar profile to that of the crystal structure with a significant correlation coefficient (Fig. 2(a), Table 1, and S1 Fig). Therefore, the AAENM seemed to reproduce the overall structure of the temperature factor profile of each crystal structure well, although the detailed profiles among atoms showed some deviations.
Furthermore, we focused on the average temperature factor for the motifs, , as recently discussed (35). We found that the obtained with the AAENMs could accurately reproduce those of the corresponding crystal structures with high correlation coefficients , where was obtained for most cases (Fig. 2(b), Table 1, and S2 Fig). These results demonstrate that the AAENM is a suitable model for describing the sequence-dependent fluctuations of the nucleotide motifs of several double-stranded DNA sequences, despite the simplicity of model construction and its easy implementation. Moreover, these results show that the sequence-dependent forms of DNA have a dominant contribution to their overall flexibilities and fluctuations.
Nevertheless, the present AAENM could not accurately reproduce well for some crystal structures. These deviations are considered to have arisen from the following primary assumption: we only considered the effects of the restrictions by the packing of DNAs in crystal form for two edge base pairs and water, whereas such effects have an influence on all atoms. Therefore, obtaining and incorporating more detailed knowledge of the restrictions of bulk sequences caused by crystal packing should help to achieve a more accurate reproduction of the molecular fluctuations for all cases.
iii.2 Comparison between the AAENMs and CGENMs of DNA
The main objective of this study was to unveil the sequence-dependent dynamic correlated motions of several long DNA sequences. We next describe these dynamics of DNA sequences with longer lengths than considered in the previous subsection. For this purpose, we also constructed a coarse-grained model, which is often useful for focusing on the slow and large-scale dynamics of molecules that essentially influence their function. Thus, to propose a coarse-grained model of long double-stranded DNA, we evaluated whether the CGENM proposed provides an appropriate coarse-grained model of the present AAENM.
We performed normal-mode analysis of the AAENM and corresponding CGENM for randomly chosen -bp sequences, and compared the mean square fluctuations of the -th nucleotide ( and ) in the directions parallel to the base pair axis ( and ), parallel to the helix axis ( and ), and in the torsional direction ( and ). Here, is employed, as in the previous study(1), and is assumed, which was manually found to provide the best fit of fluctuation profiles between AAENM and CGENM. Here, the overall fluctuation profiles of CGENM are independent of the value of since influences only on the absolute values of fluctuations. Independent of the sequences and employed helical parameters, we found that the fluctuations of each nucleotide in several directions were very similar between the AAENMs and CGENMs when these models are constructed with the same helical parameters, with average correlation coefficients (Fig. 3, Table 3, S2 Table, S3 Fig, and S4 Fig).
It is noted that the present CGENM contains only one node per nucleotide, whereas the AAENM contains nodes (atoms) per nucleotide. This fact demonstrates that the computational costs of the CGENMs are much lower than those of the AAENMs, although the accuracies of the obtained statistical aspects are almost identical between these two models. Thus, this CGENM could be used for exhaustive analysis and comparisons of the dynamic features of several sequence-dependent DNAs related to protein binding affinities, functions of transcription regulation sequences, and nucleosome positioning(16); (26); (27); (17); (18); (19); (20); (21); (22); (23); (24); (46); (55); (56); (57); (58); (59); (60). In the next subsection, we provide an example of such an analysis to determine the relationships between the nucleosome-forming abilities of several double-stranded DNA sequences and their inter-strand dynamic features.
iii.3 Exhaustive Analysis of the Sequence-dependent Behaviors of -bp DNAs with the CGENM
Nucleosome positioning is important not only for compacting DNA but also for appropriate gene regulation. Several recent studies have been performed for genome-wide nucleosome mapping and the identifications and predictions of nucleosome-forming and -inhibiting sequences for some model organisms(57); (59); (76); (77); (78); (79); (80).
As an example of the applications of the CGENM to an exhaustive analysis of the sequence-dependent behavior of DNA, we compared the dynamic features of DNA sequences of bp that were predicted as nucleosome-forming or nucleosome-inhibiting sequences in the genome of budding yeast Saccharomyces cerevisiae (5000 forming sequences and 5000 inhibiting sequences of bp)(57), nematode Caenorhabditis elegans (2567 forming sequences and 2608 inhibiting sequences of bp), Drosophila melanogaster (2900 forming sequences and 2850 inhibiting sequences of bp), and Homo sapiens (2273 forming sequences and 2300 inhibiting sequences of bp)(59). The histograms of the average relative fluctuations of DNAs for the three directions , , , , , , , and ( indicates the average for all s.) for the nucleosome-forming sequences and the nucleosome-inhibiting sequences are shown in Fig. 4 and S5 Fig S7 Fig. Here, we employed the helical parameter set (i) (Table 2) that was used in the coarse-grained molecular dynamics simulations by Freeman et al., which exhibited consistent results to some experiments(48); (61); (62); (63).
The histograms for budding yeast showed that that nucleosome-forming sequences tend to exhibit larger fluctuations in several directions compared to the inhibiting sequences. In particular, the histogram of for the nucleosome-forming sequences showed a clear shift in the direction toward larger values compared to that for the nucleosome-inhibiting sequences (Fig. 4). For the other organisms, most of the histograms showed few differences between the nucleosome- forming and -inhibiting sequences. However, similar to the case of yeast, the distribution of for the nucleosome-forming sequences shifted largely in the direction of larger values compared to that for the nucleosome-inhibiting sequences in these organisms (S5 Fig S7 Fig). These results indicate that the nucleosome-forming ability is highly correlated to the fluctuations of the inter-strand distances of DNAs, in which sequences with larger fluctuations tend to form the nucleosome.
Similar to the previous arguments, we compared the dynamical features among the random -bp DNA sequences varying in average GC-contents; the histograms, averages, and standard deviations of , , , , and were measured from 10,000 random sequences for each GC content (Fig. 5 and S8 Fig). In this case, exhibited a minimum at a GC content of , which indicates that the AT-rich sequences tend to be more rigid than the GC-rich sequences. However, the fluctuations of sequences consisting of only A or T were as large as those of the GC-rich sequences. Moreover, the GC content dependencies of , , , and showed different characteristics for GC contents larger or smaller than . In particular, the results for were similar for cases with a large GC content () but monotonically decreased with a GC content with little variance. A recent experimental study showed that the probability of nucleosome formation tends to increase with increases in the GC content ratio (81). Thus, the present results indicate that sequences with larger fluctuations of inter-strand distances tend to form the nucleosome, which is consistent with the results described above from the analysis of the genomes of the four model organisms.
Finally, we focus on the relationships between the overall geometries and fluctuations of the considered DNA sequences. The overall geometries of several DNA sequences, the nucleosome-forming and -inhibiting DNA sequences for four model organisms and random DNA sequences with different GC contents, were evaluated using scatter plots of (linearity) and (line symmetry) (Fig. 6 and S9 Fig). It is noted that and are highly correlated. The average, standard deviation, and distributions of and exhibited slight but not significant deviations between the nucleosome-forming and -inhibiting sequences.
For the random sequences, the average value of exhibited similar variations to with an increase in GC content. In particular, both values decreased with increases in GC content for GC contents , whereas they increased with increases in GC content for GC contents (Fig. 5(c) and Fig. 6(c)). In fact, and were highly correlated for random DNA sequences, regardless of their CG content. The Pearson correlation coefficient for the 110,000 sequences analyzed above with GC contents was . It is noted that showed large dispersion for each GC content, and there were significant overlaps among distributions with different GC contents (Fig 6(b) and 6(c)). This fact indicates that different DNA sequences can often show similar geometries, and such sequences also tend to show similar overall fluctuations. On the other hand, the fluctuations of inter-strand distances that may correlate to the nucleosome-forming ability did not correlate significantly to either or , with Pearson correlation coefficients of and . This fact indicates that the nucleosome-forming ability of DNA sequences are not only determined by the overall DNA geometries but also by their dynamic properties.
Iv Summary and Conclusion
In this study, simple elastic network models of double-stranded DNAs were developed in order to perform an exhaustive analysis of several sequence-dependent dynamical features. First, we constructed a simple all-atom elastic network model that could reproduce the fluctuations of the motifs of each nucleotide (sugar, phosphoric acid, and bases) of several crystal structures of short DNA sequences. Second, we proposed a simple course-grained elastic network model that could reproduce the dynamic features of the long DNA sequences obtained by the all-atom elastic network model. Through exhaustive analysis of the dynamic features of several DNA sequences with normal-mode analysis of the presented coarse-grained elastic network model, we found that the dynamic aspects of DNA are highly influenced by the properties of nucleotide sequence such as GC content. We also found that the nucleosome-forming abilities of double-stranded DNA exhibited positive correlations with their sequence-dependent inter-strand fluctuations.
In the present study, we demonstrated the sequence-dependent dynamic features for several -bp DNA sequences to evaluate the relationships between the nucleosome-forming abilities and DNA dynamics. Of course, DNA sequences longer than bp can also be analyzed using the presented course-grained model. Moreover, course-grained molecular dynamics simulations can be performed to consider the large deformations of DNA, such as formation of a super helix and nucleosome that are the basic structures of higher-order chromosome architectures, using the presented elastic network models with the excluded effect of the volume of each atom or each nucleotide. We are currently conducting these molecular dynamics simulations, and the results will be reported in the future. We did not consider the effects of solvents such as temperature and salt concentrations in the present elastic network models. Therefore, we are also planning to attempt modifications of the models so that several solvent conditions can be incorporated in future work.
The author is grateful to S. Tate and Y. Murayama for fruitful discussions.
Conceived the study: AA, NS, HN. Designed the models: NS, AA, SI. Performed the analysis: SI, AA. Wrote the paper: AA, NS, HN.
Comparisons between CGENM and AAENM. Average and standard deviation of the correlation coefficients of 500 samples of random -bp sequences between and , and , , and , and and . Helical parameter sets (ii) and (iii) (S1 Table) were used.
Temperature factor of each atom. The temperature factor of each atom obtained by the AAENMs (black curve) and X-ray crystal structure analysis (gray (red) curve) of typical double-stranded DNAs obtained from PDB ID (a) 7BNA, (b) 9BNA, (c) 1D91, (d) 1DC0, (e) 122D, (f) 123D, (g) 181D, and (h) 330D. Parameters , , and are given in S1 Table.
Average temperature factor of motifs. Average temperature factor of motifs obtained by the AAENMs (black curve) and X-ray crystal structure analysis (gray (red) curve) of typical double-stranded DNAs obtained from PDB ID (a) 7BNA, (b) 9BNA, (c) 1D91, (d) 1DC0, (e) 122D, (f) 123D, (g) 181D, and (h) 330D. Parameters , , and are given in S1 Table.
Comparisons between AAENM and CGENM. Comparisons of the fluctuations between each nucleotide in the AAENM (black curves) and CGENM (gray (red) curves) for a typical random -bp DNA sequence (5’ - AGTGGTAAGGCATGGTTCTCGAATCTCGGTTTATTTACACTGCTGCTCCA - 3’). (a) and , (b) and , (c) and , and (d) and using helical parameter set (ii) S2 Table.
Comparisons between AAENM and CGENM. Comparisons of the fluctuations between each nucleotide in the AAENM (black curves) and CGENM (gray (red) curves) for a typical random -bp DNA sequence (5’ - ATATGCTGTAGAGCGTCCCGTCCGCGCGTTGTGGTTTTTTCGGTGCTCTA - 3’). (a) and , (b) and , (c) and , and (d) and using helical parameter set (iii) S1 Table.
Histograms of the average fluctuations inCaenorhabditis elegans. Histograms of the average fluctuations of (a) and , (b) and , (c) and , and (d) and for nucleosome-forming and nucleosome-inhibiting sequences of the nematode Caenorhabditis elegans ( bp). Helical parameter set (i) (Table 2) was used.
Histograms of the average fluctuations inDrosophila melanogaster. Histograms of the average fluctuations of (a) and , (b) and , (c) and , and (d) and for nucleosome-forming and nucleosome-inhibiting sequences of Drosophila melanogaster ( bp). Helical parameter set (i) (Table 2) was used.
Histograms of the average fluctuations inHomo sapiens. Histograms of the average fluctuations of (a) and , (b) and , (c) and , and (d) and for nucleosome-forming and nucleosome-inhibiting sequences of Homo sapiens ( bp). Helical parameter set (i) (Table 2) was used.
Histograms of the average fluctuations with different GC contents. Histograms of (a) , (b) , and (c) for 10,000 samples of random -bp sequences with different average GC contents. Helical parameter set (i) (Table 2) was used.
Overall geometries of long DNA sequences in model organisms. Scatter plots of and for nucleosome-forming and -inhibiting sequences of (a) Caenorhabditis elegans ( bp), (b) Drosophila melanogaster ( bp), and (c) Homo sapiens ( bp). Helical parameter set (i) (Table 2) was used.
- Tirion MM (1996) Large amplitude elastic motions in proteins from a single-parameter, atomic analysis. Phys Rev Lett 77: 1905-1908. doi:http://dx.doi.org/10.1103/PhysRevLett.77.1905
- Hinsen K, Petrescu AJ, Dellerue S, Bellisent-Funel MC, Kneller GR (2000) Harmonicity in slow protein dynamics. Chem Phys 261: 25-37. doi:10.1016/S0301-0104(00)00222-6
- Atilgan AR, Durell SR, Jernigan RL, Demirel MC, Keskin O, et al. (2001) Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J 80: 505-515. doi:10.1016/S0006-3495(01)76033-X
- Bahar I, Rader AJ (2005) Coarse-grained normal mode analysis in structural biology. Curr Opin Struct Biol 15: 586-592. doi:10.1016/j.sbi.2005.08.007
- Bahar I, Lezon TR, Yang LW, Eyal E (2010) Global dynamics of proteins: bridging between structure and function. Annu Rev Biophys 39: 23-42. doi: 10.1146/annurev.biophys.093008.131258
- Bahar I, Lezon TR, Bakan A, Shrivastava IH (2010) Normal mode analysis of biomolecular structures: functional mechanisms of membrane proteins. Chem Rev 110: 1463-1497.
- Dykeman EC, Sankey OF (2010) Normal mode analysis and applications in biological physics. J Phys Condens Matter 22: 423202.
- Flechsig H, Mikhailov AS (2010) Tracing entire operation cycles of molecular motor hepatitis C virus helicase in structurally resolved dynamical simulations. Proc Natl Acad Sci (USA) 107: 20875-20880.
- Atilgan C, Okan OB, Atilgan AR (2012) Network-based models as tools hinting at nonevident protein functionality. Annu Rev Biophys 41: 205-25.
- Van Wynsberghe AW, Cui Q (2006) Interpreting correlated motions using normal mode analysis. Structure 14: 1647-1653. doi:10.1016/j.str.2006.09.003
- Tama F, Brooks CL (2006) Symmetry, form, and shape: guiding principles for robustness in macromolecular machines. Annu Rev Biophys Biomol Struct 35: 115-33. 0.1146/annurev.biophys.35.040405.102010
- Moritsugu K, Smith JC (2007) Coarse-grained biomolecular simulation with REACH: realistic extension algorithm via covariance Hessian. Biophys J 93: 3460-3469. doi:10.1529/biophysj.107.111898
- Yang L, Song G, Jernigan RL (2009) Protein elastic network models and the ranges of cooperativity. Proc Natl Acad Sci (USA) 106: 12347-12352. doi: 10.1073/pnas.0902159106
- Dehouck Y, Mikhailov AS (2013) Effective harmonic potentials: insights into the internal cooperativity and sequence-specificity of protein dynamics. PLoS Comp Bio 9: e1003209
- Albert B, et al. (2008) Molecular Biology of the Cell 5th edition. Garland Science, Taylor & Francis Group, LLC.
- Latchman DS (2010) Gene Control. Garland Science, Taylor & Francis Group, LLC.
- Fukue Y, Sumida N, Nishikawa J, Ohyama T (2004) Core promoter elements of eukaryotic genes have a highly distinctive mechanical property. Nucleic Acids Res 32: 5834-5840. doi: 10.1093/nar/gkh905
- Ohyama T (2005) DNA Conformation and Transcription. New York, NY, USA: Springer, Landes Bioscience.
- Valenzuela L, Kamakaka RT (2006) Chromatin insulators. Annu Rev Genet 40: 107-138. doi: 10.1146/annurev.genet.39.073003.113546
- Bushey AM, Dorman ER, Corces VG (2008) Chromatin insulators: regulatory mechanisms and epigenetic inheritance. Mol Cell 32: 1-9. doi:10.1016/j.molcel.2008.08.017
- Geyer PK, Clark I (2002) Protecting against promiscuity: the regulatory role of insulators. Cell Mol Life Sci 59: 2112-2127. doi:1420-682X/02/122112-16
- Gaszner M, Felsenfeld G (2006) Insulators: exploiting transcriptional and epigenetic mechanisms. Nat Rev Genet 7: 703-713. doi:10.1038/nrg1925
- Raab JR, Kamakaka RT (2010) Insulators and promoters: closer than we think. Nat Rev Genet 11: 439-446.
- Takagi H, Inai Y, Watanabe S, Tatemoto S, Yajima M, Akasaka K, et al. (2012) Nucleosome exclusion from the interspecies-conserved central AT-rich region of the Ars insulator. J Biochem 151: 75-87. doi: 10.1093/jb/mvr118
- Turner BM (2001) Chromatin and Gene Regulation. Blackwell Science Ltd.
- Dekker J, Rippe K, Dekker M, and Kleckner N (2002) Capturing chromosome conformation. Science 295: 1306 1311. doi: 10.1126/science.1067799
- Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W (2002) Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell 10: 1453-1465. doi:10.1016/S1097-2765(02)00781-5
- Allis CD et al. (2007) Epigenetics. Cold Spring Harbor Laboratory Press.
- Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. (2012) Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature 485: 376-380. doi:10.1038/nature11082
- Apostolou E, Ferrari F, Walsh RM, Bar-Nur O, Stadtfeld M, Cheloufi S, et al. (2013) Genome-wide chromatin interactions of the Nanog locus in pluripotency, differentiation, and reprogramming. Cell Stem Cell 13: 1-14. doi:10.1016/j.stem.2013.04.013
- Gibcus JH, Dekker J (2013) The hierarchy of the 3D genome. Mol Cell 49: 773-782. doi:10.1016/j.molcel.2013.02.011
- Dekker J, Marti-Renom MA, Mirny LA (2013) Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data. Nat Rev Gen 14: 390-403. doi:10.1038/nrg3454
- Bickmore WA, van Steensel B (2013) Genome architecture: domain organization of interphase chromosomes. Cell 152: 1270-1284. doi:10.1016/j.cell.2013.02.001
- Lin D, Matsumoto A, Go N (1997) Normal mode analysis of a double-stranded DNA dodecamer d(CGCGAATTCGCG). J Chem Phys 107: 3684-3690. doi: 10.1063/1.474724
- Matsumoto A, Olson WK (2002) Sequence-dependent motions of DNA: a normal mode analysis at the base-pair level. Biophys J 83: 22-41. doi:10.1016/S0006-3495(02)75147-3
- Lankas F, Sponer J, Langowski J, Cheatham TE 3rd (2003) DNA basepair step deformability inferred from molecular dynamics simulations. Biophys J 85: 2872-2883. doi:10.1016/S0006-3495(03)74710-9
- Fujii S, Kono H, Takenaka S, Go N, Sarai A (2007) Sequence-dependent DNA deformability studied using molecular dynamics simulations. Nucl Acids Res 35: 6063-6074. doi: 10.1093/nar/gkm627
- Orozco M, Noy A, Prez A (2008) Recent advances in the study of nucleic acid flexibility by molecular dynamics. Curr Opin Struct Biol 18: 185-193. doi:10.1016/j.sbi.2008.01.005
- Bomble YJ, Case DA (2008) Multiscale modeling of nucleic acids: Insights into DNA flexibility. Biopolymers 89: 722-731. doi:10.1002/bip.21000
- Beveridge DL, Cheatham TE 3rd, Mezei M (2012) The ABCs of molecular dynamics simulations on B-DNA, circa 2012. J Biosci 37: 379-397. doi:10.1007/s12038-012-9222-6
- Cheatham TE 3rd, Case DA (2013) Twenty-five years of nucleic acid simulations. Biopolymers 99: 969-977. doi:10.1002/bip.22331
- Pasi M, Maddocks JH, Beveridge D, Bishop TC, Case DA, Cheathm T 3rd, et al. (2014) ABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA. Nucl Acids Res 42: 12272-12283. doi: 10.1093/nar/gku855
- Knotts IV TA, Rathore N, Schwartz DC, de Pablo JJ (2007) A coarse grain model for DNA. J Chem Phys 126: 084901. doi: 10.1063/1.2431804
- Morriss-Andrews A, Rottler J, Plotkin, SS (2010). A systematically coarse-grained model for DNA and its predictions for persistence length, stacking, twist, and chirality. J Chem Phys 132: 035105. doi: 10.1063/1.326999
- Dans PD, Zida A, Machado MR, Pantano S (2010) A coarse grained model for atomic-detailed DNA simulations with explicit electrostatics. J Chem Theor Comp 6: 1711-1725. doi:10.1021/ct900653p
- Freeman GS, Hinckley DM, de Pablo JJ (2011) A coarse-grain three-site-per-nucleotide model for DNA with explicit ions. J Chem Phys 135: 165104. doi: 10.1063/1.3652956
- Freeman GS, Lequieu JP, Hinckley DM, Whitmer JK, de Pablo JJ (2014) DNA shape dominates sequence affinity in nucleosome formation. Phys Rev Lett 113: 168101. doi: http://dx.doi.org/10.1103/PhysRevLett.113.168101
- Freeman GS, Hinckley DM, Lequieu JP, Whitmer JK, Pablo JJ (2014) Coarse-grained modeling of DNA curvature. J Chem Phys 141: 165103. doi: 10.1063/1.4897649
- Pinamonti G, Bottaro S, Micheletti C, Bussi G (2015) Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments. Nucleic acids res. gkv708. doi: 10.1093/nar/gkv708
- Savelyev A, Papoian GA (2010) Chemically accurate coarse graining of double-stranded DNA. Proc Nat Acad Sci 107: 20340-20345. doi: 10.1073/pnas.1001163107
- Setny P, Zacharias M (2013) Elastic network models of nucleic acids flexibility. J Chem Theor Comp 9: 5460-5470. doi: 10.1021/ct400814n
- Gonzalez O, Petkeviciute D, Maddocks JHJ (2013) A sequence-dependent rigid-base model of DNA. J Chem Phys 138: 055102. doi: 10.1063/1.4789411
- Naome A, Laaksonen A, Vercauteren DP (2014) A solvent-mediated coarse-grained model of DNA derived with the systematic Newton inversion method. J Chem Theor Comp 10: 3541-3549. doi:10.1021/ct500222s
- Korolev N, Luo D, Lyubartsev AP, Nordensklold L (2014) A coarse-grained DNA model parameterized from atomistic simulations by inverse Monte Carlo. Polymers 6: 1655-1675. doi:10.3390/polym6061655
- Goi JR, Fenollosa C, Perez A, Torrents D, Orozco M (2008) DNAlive: a tool for the physical analysis of DNA at the genomic scale. Bioinformatics 24: 1731-1732. doi: 10.1093/bioinformatics
- Stolz RC, Bishop TC (2010) ICM Web: the interactive chromatin modeling web server. Nucleic Acids Res 38: W254-W261. doi: 10.1093/nar/gkq496
- Chen W, Lin H, Feng PM, Ding C, Zuo YC, Chou KC (2012) iNuc-PhysChem: A sequence-based predictor for identifying nucleosomes via physicochemical properties. PLos One 7: e47843. doi: 10.1371/journal.pone.0047843
- Hospital A, Faustino I, Collepardo-Guevara R, Gonzalez C, Gelpi JL, Orozco M (2013) NAFlex: a web server for the study of nucleic acid flexibility. Nucleic Acids Res 41: W47-W55. doi: 10.1093/nar/gkt378
- Guo SH, Deng EZ, Xu LQ, Ding H, Lin H, Chen W (2014) iNuc-PseKNC: a sequence-based predictor for predicting nucleosome positioning in genomes with pseudo k-tuple nucleotide composition. Bioinformatics 30: 1522-1529. doi: 10.1093/bioinformatics/btu083
- Petkevicit D, Pasi M, Gonzalez O, Maddocks JH (2014) cgDNA: a software package for the prediction of sequence-dependent coarse-grain free energies of B-form DNA. Nucleic Acids Res 42: e153. doi: 10.1093/nar/gku825
- Olson WK, Gorin AA, Lu XJ, Hock LM, Zhurkin VB (1998) DNA sequence-dependent deformability deduced from protein DNA crystal complexes. Proc Nat Acad Sci 95: 11163-11168. doi: 10.1073/pnas.95.19.11163
- Olson WK, Colasanti AV, Li Y, Ge W, Zheng G, Zhurkin VB (2006) DNA simulation benchmarks as revealed by X-ray structures. In: Computational studies of RNA and DNA. Springer Netherlands: 235-257.
- Morozov AV, Fortney K, Gaykalova DA, Studitsky VM, Widom J, Siggia ED (2009) Using DNA mechanics to predict in vitro nucleosome positions and formation energies. Nucleic Acids Res 37: 4707-4722. doi: 10.1093/nar/gkp475
- Perez A, Lankas F, Luque FJ, Orozco M (2008) Towards a molecular dynamics consensus view of B-DNA flexibility. Nucleic Acids Res 36: 2379-2394. doi: 10.1093/nar/gkn082
- Lavery R, Zakrzewska K, Beveridge D, Bishop TC, Case DA, Cheatham T, et al. (2010) A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA. Nucleic Acids Res 38: 299-313. doi: 10.1093/nar/gkp834
- Dans PD, Perez A, Faustino I, Lavery R, Orozco M (2012). Exploring polymorphisms in B-DNA helical conformations. Nucleic Acids Res 40: 10668-10678. doi: 10.1093/nar/gks884
- Lu XJ, Olson WK (2003) 3DNA: a software package for the analysis, rebuilding and visualization of three dimensional nucleic acid structures. Nucleic Acids Res 31:5108-5121. doi: 10.1093/nar/gkg680
- Drew HR, Wing RM, Takano T, Broka C, Tanaka S, Itakura K, et al. (1981) Structure of a B-DNA dodecamer: conformation and dynamics. Proc Natl Acad Sci USA 78: 2179-2183.
- Holbrook SR, Dickerson RE, Kim SH (1985). Anisotropic thermal-parameter refinement of the DNA dodecamer CGCGAATTCGCG by the segmented rigid-body method. Acta Crystallogr B 41:255-262. doi:10.1107/S0108768185002087
- Westhof E (1987). Re-refinement of the B-dodecamer d(CGCGAATTCGCG) with a comparative analysis of the solvent in it and in the Z-hexamer d(5BrCG5BrCG5BrCG). J Biomol Struct Dyn 5:581-600. doi:10.1080/07391102.1987.10506414
- Kneale G,Brown T, Kennard O, Rabinovich D (1985). G.T base-pairs in a DNA helix: the crystal structure of d(G-G-G-G-T-C-C-C). J Mol Biol 186: 805-814. doi:10.1016/0022-2836(85)90398-5
- Ng HL, Kopka ML, Dickerson RE (2000). The structure of a stable intermediate in the A ;-¿ B DNA helix transition. Proc Natl Acad Sci USA 97:2035-2039. doi: 10.1073/pnas.040571197
- Hahn M, Heinemann U (1993). DNA helix structure and refinement algorithm: comparison of models for d(CCAGGCm5CTGG) derived from NUCLSQ, TNT and X-PLOR. Acta Crystallogr D 49: 468-477. doi:10.1107/S0907444993004858
- Sadasivan C, Gautham N (1995). Sequence-dependent microheterogeneity of Z-DNA: the crystal and molecular structures of d(CACGCG).d(CGCGTG) and d(CGCACG).d(CGTGCG). J Mol Biol 248: 918-930. doi:10.1006/jmbi.1995.9894
- Timsit Y, Vilbois E, Moras D (1991). Base-pairing shift in the major groove of (CA)n tracts by B-DNA crystal structures. Nature 354: 167-170. doi:10.1038/354167a0
- Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom AC, Field Y, Moore IK, et al. (2006). A genomic code for nucleosome positioning. Nature 442: 772-778. doi:10.1038/nature04979
- Lee W, Tillo D, Bray N, Morse RH, Davis RW, et al. (2007) A high-resolution atlas of nucleosome occupancy in yeast. Nat Genet 39: 1235-1244. doi:10.1038/ng2117
- Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, et al.(2008) Dynamic regulation of nucleosome positioning in the human genome. Cell 132: 887-898. doi:10.1016/j.cell.2008.02.022
- Mavrich TN, Ioshikhes IP, Venters BJ, Jiang C, Tomsho LP, Qi J, et al. (2008) A barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome. Genome Res 18: 1073-1083. doi: 10.1101/gr.078261.108
- Mavrich TN, Jiang C, Ioshikhes IP, Li X, Venters BJ, Zanton SJ, et al. (2008) Nucleosome organization in the Drosophila genome. Nature 453: 358-362. doi:10.1038/nature06929
- Tillo D, Hughes TR (2009). G+ C content dominates intrinsic nucleosome occupancy. BMC Bioinform 1: 442. doi:10.1186/1471-2105-10-442