Explicit ion modeling predicts physicochemical interactions for chromatin organization

  1. Xingcheng Lin
  2. Bin Zhang  Is a corresponding author
  1. Department of Chemistry, Massachusetts Institute of Technology, United States
12 figures, 3 tables and 1 additional file

Figures

Illustration of the residue-level coarse-grained explicit ion model for chromatin simulations.

The left panel presents a snapshot for the simulation box of a 147 bp nucleosome in a solution of 100 mM NaCl and 0.5 mM MgCl2. The nucleosomal DNA and histone proteins are colored in red and white, …

Figure 2 with 1 supplement
Explicit ion modeling reproduces the energetics of nucleosomal DNA unwrapping.

(A) Illustration of the umbrella simulation setup using the end-to-end distance between two DNA termini as the collective variable. The same color scheme as in Figure 1 is adopted. Only ions close …

Figure 2—figure supplement 1
The explicit ion model predicts the binding affinities of protein-DNA complexes well, related to Figure 1 of the main text.

Experimental and simulated binding free energies are compared for nine protein-DNA complexes (Privalov et al., 2011), with a Pearson correlation coefficient of 0.6. The PDB ID for each complex is …

Figure 3 with 1 supplement
Explicit ion modeling predicts salt-dependent conformations of a 12-mer nucleosome array.

(A) Top: Comparison of simulated and experimental (Correll et al., 2012) sedimentation coefficients of chromatin at different salt concentrations. Bottom: Number of DNA charges neutralized by bound …

Figure 3—figure supplement 1
Probability distributions used to compute means and standard deviations of the quantities presented in Figure 3 of the main text.

(A) Probability distribution of sedimentation coefficients calculated from the simulation with Na+ ions. (B) Probability distribution of sedimentation coefficients calculated from the simulation …

Figure 4 with 4 supplements
Close contacts give rise to strong inter-nucleosomal interactions.

(A) Illustration of the simulation protocol employed to mimic the nucleosome unbinding pathway dictated by the DNA-origami device (Funke et al., 2016). The three configurations, A1, A2, and A3, …

Figure 4—figure supplement 1
Illustration of the restrained two nucleosome simulations setup, related to Figure 4 of the main text.

(A) Schematics of the DNA-origami-based force spectrometer, reproduced from Figure 1 of Funke et al., 2016. (B) Schematics for the spatial restraints imposed on nucleosomes in our simulations to …

Figure 4—figure supplement 2
Explicit ion modeling reproduces the experimental free energy profiles of nucleosome binding.

(A) Comparison between the simulated (black) and experimental (red) free energy profile as a function of the inter-nucleosome distance. Error bars were computed as the standard deviation of three …

Figure 4—figure supplement 3
Compared with DNA-origami-restrained simulations, the unrestrained simulations produce more histone-DNA contacts across nucleosomes, related to Figure 4 of the main text.

The average number of inter-nucleosome contacts between DNA and histone tails (A) or histone cores (B) is plotted as a function of the distance r. The error bars were estimated as the standard …

Figure 4—figure supplement 4
The unrestricted simulations favor a smaller angle θ between two nucleosomal planes compared to the DNA-origami-restrained simulations, related to Figure 4 of the main text.

(A) Illustration of the collective variables used in the umbrella-sampling simulation. θ is the angle between two nucleosomal planes, and r is the distance between the geometric centers of two …

Figure 5 with 3 supplements
Simulations predict significant inter-nucleosome interactions at physiological conditions.

(A) Illustration of the collective variable, θ, defined as the angle between two nucleosomal planes, and r defined as the distance between the nucleosome geometric centers. w1 and w2 represent the …

Figure 5—figure supplement 1
Dependence of inter-nucleosome interactions on the DNA sequence, related to Figure 5 of the main text.

See text section ‘Simulations at the physiological salt concentration’ for further discussions on simulation details. (A) Illustration of the collective variables used in umbrella-sampling …

Figure 5—figure supplement 2
The poly-dA:dT sequence produces a higher number of cross-nucleosome histone-DNA contacts compared to the poly-dG:dC sequence, related to Figure 5 of the main text.

(A) The average number of inter-nucleosome contacts between histone proteins and nucleosomal DNA is plotted as a function of the distance r between the geometric centers of two nucleosomes. The …

Figure 5—figure supplement 3
Free energy profiles for the interactions between a pair of nucleosomes at different nucleosome repeat lengths (NRL) and in the presence of the linker histone H1.0, related to Figure 5 of the main text.

See text section ‘Simulations at the physiological salt concentration’ for further discussions on simulation details. (A) Illustration of the collective variables used in the umbrella-sampling …

Appendix 1—figure 1
A cutoff value of 8.0 Å produces more accurate values for the number of unwrapped DNA base pairs as determined from visual inspection of representative configurations, related to Figure 2 of the main text.

See text section ‘Number of unwrapped DNA base pairs’ for additional discussions. A typical nucleosome structure with most of the outer layer DNA unwrapped was used to examine the impact of …

Author response image 1
Comparison between the radial distribution functions of Na+ (left) and Mg2+ (right) ions around the DNA phosphate groups computed from all-atom (black) and coarse-grained (red) simulations.

Figure reproduced from Figure 4 of [1]. The coarse-grained explicit ion model used in producing the red curves is identical to the one presented in the current manuscript.

