- Research article
- Open Access
Role of protein interactions in stabilizing canonical DNA features in simulations of DNA in crowded environments
BMC Biophysics volume 11, Article number: 8 (2018)
Cellular environments are highly crowded with biological macromolecules resulting in frequent non-specific interactions. While the effect of such crowding on protein structure and dynamics has been studied extensively, very little is known how cellular crowding affects the conformational sampling of nucleic acids.
The effect of protein crowding on the conformational preferences of DNA (deoxyribonucleic acid) is described from fully atomistic molecular dynamics simulations of systems containing a DNA dodecamer surrounded by protein crowders. From the simulations, it was found that DNA structures prefer to stay in B-like conformations in the presence of the crowders. The preference for B-like conformations results from non-specific interactions of crowder proteins with the DNA sugar-phosphate backbone. Moreover, the simulations suggest that the crowder interactions narrow the conformational sampling to canonical regions of the conformational space.
The overall conclusion is that crowding effects may stabilize the canonical features of DNA that are most important for biological function. The results are complementary to a previous study of DNA in reduced dielectric environments where reduced dielectric environments alone led to a conformational shift towards A-DNA. Such a shift was not observed here suggested that the reduced dielectric response of cellular environments is counteracted by non-specific interactions with protein crowders under in vivo conditions.
Biological cells are highly crowded environments due to the presence of various macromolecules. The macromolecular crowding in cells plays a crucial role in biological processes as it may alter the structure and dynamics of biomolecules . A typical biological cell has a concentration of biomolecules in the range of 300–400 mg/ml , corresponding to a macromolecular volume fraction of 20–30% . Such an environment is substantially different from dilute solutions, the frequently considered environment in most biological experiments. Recent studies have begun to consider the effects of cellular crowding and have shed light on its effects on the structure and function of biomolecules [4,5,6,7,8,9]. Three essential crowding effects have been reported from experiments  and simulations : (1) the volume exclusion effect has been suggested to favor more compact conformations based on entropic arguments, thereby generally stabilizing more compact states [12, 13]; (2) non-specific interactions between biomolecules and surrounding protein crowders have led to the destabilization of native states [14,15,16] as well as reduced diffusion ; and (3) altered solvation properties including reduced dynamic and dielectric properties  have implied a reduced hydrophobic effect [19, 20].
While much attention so far has been on proteins, nucleic acids are also affected by macromolecular crowding [21, 22]. G-quadruplex DNA structure assumes a parallel-G quadruplex form under crowded environments due to the excluded volume effect as well as alterations in the hydration of DNA [23,24,25]. Long DNA duplexes undergo a collapsing transition in the presence of polyethylene glycol (PEG) in solution, which can also be explained by the volume exclusion effect favoring states that are more compact [26, 27]. The negatively charged protein bovine serum albumin (BSA) similarly causes a compaction of large DNA molecules due to the volume exclusion effect and repulsive electrostatic interactions . Short DNA duplexes, on the other hand, have been extensively investigated by both experimental techniques and computer simulations in terms of co-solvent, salt effects, and crystallization [29,30,31,32,33,34,35,36,37,38,39]. The DNA duplex is well-known to be most stable in the B-form  in aqueous solution and in A-form in environments with depleted water and for certain sequences . Very high concentrations of salt can also induce the B- to A- form transition by bringing the negatively charged phosphate groups of DNA closer [34, 38, 42,43,44] while the addition of ethanol favors the A-form due to reduced electrostatics [32, 33, 35, 37, 45,46,47]. More recently, the effect of reduced dielectric environments on DNA as one aspect of cellular crowding was investigated and has also been shown to favor non-canonical A-form structures in implicit solvent simulations . On the other hand, a study based on a coarse-grained model has suggested that even in the presence of significant crowding there may be a solvent-rich region around DNA that is depleted in crowder molecules, which was found to have an impact on the kinetics of proteins diffusing along DNA . However, to the extent that crowder proteins do interact non-specifically with DNA, the effect of explicit protein crowder molecules on DNA duplex structures is not well understood.
Here, we describe fully atomistic molecular dynamics (MD) simulations of DNA dodecamers in the presence of explicit protein crowders in order to investigate how DNA structure and stability may be affected under such conditions. We find a general tendency of the DNA to favor the B-form in crowded environments, which is in contrast to the shift towards A-form DNA observed in the simpler reduced dielectric environments . The stabilization of B-DNA appears to be due to non-specific protein-DNA interactions. We also observe, some alterations in the hydration structure and ion distributions around DNA under crowded conditions. The results are described in detail and discussed in the following after outlining the computational methods used in this study.
MD simulations of Drew-Dickerson ((CGCGAATTCGCG)2) and GC-rich (CGCCCCGCGGGCG)2) dodecamers in crowded protein environments were carried out using the CHARMM (Chemistry at Harvard Molecular Mechanics) program package (version 41a1)  with the CHARMM36 force field [50, 51]. The initial Drew-Dickerson dodecamer structure was obtained from the X-ray structure (PDB: 1BNA) , and the initial GC-rich dodecamer structure was obtained by mutating the base sequence in the X-ray structure of the Drew-Dickerson dodecamer using the MMTSB (Molecular Modeling Tools in Structural Biology) Tool Set . In experiment, the Drew-Dickerson dodecamer is very stable in B-form, in the crystal as well as in solution . The crystal structure of the GC-rich dodecamer is in A-form , but there is less known about its conformation in solution. Generally there is little evidence for A-DNA conformations in solution unless the salt concentration is much higher than typical physiological conditions [34, 38, 42,43,44] and/or when co-solvents such as ethanol are present in significant fractions [32, 33, 35, 37, 45,46,47]. Therefore, we setup both systems in B-form as the likely conformation of both sequences in dilute aqueous solvent. The choice of sequences and initial structures also allows a direct comparison with our previous continuum dielectric study .
For each dodecamer, a dilute system without crowders (0% crowder fraction) and three systems with different protein crowder volume fractions (20, 30, 40%) were prepared. Protein G (PDB: 1PGB)  was selected as the crowder protein due to its small size and stability in computer simulations . We used neutral protein G models molecules introduced in previous work , where D36, D40, E19 and E42 are protonated. In the previous study, both, the charged and neutralized variants of protein G were studied under crowded conditions similar to the systems studied here but without DNA and both were found to be stable in simulations . Protein G is not known to specifically interact with DNA and we chose the net-neutral form to reduce electrostatic interactions with the highly charged DNA to focus on more general crowding effects while still maintaining protein-like crowders. The crowded systems (20, 30, 40%) consisted of one dodecamer and 8 protein G molecules, whereas the dilute systems only contained one dodecamer. Simulation box sizes were varied between 53.2–61.3 Å to obtain the abovementioned crowder volume fractions. The box sizes were varied instead of the number of protein copies to achieve exactly the target crowder fractions and minimize computational costs at the higher concentrations as in previous work [15, 57]. Simulation conditions of the systems are given in Table 1. There is no experimental evidence for a specific DNA-protein G complex that is stable over long time and consequently the system is assumed to be fully dynamic in the liquid state with molecular interactions varying transiently. To avoid biasing towards any specific initial protein G-DNA interaction, the initial crowded systems were set up by randomly rotating and placing the DNA dodecamer and the crowder proteins in the simulation box using a protocol developed previously . Different replicates of each system had different initial orientations and placements of the DNA and the surrounding crowders. All systems were solvated with explicit TIP3P (three-site transferable intermolecular potential)  water molecules. To neutralize the DNA dodecamer, 22 sodium ions were added to the systems. In order to keep the ion molality of all systems the same, 6 and 12 additional pairs of sodium and chloride ions were added to 30% and 0/20% systems, respectively. Therefore, all systems had 0.45 mol/kg ion molality.
The initial systems were minimized for 1000 steps using the adopted bases Newton Raphson (ABNR) algorithm and were subsequently heated by running simulations without using any restraints at 50 K, 100 K, 150 K, 200 K, 250 K for 4 ps and at 298 K for 10 ps. Productions runs were carried out at 298 K in the NVT ensemble for 1 μs with a 2 fs time step. The SHAKE algorithm  was used to constrain bond lengths involving hydrogen atoms. Temperature control was obtained by a Langevin thermostat with a 0.01 ps− 1 friction coefficient. Lennard-Jones and direct electrostatic interactions were cut off at 12 Å with a switching function becoming effective at 10 Å. Electrostatic interactions were calculated from particle-mesh Ewald  summation using 1 Å grid spacing. All simulations were performed using periodic boundary conditions. For the crowded systems, five independent simulations were carried out starting from different initial orientations. For the dilute systems without protein crowders, simulations were replicated three times starting from different initial velocities for the atoms.
The analysis of the helicoidal and backbone parameters of the dodecamers (see Additional file 1: Figure S1) were performed by using the 3DNA program package . The reported values are averages over snapshots. Radial distribution functions and 3D volume densities were analyzed by using in-house scripts. All the other analysis was carried out using the MMTSB Tool Set  in combination with CHARMM . Clustering analysis was performed by applying the k-means clustering algorithm by using the kclust program in MMTSB . For each dodecamer, all snapshots from the simulations with different protein concentrations were aggregated and clustered by using a 3 Å clustering radius. Only the last 700 ns of the simulations were analyzed because of larger variations in the helicoidal parameters during the first 300 ns (see Additional file 1: Figure S2). Only the inner eight base-pairs were taken into consideration to ignore structural distortions due to base fraying. VMD (visual molecular dynamics)  and PyMOL  were used for visualization.
Microsecond-scale molecular dynamics simulations of DNA dodecamers with and without protein crowders were carried to study the effect of crowding on DNA structure. We focused our analysis on helical properties including base geometries, groove widths and DNA bending, backbone torsions, interactions with crowder proteins, correlations between protein contacts and helical properties, and water and ion distributions around DNA.
Snapshots from the simulations were clustered to identify major conformations. Representative structures for each of the major clusters (with more than 5% population) are depicted in Fig. 1. Generally, the helices stayed intact with base fraying at the termini, which is common in simulations of short DNA fragments . The structures generally resemble B-DNA structures for both sequences. Average root mean square deviation (RMSD) values of different clusters from the initial canonical B-DNA structures vary between 1.4 and 2.0 Å for the Drew-Dickerson dodecamer and between 1.6 to 2.6 Å except for one cluster with an RMSD of 3.7 Å for the GC-rich dodecamer (see Additional file 1: Table S1). There is no clear pattern of increasing or decreasing RMSD values for the clusters most populated at different crowder concentrations.
Averages over all base-pairs excluding the first and last two terminal base-pairs with errors given in parentheses based on the variations in the independent simulations. Canonical values were averaged over the A-form structures 3V9D, 3QK4, 2B1B, 1ZEX, 1ZEY, 1ZF1, 1ZF8, 1ZF9, 1ZFA and the B-form structures 2M2C, 4AGZ, 4H0, 4AH1, 3 U05, 3 U08, 1VTJ,3U2N, 3OIE, 3BSE.
see Table 2.
Helicoidal parameters for both, Drew-Dickerson and GC-rich dodecamers were averaged from the simulations. They are summarized in Tables 2 and 3, respectively. Helicoidal parameters for crystal structures of the respective dodecamers as well as canonical A- and B-forms of DNA, averaged over ten A- and B-form crystal structures each, are provided for comparison. Average properties for each of the clusters shown in Fig. 1 are given in Additional file 1: Tables S1 and S2. The more detailed analysis of the base geometries also indicates that both dodecamers remained close to B-DNA. The Drew-Dickerson dodecamer also remained reasonably close to the respective crystal structure (1BNA), but there are larger deviations between the simulation results and the crystal structure of the GC-rich dodecamer. The crystal structure for the GC-rich dodecamer is predominantly in A-form, presumably as a result of salt concentrations above 1 M and/or the crystal environment . As mentioned above, although the crystal structure of the GC-rich dodecamer has been reported in A-form, there is no evidence that this sequence (or any other DNA sequence) assumes an A-DNA conformation in solution at sub-molar salt concentrations and in the absence of co-solvents. Therefore we expected the GC-dodecamer to remain in B-form. In the presence of the protein crowders, the helical parameters generally did not change much. We found increased X-displacement (p-values: 0.91 (Drew-Dickerson 20%), 0.16 (Drew-Dickerson 30%), 0.05 (Drew-Dickerson 40%), 0.05 (GC-rich 20%), 0.04 (GC-rich 30%), 0.07 (GC-rich 40%)) and base inclination (p-values: 0.33 (Drew-Dickerson 20%), 0.38 (Drew-Dickerson 30%), 0.01 (Drew-Dickerson 40%), 0.05 (GC-rich 20%), 0.01 (GC-rich 30%), 0.08 (GC-rich 40%)) for both dodecamers as a function of crowding. The increased x-displacement and base inclination point towards A-DNA but the values upon crowding still remained much closer to canonical B-DNA than A-DNA.
We further analyzed the displacement of phosphorus atoms relative to the horizontal plane passing between base-pairs in a base-pair step (zp) and major/minor grooves (Tables 2 & 3). The zp parameter is very different between the two forms of DNA. While B-DNA has values near − 0.3 Å, the parameter is mostly larger than 2.0 Å for A-DNA. This parameter does not show a trend upon crowding for the GC-rich dodecamer, while the Drew-Dickerson dodecamer had larger values in crowded environments (p-values: 0.0011 (20%), 0.0007 (30%), 0.0001 (40%)). Again this indicates a slight tendency towards A-DNA geometries while still remaining much closer to canonical B-DNA values. Minor and major groove widths also did not change significantly upon crowding, but we note that minor groove widths were generally overestimated compared to canonical B-DNA values. This is a general feature of the CHARMM force field that was used here . Finally, we analyzed the helical bending angles (see Additional file 1: Table S3) which also did not show a significant change upon crowding.
Sugar conformations and backbone torsions
A key feature of nucleic acid backbone is the ribose pucker conformation. A-form DNA is known to prefer C3’-endo and C2’-exo conformations whereas B-form DNA is characterized by C3’-exo and C2’-endo conformations. As shown in Fig. 2, the sugars of both dodecamers generally remained in C3’-exo and C2’-endo conformations. As expected, C3’-endo and C2’-exo sugar conformations are more prominent for the GC-rich dodecamer (see Fig. 2b). Again, there is no major change upon crowding, but in the GC-rich dodecamer, sugars shift slightly to C3’-exo and C2’-endo sugar conformations up to 30% crowding, but then revert back to more A-form conformations at 40% crowder concentrations.
We further analyzed torsion angles along the phosphate backbone. χ and δ angles are the most distinctive backbone angles to distinguish between A- and B- form DNA. We constructed potentials of mean force (PMF) as a function of δ and χ from the simulations (Fig. 3). The separation between A- and B-DNA torsion angles is readily apparent. Consistent with the ribose puckers and helical geometries, there is more sampling of B-DNA torsion angles for both dodecamers. While there is little change in the sampling of the major A- and B-form, the presence of crowders appears to affect the sampling of minor conformations with A-like δ values around 80 degrees and B-like χ values around − 100 degrees. Sampling in this region is significantly reduced in both dodecamers upon crowding (see Fig. 3). This region corresponds to a conformation where bases stay in the same orientation relative to the sugar as in B-form, but they are slightly more exposed to the environment, and apparently, this conformation is largely prevented by crowder proteins. The sampling of ε and ζ torsion angles distinguishes between BI/BII forms. A similar trend is observed where crowding reduces the sampling of minor states outside the major BI/BII basins (see Additional file 1: Figure S3). Based on this analysis, it appears that one effect of protein crowders may be to focus the sampling of DNA conformations on the major conformations.
Protein crowder conformations
In previous simulation studies involving protein G under crowded conditions, the protein G conformations remained close to the experimental structure and were not affected strongly by the concentrated environment . In the systems studied here, protein G also remains highly stable and close to the experimental structure (see Additional file 1: Figure S4). The overall average Cα RMSD value is 0.91 Å with a standard deviation of 0.24 Å between individual protein G molecules and the experimental structure (PDB ID: 3GB1 [55, 67]). The average radius of gyration is 10.76 Å with a standard deviation of 0.1 A, compared to a value of 10.65 Å for the experimental structure. A few conformations deviated slightly further from the native (as much as 2.5 Å Cα RMSD) and with slightly increased radii of gyration, especially at the highest crowder concentration (Additional file 1: Figure S4B). Further analysis via clustering revealed minor substates with slightly increased RMSD values that correlate with closer contacts to the DNA (see Additional file 1: Table S4). This suggests that the conformational sampling of protein G may be affected slightly when interacting with the DNA. Almost all of the variations are in the flexible loop involving residues 9 to 13 (see Additional file 1: Figure S4D) where root mean square fluctuations (RMSF) are largest (see Additional file 1: Figure S4C).
Protein G is not known to interact specifically with DNA but under highly crowded conditions, interactions are unavoidable. Figure 4 shows where contacts between protein G and DNA occur based on minimum distances between the major/minor grooves and sugar/phosphate groups of the DNA with different residues of protein G. More detailed contact analysis between individual base-pairs and protein G residues is shown in Additional file 1: Figures S5 and S6 for the Drew-Dickerson and GC-rich dodecamers, respectively. Most of the contacts are between the DNA sugar-phosphate backbone and protein residues 15–30, mostly in the α-helix of protein G, as well as residues at the N-terminus and near the C-terminus. Contacts involving the DNA grooves, a typical mode of interaction for DNA-binding proteins were not common with protein G. The interactions partially involve electrostatic attraction between the DNA phosphate and certain lysine residues (K4, K28, K31, and K50), but sugar oxygens O3’ and O4’ as well as phosphate oxygens also form hydrogen bonds with other polar protein residues. Representative snapshots of protein G-DNA interactions are shown in Fig. 5. As would be expected, the contacts between the proteins and DNA increase with crowder concentration and crowding seems to increase sugar-phosphate-protein contacts more for the Drew-Dickerson dodecamer than for the GC-rich dodecamer.
Correlations between DNA-protein contacts and DNA helix properties
To investigate in more detail whether the close contacts of the crowder proteins with the DNA have the potential to perturb DNA structure, we analyzed correlations between DNA-protein contacts and helicoidal properties of DNA as well as backbone torsion and pseudorotation phase angles. First, we examined the effect of close contacts on the helicoidal parameters listed in Tables 2 and 3. We found that a higher number of close protein contacts corresponded to a more narrow range of sampled values for all of the helicoidal parameters (Figs. 6-9, Additional file 1: Figures S7-S10). Among these parameters, slide (Fig. 6), x-displacement (Fig. 7), helical rise (Fig. 8) and zp (Fig. 9) values showed a clear shift towards B-form values with increasing number of contacts. These parameters focus on the displacement of bases along the x- (x-displacement) and y- (slide) axes and of phosphates along the base-pair axis (zp). All of the values approach zero with crowding. This suggests that DNA bases and phosphates undergo less displacement as a result of crowding. On the other hand, rotations of base-pairs about helical (twist) or base-pair axes (inclination) do not show a distinct shift towards any canonical values (Additional file 1: Figures S7, S8). Major and minor groove widths do not seem to be affected by contacts except for the GC-rich dodecamer, where there appears to be a clear tendency towards larger minor groove values, i.e. values more similar to A-DNA (Additional file 1: Figures S9, S10).
Similar to the helicoidal parameters, backbone torsion angles also fluctuate in a more narrow range upon crowding (Additional file 1: Figures S11–S17). This suggests that non-specific protein-DNA interactions may limit the conformational fluctuations of the DNA backbone. Particularly, δ and χ angles shift towards B-form values upon higher number of protein contacts, explaining a decrease in the sampling of non-canonical conformations shown in Fig. 3. Finally, pseudorotation angles move to B-form values with protein contacts which lead to C3’-exo and C2’-endo sugar pucker conformations (Additional file 1: Figure S18).
The results discussed here are most pronounced for the Drew-Dickerson dodecamer. In the GC-rich dodecamer, the fluctuations of helicoidal parameters and backbone angles are reduced less and a tendency to sample A-form values further complicates the picture. Overall, our results suggest that the interactions of protein crowders with DNA sugar/phosphate backbone shown in the previous section result in a stiffer DNA backbone. The stiffer backbone also prevents larger base/base-pair displacements and, therefore, restricts the conformational space of DNA. Although it appears that there is not a specific tendency towards one of the major forms of DNA upon crowding, there is a distinct effect of protein crowders on DNA structure by narrowing the conformational sampling to canonical structures.
Hydration and ion distributions around DNA
Water and ions are integral parts of DNA structures. We analyzed hydration patterns and sodium ion distributions around DNA as a function of crowding. Conditional water radial distribution functions (RDF) were obtained for water oxygen distances to the closest heavy atoms in DNA, normalized by the corresponding accessible volume at each distance and the bulk water density (0.034 Å− 3) (see Fig. 10a). The analysis shows that the first hydration shell is almost unaffected by the level of crowding, but the RDF decreases beyond the hydration shell significantly as a function of crowding. This observation is similar to what has been reported previously for the hydration around proteins under crowded environments .
Sodium RDFs were calculated in the same way as the water RDFs but normalized by the ion density of the system (0.002 Å− 3). There are two peaks in the sodium RDFs corresponding to ions in direct contact with the DNA (around 2.5 Å and largely in the minor groove) and ions interacting with the DNA through water (around 4.5 Å) [68,69,70]. While the direct contact peak is not affected significantly by crowding, the second peak shows a greater dependence on crowding. At the highest crowder fractions, the second peak is significantly reduced in both dodecamers (see Fig. 10b) and the ion density is reduced further at larger distances similar to the reduction in hydration upon crowding. The effect of crowding on the ion distributions also impacts the DNA neutralization as a function of distance (Fig. 10c). 76% of the DNA phosphate groups are neutralized as suggested by counterion condensation theory at around 9 Å for the dilute system, however, it takes up to 11–12 Å to reach 76% DNA neutralization under crowding conditions. It is interesting, that despite the impact of crowding on the second peak of the ion distribution, the counterion condensation is affected less for distances less than 6 Å. As this may seem counterintuitive, the reader is reminded that the RDF is normalized by the available volume and the overall ion density, at constant ion molality, whereas Fig. 10c simply describes the net neutralization of the DNA by the ions. The extended distance to reach 76% charge neutralization upon crowding may seem to challenge counterion condensation theory. However, the protein crowders, despite being net neutral, can provide additional charge neutralization by orienting basic lysines near the DNA surface as described above to compensate for the reduced neutralization by the sodium ions.
Finally, the 3D distributions of sodium ions around the Drew-Dickerson and GC-rich dodecamers are compared in Fig. 11. The sodium ion networks in the major and minor grooves of DNA are largely preserved for both dodecamers with little changes upon crowding. However, additional densities become apparent further away from the DNA at different locations upon crowding. Additional ordering of ions could be a result of crowder proteins interacting with the DNA and coordinating ions near the DNA. A snapshot showing a crowder protein interacting with the DNA and orienting a sodium ion at the same time is shown in Additional file 1: Figure S19.
In this study, we investigated the effect of protein crowding on the conformational preferences of DNA duplexes. In a previous study, we examined one aspect of cellular crowding, namely a reduced dielectric response of the environment due to the less available water and its slowed dynamics. Using continuum models, we found an overall shift towards A-like conformations for DNA as a result of a reduced dielectric response of its environment . Here, we included protein crowders and solvent explicitly to test whether the same conclusions would be found. In the earlier work, we compared environments with ε = 20, ε = 40, and ε = 80. Past work suggests that water under crowded conditions exhibits a reduced dielectric response of about 40 (with uncertainties) at a protein crowder volume fraction of 0.3 . If one makes the further assumption that proteins have an interior dielectric of around 10, one can estimate an average effective dielectric for the entire medium surrounding the DNA at this crowder fraction as εeff = 0.3*10 + 0.7*40 = 31 with even lower values at 40% crowder fraction. However, although some of the base parameters moved slightly towards A-like values upon crowding, B-DNA was largely maintained with the explicit crowder environment in contrast to our previous findings. This suggests that a reduced dielectric response of crowded environments and interactions with crowder proteins have different effects on DNA conformational preferences with a net effect of not altering canonical B-DNA structures much.
We found that the crowder proteins mostly interact with DNA via its phosphate-sugar backbone as previously observed in non-specific binding of proteins to DNA . These interactions arise from the electrostatic interactions between negatively charged phosphate oxygens and positively charged amino acid residues as well as the polar interactions between phosphate and/or sugar oxygens and side chains of polar amino acid residues. Previous studies have shown that DNA can undergo structural deformations from its B-form towards A-type helix as a result of forming complexes with specific DNA binding proteins [72,73,74,75,76], but we did not see such an effect here. It does appear, however, that for the system studied here, the presence of the protein crowders limits the conformational space of DNA to more canonical structures, mostly in B-form, both for the backbone torsions and the helical parameters. However, the narrowed conformational sampling appears to have little effect on the overall structural averages. Such a crowding effect on DNA structure may be understood in similar ways as protein native state stabilization due to the volume exclusion effect [1, 13, 77, 78], where the reduced space due to crowders limits the ability to widely sample conformational space. This would mean that protein crowding in vivo helps stabilize the biologically most relevant forms of DNA.
We also studied hydration patterns and ion densities around DNA in protein crowding. The first hydration shell around DNA is largely unaffected by crowding, while the water densities beyond the first solvation shell significantly reduced compared to the bulk water density under crowding effect. This result is very similar to the hydration shell around proteins upon crowding . This further confirms that, protein crowding in cells generally does not alter the first hydration shell around biomolecules. However, sodium densities around DNA are affected already when interacting with DNA through water. Only the direct-contact first peak in the sodium-DNA RDF appears to be unaffected by crowding. Moreover, the charge neutralization by ions is altered upon crowding with the classical counter-ion condensation threshold reached at larger distances from the DNA than under dilute conditions. This suggests that proteins have to play an increasing role in neutralizing DNA under highly crowded conditions.
The results obtained here shed light on the effect of protein crowding on DNA structure. We found that the crowder proteins mostly assist DNA to stay in canonical B-like conformations, limiting excursions to non-canonical conformations rather than a clear shift in the overall, average structure as suggested by a simple dielectric model of cellular environments. We hope that this hypothesis will motivate new experimental efforts to characterize DNA structure under crowded conditions. We expect that reduced conformational dynamics upon crowding would be observable via NMR spectroscopy. Another testable hypothesis is the altered ion distribution predicted by our simulations, which could be amenable to the ion-counting experiments recently carried out by the Herschlag group [79,80,81,82].
Bovine serum albumin
Chemistry at Harvard Molecular Mechanics
Molecular Modeling Tools in Structural Biology
Potential of mean force
Radial distribution function
Root mean square deviation
Root mean square fluctuations
Three-site transferable intermolecular potential
Visual molecular dynamics
Zhou H-X, Rivas G, Minton AP. Macromolecular crowding and confinement: biochemical, biophysical, and potential physiological consequences. Annu Rev Biophys. 2008;37:375–97.
Zimmerman SB, Trach SO. Estimation of macromolecule concentrations and excluded volume effects for the cytoplasm of Escherichia coli. J Mol Biol. 1991;222:599–620.
Ellis RJ. Macromolecular crowding: an important but neglected aspect of the intracellular environment. Curr Opin Struct Biol. 2001;11:114–9.
Minton AP, Wilf J. Effect of macromolecular crowding upon the structure and function of an enzyme - Glyceraldehyde-3-phosphate dehydrogenase. Biochem. 1981;20:4821–6.
McPhie P, Ni YS, Minton AP. Macromolecular crowding stabilizes the molten globule form of Apomyoglobin with respect to both cold and heat unfolding. J Mol Biol. 2006;361:7–10.
Dix JA, Verkman AS. Crowding effects on diffusion in solutions and cells. Annu Rev Biophys. 2008;37:247–63.
Minton AP. Macromolecular crowding and molecular recognition. J Mol Recognit. 1993;6:211–4.
Zimmerman SB, Minton AP. Macromolecular crowding - biochemical, biophysical, and physiological consequences. Annu Rev Biophys Biomol Struct. 1993;22:27–65.
Elcock AH. Models of macromolecular crowding effects and the need for quantitative comparisons with experiment. Curr Opin Struct Biol. 2010;20:196–206.
Rivas G, Minton AP. Macromolecular crowding in vitro, in vivo, and in between. Trends Biochem Sci. 2016;41:970–81.
Feig M, Yu I, Wang P-h, Nawrocki G, Sugita Y. Crowding in cellular environments at an atomistic level from computer simulations. J Phys Chem B. 2017;121:8009–25.
Zhou HX. Protein folding and binding in confined spaces and in crowded solutions. J Mol Recognit. 2004;17:368–75.
Cheung MS, Klimov D, Thirumalai D. Molecular crowding enhances native state stability and refolding rates of globular proteins. Proc Natl Acad Sci U S A. 2005;102:4753–8.
Senske M, Tork L, Born B, Havenith M, Herrmann C, Ebbinghaus S. Protein stabilization by macromolecular crowding through enthalpy rather than entropy. J Am Chem Soc. 2014;136:9036–41.
Harada R, Tochio N, Kigawa T, Sugita Y, Feig M. Reduced native state stability in crowded cellular environment due to protein-protein interactions. J Am Chem Soc. 2013;135:3696–701.
Feig M, Sugita Y. Variable interactions between protein Crowders and biomolecular solutes are important in understanding cellular crowding. J Phys Chem B. 2012;116:599–605.
Wang Y, Li C, Pielak GJ. Effects of proteins on protein diffusion. J Am Chem Soc. 2010;132:9392–7.
Harada R, Sugita Y, Feig M. Protein crowding affects hydration structure and dynamics. J Am Chem Soc. 2012;134:4842–9.
Tanizaki S, Clifford J, Connelly BD, Feig M. Conformational sampling of peptides in cellular environments. Biophys J. 2008;94:747–59.
Yildirim A, Sharma M, Varner BM, Fang L, Feig M. Conformational preferences of DNA in reduced dielectric environments. J Phys Chem B. 2014;118:10874–81.
S-i N, Miyoshi D, Sugimoto N. Effects of molecular crowding on the structures, interactions, and functions of nucleic acids. Chem Rev. 2014;114:2733–58.
S-i N, Sugimoto N. Model studies of the effects of intracellular crowding on nucleic acid interactions. Mol Biosyst. 2017;13:32–41.
Xue Y, Kan ZY, Wang Q, Yao Y, Liu J, Hao YH, Tan Z. Human Telomeric DNA forms parallel-stranded intramolecular G-Quadruplex in K+ solution under molecular crowding condition. J Am Chem Soc. 2007;129:11185–91.
Miyoshi D, Nakao A, Sugimoto N. Molecular Crowding Regulates the Structural Switch of the DNA G-Quadruplex. Biochem. 2002;41:15017–24.
Heddi B, Phan AT. Structure of human Telomeric DNA in crowded solution. J Am Chem Soc. 2011;133:9824–33.
Livolant F, Leforestier A. Condensed phases of DNA: structures and phase transitions. Progr Polym Sci. 1996;21:1115–64.
Bloomfield VA. DNA condensation. Curr Opin Struct Biol. 1996;6:334–41.
Yoshikawa K, Hirota S, Makita N, Yoshikawa Y. Compaction of DNA induced by like-charge protein: opposite salt-effect against the polymer-salt-induced condensation with neutral polymer. J Phys Chem Lett. 2010;1:1763–6.
Kulkarni M, Mukherjee A. Understanding B-DNA to A-DNA transition in the right-handed DNA Helix: perspective from a local to global transition. Progr Biophys Mol Biol. 2017;128:63–73.
Ivanov VI, Minchenkova LE, Minyat EE, Frank-Kamenetskii MD, Schyolkina AK. The B to a transition of DNA in solution. J Mol Biol. 1974;87:817–33.
Jose D, Porschke D. Dynamics of the B-A Transition of DNA Double Helices. Nucleic Acids Res. 2004;32:2251–8.
Zimmerman SB, Pheiffer BH. A direct demonstration that the ethanol-induced transition of DNA is between the a and B forms: an X-ray diffraction study. J Mol Biol. 1979;135:1023–7.
Cheatham TE, Crowley MF, Fox T, Kollman PA. A molecular level picture of the stabilization of A-DNA in mixed ethanol-water solutions. Proc Natl Acad Sci U S A. 1997;94:9626–30.
Cheatham TE, Kollman PA. Insight into the stabilization of A-DNA by specific ion association: spontaneous B-DNA to A-DNA transitions observed in molecular dynamics simulations of d(ACCCGCGGGT)2 in the presence of Hexaamminecobalt(III). Structure. 1997;15:1297–311.
Noy A, Perez A, Laughton CA, Orozco M. Theoretical Study of Large Conformational Transitions in DNA: The B<−->A Conformational Change in Water and Ethanol/Water. Nucleic Acids Res. 2007;35:3330–8.
Pastor N. The B- to A-DNA transition and the reorganization of solvent at the DNA surface. Biophys J. 2005;88:3262–75.
Gu B, Zhang FS, Wang ZP, Zhou HY. Solvent-induced DNA conformational transition. Phys Rev Lett. 2008;100:088104.
Arscott PG, Ma C, Wenner JR, Bloomfield VA. DNA Condensation by Cobalt Hexaammine (III) in Alcohol-Water Mixtures: Dielectric Constant and Other Solvent Effects. Biopolymers. 1995;36:345–64.
Kuzmanic A, Dans PD, Orozco M. An in-depth look at DNA crystals through the prism of molecular dynamics simulations. BioRxiv. 2018;413336. https://doi.org/10.1101/413336.
Watson JD, Crick FH. The structure of DNA Cold Spring Harbor Symp. Quant Biol. 1953;18:123–31.
Franklin RE, Gosling RG. Molecular Configuration in Sodium Thymonucleate. Nature. 1953;171:740–1.
Gao Y-G, Robinson H, Boom JHV, Wang AH-J. Influence of counter-ions on the crystal structures of DNA Decamers: binding of [co(NH3)6]3+ and Ba2+ to A-DNA. Biophys J. 1995;69:559–68.
Robinson H, Wang AH. Neomycin, Spermine and Hexaamminecobalt (III) Share Common Structural Motifs in Converting B- to A-DNA. Nucleic Acids Res. 1996;24:676–82.
Xu QW, Shoemaker RK, Braunlin WH. Induction of B-A transitions of Deoxyoligonucleotides by multivalent cations in dilute aqueous-solution. Biophys J. 1993;65:1039–49.
Fang Y, Spisz TS, Hoh JH. Ethanol-Induced Structural Transitions of DNA on Mica. Nucleic Acids Res. 1999;27:1943–9.
Young MA, Beveridge DL. Molecular dynamics simulations of an oligonucleotide duplex with adenine tracts phased by a full Helix turn. J Mol Biol. 1998;281:675–87.
Srinivasan J, Cheatham TE, Cieplak P, Kollman PA, Case DA. Continuum solvent studies of the stability of DNA, RNA, and Phosphoramidate-DNA helices. J Am Chem Soc. 1998;120:9401–9.
Dey P, Bhattacherjee A. Role of Macromolecular Crowding on the Intracellular Diffusion of DNA Binding Proteins. Sci Rep Uk. 2018;8:844.
Brooks BR, Brooks CL, Mackerell AD, Nilsson L, Petrella RJ, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner AR, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor RW, Post CB, Pu JZ, Schaefer M, Tidor B, Venable RM, Woodcock HL, Wu X, Yang W, York DM, Karplus M. CHARMM: the biomolecular simulation program. J Comput Chem. 2009;30:1545–614.
Hart K, Foloppe N, Baker CM, Denning EJ, Nilsson L, Mackerell AD Jr. Optimization of the CHARMM Additive Force Field for DNA: Improved Treatment of the BI/BII Conformational Equilibrium. J Chem Theory Comput. 2012;8:348–62.
Best RB, Zhu X, Shim J, Lopes P, Mittal J, Feig M, MacKerell AD Jr. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone ϕ, ψ and side-chain χ1 and χ2 dihedral angle. J Chem Theory Comput. 2012;8:3257–73.
Drew HR, Wing RM, Takano T, Broka C, Tanaka S, Itakura K, Dickerson RE. Structure of a B-DNA Dodecamer: conformation and dynamics. Proc Natl Acad Sci U S A. 1981;78:2179–83.
Feig M, Karanicolas J, Brooks CL III. MMTSB tool set: enhanced sampling and multiscale modeling methods for applications in structural biology. J Mol Graph Modell. 2004;22:377–95.
Malinina L, Fernandez LG, Huynh-Dinh T, Subirana JA. Structure of the d(CGCCCGCGGGCG) Dodecamer: a kinked A-DNA molecule showing some B-DNA features. J Mol Biol. 1999;285:1679–90.
Gronenborn AM, Filpula DR, Essig NZ, Achari A, Whitlow M, Wingfield PT, Clore GM. A novel Highly Stable Fold of the Immunoglobulin Binding Domain of Streptococcal Protein G. Science. 1991;253:657–61.
Wang P-H, Yu I, Feig M, Sugita Y. Influence of protein Crowder size on hydration structure and dynamics in macromolecular crowding. Chem Phys Lett. 2017;671:63–70.
Nawrocki G, Wang PH, Yu I, Sugita Y, Feig M. Slow-down in crowded protein solutions correlates with transient oligomer formation. J Phys Chem B. 2017;121:11072–84.
Feig M, Harada R, Mori T, Yu I, Takahashi K, Sugita Y. Complete atomistic model of a bacterial cytoplasm integrates physics, biochemistry, and systems biology. J Mol Graph Modell. 2015;58:1–9.
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983;79:926–35.
Ryckaert JP, Ciccotti G, Berendsen HJC. Numerical-integration of Cartesian equations of motion of a system with constraints - molecular-dynamics of N-alkanes. J Comput Phys. 1977;23:327–41.
Darden TA, York D, Pedersen LG. Particle-mesh Ewald: an N log(N) method for Ewald sums in large systems. J Chem Phys. 1993;98:10089–92.
Lu XJ, Olson WK. 3DNA: A Software Package for the Analysis, Rebuilding and Visualization of Three-Dimensional Nucleic Acid Structures. Nucleic Acids Res. 2003;31:5108–21.
Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14:33–8.
Amadei A, Apol MEF, Berendsen HJC. On the use of the quasi-Gaussian entropy theory in noncanonical ensembles. I. Prediction of temperature dependence of thermodynamic properties. J Chem Phys. 1998;109:3004–16.
Zgarbová M, Otyepka M, Šponer J, Lankaš F, Jurečka P. Base Pair Fraying in Molecular Dynamics Simulations of DNA and RNA. J Chem Theory Comput. 2014;10:3177–89.
Cheatham TE, Young MA. Molecular dynamics simulation of nucleic acids: successes Limitations, and Promise. Biopolymers. 2001;56:232–56.
Kuszewski J, Gronenborn AM, Clore GM. Improving the packing and accuracy of NMR structures with a pseudopotential for the radius of gyration. J Am Chem Soc. 1999;121:2337–8.
Feig M, Pettitt BM. Sodium and chlorine ions as part of the DNA solvation Shell. Biophys J. 1999;77:1769–81.
Korolev N, Lyubartsev AP, Laaksonen A, Nordenskiold L. On the competition between water, sodium ions, and Spermine in binding to DNA: a molecular dynamics computer simulation study. Biophys J. 2002;82:2860–75.
Bonvin AM. Localisation and dynamics of sodium Counterions around DNA in solution from molecular dynamics simulation. Eur Biophys J. 2000;29:57–60.
Kalodimos CG, Biris N, Bonvin AM, Levandoski MM, Guennuegues M, Boelens R, Kaptein R. Structure and Flexibility Adaptation in Nonspecific and Specific Protein-DNA Complexes. Science. 2004;305:386–9.
Nekludova L, Pabo CO. Distinctive DNA conformation with enlarged major groove is found in Zn-finger-DNA and other protein-DNA complexes. Proc Natl Acad Sci U S A. 1994;91:6948–52.
Shakked Z, Guzikevich-Guerstein G, Frolow F, Rabinovich D, Joachimiak A, Sigler PB. Determinants of Repressor/Operator Recognition from the Structure of the Trp Operator Binding Site. Nature. 1994;368:469–73.
Guzikevich-Guerstein G, Shakked Z. A novel form of the DNA double Helix imposed on the TATA-box by the TATA-binding protein. Nat Struct Biol. 1996;3:32–7.
Olson WK, Gorin AA, Lu XJ, Hock LM, Zhurkin VB. DNA sequence-dependent deformability deduced from protein-DNA crystal complexes. Proc Natl Acad Sci U S A. 1998;95:11163–8.
Lu XJ, Shakked Z, Olson WK. A-form conformational motifs in ligand-bound DNA structures. J Mol Biol. 2000;300:819–40.
Stagg L, Zhang SQ, Cheung MS, Wittung-Stafshede P. Molecular crowding enhances native structure and stability of α/β protein Flavodoxin. Proc Natl Acad Sci U S A. 2007;104:18976–81.
Minton AP. Implications of macromolecular crowding for protein assembly. Curr Opin Struct Biol. 2000;10:34–9.
Gebala M, Bonilla S, Bisaria N, Herschlag D. Does cation size affect occupancy and electrostatic screening of the nucleic acid ion atmosphere? J Am Chem Soc. 2016;138:10925–34.
Gebala M, Giambasu GM, Lipfert J, Bisaria N, Bonilla S, Li G, York DM, Herschlag D. Cation-anion interactions within the nucleic acid ion atmosphere revealed by ion counting. J Am Chem Soc. 2015;137:14705–15.
Giambasu GM, Luchko T, Herschlag D, York DM, Case DA. Ion counting from explicit-solvent simulations and 3D-RISM. Biophys J. 2014;106:883–94.
Allred BE, Gebala M, Herschlag D. Determination of ion atmosphere effects on the nucleic acid electrostatic potential and ligand association using AH+.C wobble formation in double-stranded DNA. J Am Chem Soc. 2017;139:7540–8.
We thank Prof. Charles Hoogstraten for discussions and exploratory analyses of the GC-rich dodecamer under different solution conditions.
Funding from NSF (MCB 1330560 and MCB 1817307) and NIH (R01 GM092949 and R35 GM126948) to carry out all aspects of the research, analysis, and writing is acknowledged. The funding bodies did not play any direct role in the design of the study, the collection, analysis, or interpretation of data, or in writing the manuscript.
Availability of data and materials
Initial conformations, simulation data, and results are available from the authors upon request.
Corresponding author: Michael Feig; 603 Wilson Road, Room BCH 218; East Lansing, MI 48824; USA; Phone: 517–432-7439; E-mail: email@example.com
Ethics approval and consent to participate
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Definition of backbone torsions and helicoidal parameters (Figure S1). Time series of helicoidal parameters (Figure S2). Potential of mean force (kcal/mol) as a function of backbone angles (Figure S3). Analysis of the protein G crowder conformational sampling (Figure S4). Average minimum distances between the crowder protein residues and DNA (Figures S5-S6). Potentials of mean force as a function of helicoidal parameters and protein contacts (Figures S7-S18). A snapshot for the crowder protein interacting with the DNA and sodium (Figure S19). Helicoidal parameters for the clusters (Table S1 and S2), bending angles for the dodecamers (Table S3), clustering analysis of protein G crowders (Table S4), and PMF error analysis (Table S5). (PDF 5597 kb)
About this article
Cite this article
Yildirim, A., Brenner, N., Sutherland, R. et al. Role of protein interactions in stabilizing canonical DNA features in simulations of DNA in crowded environments. BMC Biophys 11, 8 (2018). https://doi.org/10.1186/s13628-018-0048-y
- Molecular dynamics
- Protein G
- Conformational sampling
- Solvent interactions