Global dynamics of newly constructed oligonucleosomes of conventional and variant H2A.Z histone
© Ramaswamy and Ioshikhes; licensee BioMed Central Ltd. 2007
Received: 23 April 2007
Accepted: 08 November 2007
Published: 08 November 2007
Complexes of nucleosomes, which often occur in the gene promoter areas, are one of the fundamental levels of chromatin organization and thus are important for transcription regulation. Investigating the dynamic structure of a single nucleosome as well as nucleosome complexes is important for understanding transcription within chromatin. In a previous work, we highlighted the influence of histone variants on the functional dynamics of a single nucleosome using normal mode analysis developed by Bahar et al. The present work further analyzes the dynamics of nucleosome complexes (nucleosome oligomers or oligonucleosomes) such as dimer, trimer and tetramer (beads on a string model) with conventional core histones as well as with the H2A.Z histone variant using normal mode analysis.
The global dynamics of oligonucleosomes reveal larger amplitude of motion within the nucleosomes that contain the H2A.Z variant with in-planar and out-of-planar fluctuations as the common mode of relaxation. The docking region of H2A.Z and the L1:L1 interactions between H2A.Z monomers of nucleosome (that are responsible for the highly stable nucleosome containing variant H2A.Z-histone) are highly dynamic throughout the first two dynamic modes.
Dissection of the dynamics of oligonucleosomes discloses in-plane as well as out-of-plane fluctuations as the common mode of relaxation throughout the global motions. The dynamics of individual nucleosomes and the combination of the relaxation mechanisms expressed by the individual nucleosome are quite interesting and highly dependent on the number of nucleosome fragments present in the complexes. Distortions generated by the non-planar dynamics influence the DNA conformation, and hence the histone-DNA interactions significantly alter the dynamics of the DNA. The variant H2A.Z histone is a major source of weaker intra- and inter-molecular correlations resulting in more disordered motions.
In a previous work, we described the influence of histone variants on the functional dynamics of a single nucleosome using normal mode analysis . Knowledge of the dynamics of nucleosome complexes with conventional histones or histone variants would enrich our understanding of the role of chromatin in regulating or modifying gene transcription [2–5]. Significant advances have been made to elucidate the structure of compact nucleosome complexes in specific promoter regions where chromatin structure should be remodeled during transcription.
Variants have been described for several classes of histones, and one of the best-studied examples is the Z variant of H2A (i.e. H2A.Z). H2A.Z is a H2A subtype that has been identified in organisms as diverse as Saccharomyces cerevisiae, Tetrahymena, Drosophila, and Homo sapiens[23–26]. This protein displays 60% homology with H2A and 90% homology between species. Mutagenic assays have demonstrated that H2A.Z is essential for development in yeast  and for survival in Tetrahymena and Drosophila[29, 30]. The nucleosomes containing H2A.Z impact transcription activation within chromatin. Several pioneering studies, a few of which are mentioned here, have reported the importance and the functional diversity of the nucleosomes with H2A.Z. Suto et al. have characterized the histone-DNA and histone-histone interactions within the nucleosome core particle containing the variant H2A.Z histone . Jackson and Gorovsky demonstrated that H2A.Z genes are functionally conserved and distinct from the genes of the major H2A sequence variants . H2A.Z has been observed to be located at yeast promoters and to display a redundant role with ATP-dependent nucleosome remodeling complexes . A review on the functional diversity of histone variants was reported by Brown . Evolutionarily conserved variant histone H2A.Z has recently shown to regulate gene transcription in Saccharomyces cerevisiae. H2A.Z modulates functional interactions with transcription regulatory components, thus helping to poise chromatin for gene activation . H2A.Z alters the nucleosome surface to promote HP1α-mediated chromatin fiber folding . The H2A.Z/H2B dimer is unstable compared to the dimer containing the major H2A isoform . Studies on H2A.Z are expanding due to the importance of the H2A.Z variant histone in both positive and negative gene transcription.
Recent publications show controversial observations of the role of H2A.Z variants. Ausio group showed that histone variant H2A.Z has a destabilizing effect on both intranucleosomal histone-histone interactions and the internuclosomal structure . In contrast, Fan et al. showed that H2A.Z facilitates the intramolecular folding of nucleosome arrays while simultaneously inhibiting the formation of highly condensed structures that are a result of intermolecular association . In addition to these puzzles, a new Fluorescence Resonance Energy Transfer (FRET) approach revealed that the H2AZ variant histone stabilizes the histone octamer within the nucleosome . Continuing in the same line, Placek et al. declared that the H2A.Z/H2B dimer is the least stable histone fold characterized to date, while H2A/H2B appears to be the most stable . On the contrary, a more recent study by the Ausio group reinforced the role of H2A.Z in the stabilization of nucleosome . Recently, Greaves et al. showed that H2AZ is a core component of heterochromatin needed for its integrity and proper sister chromatid cohesion . A few more articles also explain the contrasting functions of H2A.Z and the nucleosome stability [42–45]. In order to identify the structural dynamics of nucleosome with variant histone H2A.Z within the above context, we dissected the global dynamics of the nucleosome complexes with conventional as well as variant histone H2A.Z to rationalize the stabilizing/destabilizing role of H2A.Z histone over the global dynamics of the nucleosome in normal mode analysis.
Traditional molecular simulation by full-atomic force fields [46, 47] is a very time-consuming method for exploring dynamics of the nucleosome (~50,000 atoms). Recently, Sharma et al. described a multiscale modeling of nucleosome dynamics using discrete molecular dynamics . The normal mode analysis (NMA) [49, 50] is an efficient yet physically meaningful approach. It assumes that the molecular potential energy has quadratic shape and that the generalized eigenvalue problem is solved to obtain the analytic description of the motion. The eigenvalues give the frequencies of the vibration whereas the eigenvectors represent the corresponding motion. In addition, positional fluctuations or time correlation functions can be calculated. Because of the difficulty of handling larger biomolecular systems, the simplified NMAs based on the elastic network model [51–56] were developed and successfully applied to some molecular systems [55–58]. The Gaussian Network Model (GNM) [52, 53] and Anisotropic Network Model (ANM) , proposed by Bahar and collaborators, have been extensively used to predict the structure and direction of residues as explained in the Methodology section [59–62].
This article addresses the global dynamics of the nucleosome complexes (oligonucleosomes) with conventional histones as well as variant H2A.Z and is an extension of our previous work on the structural dynamics of the single nucleosome .
Results and Discussion
The global dynamics of the nucleosomal complexes have been studied systematically starting with the analysis of the dimer and extending up to the tetramer. The specified larger cutoff distances account for the complete protein-DNA interactions. Interpretations of global modes are made by comparing the theoretically-generated thermal fluctuations with the experimentally-observed thermal fluctuations of the residues.
Validation of the constructed model of 1EQZ, 1F66 dimers with 1ZBB dimers
The correlation map of dimers H3-H4 and H2A.1-H2B.1 would be different for the two nucleosomes present in 1ZBB if the second nucleosome had influence over the dynamics of the first nucleosome and vice versa. However, it is clear from Fig. 2 and Fig. 3 that the dynamics of the first and second nucleosomes possess similar domain interactions in 1ZBB. The comparative analysis of cross-correlations shows that the constructed dimers 1F66 and 1EQZ have a similar correlation pattern as the nucleosome dimer 1ZBB. This observation reinforces that the model dimer and the experimentally proven conformation express similar inter- as well as intra- domain motions. Therefore, the modeled conformation is similar to the experimentally observed conformation, thus supporting our method of modeling the nucleosome complexes. This approach enabled us to model the complexes such as dimer, trimer and tetramer. Further extension of this approach would provide additional insights about chromatin structure/biology.
Inter- and intra-domain motions of the handshake motif in nucleosome dimer
Fig. 2 and 3 map the correlations between the residue motions of dimers H3-H4 and H2A (H2A.Z in 1F66 and H2A.1 in 1ZBB)-H2B (H2B.1 in 1ZBB) for the three examined structures, labeled I-III. In H3, the α2 helix acts as a global hinge-bending center for the anti-correlated motions expressed between the individual DNA binding loops (L1, L2), where short helices are correlated with the neighboring loops (L1 or L2 respectively). Similarly, H4 possesses two highly anti-correlated domains from N- to C- terminal, where the central residues of the helix α2 play a role of hinge domain. In this handshake motif, the two regions: (i) the loop L2 and α3 of H3 and the loop L1 and α1 of H4, (ii) L1 and α1 of H3 and L2 and α3 of H4, possess a highly correlative interaction within themselves and move in opposite directions. The dynamic N-terminus of H4 plays a vital role in gene silencing . In the presence of variant histone H2A.Z, the interactions in the cooperative and anti-cooperative regions of H3-H4 motif are weakened.
Nearly the entire H2A-H2B motif possesses correlated motions except for the C-terminal of H2A (as it inserts into H3-H4 motif) and the N-terminal of H2B. The presence of H2A.Z variant disrupts the coherent motions between the residues resulting in the anti-correlation of the uncorrelated residues. The presence of H2A.Z introduces significant differences in intra- and inter-molecular correlations. In particular, the presence of H2A.Z decreases the coupling between (i) its docking domain and the helix-loop α1L1, and (ii) the loop-helix L2α3 of H2B (regions highlighted in the magenta boxes). Such disruptive behavior of H2A.Z supports the experimentally observed chromatin-destabilizing role of H2A.Z . The histone H2B also exhibits H2A dependent inter- and intra-correlated domain motions as the cooperation of the motions decreases with the increasing H2A.Z mobility.
In general, the correlations between the motions of the handshake H2A-H2B (or H2A.Z-H2B) motives of the nucleosome dimer are similar to the correlation of H2A-H2B (or H2A.Z-H2B) motives of a single nucleosome and are also the same for H3-H4. With this argument, we have not extended the same analysis on the variant trimer as well as the tetramer. Because our main focus is on the analysis of the influence of H2A.Z over the dynamics, the interactions generated in H2A.1-H2B.1 motif of 1ZBB conformation have not been considered.
Global dynamics of nucleosome complexes
Analysis of the rigid domains of the nucleosome dimer with conventional histones (Fig. 4IA) reveals that almost half of the nucleosome encompassing the regions from the α3 helix and the L2 of H3 that interact with the DNA, linker DNA, and linker histone (regions colored from black to blue) are constrained by dynamics. The N terminal domain of H3, one of the major as well as important sites for transcription modification, is highly dynamic. In the case of the nucleosome with variant H2A.Z (Fig. 4IIA), the rigid domains are defined only by the linker DNA and linker histones. It is observed that the conventional histones restricted to move in the global modes are released for free motion in the presence of the variant H2A.Z histone. Although the dimer formed by H2A-H2B shows pronounced dynamics (colored from yellow to orange), the amplitude of dynamics for the dimer H2A.Z-H2B is even stronger (colored from orange to red) revealing the enhancement in the amplitude of dynamics from the presence of variant H2A.Z histone . This supports the experimental observation of unstable H2A.Z/H2B heterodimer compared to the dimer containing the major H2A isoform . Markedly stronger dynamics due to the presence of H2A.Z may drive the histone proteins to actively participate in the dimer global mode. Therefore, the constrained or frozen regions in the conventional nucleosome are dynamic in the nucleosome with H2A.Z variant.
Consequently, it is clear that the nucleosome dimer within the global mode, undergoes a highly dynamic motion within the regions encompassed by the linker DNA and the linker histone H1 as the hinge domains for the free dynamics of the dimer (Fig. 4). In our previous analysis, the monomeric nucleosome expressed highly symmetric dynamics with respect to the dyad axis in the global mode. In the second slowest mode, the dynamics have been identified as conjugate to that of the first slowest mode of motion, i.e., the central region of the nucleosome perpendicular to the dyad axis is highly constrained (rigid). When we compare this observation with the current observation of the nucleosomes in the dimer, it is very interesting to note that the global dynamics of the two single nucleosomes does not reflect the dynamics pattern observed for the monomeric nucleosome in its global mode. Here, in the global mode, the nucleosomes are highly dynamic while keeping the linker DNA/histone as the rigid domains.
The dynamics of the nucleosome dimer in the second slowest mode is comparatively constrained along the axis conjugate to the nucleosome dyad (Fig. 4B). Unlike the first global mode, the linker DNA and the linker histone are highly dynamic (colored red), and the regions of the H2A-H2B dimer interacting with the peripheral DNA sequences (Fig. 4IIB) are also in motion. In the case of nucleosome with H2A.Z variant, the DNA-(H2A.Z-H2B) interacting regions (α2-L2-α3-αC regions of H2B and α2-L1-α1 of H2A.Z) express stronger dynamics along with the interacting DNA segment than that of histone H2A. Furthermore, the dynamics of linker DNA/histone are more pronounced in the nucleosome dimer with H2A.Z. The vital observation is, in the second slowest mode, the nucleosome dimer reinforces the conjugate dynamics as that of the second global mode of the single nucleosome established by Ramaswamy et. al. . When compared to the relaxation of the nucleosome dimer in the slowest global mode, the nucleosome dimer with conventional histones is almost stationary in the second slowest mode.
The dynamics of a nucleosome with conventional and variant histone are comparable, but the amplitude of dynamics is significantly higher for the dimer of variant nucleosome than that of the conventional nucleosome dimer (Fig. 4). The rigid domains observed in the second slowest mode are wide spread compared to the second slowest mode of nucleosome with variant H2A.Z histone (Fig. 4B). The variation in the amplitude of nucleosomal dynamics with the presence of H2A.Z variant histone is elaborated in our previous study . The distinct behavior of linker histone in both slow modes may reinforce the observation of its ability to control the mobility of nucleosomes on DNA .
Whereas residues from N-terminal to the middle of α2-helix exist in a highly dynamic domain in the second slowest mode (Fig. 5.II), the rest of the histone fold expresses moderate dynamics (colored green). In comparison, the residues from N-terminal to the middle of α2-helix and the αC-helix of H2A.Z possess larger amplitude of dynamics than that of H2A (changed from orange to red color). In general, the docking region of H2A.Z (and H2A) is highly dynamic throughout the first two slow modes. Concerning the stability of the octameric histone core, it has been suggested that the L1:L1 interactions between H2A.Z monomers in the nucleosome may be responsible for the greater stability of H2A.Z-containing nucleosomes [14, 39]. From Fig. 5, it is obvious that the L1:L1 interactions between the two H2A.Z monomers would be highly fluctuating due to the participation of these loops in the global dynamics. The active participation of H2A.Z in the global dynamics is assumed to play a role in the transcription activation within chromatin and also reinforces the fact that H2A.Z is preferentially associated with actively transcribed chromatin .
The scale of our analysis is limited to only the small nucleosome complexes (up to tetramers) and is too short to explain whether H2A.Z facilitates intra nucleosome-nucleosome interactions to promote the formation of the 30nm fiber. From all the present observations, it can be stated that the incorporation of H2A.Z histone releases inter-residual interactions and thereby enhances the amplitude of dynamics in the global mode. We also reported the enhancement of histone-DNA interactions due to the incorporation of H2A.Z in our previous study. Such pronounced dynamics might be needed for nucleosome positioning  required to establish a specific active chromatin architecture at the promoter regions of certain genes . It doesn't answer the question whether the H2A.Z generates a completely dynamic octameric core (leading to the local destabilization ) or forces the nucleosome into the equilibrium position (leading to the global stabilization ) via the incorporation of H2A.Z.
Directionality in the global dynamics of nucleosome complexes
Figure 8IA shows the conformational deviations observed during the global mode of the conventional nucleosome dimer. The r.m.s. value between the crystal structure and new ANM conformation were calculated by using the standard command in Chimera by fitting all atoms. The calculated r.m.s. value is 1.738 Å, meaning that both conformations are nearly identical. Superimposition of the crystal and ANM conformations showed that the nucleosome dimer expresses a conformational fluctuation perpendicular to the dyadic plane. Almost half of the nucleosomes encompassing the regions such as H3-H4 dimer and the linker DNA/histone is pushed below the crystal conformation (and hence we could see the crystal conformation colored as grey) and the other half of the nucleosome is pushed above the crystal conformation (and hence we could see the magenta ANM conformation). This leads to the bending of the nucleosome dimer perpendicular to the dyad, which is the binding site for the transcription factors . Even though the GNM derived conformations revealed a highly dynamic conformation in the global mode, the amplitude of fluctuation is considerably lower for the conventional nucleosome and for this reason, no indications are given for the residual fluctuations.
Figure 8IB shows the conformational deviations observed in the global mode of variant nucleosome dimer. The H-bonding network in the dynamics of H2A.Z has been compared to the crystal structure and the ANM conformations. During the dynamics between the two ANM conformations, the two hydrogen bonds; (i) Arg 31--- G267 and (ii) Lys79 --- C278, that exist in the crystal strutcure of H2A.Z, disappear and additional H-bonds are formed. The r.m.s value between the crystal structure and new ANM conformation is 3.285 Å indicating larger deviations in the variant nucleosome and, hence, is more dynamic than the conventional one. The variant nucleosome dimer also reveals out-of-plane dynamics in the dyadic plane, yet the axis of bending is not exactly perpendicular to the dyadic axis, but is inclined ~45° to the dyad. In addition to the bending, the variant dimer also expresses notable in-plane distortions as indicated by the arrows. Such in-plane and out-of-plane motions might generate distortions in histone-DNA interactions which would weaken the stability of the nucleosome and result in nucleosome dissociation . Altogether, the in-plane as well as out-of-plane motion of the conventional/variant nucleosome reveals vibrant dynamics for the variant nucleosome and is shown with linear and curved arrows. Further analysis on nucleosome complexes such as trimer and tetramer has been focused only on the variant nucleosome since the dimer has already proven the effective dynamics/influence of variant histone H2A.Z like for the variant monomeric nucleosome.
During the dynamics of all these complexes, the interaction between the histone octamer and duplex DNA is highly altered and is expected to accompany a change in the position of the octameric histone, with respect to the DNA, playing a main role in nucleosome sliding [21, 69, 70]. Through the dynamics, the region encompassed by the linker histone and DNA are not very flexible – possibly due to the binding of the linker histone and linker DNA region. Incorporation of histone H1 into nucleosomal complexes results in essentially inactive chromatin templates with respect to both transcription elongation and initiation .
The normal mode analysis of the global dynamics of nucleosomal complexes such as dimer, trimer and tetramer with conventional as well as variant H2A.Z histones has been studied using normal mode analysis. Dissection of the dynamics of individual nucleosomes in the complex discloses in-planar and out-of-planar fluctuations as the common mode of relaxation throughout the global motions, but the loci of the axis characterizing out-of-plane motions deviate irrespective of the number of nucleosomes present in the complex. Distortions generated by the non-planar dynamics also influence the DNA sequences and hence the histone-DNA interactions significantly alter the dynamics of the DNA, conversely influencing the transcriptional activity of the nucleosome as well as chromatin .
The presence of variant histones highly influences the dynamics of the nucleosome complexes. Comparison of the dynamics of variant H2A.Z with conventional H2A histone reveals that the complete histone fold of H2A and H2A.Z is dynamic in the global mode with enhanced dynamics at (i) the αC-helix of H2A.Z, and (ii) the residues encompassing the L1:L1 interaction regions between H2A.Z monomers. However, only half of the histone fold (the residues from N-terminal to the middle of α2-helix) exists as the highly dynamic domain in the second slowest mode. The highly dynamic/active variant H2A.Z is assumed to play a part in the transcription activation within chromatin. Hence, it can be concluded that the highly dynamic H2A.Z histone is a major source of weaker intramolecular and intermolecular couplings/correlations [1, 36].
The pronounced fluctuation of H2A.Z that leads to the enhanced dynamics of the nucleosome complexes is also expected to influence the folding of H2A.Z reconstituted chromatin unlike that of the H2A chromatin . Experimentally, the replacement of H2A.Z for H2A histone in cells induces a shift in the positioning of a nucleosome located over the transcription start site . However, it is essential to determine whether the simple incorporation of H2A.Z into the nucleosome is sufficient to facilitate gene activity. The higher dynamics of the nucleosome complex is assumed to participate in the disruption of the nucleosome due to the higher mobility and the disorder of the variants. Hence, in vitro-assembled chromatin complexes containing H2A.Z are more resistant to condensation when compared to H2A-containing nucleosome complexes [43, 44]. The enhanced dynamics of H2A.Z single and oligo-nucleosomes is in agreement with their positioning in the gene promoter area and their related role in gene regulation [44, 72]. Nucleosome positioning correlates with the positioning of other promoter elements (core promoter elements and binding sites of specific transcription factors) [73–76]. Promoter chromatin structure should be remodeled to facilitate the transcription. The dynamic nature of the H2A.Z nucleosomes makes them more susceptible to this remodeling. The consistency in the dynamics pattern of the nucleosome complexes observed in the present study highlights the existence of nucleosome in various relaxation modes and the role of variant histones over the global dynamics.
Even though the nucleosome expresses the two major relaxation mechanisms [(i) in-planar motions and (ii) out-of-planar bending motions], it is very interesting to show how the individual nucleosome selects/adopts the relaxation mechanism when the dynamics of the whole nucleosome complex is considered. The choices/combination of the relaxation mechanisms adopted by the individual nucleosomes in the complex are also rather new and these observations could enrich the knowledge of the relaxation modes of nucleosome arrays (nucleosome complexes having more than six nucleosomes) when the analysis is extended beyond the order of the nucleosome complexes.
Coordinates of the linker histone H1 were taken from the monomer structure of the PDB: 1GHC. A linear B-DNA fragment of length 19bp was constructed using the Biopolymer module of the Insight II package. A complete nucleosome monomer with which the nucleosome complex was modeled was built using a nucleosome of the PDB database: 1F66 for variant nucleosome complex and 1EQZ for normal nucleosome complex. At first, two B-DNA fragments of 10 bp (base pairs) with the same sequence on both terminus of the nucleosome were modeled using the biopolymer module of Insight II and superimposed on each DNA terminus of the nucleosome. After superimposing, the newly built 10 bp were connected at the arms of the nucleosome crystal structures. The existing 10mers were removed, as these sequences possess deformed base stacking and may not extend without steric interaction for any order of nucleosome polymerization. Secondly, a B-DNA of sequence 5'-CTGCAGATTCTACCAAAAG-3' was modeled and extended at both ends to form the linker DNA region. At this stage, the nucleosome can interact with the linker histone H1 with its helix bundle. The linker histone was then located, so it would make the hydrogen bonds with the duplex of the DNA [16, 17]. After constructing the monomer, the latter was allowed for structural minimization by using the steepest descent method to attain the globally stable conformation. After optimization, the monomer was used to model the nucleosomal complexes as the "beads on a string" model. The configuration of nucleosome complexes of regular H2A histone and the variant H2A.Z histone are modeled using the crystal structures 1EQZ and 1F66 respectively. The nucleosome 1F66 corresponds to the recombinant mouse H2A.Z and recombinant Xenopus leavis H2B, H3 and H4, whereas 1EQZ refers to the chicken (Gallus gallus) histone octamer. Despite the differences in the originating organisms, the nucleosomes possess high (>95%) sequence identity with the exception of the variant histone H2A.Z with the r.m.s. deviations between α-carbon coordinates of ~0.50 Å (see  for more details). All the nucleosome complexes were modeled using the molecular modeling package Insight II. Recently, Schalch et al.  reported the crystal structure of a tetranucleosome of Xenopus leavis origin at 9 Å resolution as solved by molecular replacement using the nucleosome core structure (PDB id: 1ZBB). We used this conformation to validate the modeled complexes and to exhibit a comparison between inter- as well as intra- domain motions of the model dimer (see Results and Discussion for examples).
In the Gaussian network model (GNM), the elastic network is constructed by considering the α-Carbon atoms of each residue as nodes connected by the elastic spring constant γ. The interaction sites within a spherical shell of r c are considered coupled harmonic oscillators subjected to Gaussian fluctuations around their mean position. The Hamiltonian of the conformation is derived from the Kirchhoff matrix of contacts (Γ) by considering local packing density (pairs of amino acid or nucleotide residues within the cutoff distance, r c ) and coordination order (the separation of the contacting residues along the backbone). In proteins, the Cα network is considered with r c ≤ 10 Å. DNA, as well as RNA, has a slightly increased cutoff distance ~14 ± 2 Å to incorporate inter-strand interactions [79–81]. The Kirchhoff matrix, Γ, is an N × N symmetric matrix for a sequence of N residues with values of -1 for contacting i and j residues, and 0 for non-contacting residues. The diagonal elements are obtained using the expression where j ≠ i. In GNM, motions are expressed as a superposition of modes of different frequencies. The coarse grained dynamics of the nucleosome complex is decomposed into a collection of N-1 internal modes with frequencies λk, 2 ≤ k ≤ N.
An effective elastic network for the nucleosome dimer is generated using a cutoff distance r c ≤ 10 Å for protein, 16 Å for the DNA and the interaction between DNA and protein. In the case of the trimer, the interaction between each pair of protein residues uses the distance r c ≤ 14 Å and 20 Å for the DNA and protein-DNA interactions. For the tetramers, the cutoff distances 40 Å and 60 Å are used for proteins, DNA, and protein-DNA interactions.
ANM is an extension of GNM, where the anisotropic fluctuations from the second derivative of harmonic potential of residues are considered. Fluctuation vectors (ΔRi) in the Cartesian coordinate frame  are generated. The three components of the fluctuation vector help to find the action of individual modes in alternative conformations of the biomolecular complex. This is done by incorporating the vector ± ΔRi with the equilibrium vectors. For this analysis, the networks for 1F66 dimer and trimer were constructed using the cutoff distances 30 Å, 35 Å and 40 Å for the protein, DNA, and protein-DNA respectively. For the 1EQZ dimer, the cutoff distance is 35 Å for protein, 40 Å for DNA, and 45 for DNA-protein interaction. The network for 1F66 tetramer was built using a cutoff distance of 40 Å for protein, 45 Å for DNA, and 50 Å to include the DNA-protein interaction. All the resulting matrices (Hessian matrix) of size 3N × 3N were solved for the eigenvalue and eigenvectors for analysis. Given that the slowest modes refer to the most cooperative motions between entire subunits or subdomains [82–84], focus is only paid to the slowest modes for both GNM and ANM analysis.
The calculations were performed using a desktop PC for ANM/GNM modeling of single nucleosomes and dinucleosomes, and a Cray X1 of the Ohio Supercomputer Center for tri- and tetranucleosomes. Supplemental modeling performed using SGI Octane 2 workstation. From computational time point of view, the eigenvalue decomposition of the connectivity matrix Γ is the most expensive task in GNM calculations. The singular value decomposition (SVD) subroutine  used in GNM scales up the computing time with N3 for a network of N residues. For N < 2000, the computations are performed within minutes.
The authors thank Dr. Li Wang for assisting with the calculations, and Molly Ray, Judith O'Donnell and Tamara Freitas for proofreading and editing the text.
- Ramaswamy A, Bahar I, Ioshikhes I: Structural Dynamics of Nucleosome Core Particle: Comparison with Nucleosomes Containing Histone Variants. Proteins 2005, 58: 683–696. 10.1002/prot.20357View ArticlePubMedGoogle Scholar
- Weintraub H: Histone H1-dependent chromatin superstructures and the suppression of gene activity. Cell 1984, 38: 17–27. 10.1016/0092-8674(84)90522-1View ArticlePubMedGoogle Scholar
- Andersson K, Bjorkroth B, Daneholt B: Packing of a specific gene into higher order structures following repression of RNA synthesis. J Cell Biol 1984, 98: 1296–1303. 10.1083/jcb.98.4.1296View ArticlePubMedGoogle Scholar
- van Holde KE: Chromatin. New York, Springer; 1989.View ArticleGoogle Scholar
- Wolffe A: Chromatin: Structure and Function. San Diego, Academic Press; 1998.Google Scholar
- Luger K, Mader AW, Richmond RK, Sargent DF, Richmond TJ: Crystal structure of the nucleosome core particle at 2.8 Å resolution. Nature 1997, 389: 251–260. 10.1038/38444View ArticlePubMedGoogle Scholar
- Luger K, Rechsteiner T, Richmond TJ: Expression and purification of recombinant histones, and nucleosome reconstitution. Methods Enzymol 1999, 304: 3–19.View ArticlePubMedGoogle Scholar
- Richmond TJ, Finch JT, Rushton B, Rhodes D, Klug A: Structure of the nucleosome core particle at 7 Å resolution. Nature 1984, 311: 532–537. 10.1038/311532a0View ArticlePubMedGoogle Scholar
- Arents G, Burlingame RW, Wang BC, Love WE, Moudrianakis EN: The nucleosomal core histone octamer at 3.1 Å resolution: a tripartite protein assembly and a left-handed super helix. Proc Natl Acad Sci USA 1991, 88: 10148–10152. 10.1073/pnas.88.22.10148PubMed CentralView ArticlePubMedGoogle Scholar
- Eickbush TH, Moudrianakis EN: The histone core complex: An Octamer assembled by two sets of protein-protein interactions. Biochemistry 1978, 17: 4955–4964. 10.1021/bi00616a016View ArticlePubMedGoogle Scholar
- Hayes JJ, Clark DJ, Wolffe AP: Histone contributions to the structure of DNA in the nucleosome. Proc Natl Acad Sci USA 1991, 88: 6829–6833. 10.1073/pnas.88.15.6829PubMed CentralView ArticlePubMedGoogle Scholar
- Hayes JJ, Tullis TD, Wolffe AP: The structure of DNA in a nucleosome. Proc Natl Acad Sci USA 1990, 87: 7405–7409. 10.1073/pnas.87.19.7405PubMed CentralView ArticlePubMedGoogle Scholar
- Davie JR: Covalent modifications of histones: expression from chromatin templates. Curr Opin Genet Dev 1998, 8: 173–178. 10.1016/S0959-437X(98)80138-XView ArticlePubMedGoogle Scholar
- Suto RK, Clarkson MJ, Tremethick DJ, Luger K: Crystal structure of a nucleosome core particle containing the variant histone H2A.Z. Nat Struct Biol 2000, 7: 1121–1124. 10.1038/81971View ArticlePubMedGoogle Scholar
- Thomas JO: Histone H1: location and role. Curr Opin Cell Biol 1999, 11: 312–317. 10.1016/S0955-0674(99)80042-8View ArticlePubMedGoogle Scholar
- Thomas JO, Rees C, Finch JT: Cooperative binding of the globular domains of histones H1 and H5 to DNA. Nucl Acids Res 1992, 20: 187–194. 10.1093/nar/20.2.187PubMed CentralView ArticlePubMedGoogle Scholar
- Draves PH, Lowary PT, Widom J: Cooperative binding of the globular domain of histone H5 to DNA. J Mol Biol 1992, 225: 1105–1121. 10.1016/0022-2836(92)90108-VView ArticlePubMedGoogle Scholar
- Richmond TJ, Widom J: Nucleosome and chromatin structure. In Chromatin Structure and Gene Expression. Edited by: Elgin SCR, Workman JL. Oxford: Oxford Univ. Press; 2000:1–23.Google Scholar
- Horn PJ, Peterson CL: Chromatin higher order folding: Wrapping up transcription. Science 2002, 297: 1824–1827. 10.1126/science.1074200View ArticlePubMedGoogle Scholar
- Mohd-Sarip A, Verrijzer CP: A higher order of silence. Science 2004, 306: 1484–1485. 10.1126/science.1106767View ArticlePubMedGoogle Scholar
- Park YJ, Chodaparambil JV, Bao Y, McBryant SJ, Luger K: Nucleosome assembly protein 1 exchanges histone H2A-H2B dimers and assists nucleosome sliding. J Biol Chem 2005, 280: 1817–1825. 10.1074/jbc.M411347200View ArticlePubMedGoogle Scholar
- Saeki H, Ohsumi K, Aihara H, Ito T, Hirose S, Ura K, Kaneda Y: Linker histone variants control chromatin dynamics during early embryogenesis. Proc Natl Acad Sci USA 2005, 102: 5697–5702. 10.1073/pnas.0409824102PubMed CentralView ArticlePubMedGoogle Scholar
- Jackson JD, Falciano VT, Gorovsky MA: A likely histone H2A.F/Z variant in Saccharomyces cerevisiae . Trends Biochem Sci 1996, 21: 466–467. 10.1016/S0968-0004(96)20028-3View ArticlePubMedGoogle Scholar
- Allis CD, Glover CV, Bowen JK, Gorovsky MA: Histone variants in the transcriptionally active, amitotically dividing nucleus of a unicellular eukaryote. Cell 1980, 20: 609–617. 10.1016/0092-8674(80)90307-4View ArticlePubMedGoogle Scholar
- van Daal A, White EM, Gorovsky MA, Elgin SC: Drosophila has a single copy of the gene encoding a highly conserved histone H2A variant of the H2A.F/2 type. Nucleic Acids Res 1988, 16: 7487–7497. 10.1093/nar/16.15.7487PubMed CentralView ArticlePubMedGoogle Scholar
- Hatch CL, Bonner WM: Sequence of cDNAs for mammalian H2A.Z, an evolutionarily diverged but highly conserved basal histone H2A isoprotein species. Nucleic Acids Res 1988, 16: 1113–1124. 10.1093/nar/16.3.1113PubMed CentralView ArticlePubMedGoogle Scholar
- Jackson JD, Gorovsky MA: Histone H2A.Z has a conserved function that is distinct from that of the major H2A sequence variants. Nucleic Acids Res 2000, 28: 3811–3816. 10.1093/nar/28.19.3811PubMed CentralView ArticlePubMedGoogle Scholar
- Liu X, Li B, Gorovsky MA: Essential and non-essentialhistone 2A variants in Tetrahymena thermophila . Mol Cell Biol 1996, 16: 4305–4311.PubMed CentralPubMedGoogle Scholar
- van Daal A, Elgin SC: A histone variant, H2AvD, is essential in Drosophila melanogaster . Mol Biol Cell 1992, 3: 593–602.PubMed CentralView ArticlePubMedGoogle Scholar
- Clarkson MJ, Wells JRE, Gibson F, Saint R, Tremethick DJ: Regions of variant histone His2AvD required for Drosophila development. Nature 1999, 399: 694–697. 10.1038/21436View ArticlePubMedGoogle Scholar
- Santisteban MS, Kalashnikova T, Smith MM: Histone H2A.Z regulates transcription and is partially redundant with nucleosome remodeling complexes. Cell 2000, 103: 411–422. 10.1016/S0092-8674(00)00133-1View ArticlePubMedGoogle Scholar
- Brown DT: Histone variants: are they functionally heterogeneous? Genome Biol 2001, 2: reviews 0006.1–0006.6.Google Scholar
- Adam M, Robert F, Larochelle M, Gaudreau L: H2A.Z is required for global chromatin integrity and for recruitment of RNA polymerase II under specific conditions. Mol Cell Biol 2001, 21: 6270–6279. 10.1128/MCB.21.18.6270-6279.2001PubMed CentralView ArticlePubMedGoogle Scholar
- Larochelle M, Gaudreau L: H2A.Z has a function reminiscent of an activator required for preferential binding to intergenic DNA. EMBO Journal 2003, 22: 4512–4522. 10.1093/emboj/cdg427PubMed CentralView ArticlePubMedGoogle Scholar
- Fan JY, Rangasamy D, Luger K, Tremethick DJ: H2A.Z Alters the Nucleosome Surface to Promote HP1alpha-Mediated Chromatin Fiber Folding. Mol Cell 2004, 16: 655–661. 10.1016/j.molcel.2004.10.023View ArticlePubMedGoogle Scholar
- Placek BJ, Harrison LN, Villers BM, Gloss LM: The H2A.Z/H2B dimer is unstable compared to the dimer containing the major H2A isoform. Protein Sci 2005, 14: 514–522. 10.1110/ps.041026405PubMed CentralView ArticlePubMedGoogle Scholar
- Abbott DW, Ivanova VS, Wang X, Bonner WM, Ausio J: Characterization of the Stability and Folding of H2A.Z Chromatin Particles, implications for transcriptional activation. J Biol Chem 2001, 276: 41945–41949. 10.1074/jbc.M108217200View ArticlePubMedGoogle Scholar
- Fan JY, Gordon F, Luger K, Hansen JC, Tremethick DJ: The essential histone variant H2A.Z regulates the equilibrium between different chromatin conformational states. Nat Struct Biol 2002, 9: 172–176. 10.1038/nsb0402-316bView ArticlePubMedGoogle Scholar
- Park YJ, Dyer PN, Tremethick DJ, Luger K: A New Fluorescence Resonance Energy Transfer Approach Demonstrates That the Histone Variant H2AZ Stabilizes the Histone Octamer within the Nucleosome. J Biol Chem 2004, 279: 24274–24282. 10.1074/jbc.M313152200View ArticlePubMedGoogle Scholar
- Thambirajah AA, Dryhurst D, Ishibashi T, Li A, Maffey AH, Ausio J: H2A.Z Stabilizes Chromatin in a Way That Is Dependent on Core Histone Acetylation. J Biol Chem 2006, 281: 20036–20044. 10.1074/jbc.M601975200View ArticlePubMedGoogle Scholar
- Greaves IK, Rangasamy D, Ridgway P, Tremethick DJ: H2A.Z contributes to the unique 3D structure of the centromere. Proc Natl Acad Sci USA 2007, 104: 525–530. 10.1073/pnas.0607870104PubMed CentralView ArticlePubMedGoogle Scholar
- Marino-Ramirez L, Kann MG, Shoemaker BA, Landsman D: Histone structure and nucleosome stability. Expert Rev Proteomics 2005, 2: 719–729. 10.1586/147894126.96.36.1999PubMed CentralView ArticlePubMedGoogle Scholar
- Guillemette B, Gaudreau L: Reuniting the contrasting functions of H2A.Z. Biochem Cell Biol 2006, 84: 528–535. 10.1139/O06-077View ArticlePubMedGoogle Scholar
- Guillemette B, Bataille AR, Ge'vry N, Adam M, Blanchette M, Robert F, Gaudreau L: Variant Histone H2A.Z Is Globally Localized to the Promoters of Inactive Yeast Genes and Regulates Nucleosome Positioning. PLoS Biol 2005, 3: 2100–2110. 10.1371/journal.pbio.0030384View ArticleGoogle Scholar
- Dryhurst D, Thambirajah AA, Ausió J: New twists on H2A.Z: a histone variant with a controversial structural and functional past. Biochem Cell Biol 2004, 82: 490–497. 10.1139/o04-043View ArticlePubMedGoogle Scholar
- Bishop TC, Zhmudsky OO: Mechanical model of the nucleosome and chromatin. J Biomol Struct Dyn 2002, 19: 877–887.View ArticlePubMedGoogle Scholar
- Bishop TC: Molecular dynamics simulations of a nucleosome and free DNA. J Biomol Struct Dyn 2005, 22: 673–686.View ArticlePubMedGoogle Scholar
- Sharma S, Ding F, Dokholyan NV: Multiscale modeling of nucleosome dynamics. Biophys J 2007, 92: 1457–1470. 10.1529/biophysj.106.094805PubMed CentralView ArticlePubMedGoogle Scholar
- Levitt M, Sander C, Stern PS: Protein normal-mode dynamics: trypsin inhibitor, crambin, ribonuclease and lysozyme. J Mol Biol 1985, 181: 423–447. 10.1016/0022-2836(85)90230-XView ArticlePubMedGoogle Scholar
- Gibrat JF, Go N: Normal mode analysis of human lysozyme: study of the relative motion of the two domains and characterization of the harmonic motion. Proteins: Struct Funct Genet 1990, 8: 258–279. 10.1002/prot.340080308View ArticleGoogle Scholar
- Tirion MM: Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. Phys Rev Lett 1996, 77: 1905–1908. 10.1103/PhysRevLett.77.1905View ArticlePubMedGoogle Scholar
- Bahar I, Atilgan AR, Erman B: Direct evaluation of thermal fluctuations in proteins using a single- parameter harmonic potential. Fold Des 1997, 2: 173–181. 10.1016/S1359-0278(97)00024-2View ArticlePubMedGoogle Scholar
- Haliloglu T, Bahar I, Erman B: Gaussian dynamics of folded proteins. Phys Rev Lett 1997, 79: 3090–3093. 10.1103/PhysRevLett.79.3090View ArticleGoogle Scholar
- Hinsen K: Analysis of domain motions by approximate normal mode calculations. Proteins 1998, 33: 417–429. 10.1002/(SICI)1097-0134(19981115)33:3<417::AID-PROT10>3.0.CO;2-8View ArticlePubMedGoogle Scholar
- Tama F, Sanejouand YH: Conformational change of proteins arising from normal mode calculations. Prot Eng 2001, 14: 1–6. 10.1093/protein/14.1.1View ArticleGoogle Scholar
- Atilgan AR, Durrell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I: Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J 2001, 80: 505–515.PubMed CentralView ArticlePubMedGoogle Scholar
- Bahar I, Wallqvist A, Covell DG, Jernigan RL: Correlation between native-state hydrogen exchange and cooperative residue fluctuations from a simple model. Biochemistry 1998, 37: 1067–1075. 10.1021/bi9720641View ArticlePubMedGoogle Scholar
- Haliloglu T, Bahar I: Structure-based analysis of protein dynamics: comparison of theoretical results for hen lysozyme with X-ray diffraction and NMR relaxation data. Proteins 1999, 37: 654–667. 10.1002/(SICI)1097-0134(19991201)37:4<654::AID-PROT15>3.0.CO;2-JView ArticlePubMedGoogle Scholar
- Kundu S, Melton JS, Sorensen DC, Phillips GN Jr: Dynamics of Proteins in Crystals: Comparison of Experiment with Simple Models. Biophys J 2002, 83: 723–732.PubMed CentralView ArticlePubMedGoogle Scholar
- Ming D, Kong Y, Lambert MA, Huang Z, Ma J: How to describe protein motion without amino acid sequence and atomic coordinates. Proc Natl Acad Sci USA 2002, 99: 8620–8625. 10.1073/pnas.082148899PubMed CentralView ArticlePubMedGoogle Scholar
- Tama F, Wriggers W, Brooks CL III: Exploring global distortions of biological macromolecules and assemblies from low-resolution structural information and continuum elastic network theory. J Mol Biol 2002, 321: 297–305. 10.1016/S0022-2836(02)00627-7View ArticlePubMedGoogle Scholar
- Tama F, Valle M, Frank J, Brooks CL III: Dynamic reorganization of the functionally active ribosome explored by normal mode analysis and cryo-electron microscopy. Proc Natl Acad Sci USA 2003, 100: 9319–9323. 10.1073/pnas.1632476100PubMed CentralView ArticlePubMedGoogle Scholar
- Schalch T, Duda S, Sargent DF, Richmond TJ: X-ray structure of a tetranucleosome and its implications for the chromatin fibre. Nature 2005, 436: 138–141. 10.1038/nature03686View ArticlePubMedGoogle Scholar
- Johnson LM, Fisher Adams G, Grunstein M: Identification of a non-basic domain in the histone N-terminus required foe repression of the yeast silent mating loci. EMBO J 1992, 11: 2201–2209.PubMed CentralPubMedGoogle Scholar
- Pennings S, Meersseman G, Bradbury EM: Linker histones H1 and H5 prevent the mobility of positioned nucleosomes. Proc Nat Acad Sci USA 1994, 91: 10275–10279. 10.1073/pnas.91.22.10275PubMed CentralView ArticlePubMedGoogle Scholar
- White CL, Luger K: Defined structural changes occur in a nucleosome upon Amt1 transcription factor binding. J Mol Biol 2004, 342: 1391–1402. 10.1016/j.jmb.2004.07.080View ArticlePubMedGoogle Scholar
- Muthurajan UM, Park YJ, Edayathumangalam RS, Suto RK, Chakravarthy S, Dyer PN, Luger K: Structure and dynamics of nucleosomal DNA. Biopolymers 2003, 68: 547–556. 10.1002/bip.10317View ArticlePubMedGoogle Scholar
- Roccatano D, Barthel A, Zacharias M: Structural Flexibility of the Nucleosome Core Particle at Atomic Resolution Studied by Molecular Dynamics Simulation. Biopolymers 2007, 85: 407–427. 10.1002/bip.20690View ArticlePubMedGoogle Scholar
- Luger K: Dynamic nucleosomes. Chromosome Res 2006, 14: 5–16. 10.1007/s10577-005-1026-1View ArticlePubMedGoogle Scholar
- Smith CL, Peterson CL: ATP-dependent chromatin remodeling. Curr Top Dev Biol 2005, 65: 115–148.View ArticlePubMedGoogle Scholar
- O'Neill TE, Meersseman G, Pennings S, Bradbury EM: Deposition of histone Hi onto reconstituted nucleosome arrays inhibits both initiation and elongation of transcripts by 17 RNA polymerase. Nucleic Acids Res 1995, 23: 1075–1082. 10.1093/nar/23.6.1075PubMed CentralView ArticlePubMedGoogle Scholar
- Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom A, Field Y, Moore IK, Wang JZ, Widom J: A genomic code for nucleosome positioning. Nature 2006, 442: 772–778. 10.1038/nature04979PubMed CentralView ArticlePubMedGoogle Scholar
- Ioshikhes I, Trifonov EN, Zhang MQ: Periodical distribution of transcription factor sites in promoter regions and connection with chromatin structure. Proc Natl Acad Sci USA 1999, 96: 2891–2895. 10.1073/pnas.96.6.2891PubMed CentralView ArticlePubMedGoogle Scholar
- Gershenzon NI, Trifonov EN, Ioshikhes IP: The features of Drosophila core promoters revealed by statistical analysis. BMC Genomics 2006, 7: 161. 10.1186/1471-2164-7-161PubMed CentralView ArticlePubMedGoogle Scholar
- Ioshikhes IP, Albert I, Zanton SJ, Pugh BF: Nucleosome positions predicted through comparative genomics. Nat Genet 2006, 38: 1210–1215. 10.1038/ng1878View ArticlePubMedGoogle Scholar
- Albert I, Mavrich TN, Tomsho LP, Qi J, Zanton SJ, Schuster SC, Pugh BF: Translational and rotational settings of H2A.Z nucleosomes across the Saccharomyces cerevisiae genome. Nature 2007, 446: 572–576. 10.1038/nature05632View ArticlePubMedGoogle Scholar
- Cerf C, Lippens G, Ramakrishnan V, Muyldermans S, Segers A, Wyns L, Wodak SJ, Hallenga K: Homo- and heteronuclear two-dimensional NMR studies of the globular domain of histone H1: full assignment, tertiary structure, and comparison with the globular domain of histone H5. Biochemistry 1994, 33: 11079–11086. 10.1021/bi00203a004View ArticlePubMedGoogle Scholar
- Harp JM, Hanson BL, Timm DE, Bunick GJ: Asymmetries in the nucleosome core particle at 2.5 Å resolution. Acta Cryst 2000, D56: 1513–1534.Google Scholar
- Bahar I, Jernigan RL: Vibrational dynamics of transfer RNAs: comparison of the free and synthetase-bound forms. J Mol Biol 1998, 281: 871–884. 10.1006/jmbi.1998.1978View ArticlePubMedGoogle Scholar
- Bahar I, Erman B, Jernigan RL, Atilgan AR, Covell DG: Collective motions in HIV-1 reverse transcriptase: examination of flexibility and enzyme function. J Mol Biol 1999, 285: 1023–1037. 10.1006/jmbi.1998.2371View ArticlePubMedGoogle Scholar
- Temiz AN, Bahar I: Inhibitor Binding Alters the Directions of Domain Motions in HIV-1 Reverse Transcriptase. Struct Funct Genet 2002, 49: 61–70. 10.1002/prot.10183View ArticleGoogle Scholar
- Bahar I, Jernigan RL: Cooperative fluctuations and subunit communication in tryptophan synthase. Biochemistry 1999, 38: 3478–3490. 10.1021/bi982697vView ArticlePubMedGoogle Scholar
- Demirel MC, Atilgan AR, Jernigan RL, Erman B, Bahar I: Identification of kinetically hot residues in proteins. Protein Sci 1998, 7: 2522–2532.PubMed CentralView ArticlePubMedGoogle Scholar
- Keskin O, Jernigan RL, Bahar I: Proteins with Similar Architecture Exhibit Similar Large-Scale Dynamic Behavior. Biophys J 2000, 78: 2093–2106.PubMed CentralView ArticlePubMedGoogle Scholar
- Press WH, Teukolsky SA, Vetterling WT, Flannery BP: Numerical Recipes in Fortran. Volume Chp 2.6. Second edition. New York, Cambridge University Press; 1992:51–62.Google Scholar