© 2011, AIP Publishing. This figure is reproduced with permission from Figure 4 in Freeman GS, Hinckley DM, de Pablo JJ (2011) A coarse-grain three-site-pernucleotide model for DNA with explicit ions. The Journal of Chemical Physics 135:165104. It is not covered by the CC-BY 4.0 license and further reproduction of this figure would need permission from the copyright holder.

Author response image 2
Three-dimensional distribution of sodium ions around the nucleosome determined from all-atom explicit solvent simulations.

Darker blue colors indicate higher sodium density and high density of sodium ions around the DNA is clearly visible. The crystallographically identified acidic patch has been highlighted as spheres …

Author response image 3
Time trace of the radius of gyration (Rg) of a nucleosome with the 601-sequence along an unbiased, equilibrium trajectory.

It is evident the Rg fluctuates around the value found in the PDB structure (3.95 nm), supporting the stability of the nucleosome in our simulation.

Author response image 4
The explicit-ion model predicts the binding affinities of protein-DNA complexes well, related to Fig.

1 of the main text. Experimental and simulated binding free energies are compared for nine protein-DNA complexes [cite], with a Pearson Correlation coefficient of 0.6. The PDB ID for each complex is …

Author response image 5
Revised main text Figure 5, with Figure 5D modified for improved visual clarity.
Author response image 6
Explicit ion modeling reproduces the experimental free energy profiles of nucleosome binding.

(A) Comparison between the simulated (black) and experimental (red) free energy profile as a function of the inter-nucleosome distance. Error bars were computed as the standard deviation of three …

Tables

Table 1
Summary of parameters used to describe interactions between ions and charged particles.

See text section ‘Coarse-grained explicit ion model’ for definitions of various parameters.

Coarse-grained pairϵ(kcal/mol)σ(Å)rmϵ(Å)σϵ(Å)H1(kcal/mol)rmh1(Å)σh1(Å)H2(kcal/mol)rmh2(Å)σh2(Å)
P-P0.183796.866.860.5
Na+-P0.025104.143.441.253.154884.10.570.478016.50.4
Na+-AA+*0.2394.0653.441.253.154884.10.57
Na+-AA0.2394.0653.441.253.154884.10.570.478016.50.4
Mg2+-P0.11954.873.751.01.290636.10.50.979928.31.2
Mg2+-AA+0.2393.5563.751.01.290636.10.5
Mg2+-AA0.2393.5563.751.01.290636.10.50.979928.31.2
Cl-P0.081215.54254.20.50.836526.71.5
Cl-AA+0.2394.87254.20.50.836526.71.50.478015.60.4
Cl-AA0.2394.87254.20.50.836526.71.5
Na+-Na+0.011212.432.70.570.179255.80.57
Na+-Mg2+0.049712.372.370.5
Na+-Cl0.083873.13523.92.065.497133.30.570.478015.60.4
Mg2+-Mg2+0.894601.4121.4120.5
Mg2+-Cl0.497374.744.480.571.099435.480.440.059758.160.35
Cl-Cl0.035854.0454.20.560.239016.20.5
  1. *

    Positive amino acids.

  2. Negative amino acids.

Table 2
Summary of parameters used to describe the WCA interactions between ions and neutral particles.

See text section ‘Coarse-grained explicit ion model’ for definitions of various parameters.

Coarse-grained pairϵ(kcal/mol)σ(Å)
Na+-S*0.2394.315
Na+-A0.2393.915
Na+-T0.2394.765
Na+-G§0.2393.665
Na+-C0.2394.415
Na+-AA**0.2394.065
Mg2+-S0.2393.806
Mg2+-A0.2393.406
Mg2+-T0.2394.256
Mg2+-G0.2393.156
Mg2+-C0.2393.906
Mg2+-AA**0.2393.556
Cl-S0.2395.1225
Cl-A0.2394.7225
Cl-T0.2395.5725
Cl-G0.2394.4725
Cl-C0.2395.2225
Cl-AA**0.2394.8725
  1. *

    Sugar.

  2. Adenine base.

  3. Thymine base.

  4. §

    Guanine base.

  5. Cytosine base.

  6. **

    Non-charged amino acids.

Table 3
Summary of simulation setups used in this study.

Additional simulation details can be found in text section ‘Molecular dynamics simulation details’.

StudiesBox size (nm3)Number of Na+Number of Mg2+Number of Cl
Single nucleosome 100 mM NaCl+0.5 mM MgCl2216,00013,0176513,003
Twelve nucleosomes 5 mM NaCl1,331,000619604006
Twelve nucleosomes 150 mM NaCl216,00021,695019,505
Twelve nucleosomes 0.6 mM MgCl23,375,000023142438
Twelve nucleosomes 1 mM MgCl23,375,000031274064
Two 147 bp 601-seq nucleosomes 35 mM NaCl+11 mM MgCl2125,00029228284290
Two 147 bp 601-seq nucleosomes 150 mM NaCl+2 mM MgCl2216,00019,50526019,737
Two 147 bp poly-dA:dT nucleosomes 150 mM NaCl+2 mM MgCl2216,00019,50526019,737
Two 147 bp poly-dG:dC nucleosomes 150 mM NaCl+2 mM MgCl2216,00019,50526019,737
Two 167 bp 601-seq nucleosomes 150 mM NaCl+2 mM MgCl2216,00019,50526019,657
Two 167 bp 601-seq nucleosomes with H1.0 150 mM NaCl+2 mM MgCl2216,00019,50526019,763

Additional files

Download links